LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
CosmicClusterTagger_module.cc
Go to the documentation of this file.
1 // Class: CosmicTagger
3 // Module Type: producer
4 // File: CosmicTagger_module.cc
5 //
6 // Generated at Mon Sep 24 18:21:00 2012 by Sarah Lockwitz using artmod
7 // from art v1_02_02.
8 // artmod -e beginJob -e reconfigure -e endJob producer trkf::CosmicTagger
10 
15 
16 #include <algorithm> // std::minmax() ...
17 #include <memory> // std::unique_ptr<>
18 #include <string>
19 #include <utility> // std::pair<>, std::move()
20 
29 
30 namespace cosmic {
31  class CosmicClusterTagger;
32 }
33 
35 public:
36  explicit CosmicClusterTagger(fhicl::ParameterSet const& p);
37 
38  void produce(art::Event& e) override;
39 
40 private:
41  std::string fClusterModuleLabel;
46 };
47 
49 {
50  auto const* geo = lar::providerFrom<geo::Geometry>();
51 
52  auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
53  auto const detp =
55  double const samplingRate = sampling_rate(clock_data);
56  fClusterModuleLabel = p.get<std::string>("ClusterModuleLabel", "cluster");
57  fTickLimit = p.get<int>("TickLimit", 0);
58  const double driftVelocity = detp.DriftVelocity(detp.Efield(), detp.Temperature()); // cm/us
59 
61  2 * geo->DetHalfWidth() / (driftVelocity * samplingRate / 1000); // ~3200 for uB
62  fMinTickDrift = p.get("MinTickDrift", 3200);
64 
65  produces<std::vector<anab::CosmicTag>>();
66  produces<art::Assns<recob::Cluster, anab::CosmicTag>>();
67 }
68 
70 {
71  std::unique_ptr<std::vector<anab::CosmicTag>> cosmicTagClusterVector(
72  new std::vector<anab::CosmicTag>);
73  std::unique_ptr<art::Assns<recob::Cluster, anab::CosmicTag>> assnOutCosmicTagCluster(
75 
77  e.getByLabel(fClusterModuleLabel, Cluster_h);
78  std::vector<art::Ptr<recob::Cluster>> ClusterVec;
79  art::fill_ptr_vector(ClusterVec, Cluster_h);
80 
82  // LOOP OVER CLUSTERS
84 
85  for (unsigned int iCluster = 0; iCluster < Cluster_h->size(); iCluster++) {
86 
87  float cosmicScore = 0;
89 
90  art::Ptr<recob::Cluster> tCluster = ClusterVec.at(iCluster);
91  // art::Ptr<recob::Track> tTrack; // unused
92 
93  std::vector<float> endPt1;
94  std::vector<float> endPt2;
95 
96  // Doing some checks on the cluster to determine if it's a cosmic
97  bool failClusterTickCheck = false;
98 
99  auto const [t0, t1] = std::minmax(tCluster->StartTick(), tCluster->EndTick());
100  if (t0 + fTickLimit < fMinTickDrift) { // This is into the pre-spill window
101  failClusterTickCheck = true;
102  }
103  if (t1 - fTickLimit > fMaxTickDrift) { // This is into the post-spill window
104  failClusterTickCheck = true;
105  }
106 
107  if (failClusterTickCheck) {
108  cosmicScore = 1.;
110  }
111 
112  if (endPt1.size() < 1) {
113  for (int s = 0; s < 3; s++) {
114  endPt1.push_back(-999);
115  endPt2.push_back(-999);
116  }
117  }
118 
119  // Making stuff to save!
120  cosmicTagClusterVector->emplace_back(endPt1, endPt2, cosmicScore, tag_id);
121 
122  util::CreateAssn(e, *cosmicTagClusterVector, tCluster, *assnOutCosmicTagCluster);
123  }
124 
126  // END OF CLUSTER LOOP
128 
129  e.put(std::move(cosmicTagClusterVector));
130  e.put(std::move(assnOutCosmicTagCluster));
131 
132 } // end of produce
134 
code to link reconstructed objects back to the MC truth information
Utilities related to art service access.
TTree * t1
Definition: plottest35.C:26
void produce(art::Event &e) override
std::pair< float, float > minmax(const float a, const float b)
minmax
Unknown view.
Definition: geo_types.h:142
enum anab::cosmic_tag_id CosmicTagID_t
EDProducer(fhicl::ParameterSet const &pset)
Definition: EDProducer.cc:6
float EndTick() const
Returns the tick coordinate of the end of the cluster.
Definition: Cluster.h:331
PutHandle< PROD > put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
Definition: Event.h:77
#define DEFINE_ART_MODULE(klass)
Definition: ModuleMacros.h:65
CosmicClusterTagger(fhicl::ParameterSet const &p)
Provides recob::Track data product.
Declaration of cluster object.
bool CreateAssn(art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t index=UINT_MAX)
Creates a single one-to-one association.
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Utility object to perform functions of association.
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:306
Float_t e
Definition: plot.C:35
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Namespace collecting geometry-related classes utilities.
float StartTick() const
Returns the tick coordinate of the start of the cluster.
Definition: Cluster.h:287
art framework interface to geometry description