add geant pdg explanation

TTU-HEP · Dec 30, 2024 · cd8349b · cd8349b
1 parent e27bde1
commit cd8349b
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/sim/src/B4bSteppingAction.cc b/sim/src/B4bSteppingAction.cc
@@ -213,6 +213,14 @@ void B4bSteppingAction::UserSteppingAction(const G4Step *step)
   aHit.y = posA.y() / cm;
   aHit.z = posA.z() / cm;
   aHit.pid = pdgcode;
+  // if (fabs(pdgcode) > 1e9)
+  //{
+  //   // https://indico.ph.tum.de/event/3955/sessions/752/attachments/2741/3099/Day3_Physics.pdf
+  //   // larger than 1e9 is nuclei, ion, or excited state of atoms and nuclei
+  //   // e.g., 1000020040 is alpha, 1000010020 is deuteron, 100ZZZAAAI* is
+  //   // *ZZZ=proton number, AAA=nucleon number, I=excitation level
+  //   std::cout << "pdgcode " << pdgcode << "  mass " << mass << "  kinEnergy " << kinEnergy << std::endl;
+  // }
   aHit.calotype = caloType;
   aHit.trackid = track->GetTrackID();
   aHit.globaltime = track->GetGlobalTime() / ns;