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: G4VEmModel.hh,v 1.42 2006/06/29 19:54:43 gunter Exp $ >> 27 // GEANT4 tag $Name: geant4-08-01 $ >> 28 // 26 // ------------------------------------------- 29 // ------------------------------------------------------------------- 27 // 30 // 28 // GEANT4 Class header file 31 // GEANT4 Class header file 29 // 32 // 30 // 33 // 31 // File name: G4VEmModel 34 // File name: G4VEmModel 32 // 35 // 33 // Author: Vladimir Ivanchenko 36 // Author: Vladimir Ivanchenko 34 // 37 // 35 // Creation date: 03.01.2002 38 // Creation date: 03.01.2002 36 // 39 // 37 // Modifications: 40 // Modifications: 38 // 41 // 39 // 23-12-02 V.Ivanchenko change interface befo 42 // 23-12-02 V.Ivanchenko change interface before move to cut per region 40 // 24-01-03 Cut per region (V.Ivanchenko) 43 // 24-01-03 Cut per region (V.Ivanchenko) 41 // 13-02-03 Add name (V.Ivanchenko) 44 // 13-02-03 Add name (V.Ivanchenko) 42 // 25-02-03 Add sample theta and displacement 45 // 25-02-03 Add sample theta and displacement (V.Ivanchenko) 43 // 23-07-03 Replace G4Material by G4MaterialCu 46 // 23-07-03 Replace G4Material by G4MaterialCutCouple in dE/dx and CrossSection 44 // calculation (V.Ivanchenko) 47 // calculation (V.Ivanchenko) 45 // 01-03-04 L.Urban signature changed in Sampl 48 // 01-03-04 L.Urban signature changed in SampleCosineTheta 46 // 23-04-04 L.urban signature of SampleCosineT 49 // 23-04-04 L.urban signature of SampleCosineTheta changed back 47 // 17-11-04 Add method CrossSectionPerAtom (V. 50 // 17-11-04 Add method CrossSectionPerAtom (V.Ivanchenko) 48 // 14-03-05 Reduce number of pure virtual meth 51 // 14-03-05 Reduce number of pure virtual methods and make inline part 49 // separate (V.Ivanchenko) 52 // separate (V.Ivanchenko) 50 // 24-03-05 Remove IsInCharge and add G4VParti 53 // 24-03-05 Remove IsInCharge and add G4VParticleChange in the constructor (VI) 51 // 08-04-05 Major optimisation of internal int 54 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko) 52 // 15-04-05 optimize internal interface for ms 55 // 15-04-05 optimize internal interface for msc (V.Ivanchenko) 53 // 08-05-05 A -> N (V.Ivanchenko) 56 // 08-05-05 A -> N (V.Ivanchenko) 54 // 25-07-05 Move constructor and destructor to 57 // 25-07-05 Move constructor and destructor to the body (V.Ivanchenko) 55 // 02-02-06 ComputeCrossSectionPerAtom: defaul 58 // 02-02-06 ComputeCrossSectionPerAtom: default value A=0. (mma) 56 // 06-02-06 add method ComputeMeanFreePath() ( 59 // 06-02-06 add method ComputeMeanFreePath() (mma) 57 // 07-03-06 Optimize msc methods (V.Ivanchenko 60 // 07-03-06 Optimize msc methods (V.Ivanchenko) 58 // 29-06-06 Add member currentElement and Get/ << 59 // 29-10-07 Added SampleScattering (V.Ivanchen << 60 // 15-07-08 Reorder class members and improve << 61 // 21-07-08 Added vector of G4ElementSelector << 62 // 12-09-08 Added methods GetParticleCharge, G << 63 // CorrectionsAlongStep, ActivateNucl << 64 // 16-02-09 Moved implementations of virtual m << 65 // 07-04-09 Moved msc methods from G4VEmModel << 66 // 13-10-10 Added G4VEmAngularDistribution (VI << 67 // 61 // 68 // Class Description: 62 // Class Description: 69 // 63 // 70 // Abstract interface to energy loss models 64 // Abstract interface to energy loss models 71 65 72 // ------------------------------------------- 66 // ------------------------------------------------------------------- 73 // 67 // 74 68 75 #ifndef G4VEmModel_h 69 #ifndef G4VEmModel_h 76 #define G4VEmModel_h 1 70 #define G4VEmModel_h 1 77 71 78 #include "globals.hh" 72 #include "globals.hh" 79 #include "G4DynamicParticle.hh" 73 #include "G4DynamicParticle.hh" 80 #include "G4ParticleDefinition.hh" 74 #include "G4ParticleDefinition.hh" 81 #include "G4MaterialCutsCouple.hh" 75 #include "G4MaterialCutsCouple.hh" 82 #include "G4Material.hh" 76 #include "G4Material.hh" 83 #include "G4Element.hh" 77 #include "G4Element.hh" 84 #include "G4ElementVector.hh" 78 #include "G4ElementVector.hh" 85 #include "G4Isotope.hh" << 86 #include "G4DataVector.hh" 79 #include "G4DataVector.hh" 87 #include "G4VEmFluctuationModel.hh" 80 #include "G4VEmFluctuationModel.hh" 88 #include "G4VEmAngularDistribution.hh" << 81 #include "Randomize.hh" 89 #include "G4EmElementSelector.hh" << 90 #include <CLHEP/Random/RandomEngine.h> << 91 #include <vector> << 92 82 93 class G4ElementData; << 94 class G4PhysicsTable; 83 class G4PhysicsTable; 95 class G4Region; 84 class G4Region; 96 class G4VParticleChange; 85 class G4VParticleChange; 97 class G4ParticleChangeForLoss; << 98 class G4ParticleChangeForGamma; << 99 class G4Track; 86 class G4Track; 100 class G4LossTableManager; << 101 87 102 class G4VEmModel 88 class G4VEmModel 103 { 89 { 104 90 105 public: 91 public: 106 92 107 explicit G4VEmModel(const G4String& nam); << 93 G4VEmModel(const G4String& nam); 108 94 109 virtual ~G4VEmModel(); 95 virtual ~G4VEmModel(); 110 96 111 //------------------------------------------ 97 //------------------------------------------------------------------------ 112 // Virtual methods to be implemented for any << 98 // Virtual methods to be implemented for the concrete model 113 //------------------------------------------ 99 //------------------------------------------------------------------------ 114 100 115 virtual void Initialise(const G4ParticleDefi 101 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&) = 0; 116 102 117 virtual void SampleSecondaries(std::vector<G << 118 const G4Mater << 119 const G4Dynam << 120 G4double tmin << 121 G4double tmax << 122 << 123 //------------------------------------------ << 124 // Methods for initialisation of MT; may be << 125 //------------------------------------------ << 126 103 127 // initialisation in local thread << 104 virtual std::vector<G4DynamicParticle*>* SampleSecondaries( 128 virtual void InitialiseLocal(const G4Particl << 105 const G4MaterialCutsCouple*, 129 G4VEmModel* mas << 106 const G4DynamicParticle*, 130 << 107 G4double tmin = 0.0, 131 // initialisation of a new material at run t << 108 G4double tmax = DBL_MAX) = 0; 132 virtual void InitialiseForMaterial(const G4P << 133 const G4M << 134 << 135 // initialisation of a new element at run ti << 136 virtual void InitialiseForElement(const G4Pa << 137 G4int Z); << 138 109 139 //------------------------------------------ 110 //------------------------------------------------------------------------ 140 // Methods with standard implementation; may 111 // Methods with standard implementation; may be overwritten if needed 141 //------------------------------------------ 112 //------------------------------------------------------------------------ 142 113 143 // main method to compute dEdx << 114 virtual G4double MinEnergyCut(const G4ParticleDefinition*, 144 virtual G4double ComputeDEDXPerVolume(const << 115 const G4MaterialCutsCouple*); 145 const << 146 G4doub << 147 G4doub << 148 << 149 // main method to compute cross section per << 150 virtual G4double CrossSectionPerVolume(const << 151 const << 152 G4dou << 153 G4dou << 154 G4dou << 155 << 156 // method to get partial cross section << 157 virtual G4double GetPartialCrossSection(cons << 158 G4in << 159 cons << 160 G4do << 161 << 162 // main method to compute cross section per << 163 virtual G4double ComputeCrossSectionPerAtom( << 164 << 165 << 166 << 167 << 168 << 169 << 170 // main method to compute cross section per << 171 virtual G4double ComputeCrossSectionPerShell << 172 << 173 << 174 << 175 << 176 116 177 // Compute effective ion charge square << 117 virtual G4double ComputeDEDX(const G4MaterialCutsCouple*, 178 virtual G4double ChargeSquareRatio(const G4T << 118 const G4ParticleDefinition*, >> 119 G4double kineticEnergy, >> 120 G4double cutEnergy = DBL_MAX); >> 121 >> 122 virtual G4double CrossSection(const G4MaterialCutsCouple*, >> 123 const G4ParticleDefinition*, >> 124 G4double kineticEnergy, >> 125 G4double cutEnergy = 0.0, >> 126 G4double maxEnergy = DBL_MAX); 179 127 180 // Compute effective ion charge square << 128 virtual G4double ComputeDEDXPerVolume( 181 virtual G4double GetChargeSquareRatio(const << 129 const G4Material*, 182 const << 130 const G4ParticleDefinition*, 183 G4doub << 131 G4double kineticEnergy, >> 132 G4double cutEnergy = DBL_MAX); 184 133 185 // Compute ion charge << 186 virtual G4double GetParticleCharge(const G4P << 187 const G4M << 188 G4double << 189 << 190 // Initialisation for a new track << 191 virtual void StartTracking(G4Track*); << 192 << 193 // add correction to energy loss and compute << 194 virtual void CorrectionsAlongStep(const G4Ma << 195 const G4Dy << 196 const G4do << 197 G4double& << 198 << 199 // value which may be tabulated (by default << 200 virtual G4double Value(const G4MaterialCutsC << 201 const G4ParticleDefin << 202 G4double kineticEnerg << 203 << 204 // threshold for zero value << 205 virtual G4double MinPrimaryEnergy(const G4Ma << 206 const G4Pa << 207 G4double c << 208 134 209 // model can define low-energy limit for the << 135 virtual G4double ComputeCrossSectionPerAtom( 210 virtual G4double MinEnergyCut(const G4Partic << 136 const G4ParticleDefinition*, 211 const G4Materi << 137 G4double kinEnergy, 212 << 138 G4double Z, 213 // initialisation at run time for a given ma << 139 G4double A = 0., 214 virtual void SetupForMaterial(const G4Partic << 140 G4double cutEnergy = 0.0, >> 141 G4double maxEnergy = DBL_MAX); >> 142 >> 143 virtual G4double ComputeMeanFreePath( >> 144 const G4ParticleDefinition*, >> 145 G4double kineticEnergy, >> 146 const G4Material*, >> 147 G4double cutEnergy = 0.0, >> 148 G4double maxEnergy = DBL_MAX); >> 149 >> 150 virtual G4double CrossSectionPerVolume( 215 const G4Materi 151 const G4Material*, 216 G4double kinet << 152 const G4ParticleDefinition*, >> 153 G4double kineticEnergy, >> 154 G4double cutEnergy = 0.0, >> 155 G4double maxEnergy = DBL_MAX); >> 156 217 157 218 // add a region for the model << 158 219 virtual void DefineForRegion(const G4Region* << 159 // Methods for msc simulation 220 160 221 // fill number of different type of secondar << 161 virtual G4double ComputeTruePathLengthLimit( 222 virtual void FillNumberOfSecondaries(G4int& << 162 const G4Track& track, 223 G4int& << 163 G4PhysicsTable* theLambdaTable, >> 164 G4double currentMinimalStep); 224 165 225 // for automatic documentation << 166 virtual G4double ComputeGeomPathLength(G4double truePathLength); 226 virtual void ModelDescription(std::ostream& << 227 167 228 protected: << 168 virtual G4double ComputeTrueStepLength(G4double geomPathLength); 229 169 230 // initialisation of the ParticleChange for << 170 virtual void DefineForRegion(const G4Region*); 231 G4ParticleChangeForLoss* GetParticleChangeFo << 232 171 233 // initialisation of the ParticleChange for << 172 protected: 234 G4ParticleChangeForGamma* GetParticleChangeF << 235 173 236 // kinematically allowed max kinetic energy << 237 virtual G4double MaxSecondaryEnergy(const G4 174 virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*, 238 G4double << 175 G4double kineticEnergy); 239 << 240 public: << 241 176 242 //------------------------------------------ 177 //------------------------------------------------------------------------ 243 // Generic methods common to all models 178 // Generic methods common to all models 244 //------------------------------------------ 179 //------------------------------------------------------------------------ 245 180 246 // should be called at initialisation to bui << 181 public: 247 void InitialiseElementSelectors(const G4Part << 248 const G4Data << 249 << 250 // should be called at initialisation to acc << 251 inline std::vector<G4EmElementSelector*>* Ge << 252 << 253 // should be called at initialisation to set << 254 inline void SetElementSelectors(std::vector< << 255 << 256 // dEdx per unit length, base material appro << 257 inline G4double ComputeDEDX( const G4Materia << 258 const G4Particl << 259 G4double kineti << 260 G4double cutEne << 261 << 262 // cross section per volume, base material a << 263 inline G4double CrossSection(const G4Materia << 264 const G4Particl << 265 G4double kineti << 266 G4double cutEne << 267 G4double maxEne << 268 << 269 // compute mean free path via cross section << 270 inline G4double ComputeMeanFreePath(const G4 << 271 G4double << 272 const G4 << 273 G4double << 274 G4double << 275 << 276 // generic cross section per element << 277 inline G4double ComputeCrossSectionPerAtom(c << 278 c << 279 G << 280 G << 281 G << 282 << 283 // atom can be selected effitiantly if eleme << 284 inline const G4Element* SelectRandomAtom(con << 285 con << 286 G4d << 287 G4d << 288 G4d << 289 // same as SelectRandomAtom above but more e << 290 inline const G4Element* SelectTargetAtom(con << 291 con << 292 G4d << 293 G4d << 294 G4d << 295 G4d << 296 << 297 // to select atom cross section per volume i << 298 const G4Element* SelectRandomAtom(const G4Ma << 299 const G4Pa << 300 G4double k << 301 G4double c << 302 G4double m << 303 << 304 // to select atom if cross section is propor << 305 const G4Element* GetCurrentElement(const G4M << 306 G4int SelectRandomAtomNumber(const G4Materia << 307 << 308 // select isotope in order to have precise m << 309 const G4Isotope* GetCurrentIsotope(const G4E << 310 G4int SelectIsotopeNumber(const G4Element*) << 311 << 312 //------------------------------------------ << 313 // Get/Set methods << 314 //------------------------------------------ << 315 << 316 void SetParticleChange(G4VParticleChange*, G << 317 << 318 void SetCrossSectionTable(G4PhysicsTable*, G << 319 << 320 inline G4ElementData* GetElementData(); << 321 << 322 inline G4PhysicsTable* GetCrossSectionTable( << 323 << 324 inline G4VEmFluctuationModel* GetModelOfFluc << 325 << 326 inline G4VEmAngularDistribution* GetAngularD << 327 << 328 inline G4VEmModel* GetTripletModel(); << 329 << 330 inline void SetTripletModel(G4VEmModel*); << 331 << 332 inline void SetAngularDistribution(G4VEmAngu << 333 << 334 inline G4double HighEnergyLimit() const; << 335 << 336 inline G4double LowEnergyLimit() const; << 337 << 338 inline G4double HighEnergyActivationLimit() << 339 << 340 inline G4double LowEnergyActivationLimit() c << 341 << 342 inline G4double PolarAngleLimit() const; << 343 << 344 inline G4double SecondaryThreshold() const; << 345 << 346 inline G4bool DeexcitationFlag() const; << 347 << 348 inline G4bool ForceBuildTableFlag() const; << 349 << 350 inline G4bool UseAngularGeneratorFlag() cons << 351 << 352 inline void SetAngularGeneratorFlag(G4bool); << 353 << 354 inline void SetHighEnergyLimit(G4double); << 355 << 356 inline void SetLowEnergyLimit(G4double); << 357 << 358 inline void SetActivationHighEnergyLimit(G4d << 359 << 360 inline void SetActivationLowEnergyLimit(G4do << 361 << 362 inline G4bool IsActive(G4double kinEnergy) c << 363 << 364 inline void SetPolarAngleLimit(G4double); << 365 << 366 inline void SetSecondaryThreshold(G4double); << 367 << 368 inline void SetDeexcitationFlag(G4bool val); << 369 << 370 inline void SetForceBuildTable(G4bool val); << 371 << 372 inline void SetFluctuationFlag(G4bool val); << 373 << 374 inline G4bool IsMaster() const; << 375 << 376 inline void SetUseBaseMaterials(G4bool val); << 377 << 378 inline G4bool UseBaseMaterials() const; << 379 << 380 inline G4double MaxSecondaryKinEnergy(const << 381 << 382 inline const G4String& GetName() const; << 383 182 384 inline void SetCurrentCouple(const G4Materia << 183 void SetParticleChange(G4VParticleChange*, G4VEmFluctuationModel*); 385 184 386 inline G4bool IsLocked() const; << 185 G4VEmFluctuationModel* GetModelOfFluctuations(); 387 186 388 inline void SetLocked(G4bool); << 187 G4double HighEnergyLimit(); 389 188 390 // obsolete methods << 189 G4double LowEnergyLimit(); 391 [[deprecated("Use G4EmParameters::Instance() << 392 void SetLPMFlag(G4bool); << 393 190 394 void SetMasterThread(G4bool); << 191 void SetHighEnergyLimit(G4double); 395 192 396 // hide assignment operator << 193 void SetLowEnergyLimit(G4double); 397 G4VEmModel & operator=(const G4VEmModel &ri << 398 G4VEmModel(const G4VEmModel&) = delete; << 399 194 400 protected: << 195 G4double MaxSecondaryKinEnergy(const G4DynamicParticle* dynParticle); 401 196 402 inline const G4MaterialCutsCouple* CurrentCo << 197 const G4Element* SelectRandomAtom(const G4Material*, >> 198 const G4ParticleDefinition*, >> 199 G4double kineticEnergy, >> 200 G4double cutEnergy = 0.0, >> 201 G4double maxEnergy = DBL_MAX); 403 202 404 inline void SetCurrentElement(const G4Elemen << 203 const G4String& GetName() const; 405 204 406 private: 205 private: 407 206 408 // ======== Parameters of the class fixed at << 207 // hide assignment operator 409 << 208 G4VEmModel & operator=(const G4VEmModel &right); 410 G4VEmFluctuationModel* flucModel = null << 209 G4VEmModel(const G4VEmModel&); 411 G4VEmAngularDistribution* anglModel = null << 412 G4VEmModel* fTripletModel = << 413 const G4MaterialCutsCouple* fCurrentCouple = << 414 const G4Element* fCurrentElement << 415 std::vector<G4EmElementSelector*>* elmSelect << 416 G4LossTableManager* fEmManager; << 417 << 418 protected: << 419 << 420 G4ElementData* fElementData = << 421 G4VParticleChange* pParticleChange << 422 G4PhysicsTable* xSectionTable = << 423 const G4Material* pBaseMaterial = << 424 const std::vector<G4double>* theDensityFacto << 425 const std::vector<G4int>* theDensityIdx = << 426 << 427 G4double inveplus; << 428 G4double pFactor = 1.0; << 429 210 430 private: << 211 G4double lowLimit; >> 212 G4double highLimit; >> 213 G4double xsec[40]; 431 214 432 G4double lowLimit; << 215 G4VEmFluctuationModel* fluc; 433 G4double highLimit; << 434 G4double eMinActive = 0.0; << 435 G4double eMaxActive = DBL_MAX; << 436 G4double secondaryThreshold = DBL_MAX; << 437 G4double polarAngleLimit; << 438 216 439 G4int nSelectors = 0; << 217 const G4String name; 440 G4int nsec = 5; << 441 218 442 protected: 219 protected: 443 220 444 std::size_t currentCoupleIndex = 0; << 221 G4VParticleChange* pParticleChange; 445 std::size_t basedCoupleIndex = 0; << 446 G4bool lossFlucFlag = true; << 447 << 448 private: << 449 << 450 G4bool flagDeexcitation = false; << 451 G4bool flagForceBuildTable = false; << 452 G4bool isMaster = true; << 453 << 454 G4bool localTable = true; << 455 G4bool localElmSelectors = true; << 456 G4bool useAngularGenerator = false; << 457 G4bool useBaseMaterials = false; << 458 G4bool isLocked = false; << 459 << 460 const G4String name; << 461 std::vector<G4double> xsec; << 462 << 463 }; 222 }; 464 223 465 // ======== Run time inline methods ========== << 224 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 225 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 466 226 467 inline void G4VEmModel::SetCurrentCouple(const << 227 inline G4double G4VEmModel::HighEnergyLimit() 468 { << 469 if(fCurrentCouple != ptr) { << 470 fCurrentCouple = ptr; << 471 basedCoupleIndex = currentCoupleIndex = pt << 472 pBaseMaterial = ptr->GetMaterial(); << 473 pFactor = 1.0; << 474 if(useBaseMaterials) { << 475 basedCoupleIndex = (*theDensityIdx)[curr << 476 if(nullptr != pBaseMaterial->GetBaseMate << 477 pBaseMaterial = pBaseMaterial->GetBaseMateri << 478 pFactor = (*theDensityFactor)[currentCou << 479 } << 480 } << 481 } << 482 << 483 //....oooOO0OOooo........oooOO0OOooo........oo << 484 << 485 inline const G4MaterialCutsCouple* G4VEmModel: << 486 { << 487 return fCurrentCouple; << 488 } << 489 << 490 //....oooOO0OOooo........oooOO0OOooo........oo << 491 << 492 inline void G4VEmModel::SetCurrentElement(cons << 493 { << 494 fCurrentElement = elm; << 495 } << 496 << 497 //....oooOO0OOooo........oooOO0OOooo........oo << 498 << 499 inline << 500 G4double G4VEmModel::MaxSecondaryKinEnergy(con << 501 { << 502 return MaxSecondaryEnergy(dynPart->GetPartic << 503 dynPart->GetKineti << 504 } << 505 << 506 //....oooOO0OOooo........oooOO0OOooo........oo << 507 << 508 inline G4double G4VEmModel::ComputeDEDX(const << 509 const << 510 G4doub << 511 G4doub << 512 { << 513 SetCurrentCouple(couple); << 514 return pFactor*ComputeDEDXPerVolume(pBaseMat << 515 } << 516 << 517 //....oooOO0OOooo........oooOO0OOooo........oo << 518 << 519 inline G4double G4VEmModel::CrossSection(const << 520 const << 521 G4dou << 522 G4dou << 523 G4dou << 524 { << 525 SetCurrentCouple(couple); << 526 return pFactor*CrossSectionPerVolume(pBaseMa << 527 cutEner << 528 } << 529 << 530 //....oooOO0OOooo........oooOO0OOooo........oo << 531 << 532 inline << 533 G4double G4VEmModel::ComputeMeanFreePath(const << 534 G4dou << 535 const << 536 G4dou << 537 G4dou << 538 { << 539 G4double cross = CrossSectionPerVolume(mater << 540 return (cross > 0.0) ? 1./cross : DBL_MAX; << 541 } << 542 << 543 //....oooOO0OOooo........oooOO0OOooo........oo << 544 << 545 inline G4double << 546 G4VEmModel::ComputeCrossSectionPerAtom(const G << 547 const G << 548 G4doubl << 549 G4doubl << 550 G4doubl << 551 { << 552 fCurrentElement = elm; << 553 return ComputeCrossSectionPerAtom(part,kinEn << 554 cutEnergy, << 555 } << 556 << 557 //....oooOO0OOooo........oooOO0OOooo........oo << 558 << 559 inline const G4Element* << 560 G4VEmModel::SelectRandomAtom(const G4MaterialC << 561 const G4ParticleD << 562 G4double kinEnerg << 563 G4double cutEnerg << 564 G4double maxEnerg << 565 { << 566 SetCurrentCouple(couple); << 567 fCurrentElement = (nSelectors > 0) ? << 568 ((*elmSelectors)[couple->GetIndex()])->Sel << 569 SelectRandomAtom(pBaseMaterial,part,kinEne << 570 return fCurrentElement; << 571 } << 572 << 573 //....oooOO0OOooo........oooOO0OOooo........oo << 574 << 575 inline const G4Element* << 576 G4VEmModel::SelectTargetAtom(const G4MaterialC << 577 const G4ParticleD << 578 G4double kinEnerg << 579 G4double logKinE, << 580 G4double cutEnerg << 581 G4double maxEnerg << 582 { << 583 SetCurrentCouple(couple); << 584 fCurrentElement = (nSelectors > 0) << 585 ? ((*elmSelectors)[couple->GetIndex()])->Se << 586 : SelectRandomAtom(pBaseMaterial,part,kinEn << 587 return fCurrentElement; << 588 } << 589 << 590 // ======== Get/Set inline methods used at ini << 591 << 592 inline G4VEmFluctuationModel* G4VEmModel::GetM << 593 { << 594 return flucModel; << 595 } << 596 << 597 //....oooOO0OOooo........oooOO0OOooo........oo << 598 << 599 inline G4VEmAngularDistribution* G4VEmModel::G << 600 { << 601 return anglModel; << 602 } << 603 << 604 //....oooOO0OOooo........oooOO0OOooo........oo << 605 << 606 inline void G4VEmModel::SetAngularDistribution << 607 { << 608 if(p != anglModel) { << 609 delete anglModel; << 610 anglModel = p; << 611 } << 612 } << 613 << 614 //....oooOO0OOooo........oooOO0OOooo........oo << 615 << 616 inline G4VEmModel* G4VEmModel::GetTripletModel << 617 { << 618 return fTripletModel; << 619 } << 620 << 621 //....oooOO0OOooo........oooOO0OOooo........oo << 622 << 623 inline void G4VEmModel::SetTripletModel(G4VEmM << 624 { << 625 if(p != fTripletModel) { << 626 delete fTripletModel; << 627 fTripletModel = p; << 628 } << 629 } << 630 << 631 //....oooOO0OOooo........oooOO0OOooo........oo << 632 << 633 inline G4double G4VEmModel::HighEnergyLimit() << 634 { 228 { 635 return highLimit; 229 return highLimit; 636 } 230 } 637 231 638 //....oooOO0OOooo........oooOO0OOooo........oo << 232 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 639 233 640 inline G4double G4VEmModel::LowEnergyLimit() c << 234 inline G4double G4VEmModel::LowEnergyLimit() 641 { 235 { 642 return lowLimit; 236 return lowLimit; 643 } 237 } 644 238 645 //....oooOO0OOooo........oooOO0OOooo........oo << 239 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 646 << 647 inline G4double G4VEmModel::HighEnergyActivati << 648 { << 649 return eMaxActive; << 650 } << 651 240 652 //....oooOO0OOooo........oooOO0OOooo........oo << 241 inline void G4VEmModel::SetHighEnergyLimit(G4double val) 653 << 654 inline G4double G4VEmModel::LowEnergyActivatio << 655 { << 656 return eMinActive; << 657 } << 658 << 659 //....oooOO0OOooo........oooOO0OOooo........oo << 660 << 661 inline G4double G4VEmModel::PolarAngleLimit() << 662 { << 663 return polarAngleLimit; << 664 } << 665 << 666 //....oooOO0OOooo........oooOO0OOooo........oo << 667 << 668 inline G4double G4VEmModel::SecondaryThreshold << 669 { << 670 return secondaryThreshold; << 671 } << 672 << 673 //....oooOO0OOooo........oooOO0OOooo........oo << 674 << 675 inline G4bool G4VEmModel::DeexcitationFlag() c << 676 { << 677 return flagDeexcitation; << 678 } << 679 << 680 //....oooOO0OOooo........oooOO0OOooo........oo << 681 << 682 inline G4bool G4VEmModel::ForceBuildTableFlag( << 683 { << 684 return flagForceBuildTable; << 685 } << 686 << 687 //....oooOO0OOooo........oooOO0OOooo........oo << 688 << 689 inline G4bool G4VEmModel::UseAngularGeneratorF << 690 { << 691 return useAngularGenerator; << 692 } << 693 << 694 //....oooOO0OOooo........oooOO0OOooo........oo << 695 << 696 inline void G4VEmModel::SetAngularGeneratorFla << 697 { 242 { 698 useAngularGenerator = val; << 243 highLimit = val; 699 } 244 } 700 245 701 //....oooOO0OOooo........oooOO0OOooo........oo << 246 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 702 247 703 inline void G4VEmModel::SetFluctuationFlag(G4b << 248 inline void G4VEmModel::SetLowEnergyLimit(G4double val) 704 { 249 { 705 lossFlucFlag = val; << 250 lowLimit = val; 706 } 251 } 707 252 708 //....oooOO0OOooo........oooOO0OOooo........oo << 253 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 709 254 710 inline G4bool G4VEmModel::IsMaster() const << 255 inline void G4VEmModel::SetParticleChange(G4VParticleChange* p, >> 256 G4VEmFluctuationModel* f = 0) 711 { 257 { 712 return isMaster; << 258 if(p && pParticleChange != p) pParticleChange = p; >> 259 fluc = f; 713 } 260 } 714 261 715 //....oooOO0OOooo........oooOO0OOooo........oo << 262 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 716 << 717 inline void G4VEmModel::SetUseBaseMaterials(G4 << 718 { << 719 useBaseMaterials = val; << 720 } << 721 263 722 //....oooOO0OOooo........oooOO0OOooo........oo << 723 264 724 inline G4bool G4VEmModel::UseBaseMaterials() c << 265 inline G4VEmFluctuationModel* G4VEmModel::GetModelOfFluctuations() 725 { << 726 return useBaseMaterials; << 727 } << 728 << 729 //....oooOO0OOooo........oooOO0OOooo........oo << 730 << 731 inline void G4VEmModel::SetHighEnergyLimit(G4d << 732 { 266 { 733 highLimit = val; << 267 return fluc; 734 } 268 } 735 269 736 //....oooOO0OOooo........oooOO0OOooo........oo << 270 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 737 271 738 inline void G4VEmModel::SetLowEnergyLimit(G4do << 272 inline G4double G4VEmModel::MinEnergyCut(const G4ParticleDefinition*, >> 273 const G4MaterialCutsCouple*) 739 { 274 { 740 lowLimit = val; << 275 return 0.0; 741 } 276 } 742 277 743 //....oooOO0OOooo........oooOO0OOooo........oo << 278 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 744 279 745 inline void G4VEmModel::SetActivationHighEnerg << 280 inline G4double G4VEmModel::ComputeDEDXPerVolume( >> 281 const G4Material*, >> 282 const G4ParticleDefinition*, >> 283 G4double, >> 284 G4double) 746 { 285 { 747 eMaxActive = val; << 286 return 0.0; 748 } 287 } 749 288 750 //....oooOO0OOooo........oooOO0OOooo........oo << 289 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 751 290 752 inline void G4VEmModel::SetActivationLowEnergy << 291 inline G4double G4VEmModel::ComputeDEDX(const G4MaterialCutsCouple* c, >> 292 const G4ParticleDefinition* p, >> 293 G4double kinEnergy, >> 294 G4double cutEnergy) 753 { 295 { 754 eMinActive = val; << 296 return ComputeDEDXPerVolume(c->GetMaterial(),p,kinEnergy,cutEnergy); 755 } 297 } 756 298 757 //....oooOO0OOooo........oooOO0OOooo........oo << 299 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 758 300 759 inline G4bool G4VEmModel::IsActive(G4double ki << 301 inline G4double G4VEmModel::CrossSection(const G4MaterialCutsCouple* c, >> 302 const G4ParticleDefinition* p, >> 303 G4double kinEnergy, >> 304 G4double cutEnergy, >> 305 G4double maxEnergy) 760 { 306 { 761 return (kinEnergy >= eMinActive && kinEnergy << 307 return CrossSectionPerVolume(c->GetMaterial(),p,kinEnergy,cutEnergy,maxEnergy); 762 } 308 } 763 309 764 //....oooOO0OOooo........oooOO0OOooo........oo << 310 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 765 311 766 inline void G4VEmModel::SetPolarAngleLimit(G4d << 312 inline G4double G4VEmModel::ComputeCrossSectionPerAtom( >> 313 const G4ParticleDefinition*, >> 314 G4double, G4double, G4double, >> 315 G4double, G4double) 767 { 316 { 768 if(!isLocked) { polarAngleLimit = val; } << 317 return 0.0; 769 } 318 } 770 319 771 //....oooOO0OOooo........oooOO0OOooo........oo << 320 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 772 321 773 inline void G4VEmModel::SetSecondaryThreshold( << 322 inline const G4Element* G4VEmModel::SelectRandomAtom( 774 { << 323 const G4Material* material, 775 secondaryThreshold = val; << 324 const G4ParticleDefinition* pd, >> 325 G4double kinEnergy, >> 326 G4double tcut, >> 327 G4double tmax) >> 328 { >> 329 const G4ElementVector* theElementVector = material->GetElementVector(); >> 330 const G4Element* elm = (*theElementVector)[0]; >> 331 G4int nelm = material->GetNumberOfElements() - 1; >> 332 if (nelm > 0) { >> 333 G4double x = G4UniformRand()* >> 334 CrossSectionPerVolume(material,pd,kinEnergy,tcut,tmax); >> 335 G4int i = -1; >> 336 do {i++;} while (x > xsec[i] && i < nelm); >> 337 elm = (*theElementVector)[i]; >> 338 } >> 339 return elm; 776 } 340 } 777 341 778 //....oooOO0OOooo........oooOO0OOooo........oo << 342 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 779 343 780 inline void G4VEmModel::SetDeexcitationFlag(G4 << 344 inline G4double G4VEmModel::MaxSecondaryKinEnergy( >> 345 const G4DynamicParticle* dynParticle) 781 { 346 { 782 flagDeexcitation = val; << 347 return MaxSecondaryEnergy(dynParticle->GetDefinition(), >> 348 dynParticle->GetKineticEnergy()); 783 } 349 } 784 350 785 //....oooOO0OOooo........oooOO0OOooo........oo << 351 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 786 352 787 inline void G4VEmModel::SetForceBuildTable(G4b << 353 inline G4double G4VEmModel::MaxSecondaryEnergy(const G4ParticleDefinition*, >> 354 G4double kineticEnergy) 788 { 355 { 789 flagForceBuildTable = val; << 356 return kineticEnergy; 790 } 357 } 791 358 792 //....oooOO0OOooo........oooOO0OOooo........oo << 359 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 793 360 794 inline const G4String& G4VEmModel::GetName() c 361 inline const G4String& G4VEmModel::GetName() const 795 { 362 { 796 return name; 363 return name; 797 } 364 } 798 365 799 //....oooOO0OOooo........oooOO0OOooo........oo << 366 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 367 // Methods for msc simulation >> 368 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 800 369 801 inline std::vector<G4EmElementSelector*>* G4VE << 370 inline G4double G4VEmModel::ComputeTruePathLengthLimit( >> 371 const G4Track&, >> 372 G4PhysicsTable*, >> 373 G4double) 802 { 374 { 803 return elmSelectors; << 375 return DBL_MAX; 804 } 376 } 805 377 806 //....oooOO0OOooo........oooOO0OOooo........oo << 378 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 807 379 808 inline void << 380 inline G4double G4VEmModel::ComputeGeomPathLength(G4double truePathLength) 809 G4VEmModel::SetElementSelectors(std::vector<G4 << 810 { 381 { 811 if(p != elmSelectors) { << 382 return truePathLength; 812 elmSelectors = p; << 813 nSelectors = (nullptr != elmSelectors) ? G << 814 localElmSelectors = false; << 815 } << 816 } 383 } 817 384 818 //....oooOO0OOooo........oooOO0OOooo........oo << 385 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 819 386 820 inline G4ElementData* G4VEmModel::GetElementDa << 387 inline G4double G4VEmModel::ComputeTrueStepLength(G4double geomPathLength) 821 { 388 { 822 return fElementData; << 389 return geomPathLength; 823 } 390 } 824 391 825 //....oooOO0OOooo........oooOO0OOooo........oo << 392 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 826 393 827 inline G4PhysicsTable* G4VEmModel::GetCrossSec << 394 inline void G4VEmModel::DefineForRegion(const G4Region*) 828 { << 395 {} 829 return xSectionTable; << 830 } << 831 396 832 //....oooOO0OOooo........oooOO0OOooo........oo << 397 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 833 << 834 inline G4bool G4VEmModel::IsLocked() const << 835 { << 836 return isLocked; << 837 } << 838 << 839 //....oooOO0OOooo........oooOO0OOooo........oo << 840 << 841 inline void G4VEmModel::SetLocked(G4bool val) << 842 { << 843 isLocked = val; << 844 } << 845 << 846 //....oooOO0OOooo........oooOO0OOooo........oo << 847 398 848 #endif 399 #endif >> 400 849 401