bug fix in input handling and relnotes v01-17-01

AIDASoft · May 17, 2018 · 16de447 · 16de447
1 parent 9a8bc58
commit 16de447
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 1 deletion.
diff --git a/DDG4/src/Geant4InputHandling.cpp b/DDG4/src/Geant4InputHandling.cpp
@@ -384,10 +384,12 @@ getRelevant(set<int>& visited,
     double proper_time_Precision = pow(10.,-DBL_DIG)*me*fmax(fabs(p->time),fabs(dp->time));
     bool isProperTimeZero = (proper_time <= proper_time_Precision);
 
+    const std::set<int> leptonPDGs{11,13,15,17};
     // -- remove original if ---
     bool rejectParticle = not p.definition() // completely unknown to geant4
       or (rejectPDGs.count(abs(p->pdgID)) != 0) // quarks, gluon, "strings", W, Z etc.
-      or (isProperTimeZero and p.definition()->GetPDGStable()); // initial state electrons, etc.
+      or (isProperTimeZero and p.definition()->GetPDGStable() ) // initial state electrons, etc.
+      or (isProperTimeZero and leptonPDGs.count(abs(p->pdgID)) != 0 ) ; // charged 'documentation' leptons, e.g. in lepton pairs w/ FSR
     if (not rejectParticle) {
       map<int,G4PrimaryParticle*>::iterator ip4 = prim.find(p->id);
       G4PrimaryParticle* p4 = (ip4 == prim.end()) ? 0 : (*ip4).second;

diff --git a/doc/ReleaseNotes.md b/doc/ReleaseNotes.md
@@ -1,3 +1,10 @@
+# v01-07-01
+
+* 2018-05-17 Frank Gaede
+   - exclude leptons with zero lifetime from Geant4
+   - fixes issue: https://github.com/AIDASoft/DD4hep/issues/387
+
+
 # v01-07
 
 * 2018-03-26 Javier Cervantes Villanueva ([PR#343](https://github.com/AIDASoft/DD4hep/pull/343))