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: G4eCoulombScatteringModel.hh,v 1.5 2006/10/19 09:44:27 vnivanch Exp $ >> 27 // GEANT4 tag $Name: geant4-08-02-patch-01-ref $ 26 // 28 // 27 // ------------------------------------------- 29 // ------------------------------------------------------------------- 28 // 30 // 29 // GEANT4 Class header file 31 // GEANT4 Class header file 30 // 32 // 31 // 33 // 32 // File name: G4eCoulombScatteringModel 34 // File name: G4eCoulombScatteringModel 33 // 35 // 34 // Author: Vladimir Ivanchenko 36 // Author: Vladimir Ivanchenko 35 // 37 // 36 // Creation date: 19.02.2006 38 // Creation date: 19.02.2006 37 // 39 // 38 // Modifications: 40 // Modifications: 39 // 01.08.06 V.Ivanchenko extend upper limit of 41 // 01.08.06 V.Ivanchenko extend upper limit of table to TeV and review the 40 // logic of building - only elements 42 // logic of building - only elements from G4ElementTable 41 // 08.08.06 V.Ivanchenko build internal table 43 // 08.08.06 V.Ivanchenko build internal table in ekin scale, introduce faclim 42 // 19.08.06 V.Ivanchenko add inline function S 44 // 19.08.06 V.Ivanchenko add inline function ScreeningParameter and 43 // make some members pro 45 // make some members protected 44 // 09.10.07 V.Ivanchenko reorganized methods, << 45 // 09.06.08 V.Ivanchenko add SelectIsotope and << 46 // 17.06.09 C.Consoalndi modified SetupTarget << 47 // 27.05.10 V.Ivanchenko added G4WentzelOKandV << 48 // compute cross sections and sam << 49 // << 50 // 46 // 51 // Class Description: 47 // Class Description: 52 // 48 // 53 // Implementation of eCoulombScattering of a c << 49 // Implementation of eCoulombScattering of pointlike charge particle 54 // on Atomic Nucleus for interval of scatterin 50 // on Atomic Nucleus for interval of scattering anles in Lab system 55 // thetaMin - ThetaMax. << 51 // thetaMin - ThetaMax, nucleus recoil is neglected. 56 // The model based on analysis of J.M.Fernan 52 // The model based on analysis of J.M.Fernandez-Varea et al. 57 // NIM B73(1993)447 originated from G.Wentzel 53 // NIM B73(1993)447 originated from G.Wentzel Z.Phys. 40(1927)590 with 58 // screening parameter from H.A.Bethe Phys. Re 54 // screening parameter from H.A.Bethe Phys. Rev. 89 (1953) 1256. 59 // 55 // 60 56 61 // ------------------------------------------- 57 // ------------------------------------------------------------------- 62 // 58 // 63 59 64 #ifndef G4eCoulombScatteringModel_h 60 #ifndef G4eCoulombScatteringModel_h 65 #define G4eCoulombScatteringModel_h 1 61 #define G4eCoulombScatteringModel_h 1 66 62 67 #include "G4VEmModel.hh" 63 #include "G4VEmModel.hh" >> 64 #include "G4PhysicsTable.hh" 68 #include "globals.hh" 65 #include "globals.hh" 69 #include "G4MaterialCutsCouple.hh" << 70 #include "G4WentzelOKandVIxSection.hh" << 71 66 72 class G4ParticleChangeForGamma; 67 class G4ParticleChangeForGamma; 73 class G4ParticleDefinition; << 74 class G4IonTable; << 75 class G4NistManager; << 76 68 77 class G4eCoulombScatteringModel : public G4VEm 69 class G4eCoulombScatteringModel : public G4VEmModel 78 { 70 { 79 71 80 public: 72 public: 81 73 82 explicit G4eCoulombScatteringModel(G4bool co << 74 G4eCoulombScatteringModel(G4double thetaMin = 0.0, G4double thetaMax = pi, >> 75 G4bool build = true, G4double tlim = TeV*TeV, >> 76 const G4String& nam = "eCoulombScattering"); 83 77 84 ~G4eCoulombScatteringModel() override; << 78 virtual ~G4eCoulombScatteringModel(); 85 79 86 void Initialise(const G4ParticleDefinition*, << 80 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&); 87 81 88 void InitialiseLocal(const G4ParticleDefinit << 82 virtual G4double ComputeCrossSectionPerAtom( 89 G4VEmModel* masterModel << 90 << 91 G4double ComputeCrossSectionPerAtom( << 92 const G4Partic 83 const G4ParticleDefinition*, 93 G4double kinEnergy, << 84 G4double kinEnergy, 94 G4double Z, << 85 G4double Z, 95 G4double A, << 86 G4double A, 96 G4double cut, << 87 G4double cut, 97 G4double emax) override; << 88 G4double emax); 98 << 89 99 void SampleSecondaries(std::vector<G4Dynamic << 90 virtual std::vector<G4DynamicParticle*>* SampleSecondaries( 100 const G4MaterialCutsCouple*, << 91 const G4MaterialCutsCouple*, 101 const G4DynamicParticle*, << 92 const G4DynamicParticle*, 102 G4double tmin, << 93 G4double tmin, 103 G4double maxEnergy) override; << 94 G4double maxEnergy); 104 << 105 G4double MinPrimaryEnergy(const G4Material*, << 106 G4double) final; << 107 << 108 // defines low energy limit of the model << 109 inline void SetLowEnergyThreshold(G4double v << 110 95 111 // user definition of low-energy threshold o << 96 protected: 112 inline void SetRecoilThreshold(G4double eth) << 113 97 114 // defines low energy limit on energy transf << 98 G4double ScreeningParameter(G4double Z, G4double q2, 115 inline void SetFixedCut(G4double); << 99 G4double mom2, G4double invbeta2); 116 100 117 // low energy limit on energy transfer to at << 101 private: 118 inline G4double GetFixedCut() const; << 102 >> 103 G4double CalculateCrossSectionPerAtom(const G4ParticleDefinition*, >> 104 G4double kinEnergy, >> 105 G4double Z); 119 106 120 // hide assignment operator 107 // hide assignment operator 121 G4eCoulombScatteringModel & operator= << 108 G4eCoulombScatteringModel & operator=(const G4eCoulombScatteringModel &right); 122 (const G4eCoulombScatteringModel &right) = d << 109 G4eCoulombScatteringModel(const G4eCoulombScatteringModel&); 123 G4eCoulombScatteringModel(const G4eCoulombS << 124 110 125 protected: 111 protected: 126 112 127 inline void DefineMaterial(const G4MaterialC << 113 G4ParticleChangeForGamma* fParticleChange; 128 << 129 inline void SetupParticle(const G4ParticleDe << 130 << 131 private: << 132 << 133 G4IonTable* theIonTable; << 134 G4ParticleChangeForGamma* fParticleChange = << 135 G4WentzelOKandVIxSection* wokvi; << 136 G4NistManager* fNistManager; << 137 << 138 const std::vector<G4double>* pCuts = nullptr << 139 << 140 const G4ParticleDefinition* particle = nullp << 141 const G4ParticleDefinition* theProton; << 142 << 143 const G4MaterialCutsCouple* currentCouple = << 144 const G4Material* currentMaterial = nullptr; << 145 G4int currentMaterialIndex = 0; << 146 << 147 G4double cosThetaMin = 1.0; << 148 G4double cosThetaMax = -1.0; << 149 G4double recoilThreshold = 0.0; << 150 G4double elecRatio = 0.0; << 151 G4double fixedCut = -1.0; << 152 G4double mass; << 153 << 154 G4bool isCombined; << 155 }; << 156 114 157 //....oooOO0OOooo........oooOO0OOooo........oo << 115 G4double coeff; >> 116 G4double cosThetaMin; >> 117 G4double cosThetaMax; >> 118 G4double q2Limit; 158 119 159 inline << 160 void G4eCoulombScatteringModel::DefineMaterial << 161 { << 162 if(cup != currentCouple) { << 163 currentCouple = cup; << 164 currentMaterial = cup->GetMaterial(); << 165 currentMaterialIndex = currentCouple->GetI << 166 } << 167 } << 168 120 169 //....oooOO0OOooo........oooOO0OOooo........oo << 121 private: 170 122 171 inline << 123 G4PhysicsTable* theCrossSectionTable; 172 void G4eCoulombScatteringModel::SetupParticle( << 173 { << 174 // Initialise mass and charge << 175 if(p != particle) { << 176 particle = p; << 177 mass = particle->GetPDGMass(); << 178 wokvi->SetupParticle(p); << 179 } << 180 } << 181 124 182 //....oooOO0OOooo........oooOO0OOooo........oo << 125 G4double a0; >> 126 G4double lowKEnergy; >> 127 G4double highKEnergy; >> 128 G4double alpha2; >> 129 G4double faclim; >> 130 >> 131 G4int nbins; >> 132 G4int nmax; >> 133 G4int index[100]; 183 134 184 inline void G4eCoulombScatteringModel::SetReco << 135 G4bool buildTable; 185 { << 136 G4bool isInitialised; 186 recoilThreshold = eth; << 137 }; 187 } << 188 138 189 //....oooOO0OOooo........oooOO0OOooo........oo << 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 190 140 191 inline void G4eCoulombScatteringModel::SetFixe << 141 inline G4double G4eCoulombScatteringModel::ComputeCrossSectionPerAtom( >> 142 const G4ParticleDefinition* p, >> 143 G4double kinEnergy, >> 144 G4double Z, G4double, >> 145 G4double, G4double) 192 { 146 { 193 fixedCut = val; << 147 G4double x; >> 148 G4bool b; >> 149 if(theCrossSectionTable) { >> 150 x = std::exp((((*theCrossSectionTable)[index[G4int(Z)]])) >> 151 ->GetValue(kinEnergy, b)); >> 152 } else x = CalculateCrossSectionPerAtom(p, kinEnergy, Z); >> 153 >> 154 // G4cout << "G4eCoulombScatteringModel: e= " << kinEnergy >> 155 // << " cs= " << x << G4endl; >> 156 return x; 194 } 157 } 195 158 196 //....oooOO0OOooo........oooOO0OOooo........oo << 159 inline G4double G4eCoulombScatteringModel::ScreeningParameter( 197 << 160 G4double Z, 198 inline G4double G4eCoulombScatteringModel::Get << 161 G4double q2, >> 162 G4double mom2, >> 163 G4double invbeta2) 199 { 164 { 200 return fixedCut; << 165 G4double a = invbeta2*Z*Z*q2*alpha2; >> 166 return std::pow(Z,0.6666667)*a0*(1.13 + 3.76*a)/mom2; 201 } 167 } 202 168 203 //....oooOO0OOooo........oooOO0OOooo........oo 169 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 204 170 205 #endif 171 #endif 206 172