69 G4double kinEnergyEnd = kinEnergyStart;
76 G4double kinEnergyHalfStep = kinEnergyStart - (kinEnergyStart-kinEnergyEnd)/2.;
80 G4cout <<
" G4ErrorEnergyLoss FWD end " << kinEnergyEnd
81 <<
" halfstep " << kinEnergyHalfStep <<
G4endl;
89 kinEnergyEnd = kinEnergyStart - (kinEnergyHalfStep - kinEnergyEnd );
95 G4double kinEnergyHalfStep = kinEnergyStart - (kinEnergyStart-kinEnergyEnd)/2.;
98 G4cout <<
" G4ErrorEnergyLoss BCKD end " << kinEnergyEnd
99 <<
" halfstep " << kinEnergyHalfStep <<
G4endl;
107 kinEnergyEnd = kinEnergyStart - (kinEnergyHalfStep - kinEnergyEnd );
110 G4double edepo = kinEnergyEnd - kinEnergyStart;
114 G4cout <<
"AlongStepDoIt Estart= " << kinEnergyStart <<
" Eend " << kinEnergyEnd
115 <<
" Ediff " << kinEnergyStart-kinEnergyEnd <<
" step= " << step_length
116 <<
" mate= " << aMaterial->
GetName()
139 G4double kinEnergyLoss = kinEnergyStart;
144 kinEnergyLoss = - kinEnergyStart +
146 aMaterial, aParticleDef );
148 kinEnergyLoss = kinEnergyStart -
150 aMaterial, aParticleDef );
154 G4cout <<
" G4ErrorEnergyLoss: currentMinimumStep " <<currentMinimumStep
155 <<
" kinEnergyLoss " << kinEnergyLoss
156 <<
" kinEnergyStart " << kinEnergyStart <<
G4endl;
159 Step =
theStepLimit / (kinEnergyLoss / kinEnergyStart) * currentMinimumStep;
162 G4cout <<
" G4ErrorEnergyLoss: limiting Step " << Step
163 <<
" energy loss fraction " << kinEnergyLoss / kinEnergyStart