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: G4PenelopeRayleighModel.hh 75573 2013-11-04 11:48:15Z gcosmo $ 26 // 27 // 27 // Author: Luciano Pandola 28 // Author: Luciano Pandola 28 // 29 // 29 // History: 30 // History: 30 // ----------- 31 // ----------- 31 // 03 Dec 2009 L. Pandola 1st implementati 32 // 03 Dec 2009 L. Pandola 1st implementation. 32 // 25 May 2011 L. Pandola Renamed (make v2 33 // 25 May 2011 L. Pandola Renamed (make v2008 as default Penelope) 33 // 27 Sep 2013 L. Pandola Migration to MT p 34 // 27 Sep 2013 L. Pandola Migration to MT paradigm 34 // 04 Mar 2021 L. Pandola Replace maps wit << 35 // 35 // 36 // ------------------------------------------- 36 // ------------------------------------------------------------------- 37 // 37 // 38 // Class description: 38 // Class description: 39 // Low Energy Electromagnetic Physics, Rayleig 39 // Low Energy Electromagnetic Physics, Rayleigh Scattering 40 // with the model from Penelope, version 2008 40 // with the model from Penelope, version 2008 41 // ------------------------------------------- 41 // ------------------------------------------------------------------- 42 42 43 #ifndef G4PENELOPERAYLEIGHMODEL_HH 43 #ifndef G4PENELOPERAYLEIGHMODEL_HH 44 #define G4PENELOPERAYLEIGHMODEL_HH 1 44 #define G4PENELOPERAYLEIGHMODEL_HH 1 45 45 46 #include "globals.hh" 46 #include "globals.hh" 47 #include "G4VEmModel.hh" 47 #include "G4VEmModel.hh" 48 #include "G4DataVector.hh" 48 #include "G4DataVector.hh" 49 #include "G4ParticleChangeForGamma.hh" 49 #include "G4ParticleChangeForGamma.hh" 50 50 51 class G4ParticleDefinition; 51 class G4ParticleDefinition; 52 class G4DynamicParticle; 52 class G4DynamicParticle; 53 class G4MaterialCutsCouple; 53 class G4MaterialCutsCouple; 54 class G4Material; 54 class G4Material; 55 class G4PhysicsFreeVector; 55 class G4PhysicsFreeVector; 56 class G4PenelopeSamplingData; 56 class G4PenelopeSamplingData; 57 57 58 class G4PenelopeRayleighModel : public G4VEmMo 58 class G4PenelopeRayleighModel : public G4VEmModel 59 { 59 { 60 60 61 public: 61 public: 62 explicit G4PenelopeRayleighModel(const G4Par << 62 G4PenelopeRayleighModel(const G4ParticleDefinition* p=0, 63 const G4String& processName ="PenRa << 63 const G4String& processName ="PenRayleigh"); 64 64 65 virtual ~G4PenelopeRayleighModel(); 65 virtual ~G4PenelopeRayleighModel(); 66 66 67 void Initialise(const G4ParticleDefinition*, << 67 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&); 68 void InitialiseLocal(const G4ParticleDefinit << 68 virtual void InitialiseLocal(const G4ParticleDefinition*, 69 G4VEmModel *masterModel) override; << 69 G4VEmModel *masterModel); 70 << 70 71 G4double ComputeCrossSectionPerAtom(const G4 << 71 virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, 72 G4double kinEnergy, << 72 G4double kinEnergy, 73 G4double Z, << 73 G4double Z, 74 G4double A=0, << 74 G4double A=0, 75 G4double cut=0, << 75 G4double cut=0, 76 G4double emax=DBL_MAX) override; << 76 G4double emax=DBL_MAX); 77 << 77 78 void SampleSecondaries(std::vector<G4Dynamic << 78 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*, 79 const G4MaterialCutsCouple*, << 79 const G4MaterialCutsCouple*, 80 const G4DynamicParticle*, << 80 const G4DynamicParticle*, 81 G4double tmin, << 81 G4double tmin, 82 G4double maxEnergy) override; << 82 G4double maxEnergy); 83 83 84 void SetVerbosityLevel(G4int lev){fVerboseLe << 84 void SetVerbosityLevel(G4int lev){verboseLevel = lev;}; 85 G4int GetVerbosityLevel(){return fVerboseLev << 85 G4int GetVerbosityLevel(){return verboseLevel;}; 86 86 87 //Testing purposes 87 //Testing purposes 88 void DumpFormFactorTable(const G4Material*); 88 void DumpFormFactorTable(const G4Material*); 89 89 90 G4PenelopeRayleighModel& operator=(const G4P << 91 G4PenelopeRayleighModel(const G4PenelopeRayl << 92 << 93 protected: 90 protected: 94 G4ParticleChangeForGamma* fParticleChange; 91 G4ParticleChangeForGamma* fParticleChange; 95 const G4ParticleDefinition* fParticle; 92 const G4ParticleDefinition* fParticle; 96 93 97 private: 94 private: >> 95 G4PenelopeRayleighModel& operator=(const G4PenelopeRayleighModel &right); >> 96 G4PenelopeRayleighModel(const G4PenelopeRayleighModel&); >> 97 98 void SetParticle(const G4ParticleDefinition* 98 void SetParticle(const G4ParticleDefinition*); >> 99 >> 100 //Intrinsic energy limits of the model: cannot be extended by >> 101 //the parent process >> 102 G4double fIntrinsicLowEnergyLimit; >> 103 G4double fIntrinsicHighEnergyLimit; >> 104 >> 105 G4int verboseLevel; >> 106 G4bool isInitialised; >> 107 >> 108 //Internal tables and manager methods >> 109 std::map<G4int,G4PhysicsFreeVector*> *logAtomicCrossSection; >> 110 std::map<G4int,G4PhysicsFreeVector*> *atomicFormFactor; >> 111 >> 112 >> 113 G4DataVector logQSquareGrid; //log(Q^2) grid for interpolation >> 114 std::map<const G4Material*,G4PhysicsFreeVector*> *logFormFactorTable; //log(Q^2) vs. log(F^2) >> 115 >> 116 G4DataVector logEnergyGridPMax; //energy grid for PMac (and originally for the x-section) >> 117 std::map<const G4Material*,G4PhysicsFreeVector*> *pMaxTable; //E vs. Pmax >> 118 >> 119 std::map<const G4Material*,G4PenelopeSamplingData*> *samplingTable; >> 120 99 //Helper methods 121 //Helper methods 100 void ReadDataFile(G4int); 122 void ReadDataFile(G4int); 101 void ClearTables(); 123 void ClearTables(); 102 void BuildFormFactorTable(const G4Material*) 124 void BuildFormFactorTable(const G4Material*); 103 void GetPMaxTable(const G4Material*); 125 void GetPMaxTable(const G4Material*); 104 126 105 G4double GetFSquared(const G4Material*,const 127 G4double GetFSquared(const G4Material*,const G4double); 106 void InitializeSamplingAlgorithm(const G4Mat 128 void InitializeSamplingAlgorithm(const G4Material*); 107 129 108 G4DataVector fLogQSquareGrid; //log(Q^2) gri << 109 std::map<const G4Material*,G4PhysicsFreeVect << 110 << 111 G4DataVector fLogEnergyGridPMax; //energy gr << 112 std::map<const G4Material*,G4PhysicsFreeVect << 113 std::map<const G4Material*,G4PenelopeSamplin << 114 << 115 //Internal tables and manager methods << 116 static const G4int fMaxZ =99; << 117 static G4PhysicsFreeVector* fLogAtomicCrossS << 118 static G4PhysicsFreeVector* fAtomicFormFacto << 119 << 120 //Intrinsic energy limits of the model: cann << 121 //the parent process << 122 G4double fIntrinsicLowEnergyLimit; << 123 G4double fIntrinsicHighEnergyLimit; << 124 G4int fVerboseLevel; << 125 G4bool fIsInitialised; << 126 //Used only for G4EmCalculator and Unit Test 130 //Used only for G4EmCalculator and Unit Tests 127 G4bool fLocalTable; 131 G4bool fLocalTable; >> 132 128 }; 133 }; 129 134 130 #endif 135 #endif 131 136 132 137