Add local trigger time to ntuples

rat-pac · Jul 16, 2024 · b77cbfd · b77cbfd
1 parent 9778277
commit b77cbfd
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 0 deletions.
diff --git a/src/io/include/RAT/OutNtupleProc.hh b/src/io/include/RAT/OutNtupleProc.hh
@@ -127,6 +127,7 @@ class OutNtupleProc : public Processor {
   std::vector<double> hitPMTDigitizedTime;
   std::vector<double> hitPMTDigitizedCharge;
   std::vector<int> hitPMTNCrossings;
+  std::vector<double> hitPMTDigitizedLocalTriggerTime;
   // Tracking
   std::map<std::string, int> processCodeMap;
   std::vector<int> processCodeIndex;

diff --git a/src/io/src/OutNtupleProc.cc b/src/io/src/OutNtupleProc.cc
@@ -128,6 +128,7 @@ bool OutNtupleProc::OpenFile(std::string filename) {
     outputTree->Branch("hitPMTDigitizedTime", &hitPMTDigitizedTime);
     outputTree->Branch("hitPMTDigitizedCharge", &hitPMTDigitizedCharge);
     outputTree->Branch("hitPMTNCrossings", &hitPMTNCrossings);
+    outputTree->Branch("hitPMTDigitizedLocalTriggerTime", &hitPMTDigitizedLocalTriggerTime);
   }
   if (options.mchits) {
     // Save full MC PMT hit information
@@ -407,6 +408,7 @@ Processor::Result OutNtupleProc::DSEvent(DS::Root *ds) {
       hitPMTDigitizedTime.clear();
       hitPMTDigitizedCharge.clear();
       hitPMTNCrossings.clear();
+      hitPMTDigitizedLocalTriggerTime.clear();
 
       for (int pmtc = 0; pmtc < ev->GetPMTCount(); pmtc++) {
         RAT::DS::PMT *pmt = ev->GetPMT(pmtc);
@@ -419,6 +421,7 @@ Processor::Result OutNtupleProc::DSEvent(DS::Root *ds) {
         hitPMTDigitizedTime.push_back(digitpmt->GetDigitizedTime());
         hitPMTDigitizedCharge.push_back(digitpmt->GetDigitizedCharge());
         hitPMTNCrossings.push_back(digitpmt->GetNCrossings());
+        hitPMTDigitizedLocalTriggerTime.push_back(digitpmt->GetLocalTriggerTime());
       }
     }
     this->FillEvent(ds, ev);
@@ -434,6 +437,7 @@ Processor::Result OutNtupleProc::DSEvent(DS::Root *ds) {
       hitPMTDigitizedTime.clear();
       hitPMTDigitizedCharge.clear();
       hitPMTNCrossings.clear();
+      hitPMTDigitizedLocalTriggerTime.clear();
     }
     this->FillNoTriggerEvent(ds);
     outputTree->Fill();