Geant4  10.02.p03
Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
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 ()
 
G4HadronElastic * GetNeutronModel ()
 
G4HadronicProcess * GetNeutronProcess ()
 
- Public Member Functions inherited from G4VPhysicsConstructor
 G4VPhysicsConstructor (const G4String &="")
 
 G4VPhysicsConstructor (const G4String &name, G4int physics_type)
 
virtual ~G4VPhysicsConstructor ()
 
void SetPhysicsName (const G4String &="")
 
const G4String & GetPhysicsName () const
 
void SetPhysicsType (G4int)
 
G4int GetPhysicsType () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
G4int GetInstanceID () const
 

Private Member Functions

 ExExChHadronElasticPhysics (ExExChHadronElasticPhysics &)
 
ExExChHadronElasticPhysics & operator= (const ExExChHadronElasticPhysics &right)
 

Private Attributes

G4int verbose
 

Static Private Attributes

static G4ThreadLocal G4bool wasActivated = false
 
static G4ThreadLocal G4HadronElastic * neutronModel = 0
 
static G4ThreadLocal G4HadronicProcess * neutronProcess = 0
 

Additional Inherited Members

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

Detailed Description

Definition at line 36 of file ExExChHadronElasticPhysics.hh.

Constructor & Destructor Documentation

◆ ExExChHadronElasticPhysics() [1/2]

ExExChHadronElasticPhysics::ExExChHadronElasticPhysics ( G4int  ver = 0)

Definition at line 72 of file ExExChHadronElasticPhysics.cc.

73 : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver)
74 {
75  if(verbose > 1) {
76  G4cout << "### ExExChHadronElasticPhysics: " << GetPhysicsName()
77  << G4endl;
78  }
79 }
G4VPhysicsConstructor::GetPhysicsName
const G4String & GetPhysicsName() const
Definition: G4VPhysicsConstructor.hh:195
G4cout
G4GLOB_DLL std::ostream G4cout
G4VPhysicsConstructor::G4VPhysicsConstructor
G4VPhysicsConstructor(const G4String &="")
Definition: G4VPhysicsConstructor.cc:48
G4endl
#define G4endl
Definition: G4ios.hh:61
Here is the call graph for this function:

◆ ~ExExChHadronElasticPhysics()

ExExChHadronElasticPhysics::~ExExChHadronElasticPhysics ( )
virtual

Definition at line 81 of file ExExChHadronElasticPhysics.cc.

82 {}

◆ ExExChHadronElasticPhysics() [2/2]

ExExChHadronElasticPhysics::ExExChHadronElasticPhysics ( ExExChHadronElasticPhysics &  )
private

Member Function Documentation

◆ ConstructParticle()

void ExExChHadronElasticPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 84 of file ExExChHadronElasticPhysics.cc.

85 {
86  G4MesonConstructor pMesonConstructor;
87  pMesonConstructor.ConstructParticle();
88 
89  G4BaryonConstructor pBaryonConstructor;
90  pBaryonConstructor.ConstructParticle();
91 
92  G4IonConstructor pConstructor;
93  pConstructor.ConstructParticle();
94 }
G4BaryonConstructor::ConstructParticle
static void ConstructParticle()
Definition: G4BaryonConstructor.cc:102
G4MesonConstructor::ConstructParticle
static void ConstructParticle()
Definition: G4MesonConstructor.cc:80
G4IonConstructor::ConstructParticle
static void ConstructParticle()
Definition: G4IonConstructor.cc:59
Here is the call graph for this function:

◆ ConstructProcess()

void ExExChHadronElasticPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 96 of file ExExChHadronElasticPhysics.cc.

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

◆ GetNeutronModel()

G4HadronElastic * ExExChHadronElasticPhysics::GetNeutronModel ( )
inline

Definition at line 72 of file ExExChHadronElasticPhysics.hh.

73 {
74  return neutronModel;
75 }
ExExChHadronElasticPhysics::neutronModel
static G4ThreadLocal G4HadronElastic * neutronModel
Definition: ExExChHadronElasticPhysics.hh:67

◆ GetNeutronProcess()

G4HadronicProcess * ExExChHadronElasticPhysics::GetNeutronProcess ( )
inline

Definition at line 77 of file ExExChHadronElasticPhysics.hh.

78 {
79  return neutronProcess;
80 }
ExExChHadronElasticPhysics::neutronProcess
static G4ThreadLocal G4HadronicProcess * neutronProcess
Definition: ExExChHadronElasticPhysics.hh:68

◆ operator=()

ExExChHadronElasticPhysics& ExExChHadronElasticPhysics::operator= ( const ExExChHadronElasticPhysics &  right)
private

Member Data Documentation

◆ neutronModel

G4ThreadLocal G4HadronElastic * ExExChHadronElasticPhysics::neutronModel = 0
staticprivate

Definition at line 67 of file ExExChHadronElasticPhysics.hh.

◆ neutronProcess

G4ThreadLocal G4HadronicProcess * ExExChHadronElasticPhysics::neutronProcess = 0
staticprivate

Definition at line 68 of file ExExChHadronElasticPhysics.hh.

◆ verbose

G4int ExExChHadronElasticPhysics::verbose
private

Definition at line 64 of file ExExChHadronElasticPhysics.hh.

◆ wasActivated

G4ThreadLocal G4bool ExExChHadronElasticPhysics::wasActivated = false
staticprivate

Definition at line 65 of file ExExChHadronElasticPhysics.hh.

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