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 // ******************************************* 26 // ********************************************************************* 27 // | 27 // | | 28 // | G4LowEPComptonModel-- Geant4 28 // | G4LowEPComptonModel-- Geant4 Monash University | 29 // | low energy Compton scat 29 // | low energy Compton scattering model. | 30 // | J. M. C. Brown, Monash Univer 30 // | J. M. C. Brown, Monash University, Australia | 31 // | ## Unpolarised photons 31 // | ## Unpolarised photons only ## | 32 // | 32 // | | 33 // | 33 // | | 34 // ******************************************* 34 // ********************************************************************* 35 // | 35 // | | 36 // | The following is a Geant4 class to simula 36 // | The following is a Geant4 class to simulate the process of | 37 // | bound electron Compton scattering. Genera 37 // | bound electron Compton scattering. General code structure is | 38 // | based on G4LowEnergyCompton.cc and G4Live 38 // | based on G4LowEnergyCompton.cc and G4LivermoreComptonModel.cc. | 39 // | Algorithms for photon energy, and ejected 39 // | Algorithms for photon energy, and ejected Compton electron | 40 // | direction taken from: 40 // | direction taken from: | 41 // | 41 // | | 42 // | J. M. C. Brown, M. R. Dimmock, J. E. Gill 42 // | J. M. C. Brown, M. R. Dimmock, J. E. Gillam and D. M. Paganin, | 43 // | "A low energy bound atomic electron Compt 43 // | "A low energy bound atomic electron Compton scattering model | 44 // | for Geant4", NIMB, Vol. 338, 77-88, 2014 << 44 // | for Geant4", IEEE Transactions on Nuclear Science, submitted. | 45 // | 45 // | | 46 // | The author acknowledges the work of the G 46 // | The author acknowledges the work of the Geant4 collaboration | 47 // | in developing the following algorithms th 47 // | in developing the following algorithms that have been employed | 48 // | or adapeted for the present software: 48 // | or adapeted for the present software: | 49 // | 49 // | | 50 // | # sampling of photon scattering angle, 50 // | # sampling of photon scattering angle, | 51 // | # target element selection in composite 51 // | # target element selection in composite materials, | 52 // | # target shell selection in element, 52 // | # target shell selection in element, | 53 // | # and sampling of bound electron momentu 53 // | # and sampling of bound electron momentum from Compton profiles. | 54 // | 54 // | | 55 // ******************************************* 55 // ********************************************************************* 56 // | 56 // | | 57 // | History: 57 // | History: | 58 // | -------- 58 // | -------- | 59 // | 59 // | | 60 // | Nov. 2011 JMCB - First version 60 // | Nov. 2011 JMCB - First version | 61 // | Feb. 2012 JMCB - Migration to Geant 61 // | Feb. 2012 JMCB - Migration to Geant4 9.5 | 62 // | Sep. 2012 JMCB - Final fixes for Ge 62 // | Sep. 2012 JMCB - Final fixes for Geant4 9.6 | 63 // | Feb. 2013 JMCB - Geant4 9.6 FPE fix << 64 // | Jan. 2015 JMCB - Migration to MT (B << 65 // | implementation) << 66 // | Feb. 2016 JMCB - Geant4 10.2 FPE fi << 67 // | 63 // | | 68 // ******************************************* 64 // ********************************************************************* 69 65 70 #ifndef G4LowEPComptonModel_h 66 #ifndef G4LowEPComptonModel_h 71 #define G4LowEPComptonModel_h 1 67 #define G4LowEPComptonModel_h 1 72 68 73 #include "G4VEmModel.hh" 69 #include "G4VEmModel.hh" 74 #include "G4PhysicsFreeVector.hh" << 70 #include "G4ShellData.hh" >> 71 #include "G4DopplerProfile.hh" 75 72 76 class G4ParticleChangeForGamma; 73 class G4ParticleChangeForGamma; >> 74 class G4VCrossSectionHandler; 77 class G4VAtomDeexcitation; 75 class G4VAtomDeexcitation; 78 class G4ShellData; << 76 class G4VEMDataSet; 79 class G4DopplerProfile; << 80 77 81 class G4LowEPComptonModel : public G4VEmModel 78 class G4LowEPComptonModel : public G4VEmModel 82 { 79 { 83 80 84 public: 81 public: 85 82 86 explicit G4LowEPComptonModel(const G4Particl << 83 G4LowEPComptonModel(const G4ParticleDefinition* p = 0, 87 const G4String& nam = "LowEPComp 84 const G4String& nam = "LowEPComptonModel"); 88 virtual ~G4LowEPComptonModel(); << 85 >> 86 virtual ~G4LowEPComptonModel(); 89 87 90 void Initialise(const G4ParticleDefinition*, << 88 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&); 91 void InitialiseLocal(const G4ParticleDefinit << 92 G4VEmModel* mas << 93 << 94 void InitialiseForElement(const G4ParticleDe << 95 << 96 G4double ComputeCrossSectionPerAtom( const G << 97 G4double kinEnergy, << 98 G4double Z, << 99 G4double A=0, << 100 G4double cut=0, << 101 G4double emax=DBL_MAX ) overrid << 102 << 103 void SampleSecondaries(std::vector<G4Dynamic << 104 const G4MaterialCutsCouple*, << 105 const G4DynamicParticle*, << 106 G4double tmin, << 107 G4double maxEnergy) override; << 108 89 109 G4LowEPComptonModel & operator=(const G4Low << 90 virtual G4double ComputeCrossSectionPerAtom( const G4ParticleDefinition*, 110 G4LowEPComptonModel(const G4LowEPComptonMod << 91 G4double kinEnergy, >> 92 G4double Z, >> 93 G4double A=0, >> 94 G4double cut=0, >> 95 G4double emax=DBL_MAX ); >> 96 >> 97 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*, >> 98 const G4MaterialCutsCouple*, >> 99 const G4DynamicParticle*, >> 100 G4double tmin, >> 101 G4double maxEnergy); 111 102 112 private: << 103 protected: 113 void ReadData(size_t Z, const char* path = 0 << 114 G4double ComputeScatteringFunction(G4double << 115 104 116 G4ParticleChangeForGamma* fParticleChange; 105 G4ParticleChangeForGamma* fParticleChange; 117 G4VAtomDeexcitation* fAtomDeexcitation; << 106 118 static G4ShellData* shellData; << 107 private: 119 static G4DopplerProfile* profileData; << 108 120 << 109 G4double lowEnergyLimit; 121 static const G4int maxZ = 99; << 110 G4double highEnergyLimit; 122 static G4PhysicsFreeVector* data[100]; << 123 static const G4double ScatFuncFitParam[101][ << 124 << 125 G4int verboseLevel; << 126 G4bool isInitialised; 111 G4bool isInitialised; >> 112 G4int verboseLevel; >> 113 >> 114 G4VEMDataSet* scatterFunctionData; >> 115 G4VCrossSectionHandler* crossSectionHandler; >> 116 >> 117 G4VAtomDeexcitation* fAtomDeexcitation; >> 118 >> 119 G4ShellData shellData; >> 120 G4DopplerProfile profileData; >> 121 >> 122 G4LowEPComptonModel & operator=(const G4LowEPComptonModel &right); >> 123 G4LowEPComptonModel(const G4LowEPComptonModel&); >> 124 >> 125 127 }; 126 }; 128 127 129 //******************************************** 128 //**************************************************************************** 130 129 131 #endif 130 #endif 132 131