Merge pull request #140 from tannerbk/mcindex

add mcid to ntuple
rat-pac · Jul 11, 2024 · a99058e · a99058e
2 parents 07d87fa + da70cd5
commit a99058e
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 0 deletions.
diff --git a/src/io/include/RAT/OutNtupleProc.hh b/src/io/include/RAT/OutNtupleProc.hh
@@ -107,6 +107,7 @@ class OutNtupleProc : public Processor {
   std::vector<double> mcpecharge;
   // MCParticles
   int mcpcount;
+  int mcid;
   std::vector<Int_t> pdgcodes;
   std::vector<double> mcKEnergies;
   std::vector<double> mcPosx;

diff --git a/src/io/src/OutNtupleProc.cc b/src/io/src/OutNtupleProc.cc
@@ -96,6 +96,7 @@ bool OutNtupleProc::OpenFile(std::string filename) {
   outputTree->Branch("nhits", &nhits);
   outputTree->Branch("triggerTime", &triggerTime);
   // MC Information
+  outputTree->Branch("mcid", &mcid);
   outputTree->Branch("mcparticlecount", &mcpcount);
   outputTree->Branch("mcpecount", &mcpecount);
   outputTree->Branch("mcnhits", &mcnhits);
@@ -187,6 +188,7 @@ Processor::Result OutNtupleProc::DSEvent(DS::Root *ds) {
   mcTime.clear();
 
   DS::MC *mc = ds->GetMC();
+  mcid = mc->GetID();
   mcpcount = mc->GetMCParticleCount();
   for (int pid = 0; pid < mcpcount; pid++) {
     DS::MCParticle *particle = mc->GetMCParticle(pid);