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 // 28 // GEANT4 Class header file 29 // GEANT4 Class header file 29 // 30 // 30 // 31 // 31 // File name: G4LindhardSorensenIonModel 32 // File name: G4LindhardSorensenIonModel 32 // 33 // 33 // Author: Alexander Bagulya & Vladimir 34 // Author: Alexander Bagulya & Vladimir Ivanchenko 34 // 35 // 35 // Creation date: 16.04.2018 36 // Creation date: 16.04.2018 36 // 37 // 37 // 38 // 38 // Class Description: 39 // Class Description: 39 // 40 // 40 // Implementation of ion ionisation energy los 41 // Implementation of ion ionisation energy loss and delta-electron 41 // production by heavy charged particles accor 42 // production by heavy charged particles according to 42 // J. Lindhard & A.H. Sorensen, Phys. Rev. A 5 43 // J. Lindhard & A.H. Sorensen, Phys. Rev. A 53 (1996) 2443-2455 43 44 44 // ------------------------------------------- 45 // ------------------------------------------------------------------- 45 // 46 // 46 47 47 #ifndef G4LindhardSorensenIonModel_h 48 #ifndef G4LindhardSorensenIonModel_h 48 #define G4LindhardSorensenIonModel_h 1 49 #define G4LindhardSorensenIonModel_h 1 49 50 50 #include <vector> << 51 #include <CLHEP/Units/SystemOfUnits.h> 51 52 52 #include "G4VEmModel.hh" 53 #include "G4VEmModel.hh" 53 #include "G4NistManager.hh" 54 #include "G4NistManager.hh" 54 55 55 class G4EmCorrections; 56 class G4EmCorrections; 56 class G4ParticleChangeForLoss; 57 class G4ParticleChangeForLoss; 57 class G4LindhardSorensenData; 58 class G4LindhardSorensenData; 58 class G4BraggModel; << 59 class G4BetheBlochModel; << 60 class G4IonICRU73Data; << 61 59 62 class G4LindhardSorensenIonModel : public G4VE 60 class G4LindhardSorensenIonModel : public G4VEmModel 63 { 61 { >> 62 64 public: 63 public: 65 64 66 explicit G4LindhardSorensenIonModel(const G4 65 explicit G4LindhardSorensenIonModel(const G4ParticleDefinition* p = nullptr, 67 const G4String& nam = "LindhardS 66 const G4String& nam = "LindhardSorensen"); 68 67 69 ~G4LindhardSorensenIonModel() override; << 68 virtual ~G4LindhardSorensenIonModel(); 70 69 71 void Initialise(const G4ParticleDefinition*, << 70 virtual void Initialise(const G4ParticleDefinition*, >> 71 const G4DataVector&) override; 72 72 73 G4double MinEnergyCut(const G4ParticleDefini << 73 virtual G4double MinEnergyCut(const G4ParticleDefinition*, 74 const G4MaterialCutsCo << 74 const G4MaterialCutsCouple* couple) override; 75 75 76 G4double ComputeCrossSectionPerElectron( << 76 virtual G4double ComputeCrossSectionPerElectron( 77 const G4ParticleDefinitio << 77 const G4ParticleDefinition*, 78 G4double kineticEnergy, 78 G4double kineticEnergy, 79 G4double cutEnergy, 79 G4double cutEnergy, 80 G4double maxEnergy); 80 G4double maxEnergy); 81 81 82 G4double ComputeCrossSectionPerAtom( << 82 virtual G4double ComputeCrossSectionPerAtom( 83 const G4ParticleDefinition*, 83 const G4ParticleDefinition*, 84 G4double kineticEnergy, 84 G4double kineticEnergy, 85 G4double Z, G4double A, 85 G4double Z, G4double A, 86 G4double cutEnergy, 86 G4double cutEnergy, 87 G4double maxEnergy) override; 87 G4double maxEnergy) override; 88 88 89 G4double CrossSectionPerVolume(const G4Mater << 89 virtual G4double CrossSectionPerVolume(const G4Material*, 90 const G4ParticleDefinition*, 90 const G4ParticleDefinition*, 91 G4double kineticEnergy, 91 G4double kineticEnergy, 92 G4double cutEnergy, 92 G4double cutEnergy, 93 G4double maxEnergy) override; 93 G4double maxEnergy) override; 94 << 94 95 G4double ComputeDEDXPerVolume(const G4Materi << 95 virtual G4double ComputeDEDXPerVolume(const G4Material*, 96 const G4Partic << 96 const G4ParticleDefinition*, 97 G4double kinet << 97 G4double kineticEnergy, 98 G4double cutEn << 98 G4double cutEnergy) override; 99 << 99 100 G4double GetChargeSquareRatio(const G4Partic << 100 virtual G4double GetChargeSquareRatio(const G4ParticleDefinition* p, 101 const G4Materi << 101 const G4Material* mat, 102 G4double kinet << 102 G4double kineticEnergy) override; 103 << 103 104 G4double GetParticleCharge(const G4ParticleD << 104 virtual G4double GetParticleCharge(const G4ParticleDefinition* p, 105 const G4Material* mat, << 105 const G4Material* mat, 106 G4double kineticE << 106 G4double kineticEnergy) override; 107 << 107 108 void CorrectionsAlongStep(const G4MaterialCu << 108 virtual void CorrectionsAlongStep(const G4MaterialCutsCouple* couple, 109 const G4DynamicPar << 109 const G4DynamicParticle* dp, 110 const G4double& le << 110 G4double& eloss, 111 G4double& eloss) o << 111 G4double&, 112 << 112 G4double length) override; 113 void SampleSecondaries(std::vector<G4Dynamic << 113 114 const G4MaterialCutsC << 114 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*, 115 const G4DynamicParticle*, << 115 const G4MaterialCutsCouple*, 116 G4double tmin, << 116 const G4DynamicParticle*, 117 G4double maxEnergy) overrid << 117 G4double tmin, 118 << 118 G4double maxEnergy) override; 119 // hide assignment operator << 120 G4LindhardSorensenIonModel & operator= << 121 (const G4LindhardSorensenIonModel &right) = << 122 G4LindhardSorensenIonModel(const G4Lindhard << 123 119 124 protected: 120 protected: 125 121 126 G4double MaxSecondaryEnergy(const G4Particle << 122 virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*, 127 G4double kinEner << 123 G4double kinEnergy) override; >> 124 >> 125 inline G4double GetChargeSquareRatio() const; 128 126 129 inline void SetChargeSquareRatio(G4double va 127 inline void SetChargeSquareRatio(G4double val); 130 128 131 private: 129 private: 132 130 133 void SetupParameters(); 131 void SetupParameters(); 134 132 135 void InitialiseLS(); << 136 << 137 G4double ComputeDEDXPerVolumeLS(const G4Mate << 138 const G4Part << 139 G4double kinEnergy, G4double cutEner << 140 << 141 inline void SetParticle(const G4ParticleDefi 133 inline void SetParticle(const G4ParticleDefinition* p); 142 134 143 // static const G4int MAXZION = 93; << 135 // hide assignment operator >> 136 G4LindhardSorensenIonModel & operator= >> 137 (const G4LindhardSorensenIonModel &right) = delete; >> 138 G4LindhardSorensenIonModel(const G4LindhardSorensenIonModel&) = delete; 144 139 145 static G4IonICRU73Data* fIonData; << 146 static G4LindhardSorensenData* lsdata; 140 static G4LindhardSorensenData* lsdata; 147 141 148 const G4ParticleDefinition* particle = nullp << 142 const G4ParticleDefinition* particle; 149 G4ParticleDefinition* theElectron; << 143 G4ParticleDefinition* theElectron; 150 G4EmCorrections* corr; << 144 G4EmCorrections* corr; 151 G4ParticleChangeForLoss* fParticleChange = n << 145 G4ParticleChangeForLoss* fParticleChange; 152 G4NistManager* nist; << 146 G4NistManager* nist; 153 G4BraggModel* fBraggModel; << 147 154 G4BetheBlochModel* fBBModel; << 148 G4int Zin; 155 << 149 G4double mass; 156 G4int Zin = 1; << 150 G4double tlimit; 157 G4double mass = 0.0; << 151 G4double spin; 158 G4double tlimit = DBL_MAX; << 152 G4double magMoment2; 159 G4double spin = 0.0; << 153 G4double chargeSquare; 160 G4double magMoment2 = 0.0; << 154 G4double charge; 161 G4double chargeSquare = 1.0; << 155 G4double ratio; 162 G4double charge = 1.0; << 156 G4double formfact; 163 G4double eRatio = 0.0; << 164 G4double pRatio = 1.0; << 165 G4double formfact = 0.0; << 166 G4double twoln10; 157 G4double twoln10; 167 G4double fElimit; << 168 G4bool isFirst = false; << 169 }; 158 }; 170 159 171 //....oooOO0OOooo........oooOO0OOooo........oo 160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 172 161 173 inline void 162 inline void 174 G4LindhardSorensenIonModel::SetParticle(const 163 G4LindhardSorensenIonModel::SetParticle(const G4ParticleDefinition* p) 175 { 164 { 176 if(particle != p) { 165 if(particle != p) { 177 particle = p; 166 particle = p; 178 SetupParameters(); 167 SetupParameters(); 179 } 168 } >> 169 } >> 170 >> 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 172 >> 173 inline G4double G4LindhardSorensenIonModel::GetChargeSquareRatio() const >> 174 { >> 175 return chargeSquare; 180 } 176 } 181 177 182 //....oooOO0OOooo........oooOO0OOooo........oo 178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 183 179 184 inline void G4LindhardSorensenIonModel::SetCha 180 inline void G4LindhardSorensenIonModel::SetChargeSquareRatio(G4double val) 185 { 181 { 186 chargeSquare = val; 182 chargeSquare = val; 187 } 183 } 188 184 189 //....oooOO0OOooo........oooOO0OOooo........oo 185 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 190 186 191 #endif 187 #endif 192 188