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

#include <G4HadronDElasticPhysics.hh>

Inheritance diagram for G4HadronDElasticPhysics:
Collaboration diagram for G4HadronDElasticPhysics:

Public Member Functions

 G4HadronDElasticPhysics (G4int ver=1)
 
virtual ~G4HadronDElasticPhysics ()
 
virtual void ConstructParticle ()
 
virtual void ConstructProcess ()
 
- 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 45 of file G4HadronDElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronDElasticPhysics::G4HadronDElasticPhysics ( G4int  ver = 1)

Definition at line 89 of file G4HadronDElasticPhysics.cc.

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

G4HadronDElasticPhysics::~G4HadronDElasticPhysics ( )
virtual

Definition at line 98 of file G4HadronDElasticPhysics.cc.

99 {}

Member Function Documentation

void G4HadronDElasticPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 101 of file G4HadronDElasticPhysics.cc.

102 {
103  // G4cout << "G4HadronDElasticPhysics::ConstructParticle" << G4endl;
104  G4MesonConstructor pMesonConstructor;
105  pMesonConstructor.ConstructParticle();
106 
107  G4BaryonConstructor pBaryonConstructor;
108  pBaryonConstructor.ConstructParticle();
109 
110  // Construct light ions
111  G4IonConstructor pConstructor;
112  pConstructor.ConstructParticle();
113 }
static void ConstructParticle()
static void ConstructParticle()
static void ConstructParticle()

Here is the call graph for this function:

void G4HadronDElasticPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 115 of file G4HadronDElasticPhysics.cc.

116 {
117  if(wasActivated) return;
118  wasActivated = true;
119 
120  const G4double elimitAntiNuc = 100.1*MeV;
121  if(verbose > 1) {
122  G4cout << "### HadronDElasticPhysics Construct Processes "
123  << " for anti-neuclei "
124  << elimitAntiNuc/GeV << " GeV" << G4endl;
125  }
126 
127  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
128  anuc->SetMinEnergy(elimitAntiNuc);
129  G4CrossSectionElastic* anucxs =
131 
132  G4HadronElastic* lhep0 = new G4HadronElastic();
133  G4HadronElastic* lhep1 = new G4HadronElastic();
134  lhep1->SetMaxEnergy(10.1*MeV);
135  G4HadronElastic* lhep2 = new G4HadronElastic();
136  lhep2->SetMaxEnergy(elimitAntiNuc);
137 
138  G4DiffuseElastic* model = 0;
139 
140  auto myParticleIterator=GetParticleIterator();
141  myParticleIterator->reset();
142  while( (*myParticleIterator)() )
143  {
144  G4ParticleDefinition* particle = myParticleIterator->value();
145  G4ProcessManager* pmanager = particle->GetProcessManager();
146  G4String pname = particle->GetParticleName();
147  if(pname == "anti_lambda" ||
148  pname == "anti_neutron" ||
149  pname == "anti_omega-" ||
150  pname == "anti_sigma-" ||
151  pname == "anti_sigma+" ||
152  pname == "anti_xi-" ||
153  pname == "anti_xi0" ||
154  pname == "lambda" ||
155  pname == "omega-" ||
156  pname == "sigma-" ||
157  pname == "sigma+" ||
158  pname == "xi-" ||
159  pname == "alpha" ||
160  pname == "deuteron" ||
161  pname == "triton"
162  ) {
163 
165  hel->RegisterMe(lhep0);
166  pmanager->AddDiscreteProcess(hel);
167  if(verbose > 1) {
168  G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
169  << " added for " << particle->GetParticleName() << G4endl;
170  }
171 
172  } else if(pname == "proton") {
173 
175  hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
176  //hel->AddDataSet(new G4CHIPSElasticXS());
177  model = new G4DiffuseElastic();
178  hel->RegisterMe(lhep1);
179  hel->RegisterMe(model);
180  pmanager->AddDiscreteProcess(hel);
181  if(verbose > 1) {
182  G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
183  << " added for " << particle->GetParticleName() << G4endl;
184  }
185 
186  } else if(pname == "neutron") {
187 
190  model = new G4DiffuseElastic();
191  hel->RegisterMe(lhep1);
192  hel->RegisterMe(model);
193  pmanager->AddDiscreteProcess(hel);
194  if(verbose > 1) {
195  G4cout << "### HadronDElasticPhysics: "
196  << hel->GetProcessName()
197  << " added for " << particle->GetParticleName() << G4endl;
198  }
199 
200  } else if (pname == "pi+" || pname == "pi-") {
201 
203  hel->AddDataSet(new G4BGGPionElasticXS(particle));
204  model = new G4DiffuseElastic();
205  hel->RegisterMe(lhep1);
206  hel->RegisterMe(model);
207  pmanager->AddDiscreteProcess(hel);
208  if(verbose > 1) {
209  G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
210  << " added for " << particle->GetParticleName() << G4endl;
211  }
212 
213  } else if(pname == "kaon-") {
214 
217  model = new G4DiffuseElastic();
218  hel->RegisterMe(lhep1);
219  hel->RegisterMe(model);
220  pmanager->AddDiscreteProcess(hel);
221  if(verbose > 1) {
222  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
223  << " added for " << particle->GetParticleName() << G4endl;
224  }
225  } else if(pname == "kaon+") {
226 
229  model = new G4DiffuseElastic();
230  hel->RegisterMe(lhep1);
231  hel->RegisterMe(model);
232  pmanager->AddDiscreteProcess(hel);
233  if(verbose > 1) {
234  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
235  << " added for " << particle->GetParticleName() << G4endl;
236  }
237  } else if(pname == "kaon0S" ||
238  pname == "kaon0L"
239  ) {
240 
243  model = new G4DiffuseElastic();
244  hel->RegisterMe(lhep1);
245  hel->RegisterMe(model);
246  pmanager->AddDiscreteProcess(hel);
247  if(verbose > 1) {
248  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
249  << " added for " << particle->GetParticleName() << G4endl;
250  }
251  } else if(
252  pname == "anti_proton" ||
253  pname == "anti_alpha" ||
254  pname == "anti_deuteron" ||
255  pname == "anti_triton" ||
256  pname == "anti_He3" ) {
257 
259  hel->AddDataSet(anucxs);
260  hel->RegisterMe(lhep2);
261  hel->RegisterMe(anuc);
262  pmanager->AddDiscreteProcess(hel);
263  }
264  }
265  if(verbose > 1) {
266  G4cout << "### HadronDElasticPhysics Construct Processes " << G4endl;
267  }
268 }
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()
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
static const char * Default_Name()
void SetMaxEnergy(const G4double anEnergy)
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
static const char * Default_Name()
double G4double
Definition: G4Types.hh:76
const XML_Char XML_Content * model
Definition: expat.h:151

Here is the call graph for this function:


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