Geant4  10.02.p01
G4HadronElasticPhysics.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4HadronElasticPhysics.cc 95780 2016-02-24 08:34:54Z gcosmo $
27 //
28 //---------------------------------------------------------------------------
29 //
30 // ClassName: G4HadronElasticPhysics
31 //
32 // Author: 23 November 2006 V. Ivanchenko
33 //
34 // Modified:
35 // 21.03.2007 V.Ivanchenko Use G4BGGNucleonElasticXS and G4BGGPionElasticXS;
36 // Reduce thresholds for HE and Q-models to zero
37 // 03.06.2010 V.Ivanchenko cleanup constructors and ConstructProcess method
38 // 29.07.2010 V.Ivanchenko rename this class from G4HadronHElasticPhysics to
39 // G4HadronElasticPhysics, old version of the class
40 // is renamed to G4HadronElasticPhysics93
41 //
42 //----------------------------------------------------------------------------
43 //
44 // CHIPS for sampling scattering for p and n
45 // Glauber model for samplimg of high energy pi+- (E > 1GeV)
46 // LHEP sampling model for the other particle
47 // BBG cross sections for p, n and pi+-
48 // LHEP cross sections for other particles
49 
51 
52 #include "G4SystemOfUnits.hh"
53 #include "G4ParticleDefinition.hh"
54 #include "G4ProcessManager.hh"
55 
56 #include "G4MesonConstructor.hh"
57 #include "G4BaryonConstructor.hh"
58 #include "G4IonConstructor.hh"
59 
61 #include "G4HadronElastic.hh"
62 #include "G4ChipsElasticModel.hh"
64 #include "G4AntiNuclElastic.hh"
65 
67 
68 #include "G4BGGNucleonElasticXS.hh"
69 #include "G4BGGPionElasticXS.hh"
70 
72 
74 
75 #include "G4NeutronElasticXS.hh"
76 
77 #include "G4CrossSectionElastic.hh"
78 
79 // factory
81 //
83 //
87 
89  : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver)
90 {
91  if(verbose > 1) {
92  G4cout << "### G4HadronElasticPhysics: " << GetPhysicsName()
93  << G4endl;
94  }
95 }
96 
98 {}
99 
101 {
102  // G4cout << "G4HadronElasticPhysics::ConstructParticle" << G4endl;
103  G4MesonConstructor pMesonConstructor;
104  pMesonConstructor.ConstructParticle();
105 
106  G4BaryonConstructor pBaryonConstructor;
107  pBaryonConstructor.ConstructParticle();
108 
109  G4IonConstructor pConstructor;
110  pConstructor.ConstructParticle();
111 }
112 
114 {
115  if(wasActivated) { return; }
116  wasActivated = true;
117 
118  const G4double elimitPi = 1.0*GeV;
119  const G4double elimitAntiNuc = 100.*MeV;
120  const G4double delta = 0.1*MeV;
121  if(verbose > 1) {
122  G4cout << "### HadronElasticPhysics::ConstructProcess: Elimit for pi "
123  << elimitPi/GeV << " GeV" << G4endl;
124  G4cout << " for anti-neuclei "
125  << elimitAntiNuc/GeV << " GeV" << G4endl;
126  }
127 
128  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
129  anuc->SetMinEnergy(elimitAntiNuc);
130  G4CrossSectionElastic* anucxs =
132 
133  G4HadronElastic* lhep0 = new G4HadronElastic();
134  G4HadronElastic* lhep1 = new G4HadronElastic();
135  G4HadronElastic* lhep2 = new G4HadronElastic();
136  lhep1->SetMaxEnergy(elimitPi+delta);
137  lhep2->SetMaxEnergy(elimitAntiNuc+delta);
138 
141 
143  he->SetMinEnergy(elimitPi);
144 
145  aParticleIterator->reset();
146  while( (*aParticleIterator)() )
147  {
148  G4ParticleDefinition* particle = aParticleIterator->value();
149  G4ProcessManager* pmanager = particle->GetProcessManager();
150  G4String pname = particle->GetParticleName();
151  if(pname == "anti_lambda" ||
152  pname == "anti_neutron" ||
153  pname == "anti_omega-" ||
154  pname == "anti_sigma-" ||
155  pname == "anti_sigma+" ||
156  pname == "anti_xi-" ||
157  pname == "anti_xi0" ||
158  pname == "lambda" ||
159  pname == "omega-" ||
160  pname == "sigma-" ||
161  pname == "sigma+" ||
162  pname == "xi-" ||
163  pname == "alpha" ||
164  pname == "deuteron" ||
165  pname == "triton"
166  ) {
167 
169  hel->RegisterMe(lhep0);
170  pmanager->AddDiscreteProcess(hel);
171  if(verbose > 1) {
172  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
173  << " added for " << particle->GetParticleName() << G4endl;
174  }
175 
176  } else if(pname == "He3") {
178  G4VCrossSectionDataSet* theComponentGGNuclNuclData =
180  hel->AddDataSet(theComponentGGNuclNuclData);
181  hel->RegisterMe(lhep0);
182  pmanager->AddDiscreteProcess(hel);
183  if(verbose > 1) {
184  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
185  << " added for " << particle->GetParticleName() << G4endl;
186  }
187 
188  } else if(pname == "proton") {
189 
191  //hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
192 
193  // hel->AddDataSet(new G4ChipsProtonElasticXS());
195 
196  hel->RegisterMe(chipsp);
197  pmanager->AddDiscreteProcess(hel);
198  if(verbose > 1) {
199  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
200  << " added for " << particle->GetParticleName() << G4endl;
201  }
202 
203  } else if(pname == "neutron") {
204 
209  if(verbose > 1) {
210  G4cout << "### HadronElasticPhysics: "
212  << " added for " << particle->GetParticleName() << G4endl;
213  }
214 
215  } else if (pname == "pi+" || pname == "pi-") {
216 
218  hel->AddDataSet(new G4BGGPionElasticXS(particle));
219  hel->RegisterMe(lhep1);
220  hel->RegisterMe(he);
221  pmanager->AddDiscreteProcess(hel);
222  if(verbose > 1) {
223  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
224  << " added for " << particle->GetParticleName() << G4endl;
225  }
226 
227  } else if(pname == "kaon-" ||
228  pname == "kaon+" ||
229  pname == "kaon0S" ||
230  pname == "kaon0L"
231  ) {
232 
234  hel->RegisterMe(lhep0);
235  pmanager->AddDiscreteProcess(hel);
236  if(verbose > 1) {
237  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
238  << " added for " << particle->GetParticleName() << G4endl;
239  }
240 
241  } else if(
242  pname == "anti_proton" ||
243  pname == "anti_alpha" ||
244  pname == "anti_deuteron" ||
245  pname == "anti_triton" ||
246  pname == "anti_He3" ) {
247 
249  hel->AddDataSet(anucxs);
250  hel->RegisterMe(lhep2);
251  hel->RegisterMe(anuc);
252  pmanager->AddDiscreteProcess(hel);
253  }
254  }
255 }
256 
257 
259 {
260  return neutronModel;
261 }
262 
264 {
265  return neutronProcess;
266 }
267 
268 
static const double MeV
Definition: G4SIunits.hh:211
G4HadronElastic * GetNeutronModel()
G4HadronicProcess * GetNeutronProcess()
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
static void ConstructParticle()
#define G4ThreadLocal
Definition: tls.hh:89
static void ConstructParticle()
G4ProcessManager * GetProcessManager() const
int G4int
Definition: G4Types.hh:78
G4ComponentAntiNuclNuclearXS * GetComponentCrossSection()
const G4String & GetParticleName() const
static void ConstructParticle()
void RegisterMe(G4HadronicInteraction *a)
void SetMinEnergy(G4double anEnergy)
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
#define aParticleIterator
static G4ThreadLocal G4bool wasActivated
G4_DECLARE_PHYSCONSTR_FACTORY(G4HadronElasticPhysics)
static G4CrossSectionDataSetRegistry * Instance()
static const double GeV
Definition: G4SIunits.hh:214
const G4String & GetPhysicsName() const
const G4String & GetProcessName() const
Definition: G4VProcess.hh:408
static G4ThreadLocal G4HadronElastic * neutronModel
static G4ThreadLocal G4HadronicProcess * neutronProcess
static const char * Default_Name()
void SetMaxEnergy(const G4double anEnergy)
static const char * Default_Name()
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76