StripClusterProcessor.h

#pragma once
 
#include <map>
#include <memory>
#include <string>
 
#include "Framework/Configure/Parameters.h"
#include "Framework/Event.h"
#include "Tracking/Digitization/SiStripConstants.h"
#include "Tracking/Digitization/StripClusterer.h"
#include "Tracking/Reco/TrackerDaqMap.h"
#include "Tracking/Reco/TrackingGeometryUser.h"
 
namespace tracking::reco {
 
class StripClusterProcessor : public TrackingGeometryUser {
 public:
  StripClusterProcessor(const std::string& name, framework::Process& process);
  virtual ~StripClusterProcessor() = default;
 
  void configure(framework::config::Parameters& parameters) override;
  void onProcessStart() override;
  void produce(framework::Event& event) override;
 
 private:
  // I/O
  std::string in_collection_{"FittedSiStripHits"};
  std::string in_pass_{""};
  std::string out_collection_{"StripMeasurements"};
 
  // Optional DAQ map: when set, the local-U centre offset for a layer uses that
  // sensor's real strip count instead of the fixed N_READOUT_STRIPS constant.
  // Empty (the default, MC) keeps the constant and leaves the MC path
  // unchanged.
  std::string daq_map_file_{""};
  std::map<int, int> layer_n_strips_;  
 
  // Clustering parameters (forwarded to StripClusterer)
  double seed_threshold_{4.0};
  double neighbor_threshold_{3.0};
  double cluster_threshold_{4.0};
  double mean_time_ns_{0.0};
  double time_window_ns_{-1.0};
  double neighbor_delta_t_ns_{-1.0};
  double max_chi2_ndf_{-1.0};
 
  std::unique_ptr<tracking::digitization::StripClusterer> clusterer_;
};
 
}  // namespace tracking::reco