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

#include <G4HadronElasticPhysics.hh>

Inheritance diagram for G4HadronElasticPhysics:
Collaboration diagram for G4HadronElasticPhysics:

Public Member Functions

 G4HadronElasticPhysics (G4int ver=0)
 
virtual ~G4HadronElasticPhysics ()
 
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 52 of file G4HadronElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronElasticPhysics::G4HadronElasticPhysics ( G4int  ver = 0)

Definition at line 88 of file G4HadronElasticPhysics.cc.

89  : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver)
90 {
91  if(verbose > 1) {
92  G4cout << "### G4HadronElasticPhysics: " << GetPhysicsName()
93  << G4endl;
94  }
95 }
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:

G4HadronElasticPhysics::~G4HadronElasticPhysics ( )
virtual

Definition at line 97 of file G4HadronElasticPhysics.cc.

98 {}

Member Function Documentation

void G4HadronElasticPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 100 of file G4HadronElasticPhysics.cc.

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 }
static void ConstructParticle()
static void ConstructParticle()
static void ConstructParticle()

Here is the call graph for this function:

Here is the caller graph for this function:

void G4HadronElasticPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Reimplemented in HadronElasticPhysicsHP, HadronElasticPhysicsHP, and HadronElasticPhysicsHP.

Definition at line 113 of file G4HadronElasticPhysics.cc.

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 
140  neutronModel = new G4ChipsElasticModel();
141 
143  he->SetMinEnergy(elimitPi);
144 
145  auto myParticleIterator=GetParticleIterator();
146  myParticleIterator->reset();
147  while( (*myParticleIterator)() )
148  {
149  G4ParticleDefinition* particle = myParticleIterator->value();
150  G4ProcessManager* pmanager = particle->GetProcessManager();
151  G4String pname = particle->GetParticleName();
152  if(pname == "anti_lambda" ||
153  pname == "anti_neutron" ||
154  pname == "anti_omega-" ||
155  pname == "anti_sigma-" ||
156  pname == "anti_sigma+" ||
157  pname == "anti_xi-" ||
158  pname == "anti_xi0" ||
159  pname == "lambda" ||
160  pname == "omega-" ||
161  pname == "sigma-" ||
162  pname == "sigma+" ||
163  pname == "xi-" ||
164  pname == "alpha" ||
165  pname == "deuteron" ||
166  pname == "triton"
167  ) {
168 
170  hel->RegisterMe(lhep0);
171  pmanager->AddDiscreteProcess(hel);
172  if(verbose > 1) {
173  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
174  << " added for " << particle->GetParticleName() << G4endl;
175  }
176 
177  } else if(pname == "He3") {
179  G4VCrossSectionDataSet* theComponentGGNuclNuclData =
181  hel->AddDataSet(theComponentGGNuclNuclData);
182  hel->RegisterMe(lhep0);
183  pmanager->AddDiscreteProcess(hel);
184  if(verbose > 1) {
185  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
186  << " added for " << particle->GetParticleName() << G4endl;
187  }
188 
189  } else if(pname == "proton") {
190 
192  //hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
193 
194  // hel->AddDataSet(new G4ChipsProtonElasticXS());
196 
197  hel->RegisterMe(chipsp);
198  pmanager->AddDiscreteProcess(hel);
199  if(verbose > 1) {
200  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
201  << " added for " << particle->GetParticleName() << G4endl;
202  }
203 
204  } else if(pname == "neutron") {
205 
206  neutronProcess = new G4HadronElasticProcess();
207  neutronProcess->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronElasticXS::Default_Name()));
208  neutronProcess->RegisterMe(neutronModel);
209  pmanager->AddDiscreteProcess(neutronProcess);
210  if(verbose > 1) {
211  G4cout << "### HadronElasticPhysics: "
212  << neutronProcess->GetProcessName()
213  << " added for " << particle->GetParticleName() << G4endl;
214  }
215 
216  } else if (pname == "pi+" || pname == "pi-") {
217 
219  hel->AddDataSet(new G4BGGPionElasticXS(particle));
220  hel->RegisterMe(lhep1);
221  hel->RegisterMe(he);
222  pmanager->AddDiscreteProcess(hel);
223  if(verbose > 1) {
224  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
225  << " added for " << particle->GetParticleName() << G4endl;
226  }
227 
228  } else if(pname == "kaon-" ||
229  pname == "kaon+" ||
230  pname == "kaon0S" ||
231  pname == "kaon0L"
232  ) {
233 
235  hel->RegisterMe(lhep0);
236  pmanager->AddDiscreteProcess(hel);
237  if(verbose > 1) {
238  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
239  << " added for " << particle->GetParticleName() << G4endl;
240  }
241 
242  } else if(
243  pname == "anti_proton" ||
244  pname == "anti_alpha" ||
245  pname == "anti_deuteron" ||
246  pname == "anti_triton" ||
247  pname == "anti_He3" ) {
248 
250  hel->AddDataSet(anucxs);
251  hel->RegisterMe(lhep2);
252  hel->RegisterMe(anuc);
253  pmanager->AddDiscreteProcess(hel);
254  }
255  }
256 }
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 G4CrossSectionDataSetRegistry * Instance()
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
static const char * Default_Name()
void SetMaxEnergy(const G4double anEnergy)
static const char * Default_Name()
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

G4HadronElastic * G4HadronElasticPhysics::GetNeutronModel ( )

Definition at line 259 of file G4HadronElasticPhysics.cc.

260 {
261  return neutronModel;
262 }

Here is the caller graph for this function:

G4HadronicProcess * G4HadronElasticPhysics::GetNeutronProcess ( )

Definition at line 264 of file G4HadronElasticPhysics.cc.

265 {
266  return neutronProcess;
267 }

Here is the caller graph for this function:


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