- Timestamp:
- Apr 6, 2009, 12:21:12 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/standard/include/G4MultipleScattering71.hh
r819 r961 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4MultipleScattering71.hh,v 1. 5 2007/05/22 17:34:36vnivanch Exp $27 // GEANT4 tag $Name: $26 // $Id: G4MultipleScattering71.hh,v 1.6 2008/07/16 11:27:41 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-02-ref-02 $ 28 28 // 29 29 // … … 58 58 //------------------------------------------------------------------------------ 59 59 // 60 // $Id: G4MultipleScattering71.hh,v 1. 5 2007/05/22 17:34:36vnivanch Exp $61 // GEANT4 tag $Name: $60 // $Id: G4MultipleScattering71.hh,v 1.6 2008/07/16 11:27:41 vnivanch Exp $ 61 // GEANT4 tag $Name: geant4-09-02-ref-02 $ 62 62 63 63 // class description … … 180 180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 181 181 182 inline G4double G4MultipleScattering71::GetContinuousStepLimit(183 const G4Track& track,184 G4double,185 G4double currentMinimalStep,186 G4double&)187 {188 DefineMaterial(track.GetMaterialCutsCouple());189 const G4MaterialCutsCouple* couple = CurrentMaterialCutsCouple();190 G4double e = track.GetKineticEnergy();191 model = dynamic_cast<G4MscModel71*>(SelectModel(e));192 const G4ParticleDefinition* p = track.GetDefinition();193 G4double lambda0 = GetLambda(p, e);194 range = G4LossTableManager::Instance()->GetRangeFromRestricteDEDX(p,e,couple);195 if(range < currentMinimalStep) currentMinimalStep = range;196 truePathLength = TruePathLengthLimit(track,lambda0,currentMinimalStep);197 // G4cout << "StepLimit: tpl= " << truePathLength << " lambda0= "198 // << lambda0 << " range= " << currentRange199 // << " currentMinStep= " << currentMinimalStep << G4endl;200 if (truePathLength < currentMinimalStep) valueGPILSelectionMSC = CandidateForSelection;201 geomPathLength = model->GeomPathLength(LambdaTable(),couple,202 p,e,lambda0,range,truePathLength);203 if(geomPathLength > lambda0) geomPathLength = lambda0;204 return geomPathLength;205 }206 207 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....208 209 182 inline G4VParticleChange* G4MultipleScattering71::AlongStepDoIt( 210 183 const G4Track&,
Note: See TracChangeset
for help on using the changeset viewer.