TestBeamClusterProducer.h

Go to the documentation of this file.

 
#ifndef TRIGSCINT_TESTBEAMCLUSTERPRODUCER_H
#define TRIGSCINT_TESTBEAMCLUSTERPRODUCER_H
 
// LDMX Framework
#include "Framework/Configure/Parameters.h"  // Needed to import parameters from configuration file
#include "Framework/Event.h"
#include "Framework/EventProcessor.h"  //Needed to declare processor
#include "TrigScint/Event/TestBeamHit.h"
#include "TrigScint/Event/TrigScintCluster.h"
 
namespace trigscint {
 
class TestBeamClusterProducer : public framework::Producer {
 public:
  TestBeamClusterProducer(const std::string& name, framework::Process& process)
      : Producer(name, process) {}
 
  virtual void configure(framework::config::Parameters& ps);
 
  virtual void produce(framework::Event& event);
 
  virtual void addHit(uint idx, trigscint::TestBeamHit hit);
 
  virtual void onProcessStart();
 
  virtual void onProcessEnd();
 
 private:
  // collection of clusters produced
  std::vector<ldmx::TrigScintCluster> clusters_;
 
  // cluster seeding threshold
  double seed_{0.};
 
  // min threshold for adding a hit to a cluster
  double min_thr_{0.};
 
  // max number of neighboring hits to combine when forming a cluster
  int max_width_{2};
 
  // max channel number (to avoid unused channels)
  int max_channel_id_{11};
 
  // specific verbosity of this producer
  int verbose_{0};
 
  // expected arrival time of hits in the pad [ns]
  double pad_time_{0.};
 
  // maximum allowed delay for hits to be considered for clustering
  double time_tolerance_{0.};
 
  // input collection (hits)
  std::string input_collection_;
 
  // output collection (clusters)
  std::string output_collection_;
 
  // specific pass name to use for track making
  std::string pass_name_{""};
 
  // cluster channel nb centroid (will be content weighted)
  float centroid_{0.};
 
  // energy (edep), PE, or sth
  float val_{0.};
 
  // edep content, only; leave val_ for PE
  float val_e_{0.};
 
  // book keep which channels have already been added to the cluster at hand
  std::vector<unsigned int> v_added_indices_;
 
  // book keep which channels have already been added to any cluster
  std::vector<unsigned int> v_used_indices_;
 
  // fraction of cluster energy deposition associated with beam electron sim
  // hits
  // -- could convert this to instead be a "cleanb frac"; fraction of cluster
  // energy coming from clean hits
  float beam_e_{0.};
 
  bool do_clean_hits_{false};
 
  // cluster time (energy weighted based on hit time)
  float time_{0.};
 
  // empty map container
  std::map<int, int> hit_channel_map_;
};
 
}  // namespace trigscint
 
#endif /* TRIGSCINT_TESTBEAMCLUSTERPRODUCER_H */