Problem 58 - In muons/src/G4MuPairProduction.cc log() from 0 calculated
Summary: In muons/src/G4MuPairProduction.cc log() from 0 calculated
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/electromagnetic (show other problems)
Version: 1.0
Hardware: PC Other
: P2 normal
Assignee: Michel.Maire
URL:
Depends on:
Blocks:
 
Reported: 2000-02-09 00:34 CET by isupov
Modified: 2012-02-15 04:30 CET (History)
0 users

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description isupov 2000-02-09 00:34:10 CET 
In the muons/src/G4MuPairProduction.cc, G4MuPairProduction::PostStepDoIt()
we have (lines 649-650):

   xc = log(CutInPairEnergy/MinPairEnergy)/log(MaxPairEnergy/MinPairEnergy) ;
   yc = log(xc) ;

which WILL calculate log(0), if CutInPairEnergy == MinPairEnergy  .
But in line 601 we see:

   if (CutInPairEnergy < MinPairEnergy) CutInPairEnergy = MinPairEnergy ;

:-)
But what we must do in such case ? May be, simply

   if (CutInPairEnergy < MinPairEnergy)
     return G4VContinuousDiscreteProcess::PostStepDoIt(trackData,stepData);

?
Comment 1 Michel.Maire 2000-03-02 10:13:59 CET 
fixed in tag em-V01-00-00