Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/highenergy/src/G4eeToHadronsMultiModel.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/highenergy/src/G4eeToHadronsMultiModel.cc (Version 11.3.0) and /processes/electromagnetic/highenergy/src/G4eeToHadronsMultiModel.cc (Version 9.4.p2)


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