Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // >> 26 // $Id: G4WentzelVIModel.hh 89893 2015-05-04 07:29:17Z gcosmo $ 26 // 27 // 27 // ------------------------------------------- 28 // ------------------------------------------------------------------- 28 // 29 // 29 // 30 // 30 // GEANT4 Class header file 31 // GEANT4 Class header file 31 // 32 // 32 // 33 // 33 // File name: G4WentzelVIModel 34 // File name: G4WentzelVIModel 34 // 35 // 35 // Author: V.Ivanchenko 36 // Author: V.Ivanchenko 36 // 37 // 37 // Creation date: 09.04.2008 from G4MuMscModel 38 // Creation date: 09.04.2008 from G4MuMscModel 38 // 39 // 39 // Modifications: 40 // Modifications: 40 // 27-05-2010 V.Ivanchenko added G4WentzelOKan 41 // 27-05-2010 V.Ivanchenko added G4WentzelOKandVIxSection class to 41 // compute cross sections and sam 42 // compute cross sections and sample scattering angle 42 // 43 // 43 // Class Description: 44 // Class Description: 44 // 45 // 45 // Implementation of the model of multiple sca 46 // Implementation of the model of multiple scattering based on 46 // G.Wentzel, Z. Phys. 40 (1927) 590. 47 // G.Wentzel, Z. Phys. 40 (1927) 590. 47 // H.W.Lewis, Phys Rev 78 (1950) 526. 48 // H.W.Lewis, Phys Rev 78 (1950) 526. 48 // J.M. Fernandez-Varea et al., NIM B73 (1993) 49 // J.M. Fernandez-Varea et al., NIM B73 (1993) 447. 49 // L.Urban, CERN-OPEN-2006-077. 50 // L.Urban, CERN-OPEN-2006-077. 50 51 51 // ------------------------------------------- 52 // ------------------------------------------------------------------- 52 // 53 // 53 54 54 #ifndef G4WentzelVIModel_h 55 #ifndef G4WentzelVIModel_h 55 #define G4WentzelVIModel_h 1 56 #define G4WentzelVIModel_h 1 56 57 57 //....oooOO0OOooo........oooOO0OOooo........oo 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 58 59 59 #include "G4VMscModel.hh" 60 #include "G4VMscModel.hh" 60 #include "G4MaterialCutsCouple.hh" 61 #include "G4MaterialCutsCouple.hh" 61 #include "G4WentzelOKandVIxSection.hh" 62 #include "G4WentzelOKandVIxSection.hh" 62 63 63 //....oooOO0OOooo........oooOO0OOooo........oo 64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 65 65 class G4WentzelVIModel : public G4VMscModel 66 class G4WentzelVIModel : public G4VMscModel 66 { 67 { 67 68 68 public: 69 public: 69 70 70 explicit G4WentzelVIModel(G4bool comb=true, << 71 G4WentzelVIModel(G4bool comb = true, const G4String& nam = "WentzelVIUni"); 71 72 72 ~G4WentzelVIModel() override; << 73 virtual ~G4WentzelVIModel(); 73 74 74 void Initialise(const G4ParticleDefinition*, << 75 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&); 75 76 76 void InitialiseLocal(const G4ParticleDefinit << 77 virtual void InitialiseLocal(const G4ParticleDefinition*, 77 G4VEmModel* masterModel) override; << 78 G4VEmModel* masterModel); 78 79 79 void StartTracking(G4Track*) override; << 80 void StartTracking(G4Track*); 80 81 81 G4double ComputeCrossSectionPerAtom(const G4 << 82 virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, 82 G4double KineticEnergy, << 83 G4double KineticEnergy, 83 G4double AtomicNumber, << 84 G4double AtomicNumber, 84 G4double AtomicWeight=0., << 85 G4double AtomicWeight=0., 85 G4double cut = DBL_MAX, << 86 G4double cut = DBL_MAX, 86 G4double emax= DBL_MAX) override << 87 G4double emax= DBL_MAX); 87 88 88 G4ThreeVector& SampleScattering(const G4Thre << 89 virtual G4ThreeVector& SampleScattering(const G4ThreeVector&, 89 G4double safety) override; << 90 G4double safety); 90 91 91 G4double << 92 virtual G4double ComputeTruePathLengthLimit(const G4Track& track, 92 ComputeTruePathLengthLimit(const G4Track& tr << 93 G4double& currentMinimalStep); 93 G4double& currentMinimalStep) overr << 94 94 95 G4double ComputeGeomPathLength(G4double true << 95 virtual G4double ComputeGeomPathLength(G4double truePathLength); 96 96 97 G4double ComputeTrueStepLength(G4double geom << 97 virtual G4double ComputeTrueStepLength(G4double geomStepLength); 98 98 99 // defines low energy limit on energy transf 99 // defines low energy limit on energy transfer to atomic electron 100 void SetFixedCut(G4double); << 100 inline void SetFixedCut(G4double); 101 101 102 // low energy limit on energy transfer to at 102 // low energy limit on energy transfer to atomic electron 103 G4double GetFixedCut() const; << 103 inline G4double GetFixedCut() const; 104 104 105 // access to cross section class 105 // access to cross section class 106 void SetWVICrossSection(G4WentzelOKandVIxSec << 106 inline G4WentzelOKandVIxSection* GetWVICrossSection(); 107 107 108 G4WentzelOKandVIxSection* GetWVICrossSection << 108 inline void SetUseSecondMoment(G4bool); 109 109 110 void SetUseSecondMoment(G4bool); << 110 inline G4bool UseSecondMoment() const; 111 111 112 G4bool UseSecondMoment() const; << 112 inline G4PhysicsTable* GetSecondMomentTable(); 113 113 114 G4PhysicsTable* GetSecondMomentTable(); << 114 inline G4double SecondMoment(const G4ParticleDefinition*, 115 << 115 const G4MaterialCutsCouple*, 116 G4double SecondMoment(const G4ParticleDefini << 116 G4double kineticEnergy); 117 const G4MaterialCutsCouple*, << 118 G4double kineticEnergy); << 119 117 120 void SetSingleScatteringFactor(G4double); 118 void SetSingleScatteringFactor(G4double); 121 119 122 void DefineMaterial(const G4MaterialCutsCoup << 123 << 124 G4WentzelVIModel & operator=(const G4Wentzel << 125 G4WentzelVIModel(const G4WentzelVIModel&) = << 126 << 127 protected: 120 protected: 128 121 129 G4double ComputeTransportXSectionPerVolume(G << 122 inline void DefineMaterial(const G4MaterialCutsCouple*); 130 << 131 inline void SetupParticle(const G4ParticleDe << 132 123 133 private: 124 private: 134 125 >> 126 G4double ComputeTransportXSectionPerVolume(G4double cosTheta); >> 127 135 G4double ComputeSecondMoment(const G4Particl 128 G4double ComputeSecondMoment(const G4ParticleDefinition*, 136 G4double kineticEnergy); 129 G4double kineticEnergy); 137 130 >> 131 inline void SetupParticle(const G4ParticleDefinition*); >> 132 >> 133 // hide assignment operator >> 134 G4WentzelVIModel & operator=(const G4WentzelVIModel &right); >> 135 G4WentzelVIModel(const G4WentzelVIModel&); >> 136 138 protected: 137 protected: 139 138 140 G4WentzelOKandVIxSection* wokvi; 139 G4WentzelOKandVIxSection* wokvi; 141 const G4MaterialCutsCouple* currentCouple = << 142 const G4Material* currentMaterial = nullptr; << 143 << 144 const G4ParticleDefinition* particle = nullp << 145 G4ParticleChangeForMSC* fParticleChange = nu << 146 const G4DataVector* currentCuts = nullptr; << 147 G4PhysicsTable* fSecondMoments = nullptr; << 148 140 149 G4double lowEnergyLimit; << 150 G4double tlimitminfix; 141 G4double tlimitminfix; 151 G4double ssFactor = 1.05; << 142 G4double ssFactor; 152 G4double invssFactor = 1.0; << 143 G4double invssFactor; 153 144 154 // cache kinematics 145 // cache kinematics 155 G4double preKinEnergy = 0.0; << 146 G4double preKinEnergy; 156 G4double tPathLength = 0.0; << 147 G4double tPathLength; 157 G4double zPathLength = 0.0; << 148 G4double zPathLength; 158 G4double lambdaeff = 0.0; << 149 G4double lambdaeff; 159 G4double currentRange = 0.0; << 150 G4double currentRange; 160 G4double cosTetMaxNuc = 0.0; << 151 G4double cosTetMaxNuc; >> 152 >> 153 // cache material >> 154 G4int currentMaterialIndex; >> 155 const G4MaterialCutsCouple* currentCouple; >> 156 const G4Material* currentMaterial; >> 157 >> 158 const G4ParticleDefinition* particle; >> 159 >> 160 // flags >> 161 G4bool singleScatteringMode; >> 162 >> 163 private: >> 164 >> 165 G4ParticleChangeForMSC* fParticleChange; >> 166 const G4DataVector* currentCuts; 161 167 162 G4double fixedCut = -1.0; << 168 G4double invsqrt12; >> 169 G4double fixedCut; 163 170 164 // cache kinematics 171 // cache kinematics 165 G4double effKinEnergy = 0.0; << 172 G4double effKinEnergy; 166 173 167 // single scattering parameters 174 // single scattering parameters 168 G4double cosThetaMin = 1.0; << 175 G4double cosThetaMin; 169 G4double cosThetaMax = -1.0; << 176 G4double cosThetaMax; 170 G4double xtsec = 0.0; << 171 177 172 G4int currentMaterialIndex = 0; << 178 G4PhysicsTable* fSecondMoments; 173 size_t idx2 = 0; << 179 size_t idx2; 174 180 175 // data for single scattering mode 181 // data for single scattering mode 176 G4int nelments = 0; << 182 G4double xtsec; >> 183 std::vector<G4double> xsecn; >> 184 std::vector<G4double> prob; >> 185 G4int nelments; >> 186 >> 187 G4double numlimit; >> 188 >> 189 // projectile >> 190 G4double lowEnergyLimit; 177 191 178 // flags 192 // flags 179 G4bool singleScatteringMode; << 180 G4bool isCombined; 193 G4bool isCombined; 181 G4bool useSecondMoment; 194 G4bool useSecondMoment; 182 << 183 std::vector<G4double> xsecn; << 184 std::vector<G4double> prob; << 185 }; 195 }; 186 196 187 //....oooOO0OOooo........oooOO0OOooo........oo 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 188 //....oooOO0OOooo........oooOO0OOooo........oo 198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 189 199 >> 200 inline >> 201 void G4WentzelVIModel::DefineMaterial(const G4MaterialCutsCouple* cup) >> 202 { >> 203 if(cup != currentCouple) { >> 204 currentCouple = cup; >> 205 SetCurrentCouple(cup); >> 206 currentMaterial = cup->GetMaterial(); >> 207 currentMaterialIndex = currentCouple->GetIndex(); >> 208 } >> 209 } >> 210 >> 211 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 212 190 inline void G4WentzelVIModel::SetupParticle(co 213 inline void G4WentzelVIModel::SetupParticle(const G4ParticleDefinition* p) 191 { 214 { 192 // Initialise mass and charge 215 // Initialise mass and charge 193 if(p != particle) { 216 if(p != particle) { 194 particle = p; 217 particle = p; 195 wokvi->SetupParticle(p); 218 wokvi->SetupParticle(p); 196 } 219 } 197 } 220 } 198 221 199 //....oooOO0OOooo........oooOO0OOooo........oo 222 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 200 223 201 inline void G4WentzelVIModel::SetFixedCut(G4do 224 inline void G4WentzelVIModel::SetFixedCut(G4double val) 202 { 225 { 203 fixedCut = val; 226 fixedCut = val; 204 } 227 } 205 228 206 //....oooOO0OOooo........oooOO0OOooo........oo 229 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 207 230 208 inline G4double G4WentzelVIModel::GetFixedCut( 231 inline G4double G4WentzelVIModel::GetFixedCut() const 209 { 232 { 210 return fixedCut; 233 return fixedCut; 211 } 234 } 212 235 213 //....oooOO0OOooo........oooOO0OOooo........oo 236 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 214 237 215 inline void G4WentzelVIModel::SetWVICrossSecti << 216 { << 217 if(ptr != wokvi) { << 218 delete wokvi; << 219 wokvi = ptr; << 220 } << 221 } << 222 << 223 //....oooOO0OOooo........oooOO0OOooo........oo << 224 << 225 inline G4WentzelOKandVIxSection* G4WentzelVIMo 238 inline G4WentzelOKandVIxSection* G4WentzelVIModel::GetWVICrossSection() 226 { 239 { 227 return wokvi; 240 return wokvi; 228 } 241 } 229 242 230 //....oooOO0OOooo........oooOO0OOooo........oo 243 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 231 244 232 inline void G4WentzelVIModel::SetUseSecondMome 245 inline void G4WentzelVIModel::SetUseSecondMoment(G4bool val) 233 { 246 { 234 useSecondMoment = val; 247 useSecondMoment = val; 235 } 248 } 236 249 237 //....oooOO0OOooo........oooOO0OOooo........oo 250 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 238 251 239 inline G4bool G4WentzelVIModel::UseSecondMomen 252 inline G4bool G4WentzelVIModel::UseSecondMoment() const 240 { 253 { 241 return useSecondMoment; 254 return useSecondMoment; 242 } 255 } 243 256 244 //....oooOO0OOooo........oooOO0OOooo........oo 257 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 245 258 246 inline G4PhysicsTable* G4WentzelVIModel::GetSe 259 inline G4PhysicsTable* G4WentzelVIModel::GetSecondMomentTable() 247 { 260 { 248 return fSecondMoments; 261 return fSecondMoments; 249 } 262 } 250 263 251 //....oooOO0OOooo........oooOO0OOooo........oo 264 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 252 265 253 inline G4double 266 inline G4double 254 G4WentzelVIModel::SecondMoment(const G4Particl 267 G4WentzelVIModel::SecondMoment(const G4ParticleDefinition* part, 255 const G4MaterialCutsCouple* coupl 268 const G4MaterialCutsCouple* couple, 256 G4double ekin) 269 G4double ekin) 257 { 270 { 258 G4double x = 0.0; 271 G4double x = 0.0; 259 if(useSecondMoment) { 272 if(useSecondMoment) { 260 DefineMaterial(couple); 273 DefineMaterial(couple); 261 x = (fSecondMoments) ? << 274 if(fSecondMoments) { 262 (*fSecondMoments)[(*theDensityIdx)[curre << 275 x = (*fSecondMoments)[(*theDensityIdx)[currentMaterialIndex]] 263 *(*theDensityFactor)[currentMaterialInde << 276 ->Value(ekin, idx2) 264 : ComputeSecondMoment(part, ekin); << 277 *(*theDensityFactor)[currentMaterialIndex]/(ekin*ekin); >> 278 } else { >> 279 x = ComputeSecondMoment(part, ekin); >> 280 } 265 } 281 } 266 return x; 282 return x; 267 } 283 } 268 284 269 //....oooOO0OOooo........oooOO0OOooo........oo 285 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 270 286 271 #endif 287 #endif 272 288 273 289