Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | List of all members
G4HadronHElasticPhysics Class Reference

#include <G4HadronHElasticPhysics.hh>

Inheritance diagram for G4HadronHElasticPhysics:
Inheritance graph
[legend]
Collaboration diagram for G4HadronHElasticPhysics:
Collaboration graph
[legend]

Public Member Functions

 G4HadronHElasticPhysics (G4int ver=0, G4bool diffraction=false)
 
virtual ~G4HadronHElasticPhysics ()
 
virtual void ConstructParticle ()
 
virtual void ConstructProcess ()
 
void SetDiffraction (G4bool val)
 
- 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
 

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 47 of file G4HadronHElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronHElasticPhysics::G4HadronHElasticPhysics ( G4int  ver = 0,
G4bool  diffraction = false 
)

Definition at line 90 of file G4HadronHElasticPhysics.cc.

91  : G4VPhysicsConstructor( "hElastic_BEST" ), verbose( ver ),
92  fDiffraction(diffraction)
93 {
94  if ( verbose > 1 ) {
95  G4cout << "### G4HadronHElasticPhysics: " << GetPhysicsName()
96  << " low-mass diffraction: " << fDiffraction << G4endl;
97  }
98 }
G4cout
G4GLOB_DLL std::ostream G4cout
G4VPhysicsConstructor::GetPhysicsName
const G4String & GetPhysicsName() const
Definition: G4VPhysicsConstructor.hh:195
G4VPhysicsConstructor::G4VPhysicsConstructor
G4VPhysicsConstructor(const G4String &="")
Definition: G4VPhysicsConstructor.cc:48
G4endl
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

G4HadronHElasticPhysics::~G4HadronHElasticPhysics ( )
virtual

Definition at line 100 of file G4HadronHElasticPhysics.cc.

100 {}

Member Function Documentation

void G4HadronHElasticPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 103 of file G4HadronHElasticPhysics.cc.

103  {
104  // G4cout << "G4HadronElasticPhysics::ConstructParticle" << G4endl;
105  G4MesonConstructor pMesonConstructor;
106  pMesonConstructor.ConstructParticle();
107 
108  G4BaryonConstructor pBaryonConstructor;
109  pBaryonConstructor.ConstructParticle();
110 
111  // Construct light ions
112  G4IonConstructor pConstructor;
113  pConstructor.ConstructParticle();
114 }
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:

void G4HadronHElasticPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 116 of file G4HadronHElasticPhysics.cc.

116  {
117 
118  const G4double elimitDiffuse = 0.0;
119  const G4double elimitAntiNuc = 100.0*MeV;
120  const G4double delta = 0.1*MeV;
121 
122  if ( verbose > 1 ) {
123  G4cout << "### HadronHElasticPhysics::ConstructProcess: lower energy limit for DiffuseElastic : "
124  << elimitDiffuse/GeV << " GeV" << G4endl
125  << " transition energy for anti-nuclei : "
126  << elimitAntiNuc/GeV << " GeV" << G4endl;
127  }
128 
129  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
130  anuc->SetMinEnergy( elimitAntiNuc );
131  G4CrossSectionElastic* anucxs =
132  new G4CrossSectionElastic( anuc->GetComponentCrossSection() );
133 
134  G4HadronElastic* lhep = new G4HadronElastic();
135  lhep->SetMaxEnergy( elimitAntiNuc + delta );
136 
137  // Three instances of Chips elastic model: one used everywhere,
138  // one used below a energy threshold, and one used only for the
139  // hydrogen element.
140  G4ChipsElasticModel* chips1 = new G4ChipsElasticModel();
141  G4ChipsElasticModel* chips2 = new G4ChipsElasticModel();
142  chips2->SetMaxEnergy( elimitAntiNuc + delta );
143  G4ChipsElasticModel* chipsH = new G4ChipsElasticModel();
144  const G4ElementTable* theElementTable = G4Element::GetElementTable();
145  for ( size_t i_ele = 0; i_ele < theElementTable->size(); i_ele++ ) {
146  G4Element* element = (*theElementTable)[ i_ele ];
147  if ( element->GetZ() > 1.0 ) chipsH->DeActivateFor( element );
148  }
149 
150  G4NuclNuclDiffuseElastic* diffuseNuclNuclElastic = new G4NuclNuclDiffuseElastic();
151  diffuseNuclNuclElastic->SetMinEnergy( elimitDiffuse );
152 
153  G4VCrossSectionDataSet* theComponentGGHadronNucleusData =
154  new G4CrossSectionElastic( new G4ComponentGGHadronNucleusXsc );
155 
156  G4VCrossSectionDataSet* theComponentGGNuclNuclData =
157  new G4CrossSectionElastic( new G4ComponentGGNuclNuclXsc() );
158 
159  G4LMsdGenerator* diffGen = 0;
160  if(fDiffraction) {
161  diffGen = new G4LMsdGenerator("LMsdDiffraction");
162  diffRatio = new G4DiffElasticRatio();
163  G4AutoDelete::Register(diffRatio);
164  }
165 
166  auto myParticleIterator=GetParticleIterator();
167  myParticleIterator->reset();
168  while( (*myParticleIterator)() ) {
169 
170  G4ParticleDefinition* particle = myParticleIterator->value();
171  G4ProcessManager* pmanager = particle->GetProcessManager();
172  G4String pname = particle->GetParticleName();
173 
174  if ( pname == "anti_lambda" ||
175  pname == "anti_sigma-" ||
176  pname == "anti_sigma0" ||
177  pname == "anti_sigma+" ||
178  pname == "anti_xi-" ||
179  pname == "anti_xi0" ||
180  pname == "anti_omega-"
181  ) {
182  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
183  hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4ChipsAntiBaryonElasticXS::Default_Name() ) );
184  hel->RegisterMe( chips1 );
185  pmanager->AddDiscreteProcess( hel );
186  if ( verbose > 1 ) {
187  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
188  << " added for " << particle->GetParticleName() << G4endl;
189  }
190 
191  } else if ( pname == "lambda" ||
192  pname == "sigma-" ||
193  pname == "sigma0" ||
194  pname == "sigma+" ||
195  pname == "xi-" ||
196  pname == "xi0" ||
197  pname == "omega-"
198  ) {
199  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
200  hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4ChipsHyperonElasticXS::Default_Name() ) );
201  hel->RegisterMe( chips1 );
202  pmanager->AddDiscreteProcess( hel );
203  if ( verbose > 1 ) {
204  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
205  << " added for " << particle->GetParticleName() << G4endl;
206  }
207 
208  } else if ( pname == "proton" ) {
209  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
210  hel->AddDataSet( new G4BGGNucleonElasticXS( particle ) );
211  // To preserve reproducibility, a different instance of
212  // G4DiffuseElastic must be used for each particle type.
213  G4DiffuseElastic* protonDiffuseElastic = new G4DiffuseElastic();
214  protonDiffuseElastic->SetMinEnergy( elimitDiffuse );
215  hel->RegisterMe( chipsH ); // Use Chips only for Hydrogen element
216  hel->RegisterMe( protonDiffuseElastic );
217  pmanager->AddDiscreteProcess( hel );
218  if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
219  if ( verbose > 1 ) {
220  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
221  << " added for " << particle->GetParticleName() << G4endl;
222  }
223 
224  } else if ( pname == "neutron" ) {
225  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
226  hel->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronElasticXS::Default_Name()) );
227  // To preserve reproducibility, a different instance of
228  // G4DiffuseElastic must be used for each particle type.
229  G4DiffuseElastic* neutronDiffuseElastic = new G4DiffuseElastic();
230  neutronDiffuseElastic->SetMinEnergy( elimitDiffuse );
231  hel->RegisterMe( chipsH ); // Use Chips only for Hydrogen element
232  hel->RegisterMe( neutronDiffuseElastic );
233  pmanager->AddDiscreteProcess( hel );
234  if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
235  if ( verbose > 1 ) {
236  G4cout << "### HadronElasticPhysics: "
237  << hel->GetProcessName()
238  << " added for " << particle->GetParticleName() << G4endl;
239  }
240 
241  } else if ( pname == "pi-" ) {
242  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
243  hel->AddDataSet( new G4BGGPionElasticXS( particle ) );
244  // To preserve reproducibility, a different instance of
245  // G4DiffuseElastic must be used for each particle type.
246  G4DiffuseElastic* pionMinusDiffuseElastic = new G4DiffuseElastic();
247  pionMinusDiffuseElastic->SetMinEnergy( elimitDiffuse );
248  hel->RegisterMe( chipsH ); // Use Chips only for Hydrogen element
249  hel->RegisterMe( pionMinusDiffuseElastic );
250  pmanager->AddDiscreteProcess( hel );
251  if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
252  if ( verbose > 1 ) {
253  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
254  << " added for " << particle->GetParticleName() << G4endl;
255  }
256 
257  } else if ( pname == "pi+" ) {
258  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
259  hel->AddDataSet( new G4BGGPionElasticXS( particle ) );
260  // To preserve reproducibility, a different instance of
261  // G4DiffuseElastic must be used for each particle type.
262  G4DiffuseElastic* pionPlusDiffuseElastic = new G4DiffuseElastic();
263  hel->RegisterMe( chipsH ); // Use Chips only for Hydrogen element
264  pionPlusDiffuseElastic->SetMinEnergy( elimitDiffuse );
265  hel->RegisterMe( pionPlusDiffuseElastic );
266  pmanager->AddDiscreteProcess( hel );
267  if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
268  if ( verbose > 1 ) {
269  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
270  << " added for " << particle->GetParticleName() << G4endl;
271  }
272 
273  } else if ( pname == "kaon-" ||
274  pname == "kaon+" ||
275  pname == "kaon0S" ||
276  pname == "kaon0L"
277  ) {
278  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
279  //AR-14Aug2017 : Replaced Chips elastic kaon cross sections with
280  // Grichine's Glauber-Gribov ones. In this way, the
281  // total (elastic + inelastic) kaon cross sections
282  // are consistent with the PDG ones.
283  // For the time being, kept Chips elastic as
284  // final-state model.
285  hel->AddDataSet( theComponentGGHadronNucleusData );
286  hel->RegisterMe( chips1 );
287  pmanager->AddDiscreteProcess( hel );
288  if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
289  if ( verbose > 1 ) {
290  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
291  << " added for " << particle->GetParticleName() << G4endl;
292  }
293 
294  } else if (
295  pname == "deuteron" ||
296  pname == "triton" ||
297  pname == "He3" ||
298  pname == "alpha"
299  ) {
300  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
301  hel->AddDataSet( theComponentGGNuclNuclData );
302  // To preserve reproducibility, replace temporarily
303  // G4NuclNuclDiffuseElastic with the Gheisha elastic model.
304  //hel->RegisterMe( diffuseNuclNuclElastic );
305  G4HadronElastic* lhepLightIon = new G4HadronElastic();
306  hel->RegisterMe( lhepLightIon );
307  pmanager->AddDiscreteProcess( hel );
308  if ( verbose > 1 ) {
309  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
310  << " added for " << particle->GetParticleName() << G4endl;
311  }
312 
313  } else if ( pname == "anti_proton" || pname == "anti_neutron" ) {
314  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
315  hel->AddDataSet( anucxs );
316  hel->RegisterMe( chips2 );
317  hel->RegisterMe( anuc );
318  pmanager->AddDiscreteProcess( hel );
319  if ( verbose > 1 ) {
320  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
321  << " added for " << particle->GetParticleName() << G4endl;
322  }
323 
324  } else if ( pname == "anti_deuteron" ||
325  pname == "anti_triton" ||
326  pname == "anti_He3" ||
327  pname == "anti_alpha"
328  ) {
329  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
330  hel->AddDataSet( anucxs );
331  hel->RegisterMe( lhep );
332  hel->RegisterMe( anuc );
333  pmanager->AddDiscreteProcess( hel );
334  if ( verbose > 1 ) {
335  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
336  << " added for " << particle->GetParticleName() << G4endl;
337  }
338 
339  } else if ( pname == "GenericIon" ) {
340  // To preserve reproducibility, disable temporarily
341  // G4NuclNuclDiffuseElastic.
342  //G4HadronElasticProcess* hel = new G4HadronElasticProcess();
343  //hel->AddDataSet( theComponentGGNuclNuclData );
344  //hel->RegisterMe( diffuseNuclNuclElastic );
345  //pmanager->AddDiscreteProcess( hel );
346  //if ( verbose > 1 ) {
347  // G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
348  // << " added for " << particle->GetParticleName() << G4endl;
349  //}
350 
351  }
352 
353  }
354 
355 }
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4HadronElasticProcess::SetDiffraction
void SetDiffraction(G4HadronicInteraction *, G4VCrossSectionRatio *)
Definition: G4HadronElasticProcess.cc:315
G4AntiNuclElastic::GetComponentCrossSection
G4ComponentAntiNuclNuclearXS * GetComponentCrossSection()
Definition: G4AntiNuclElastic.hh:123
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4HadronicProcess::RegisterMe
void RegisterMe(G4HadronicInteraction *a)
Definition: G4HadronicProcess.cc:153
G4ChipsHyperonElasticXS::Default_Name
static const char * Default_Name()
Definition: G4ChipsHyperonElasticXS.hh:54
G4VPhysicsConstructor::GetParticleIterator
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
Definition: G4VPhysicsConstructor.cc:78
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4HadronicProcess::AddDataSet
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
Definition: G4HadronicProcess.hh:111
G4AutoDelete::Register
void Register(T *inst)
Definition: G4AutoDelete.hh:65
G4cout
G4GLOB_DLL std::ostream G4cout
G4CrossSectionDataSetRegistry::Instance
static G4CrossSectionDataSetRegistry * Instance()
Definition: G4CrossSectionDataSetRegistry.cc:79
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:408
eplot.pname
string pname
Definition: eplot.py:33
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:247
GeV
static constexpr double GeV
Definition: G4SIunits.hh:217
G4NeutronElasticXS::Default_Name
static const char * Default_Name()
Definition: G4NeutronElasticXS.hh:69
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4endl
#define G4endl
Definition: G4ios.hh:61
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4ElementTable
std::vector< G4Element * > G4ElementTable
Definition: G4ElementTable.hh:41
G4double
double G4double
Definition: G4Types.hh:76
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:398
G4HadronicInteraction::DeActivateFor
void DeActivateFor(const G4Material *aMaterial)
Definition: G4HadronicInteraction.cc:216

Here is the call graph for this function:

void G4HadronHElasticPhysics::SetDiffraction ( G4bool  val)
inline

Definition at line 64 of file G4HadronHElasticPhysics.hh.

64 {fDiffraction = val;};
The documentation for this class was generated from the following files: