Skip to content
Snippets Groups Projects
Commit 039aaadd authored by Andre Sailer's avatar Andre Sailer
Browse files

Geant4InputHandling: protect propertime against negative masses

parent f6f3f7fc
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
DDG4/src/Geant4InputHandling.cpp
+ 2
− 2
View file @ 039aaadd
...@@ -385,8 +385,8 @@ getRelevant(set<int>& visited, ...@@ -385,8 +385,8 @@ getRelevant(set<int>& visited,
double me = en > std::numeric_limits<double>::epsilon() ? p->mass / en : 0.0; double me = en > std::numeric_limits<double>::epsilon() ? p->mass / en : 0.0;
// fix by S.Morozov for real != 0 // fix by S.Morozov for real != 0
double proper_time = fabs(dp->time-p->time) * me; double proper_time = fabs(dp->time-p->time) * me;
double proper_time_Precision = pow(10.,-DBL_DIG)*me*fmax(fabs(p->time),fabs(dp->time)); double proper_time_Precision = pow(10.,-DBL_DIG)*fabs(me)*fmax(fabs(p->time),fabs(dp->time));
bool isProperTimeZero = (proper_time <= proper_time_Precision); bool isProperTimeZero = (fabs(proper_time) <= fabs(proper_time_Precision));
// -- remove original if --- // -- remove original if ---
bool rejectParticle = not p.definition() // completely unknown to geant4 bool rejectParticle = not p.definition() // completely unknown to geant4
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment