Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
ExExChHadronElasticPhysics Class Reference

#include <ExExChHadronElasticPhysics.hh>

Inheritance diagram for ExExChHadronElasticPhysics:
Collaboration diagram for ExExChHadronElasticPhysics:

Public Member Functions

 ExExChHadronElasticPhysics (G4int ver=0)
 
virtual ~ExExChHadronElasticPhysics ()
 
virtual void ConstructParticle ()
 
virtual void ConstructProcess ()
 
G4HadronElasticGetNeutronModel ()
 
G4HadronicProcessGetNeutronProcess ()
 
- Public Member Functions inherited from G4VPhysicsConstructor
 G4VPhysicsConstructor (const G4String &="")
 
 G4VPhysicsConstructor (const G4String &name, G4int physics_type)
 
virtual ~G4VPhysicsConstructor ()
 
void SetPhysicsName (const G4String &="")
 
const G4StringGetPhysicsName () const
 
void SetPhysicsType (G4int)
 
G4int GetPhysicsType () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
G4int GetInstanceID () const
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManagerGetSubInstanceManager ()
 
- Protected Member Functions inherited from G4VPhysicsConstructor
G4bool RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
 
G4ParticleTable::G4PTblDicIteratorGetParticleIterator () const
 
- Protected Attributes inherited from G4VPhysicsConstructor
G4int verboseLevel
 
G4String namePhysics
 
G4int typePhysics
 
G4ParticleTabletheParticleTable
 
G4int g4vpcInstanceID
 
- Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager subInstanceManager
 

Detailed Description

Definition at line 39 of file ExExChHadronElasticPhysics.hh.

Constructor & Destructor Documentation

ExExChHadronElasticPhysics::ExExChHadronElasticPhysics ( G4int  ver = 0)

Definition at line 76 of file ExExChHadronElasticPhysics.cc.

77 : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver)
78 {
79  if(verbose > 1) {
80  G4cout << "### ExExChHadronElasticPhysics: " << GetPhysicsName()
81  << G4endl;
82  }
83 }
G4GLOB_DLL std::ostream G4cout
const G4String & GetPhysicsName() const
G4VPhysicsConstructor(const G4String &="")
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

ExExChHadronElasticPhysics::~ExExChHadronElasticPhysics ( )
virtual

Definition at line 87 of file ExExChHadronElasticPhysics.cc.

88 {}

Member Function Documentation

void ExExChHadronElasticPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 92 of file ExExChHadronElasticPhysics.cc.

93 {
94  G4MesonConstructor pMesonConstructor;
95  pMesonConstructor.ConstructParticle();
96 
97  G4BaryonConstructor pBaryonConstructor;
98  pBaryonConstructor.ConstructParticle();
99 
100  G4IonConstructor pConstructor;
101  pConstructor.ConstructParticle();
102 }
static void ConstructParticle()
static void ConstructParticle()
static void ConstructParticle()

Here is the call graph for this function:

void ExExChHadronElasticPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 106 of file ExExChHadronElasticPhysics.cc.

107 {
108  if(wasActivated) { return; }
109  wasActivated = true;
110 
111  const G4double elimitPi = 1.0*GeV;
112  const G4double elimitAntiNuc = 100.*MeV;
113  const G4double delta = 0.1*MeV;
114  if(verbose > 1) {
115  G4cout << "### HadronElasticPhysics::ConstructProcess: Elimit for pi "
116  << elimitPi/GeV << " GeV" << G4endl;
117  G4cout << " for anti-neuclei "
118  << elimitAntiNuc/GeV << " GeV" << G4endl;
119  }
120 
121  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
122  anuc->SetMinEnergy(elimitAntiNuc);
123  G4CrossSectionElastic* anucxs =
125 
126  G4HadronElastic* lhep0 = new G4HadronElastic();
127  G4HadronElastic* lhep1 = new G4HadronElastic();
128  G4HadronElastic* lhep2 = new G4HadronElastic();
129  lhep1->SetMaxEnergy(elimitPi+delta);
130  lhep2->SetMaxEnergy(elimitAntiNuc+delta);
131 
133  neutronModel = new G4ChipsElasticModel();
134 
136  he->SetMinEnergy(elimitPi);
137 
139  particleIterator->reset();
140  while( (*particleIterator)() )
141  {
142  G4ParticleDefinition* particle = particleIterator->value();
143  G4ProcessManager* pmanager = particle->GetProcessManager();
144  G4String pname = particle->GetParticleName();
145  if(pname == "anti_lambda" ||
146  pname == "anti_neutron" ||
147  pname == "anti_omega-" ||
148  pname == "anti_sigma-" ||
149  pname == "anti_sigma+" ||
150  pname == "anti_xi-" ||
151  pname == "anti_xi0" ||
152  pname == "lambda" ||
153  pname == "omega-" ||
154  pname == "sigma-" ||
155  pname == "sigma+" ||
156  pname == "xi-" ||
157  pname == "alpha" ||
158  pname == "deuteron" ||
159  pname == "triton"
160  ) {
161 
163  hel->RegisterMe(lhep0);
164 
165  XWrapperDiscreteProcess* hel_wrapper =
167  hel_wrapper->RegisterProcess(hel,1);
168  pmanager->AddDiscreteProcess(hel_wrapper);
169 
170  if(verbose > 1) {
171  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
172  << " added for " << particle->GetParticleName() << G4endl;
173  }
174 
175  } else if(pname == "proton") {
176 
178 
180  Instance()->GetCrossSectionDataSet(
182 
183  hel->RegisterMe(chipsp);
184 
185  XWrapperDiscreteProcess* hel_wrapper =
187  hel_wrapper->RegisterProcess(hel,1);
188  pmanager->AddDiscreteProcess(hel_wrapper);
189  if(verbose > 1) {
190  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
191  << " added for " << particle->GetParticleName() << G4endl;
192  }
193 
194  } else if(pname == "neutron") {
195 
196  neutronProcess = new G4HadronElasticProcess();
197  //neutronProcess->AddDataSet(new G4BGGNucleonElasticXS(particle));
198  neutronProcess->AddDataSet(
200  Instance()->GetCrossSectionDataSet(
202  neutronProcess->RegisterMe(neutronModel);
203  pmanager->AddDiscreteProcess(neutronProcess);
204  if(verbose > 1) {
205  G4cout << "### HadronElasticPhysics: "
206  << neutronProcess->GetProcessName()
207  << " added for " << particle->GetParticleName() << G4endl;
208  }
209 
210  } else if (pname == "pi+" || pname == "pi-") {
211 
213  hel->AddDataSet(new G4BGGPionElasticXS(particle));
214  hel->RegisterMe(lhep1);
215  hel->RegisterMe(he);
216  XWrapperDiscreteProcess* hel_wrapper =
218  hel_wrapper->RegisterProcess(hel,1);
219  pmanager->AddDiscreteProcess(hel_wrapper);
220 
221  if(verbose > 1) {
222  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
223  << " added for " << particle->GetParticleName() << G4endl;
224  }
225 
226  } else if(pname == "kaon-" ||
227  pname == "kaon+" ||
228  pname == "kaon0S" ||
229  pname == "kaon0L"
230  ) {
231 
233  hel->RegisterMe(lhep0);
234  XWrapperDiscreteProcess* hel_wrapper =
236  hel_wrapper->RegisterProcess(hel,1);
237  pmanager->AddDiscreteProcess(hel_wrapper);
238  if(verbose > 1) {
239  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
240  << " added for " << particle->GetParticleName() << G4endl;
241  }
242 
243  } else if(
244  pname == "anti_proton" ||
245  pname == "anti_alpha" ||
246  pname == "anti_deuteron" ||
247  pname == "anti_triton" ||
248  pname == "anti_He3" ) {
249 
251  hel->AddDataSet(anucxs);
252  hel->RegisterMe(lhep2);
253  hel->RegisterMe(anuc);
254  XWrapperDiscreteProcess* hel_wrapper =
256  hel_wrapper->RegisterProcess(hel,1);
257  pmanager->AddDiscreteProcess(hel_wrapper);
258  }
259  }
260 }
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4ComponentAntiNuclNuclearXS * GetComponentCrossSection()
const G4String & GetParticleName() const
void RegisterMe(G4HadronicInteraction *a)
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
void SetMinEnergy(G4double anEnergy)
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
G4GLOB_DLL std::ostream G4cout
static const char * Default_Name()
const G4String & GetProcessName() const
Definition: G4VProcess.hh:408
string pname
Definition: eplot.py:33
G4ProcessManager * GetProcessManager() const
static constexpr double GeV
Definition: G4SIunits.hh:217
void SetMaxEnergy(const G4double anEnergy)
std::vector< G4InuclElementaryParticle >::iterator particleIterator
Definition: G4BigBanger.cc:65
static const char * Default_Name()
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
void RegisterProcess(G4VDiscreteProcess *)
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4HadronElastic * ExExChHadronElasticPhysics::GetNeutronModel ( )
inline

Definition at line 75 of file ExExChHadronElasticPhysics.hh.

76 {
77  return neutronModel;
78 }
G4HadronicProcess * ExExChHadronElasticPhysics::GetNeutronProcess ( )
inline

Definition at line 80 of file ExExChHadronElasticPhysics.hh.

81 {
82  return neutronProcess;
83 }

The documentation for this class was generated from the following files: