Add histogram of number of triggered channels in muon candidate events

sfegan · Oct 8, 2024 · 249a947 · 249a947
1 parent 5c436c6
commit 249a947
Show file tree

Hide file tree

Showing 3 changed files with 20 additions and 1 deletion.
diff --git a/include/diagnostics/simple_charge_stats.hpp b/include/diagnostics/simple_charge_stats.hpp
@@ -160,14 +160,16 @@ class SimpleChargeStatsParallelEventVisitor:
       num_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0)),
       num_contiguous_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0)),
       phys_trig_num_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0)),
-      phys_trig_num_contiguous_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0))
+      phys_trig_num_contiguous_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0)),
+      muon_candidate_num_channel_triggered_hist(new calin::math::histogram::Histogram1D(1.0))
     { /* nothing to see here */ }
 
     ~CameraHists() {
       delete num_channel_triggered_hist;
       delete num_contiguous_channel_triggered_hist;
       delete phys_trig_num_channel_triggered_hist;
       delete phys_trig_num_contiguous_channel_triggered_hist;
+      delete muon_candidate_num_channel_triggered_hist;
       delete high_gain;
       delete low_gain;
     }
@@ -179,6 +181,7 @@ class SimpleChargeStatsParallelEventVisitor:
     calin::math::histogram::Histogram1D* num_contiguous_channel_triggered_hist = nullptr;
     calin::math::histogram::Histogram1D* phys_trig_num_channel_triggered_hist = nullptr;
     calin::math::histogram::Histogram1D* phys_trig_num_contiguous_channel_triggered_hist = nullptr;
+    calin::math::histogram::Histogram1D* muon_candidate_num_channel_triggered_hist = nullptr;
   };
 
   void integrate_one_gain_partials(

diff --git a/proto/diagnostics/simple_charge_stats.proto b/proto/diagnostics/simple_charge_stats.proto
@@ -237,6 +237,9 @@ message SimpleChargeStats {
   calin.ix.math.histogram.Histogram1DData phy_trigger_num_contiguous_channel_triggered_hist = 35 [
     (CFO).desc = "Histogram of number of channels triggered in physics events.",
     (CFO).units = "channels (bin-units), events (weight-units)" ];
+  calin.ix.math.histogram.Histogram1DData muon_candidate_num_channel_triggered_hist = 36 [
+    (CFO).desc = "Histogram of number of channels triggered in muon-candidate events.",
+    (CFO).units = "channels (bin-units), events (weight-units)" ];
 };
 
 message PartialOneGainChannelSimpleChargeStats {
@@ -432,6 +435,10 @@ message PartialCameraSimpleChargeStats {
   calin.ix.math.histogram.Histogram1DData phys_trig_num_contiguous_channel_triggered_hist = 25 [
     (CFO).desc = "Histogram of number of channels triggered in physics-trigger events.",
     (CFO).units = "channels (bin-units), events (weight-units)" ];
+  calin.ix.math.histogram.Histogram1DData muon_candidate_num_channel_triggered_hist = 26 [
+    (CFO).desc = "Histogram of number of channels triggered in muon-candidate events.",
+    (CFO).units = "channels (bin-units), events (weight-units)" ];
+
 };
 
 message PartialSimpleChargeStats {

diff --git a/src/diagnostics/simple_charge_stats.cpp b/src/diagnostics/simple_charge_stats.cpp
@@ -515,6 +515,10 @@ bool SimpleChargeStatsParallelEventVisitor::leave_telescope_run(
   partials_.mutable_camera()->mutable_phys_trig_num_contiguous_channel_triggered_hist()->IntegrateFrom(*hp);
   delete hp;
 
+  hp = camera_hists_->muon_candidate_num_channel_triggered_hist->dump_as_proto();
+  partials_.mutable_camera()->mutable_muon_candidate_num_channel_triggered_hist()->IntegrateFrom(*hp);
+  delete hp;
+
   if(parent_)return true;
 
   for(int ichan = 0; ichan<partials_.channel_size(); ichan++) {
@@ -541,6 +545,8 @@ bool SimpleChargeStatsParallelEventVisitor::leave_telescope_run(
     partials_.camera().phys_trig_num_channel_triggered_hist());
   results_.mutable_phy_trigger_num_contiguous_channel_triggered_hist()->IntegrateFrom(
     partials_.camera().phys_trig_num_contiguous_channel_triggered_hist());
+  results_.mutable_muon_candidate_num_channel_triggered_hist()->IntegrateFrom(
+    partials_.camera().muon_candidate_num_channel_triggered_hist());
 
   partials_.Clear();
   return true;
@@ -775,6 +781,9 @@ bool SimpleChargeStatsParallelEventVisitor::visit_telescope_event(uint64_t seq_i
         }
       }
     }
+    if(event->is_muon_candidate()) {
+      camera_hists_->muon_candidate_num_channel_triggered_hist->insert(event->trigger_map().hit_channel_id_size());
+    }
   }
   return true;
 }