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
G4double GetStepLength() const
const G4DynamicParticle * GetDynamicParticle() const
G4VParticleChange * AlongStepDoIt(const G4Track &aTrack, const G4Step &aStep)
const G4String & GetName() const
G4EnergyLossForExtrapolator * theELossForExtrapolator
G4ErrorEnergyLoss(const G4String &processName="G4ErrorEnergyLoss", G4ProcessType type=fElectromagnetic)
G4ParticleDefinition * GetDefinition() const
virtual ~G4ErrorEnergyLoss()
const G4String & GetParticleName() const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
G4ErrorMode GetMode() const
void InstantiateEforExtrapolator()
const G4String & GetProcessName() const
G4Material * GetMaterial() const
virtual void Initialize(const G4Track &)
void SetNumberOfSecondaries(G4int totSecondaries)
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange aParticleChange
static G4ErrorPropagatorData * GetErrorPropagatorData()
G4double GetContinuousStepLimit(const G4Track &aTrack, G4double, G4double currentMinimumStep, G4double &)