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 file 30 // 31 // 32 // File name: G4ePairProduction 33 // 34 // Author: Vladimir Ivanchenko 35 // 36 // Creation date: 17.03.2016 37 // 38 // Modifications: 39 // 40 // ------------------------------------------- 41 // 42 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oo 44 45 #include "G4ePairProduction.hh" 46 #include "G4SystemOfUnits.hh" 47 #include "G4Electron.hh" 48 #include "G4Positron.hh" 49 #include "G4VEmModel.hh" 50 #include "G4VEmFluctuationModel.hh" 51 #include "G4MuPairProductionModel.hh" 52 #include "G4ElementData.hh" 53 #include "G4EmParameters.hh" 54 55 //....oooOO0OOooo........oooOO0OOooo........oo 56 57 G4ePairProduction::G4ePairProduction(const G4S 58 : G4VEnergyLossProcess(name), 59 lowestKinEnergy(100.*CLHEP::MeV) 60 { 61 SetProcessSubType(fPairProdByCharged); 62 SetSecondaryParticle(G4Positron::Positron()) 63 SetIonisation(false); 64 SetSpline(false); 65 } 66 67 //....oooOO0OOooo........oooOO0OOooo........oo 68 69 G4bool G4ePairProduction::IsApplicable(const G 70 { 71 return (G4Electron::Electron() == &p || G4Po 72 } 73 74 //....oooOO0OOooo........oooOO0OOooo........oo 75 76 G4double G4ePairProduction::MinPrimaryEnergy(c 77 const G4Material*, 78 G4double) 79 { 80 return lowestKinEnergy; 81 } 82 83 //....oooOO0OOooo........oooOO0OOooo........oo 84 85 void G4ePairProduction::InitialiseEnergyLossPr 86 const G4ParticleDefin 87 const G4ParticleDefinition*) 88 { 89 if (!isInitialised) { 90 isInitialised = true; 91 92 theParticle = part; 93 94 G4MuPairProductionModel* mod = new G4MuPai 95 SetEmModel(mod); 96 97 lowestKinEnergy = std::max(lowestKinEnergy 98 mod->SetLowestKineticEnergy(lowestKinEnerg 99 100 G4VEmFluctuationModel* fm = nullptr; 101 G4EmParameters* param = G4EmParameters::In 102 mod->SetLowEnergyLimit(param->MinKinEnergy 103 mod->SetHighEnergyLimit(param->MaxKinEnerg 104 AddEmModel(1, mod, fm); 105 } 106 } 107 108 //....oooOO0OOooo........oooOO0OOooo........oo 109 110 void G4ePairProduction::StreamProcessInfo(std: 111 { 112 G4ElementData* ed = EmModel(0)->GetElementDa 113 if(ed) { 114 for(G4int Z=1; Z<93; ++Z) { 115 G4Physics2DVector* pv = ed->GetElement2D 116 if(pv) { 117 out << " Sampling table " << pv-> 118 << "x" << pv->GetLengthX() << "; from " 119 << G4Exp(pv->GetY(0))/GeV << " GeV to " 120 << G4Exp(pv->GetY(pv->GetLengthY()-1))/T 121 << " TeV " << G4endl; 122 break; 123 } 124 } 125 } 126 } 127 128 //....oooOO0OOooo........oooOO0OOooo........oo 129 130 void G4ePairProduction::ProcessDescription(std 131 { 132 out << " Electron-positron pair production 133 G4VEnergyLossProcess::ProcessDescription(out 134 } 135 136 //....oooOO0OOooo........oooOO0OOooo........oo 137