Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // ------------------------------------------- 28 // 29 // GEANT4 Class header file 30 // 31 // 32 // File name: G4PAIPhotModel 33 // 34 // Author: V. Grichine based on G4PAIMo 35 // 36 // Creation date: 07.10.2013 37 // 38 // 39 // Class Description: 40 // 41 // Implementation of PAI photon model of energ 42 // delta-electron or X-ray production by charg 43 // 44 // ------------------------------------------- 45 // 46 47 #ifndef G4PAIPhotModel_h 48 #define G4PAIPhotModel_h 1 49 50 #include <CLHEP/Units/PhysicalConstants.h> 51 52 #include "G4VEmModel.hh" 53 #include "G4VEmFluctuationModel.hh" 54 #include "globals.hh" 55 #include <vector> 56 57 class G4Region; 58 class G4MaterialCutsCouple; 59 class G4ParticleChangeForLoss; 60 61 class G4PAIPhotData; 62 63 class G4PAIPhotModel final : public G4VEmModel 64 { 65 66 public: 67 68 explicit G4PAIPhotModel(const G4ParticleDefi 69 const G4String& nam = "PAI"); 70 71 ~G4PAIPhotModel() final; 72 73 void Initialise(const G4ParticleDefinition*, 74 75 void InitialiseLocal(const G4ParticleDefinit 76 G4VEmModel* masterModel 77 78 G4double MinEnergyCut(const G4ParticleDefini 79 const G4MaterialCutsCo 80 81 G4double ComputeDEDXPerVolume(const G4Materi 82 const G4ParticleDefinition 83 G4double kineticEnergy, 84 G4double cutEnergy) final; 85 86 G4double CrossSectionPerVolume(const G4Mater 87 const G4ParticleDefinition*, 88 G4double kineticEnergy, 89 G4double cutEnergy, 90 G4double maxEnergy) final; 91 92 void SampleSecondaries(std::vector<G4Dynamic 93 const G4MaterialCutsCouple*, 94 const G4DynamicParticle*, 95 G4double tmin, 96 G4double maxEnergy) final; 97 98 G4double SampleFluctuations(const G4Material 99 const G4DynamicParticle*, 100 const G4double, const G4double, 101 const G4double, 102 103 G4double Dispersion(const G4Material*, const 104 const G4double, const G4double, cons 105 106 void DefineForRegion(const G4Region* r) fina 107 108 inline G4PAIPhotData* GetPAIPhotData(); 109 110 inline const std::vector<const G4MaterialCut 111 112 inline G4double ComputeMaxEnergy(G4double sc 113 114 inline void SetVerboseLevel(G4int verbose); 115 116 // hide assignment operator 117 G4PAIPhotModel & operator=(const G4PAIPhotM 118 G4PAIPhotModel(const G4PAIPhotModel&) = del 119 120 protected: 121 122 G4double MaxSecondaryEnergy(const G4Particle 123 G4double kinEnergy) final; 124 125 private: 126 127 inline G4int FindCoupleIndex(const G4Materia 128 129 inline void SetParticle(const G4ParticleDefi 130 131 G4int fVerbose; 132 133 G4PAIPhotData* fModelData; 134 135 std::vector<const G4MaterialCutsCouple*> fMa 136 std::vector<const G4Region*> fPAIRegion 137 138 const G4ParticleDefinition* fParticle; 139 const G4ParticleDefinition* fElectron; 140 const G4ParticleDefinition* fPositron; 141 G4ParticleChangeForLoss* fParticleChange; 142 143 G4double fMass; 144 G4double fRatio; 145 G4double fChargeSquare; 146 G4double fLowestTcut; 147 }; 148 149 inline G4PAIPhotData* G4PAIPhotModel::GetPAIPh 150 { 151 return fModelData; 152 } 153 154 inline const std::vector<const G4MaterialCutsC 155 G4PAIPhotModel::GetVectorOfCouples() 156 { 157 return fMaterialCutsCoupleVector; 158 } 159 160 inline G4double G4PAIPhotModel::ComputeMaxEner 161 { 162 return MaxSecondaryEnergy(fParticle, scaledE 163 } 164 165 inline void G4PAIPhotModel::SetVerboseLevel(G4 166 { 167 fVerbose=verbose; 168 } 169 170 inline G4int G4PAIPhotModel::FindCoupleIndex(c 171 { 172 G4int idx = -1; 173 G4int jMatMax = (G4int)fMaterialCutsCoupleVe 174 for(G4int jMat = 0;jMat < jMatMax; ++jMat) { 175 if(couple == fMaterialCutsCoupleVector[jMa 176 idx = jMat; 177 break; 178 } 179 } 180 return idx; 181 } 182 183 inline void G4PAIPhotModel::SetParticle(const 184 { 185 if(fParticle != p) { 186 fParticle = p; 187 fMass = fParticle->GetPDGMass(); 188 fRatio = CLHEP::proton_mass_c2/fMass; 189 G4double q = fParticle->GetPDGCharge()/CLH 190 fChargeSquare = q*q; 191 } 192 } 193 194 #endif 195