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: G4eeToHadronsMultiModel 33 // 34 // Author: Vladimir Ivanchenko 35 // 36 // Creation date: 02.08.2004 37 // 38 // Modifications: 39 // 08-11-04 Migration to new interface of Stor 40 // 08-04-05 Major optimisation of internal int 41 // 42 43 // 44 // ------------------------------------------- 45 // 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 //....oooOO0OOooo........oooOO0OOooo........oo 48 49 #include "G4eeToHadronsMultiModel.hh" 50 #include "G4PhysicalConstants.hh" 51 #include "G4SystemOfUnits.hh" 52 #include "G4eeToTwoPiModel.hh" 53 #include "G4eeTo3PiModel.hh" 54 #include "G4eeToPGammaModel.hh" 55 #include "G4ee2KNeutralModel.hh" 56 #include "G4ee2KChargedModel.hh" 57 #include "G4eeCrossSections.hh" 58 #include "G4Vee2hadrons.hh" 59 #include "G4Positron.hh" 60 61 //....oooOO0OOooo........oooOO0OOooo........oo 62 63 using namespace std; 64 65 G4eeToHadronsMultiModel::G4eeToHadronsMultiMod 66 const G4String& mname) : G4VEmModel(mname), 67 { 68 maxKineticEnergy = 4.521*CLHEP::GeV; //crre 69 delta = 1.0*CLHEP::MeV; //for bin width 70 } 71 72 //....oooOO0OOooo........oooOO0OOooo........oo 73 74 G4eeToHadronsMultiModel::~G4eeToHadronsMultiMo 75 { 76 delete cross; 77 } 78 79 //....oooOO0OOooo........oooOO0OOooo........oo 80 81 void G4eeToHadronsMultiModel::Initialise(const 82 const G4DataVector& cuts) 83 { 84 if(!isInitialised) { 85 isInitialised = true; 86 87 //G4cout<<"###Initialise in HadronMultiMod 88 89 cross = new G4eeCrossSections(); 90 91 G4eeToTwoPiModel* m2pi = 92 new G4eeToTwoPiModel(cross,maxKineticEne 93 AddEEModel(m2pi,cuts); 94 95 G4eeTo3PiModel* m3pi = 96 new G4eeTo3PiModel(cross,maxKineticEnerg 97 AddEEModel(m3pi,cuts); 98 99 G4ee2KChargedModel* m2kc = 100 new G4ee2KChargedModel(cross,maxKineticE 101 AddEEModel(m2kc,cuts); 102 103 G4ee2KNeutralModel* m2kn = 104 new G4ee2KNeutralModel(cross,maxKineticE 105 AddEEModel(m2kn,cuts); 106 107 G4eeToPGammaModel* mpg1 = 108 new G4eeToPGammaModel(cross,"pi0",maxKin 109 AddEEModel(mpg1,cuts); 110 111 G4eeToPGammaModel* mpg2 = 112 new G4eeToPGammaModel(cross,"eta",maxKin 113 AddEEModel(mpg2,cuts); 114 115 nModels = (G4int)models.size(); 116 117 fParticleChange = GetParticleChangeForGamm 118 } 119 } 120 121 //....oooOO0OOooo........oooOO0OOooo........oo 122 123 void G4eeToHadronsMultiModel::AddEEModel(G4Vee 124 const G4DataVector& cuts) 125 { 126 G4eeToHadronsModel* model = new G4eeToHadron 127 models.push_back(model); 128 G4double elow = mod->LowEnergy(); 129 ekinMin.push_back(elow); 130 if(thKineticEnergy > elow) { thKineticEnergy 131 ekinMax.push_back(mod->HighEnergy()); 132 ekinPeak.push_back(mod->PeakEnergy()); 133 cumSum.push_back(0.0); 134 135 const G4ParticleDefinition* positron = G4Pos 136 model->Initialise(positron,cuts); 137 } 138 139 //....oooOO0OOooo........oooOO0OOooo........oo 140 141 G4double G4eeToHadronsMultiModel::CrossSection 142 const G4Material* mat, 143 const G4ParticleDefinition* p, 144 G4double kineticEnergy, 145 G4double, G4double) 146 { 147 return mat->GetElectronDensity()* 148 ComputeCrossSectionPerElectron(p, kineticE 149 } 150 151 //....oooOO0OOooo........oooOO0OOooo........oo 152 153 G4double G4eeToHadronsMultiModel::ComputeCross 154 const G4 155 G4double kineticEnergy, 156 G4double Z, G4double, 157 G4double, G4double) 158 { 159 return Z*ComputeCrossSectionPerElectron(p, k 160 } 161 162 //....oooOO0OOooo........oooOO0OOooo........oo 163 164 G4double G4eeToHadronsMultiModel::ComputeCross 165 G4double kineticEnergy, 166 G4double, G4double) 167 { 168 G4double res = 0.0; 169 170 G4double energy = LabToCM(kineticEnergy); 171 172 if (energy > thKineticEnergy) { 173 for(G4int i=0; i<nModels; i++) { 174 if(energy >= ekinMin[i] && energy <= eki 175 res += (models[i])->ComputeCrossSectio 176 } 177 cumSum[i] = res; 178 } 179 } 180 return res*csFactor; 181 } 182 183 //....oooOO0OOooo........oooOO0OOooo........oo 184 185 void G4eeToHadronsMultiModel::SampleSecondarie 186 std::vector<G4DynamicParticle*>* ne 187 const G4MaterialCutsCouple* couple, 188 const G4DynamicParticle* dp, 189 G4double, G4double) 190 { 191 G4double kinEnergy = dp->GetKineticEnergy(); 192 G4double energy = LabToCM(kinEnergy); 193 if (energy > thKineticEnergy) { 194 G4double q = cumSum[nModels-1]*G4UniformRa 195 for(G4int i=0; i<nModels; ++i) { 196 if(q <= cumSum[i]) { 197 (models[i])->SampleSecondaries(newp, c 198 if(newp->size() > 0) { 199 fParticleChange->ProposeTrackStatus(fStopA 200 } 201 break; 202 } 203 } 204 } 205 } 206 207 //....oooOO0OOooo........oooOO0OOooo........oo 208 209 void G4eeToHadronsMultiModel::ModelDescription 210 { 211 if(verbose > 0) { 212 G4double e1 = 0.5*thKineticEnergy*thKineti 213 - 2.0*electron_mass_c2; 214 G4double e2 = 0.5*maxKineticEnergy*maxKine 215 - 2.0*electron_mass_c2; 216 outFile << " e+ annihilation into had 217 << e1/GeV << " GeV to " << e2/GeV << " G 218 } 219 } 220 221 //....oooOO0OOooo........oooOO0OOooo........oo 222 223 void G4eeToHadronsMultiModel::SetCrossSecFacto 224 { 225 if(fac > 1.0) { 226 csFactor = fac; 227 if(verbose > 0) { 228 G4cout << "### G4eeToHadronsMultiModel: 229 << "G4eeToHadronsMultiModel is increase 230 << csFactor << " times" << G4endl; 231 } 232 } 233 } 234 235 //....oooOO0OOooo........oooOO0OOooo........oo 236