Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // ------------------------------------------- 28 // 29 // GEANT4 Class file 30 // 31 // 32 // File name: G4LowEWentzelVIModel 33 // 34 // Author: V.Ivanchenko 35 // 36 // Creation date: 11.02.2014 from G4WentzelVIM 37 // 38 // Modifications: 39 // 40 // Class Description: 41 // 42 43 // ------------------------------------------- 44 // 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 //....oooOO0OOooo........oooOO0OOooo........oo 48 49 #include "G4LowEWentzelVIModel.hh" 50 #include "G4PhysicalConstants.hh" 51 #include "G4SystemOfUnits.hh" 52 53 //....oooOO0OOooo........oooOO0OOooo........oo 54 55 G4LowEWentzelVIModel::G4LowEWentzelVIModel() : 56 G4WentzelVIModel(true, "LowEnWentzelVI") 57 { 58 SetSingleScatteringFactor(0.5); 59 } 60 61 //....oooOO0OOooo........oooOO0OOooo........oo 62 63 G4double G4LowEWentzelVIModel::ComputeTruePath 64 const G4Track& tr 65 G4double& currentMinimalStep) 66 { 67 G4double tlimit = currentMinimalStep; 68 const G4DynamicParticle* dp = track.GetDynam 69 G4StepPoint* sp = track.GetStep()->GetPreSte 70 G4StepStatus stepStatus = sp->GetStepStatus( 71 singleScatteringMode = false; 72 73 // initialisation for each step, lambda may 74 preKinEnergy = dp->GetKineticEnergy(); 75 DefineMaterial(track.GetMaterialCutsCouple() 76 lambdaeff = GetTransportMeanFreePath(particl 77 currentRange = GetRange(particle,preKinEnerg 78 cosTetMaxNuc = wokvi->SetupKinematic(preKinE 79 80 // extra check for abnormal situation 81 // this check needed to run MSC with eIoni a 82 tlimit = std::min(tlimit, currentRange); 83 84 // stop here if small range particle 85 if(tlimit < tlimitminfix) { 86 return ConvertTrueToGeom(tlimit, currentMi 87 } 88 89 // pre step 90 G4double presafety = sp->GetSafety(); 91 // far from geometry boundary 92 if(currentRange < presafety) { 93 return ConvertTrueToGeom(tlimit, currentMi 94 } 95 96 // compute presafety again if presafety <= 0 97 // i.e. when it is needed for optimization p 98 if(stepStatus != fGeomBoundary && presafety 99 presafety = ComputeSafety(sp->GetPosition( 100 if(currentRange < presafety) { 101 return ConvertTrueToGeom(tlimit, current 102 } 103 } 104 /* 105 G4cout << "e(MeV)= " << preKinEnergy/MeV 106 << " " << particle->GetParticleName() 107 << " CurLimit(mm)= " << tlimit/mm <<" safet 108 << " R(mm)= " <<currentRange/mm 109 << " L0(mm^-1)= " << lambdaeff*mm 110 <<G4endl; 111 */ 112 // natural limit for high energy 113 G4double rlimit = std::max(facrange*currentR 114 115 // low-energy e- 116 rlimit = std::max(rlimit, facsafety*presafet 117 118 tlimit = std::min(tlimit, rlimit); 119 tlimit = std::max(tlimit, tlimitminfix); 120 121 // step limit in infinite media 122 tlimit = std::min(tlimit, 50*currentMaterial 123 124 //compute geomlimit and force few steps with 125 if (steppingAlgorithm == fUseDistanceToBound 126 && stepStatus == fGeomBoundary) { 127 128 G4double geomlimit = ComputeGeomLimit(trac 129 tlimit = std::min(tlimit, geomlimit/facgeo 130 } 131 /* 132 G4cout << particle->GetParticleName() << " E 133 << " L0= " << lambdaeff << " R= " << curren 134 << " tlimit= " << tlimit 135 << " currentMinimalStep= " << currentMini 136 */ 137 return ConvertTrueToGeom(tlimit, currentMini 138 } 139 140 //....oooOO0OOooo........oooOO0OOooo........oo 141