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

#include <G4HadronInelasticQBBC.hh>

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

Public Member Functions

 G4HadronInelasticQBBC (G4int ver=1)
 
 G4HadronInelasticQBBC (const G4String &name, G4int ver=1, G4bool ftf=false, G4bool bert=false, G4bool chips=false, G4bool hp=false, G4bool glauber=false)
 
virtual ~G4HadronInelasticQBBC ()
 
virtual void ConstructProcess ()
 
- Public Member Functions inherited from G4VHadronPhysics
 G4VHadronPhysics (const G4String &name="hInelastic", G4int verbose=0)
 
virtual ~G4VHadronPhysics ()
 
virtual void ConstructParticle ()
 
G4HadronicInteraction * BuildModel (G4VHadronModelBuilder *, G4double emin, G4double emax)
 
G4HadronicInteraction * NewModel (G4HadronicInteraction *, G4double emin, G4double emax)
 
void AddInelasticCrossSection (const G4String &, G4VCrossSectionDataSet *)
 
void AddInelasticCrossSection (const G4ParticleDefinition *, G4VCrossSectionDataSet *)
 
void AddElasticCrossSection (const G4String &, G4VCrossSectionDataSet *)
 
void AddElasticCrossSection (const G4ParticleDefinition *, G4VCrossSectionDataSet *)
 
void AddCaptureCrossSection (G4VCrossSectionDataSet *)
 
void AddFissionCrossSection (G4VCrossSectionDataSet *)
 
- 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 G4VHadronPhysics
G4HadronicProcess * FindInelasticProcess (const G4String &)
 
G4HadronicProcess * FindInelasticProcess (const G4ParticleDefinition *)
 
G4HadronicProcess * FindElasticProcess (const G4String &)
 
G4HadronicProcess * FindElasticProcess (const G4ParticleDefinition *)
 
G4HadronicProcess * FindCaptureProcess ()
 
G4HadronicProcess * FindFissionProcess ()
 
- 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 48 of file G4HadronInelasticQBBC.hh.

Constructor & Destructor Documentation

G4HadronInelasticQBBC::G4HadronInelasticQBBC ( G4int  ver = 1)

Definition at line 85 of file G4HadronInelasticQBBC.cc.

86  : G4VHadronPhysics("hInelasticQBBC"),verbose(ver)
87 {
88  theAntiNuclXS = nullptr;
89  theKaonXS = nullptr;
90 }
G4VHadronPhysics::G4VHadronPhysics
G4VHadronPhysics(const G4String &name="hInelastic", G4int verbose=0)
Definition: G4VHadronPhysics.cc:59
G4HadronInelasticQBBC::G4HadronInelasticQBBC ( const G4String &  name,
G4int  ver = 1,
G4bool  ftf = false,
G4bool  bert = false,
G4bool  chips = false,
G4bool  hp = false,
G4bool  glauber = false 
)

Definition at line 92 of file G4HadronInelasticQBBC.cc.

93  : G4HadronInelasticQBBC(ver)
94 {}
G4HadronInelasticQBBC::~G4HadronInelasticQBBC ( )
virtual

Definition at line 96 of file G4HadronInelasticQBBC.cc.

97 {
98  delete theAntiNuclXS; theAntiNuclXS=nullptr;
99  delete theKaonXS; theKaonXS=nullptr;
100 }

Member Function Documentation

void G4HadronInelasticQBBC::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 102 of file G4HadronInelasticQBBC.cc.

103 {
104  if(verbose > 1) {
105  G4cout << "### HadronInelasticQBBC Construct Process " << G4endl;
106  }
107 
108  G4double emax = 100.*TeV;
109 
110  //G4cout << "G4HadronInelasticQBBC::ConstructProcess new PRECO"<< G4endl;
111 
112  // PreCompound and Evaporation models are instantiated here
113  G4PreCompoundModel* thePreCompound = nullptr;
114  G4HadronicInteraction* p =
115  G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
116  thePreCompound = static_cast<G4PreCompoundModel*>(p);
117  if(!thePreCompound) { thePreCompound = new G4PreCompoundModel(); }
118 
119  // configure models
120  //G4HadronicInteraction* theQGSP =
121  // BuildModel(new G4QGSBuilder("QGSP",thePreCompound,true,false),12.5*GeV,emax);
122  G4HadronicInteraction* theFTFP =
123  BuildModel(new G4FTFBuilder("FTFP",thePreCompound),3.0*GeV,emax);
124  G4HadronicInteraction* theFTFP1 =
125  BuildModel(new G4FTFBuilder("FTFP",thePreCompound),3.0*GeV,emax);
126  G4HadronicInteraction* theFTFP2 =
127  BuildModel(new G4FTFBuilder("FTFP",thePreCompound),0.0,emax);
128 
129  G4CascadeInterface* casc = new G4CascadeInterface();
130  casc->usePreCompoundDeexcitation();
131  G4HadronicInteraction* theBERT = NewModel(casc,1.0*GeV,5.0*GeV);
132 
133  casc = new G4CascadeInterface();
134  casc->usePreCompoundDeexcitation();
135  G4HadronicInteraction* theBERT1 = NewModel(casc,0.0*GeV,5.0*GeV);
136 
137  //G4cout << "G4HadronInelasticQBBC::ConstructProcess new Binary"<< G4endl;
138  G4BinaryCascade* bic = new G4BinaryCascade(thePreCompound);
139  G4HadronicInteraction* theBIC = NewModel(bic,0.0,1.5*GeV);
140 
141  // cross sections
142  theAntiNuclXS = new G4ComponentAntiNuclNuclearXS();
143  G4CrossSectionInelastic* anucxs =
144  new G4CrossSectionInelastic(theAntiNuclXS);
145  theKaonXS = new G4ComponentGGHadronNucleusXsc();
146  G4CrossSectionInelastic* kaonxs = new G4CrossSectionInelastic(theKaonXS);
147 
148  // loop over particles
149  auto myParticleIterator=GetParticleIterator();
150  myParticleIterator->reset();
151  while( (*myParticleIterator)() ) {
152  G4ParticleDefinition* particle = myParticleIterator->value();
153  G4String pname = particle->GetParticleName();
154  if(verbose > 1) {
155  G4cout << "### HadronInelasticQBBC: " << pname << G4endl;
156  }
157 
158  //
159  // model and X-section configuration per particle type
160  //
161  if(pname == "proton") {
162  G4HadronicProcess* hp = FindInelasticProcess(particle);
163  hp->AddDataSet(new G4BGGNucleonInelasticXS(particle));
164 
165  //hp->RegisterMe(theQGSP);
166  hp->RegisterMe(theFTFP);
167  hp->RegisterMe(theBERT);
168  hp->RegisterMe(theBIC);
169 
170  } else if(pname == "neutron") {
171  G4HadronicProcess* hp = FindInelasticProcess(particle);
172  hp->AddDataSet((G4NeutronInelasticXS*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronInelasticXS::Default_Name()));
173  //hp->RegisterMe(theQGSP);
174  hp->RegisterMe(theFTFP);
175 
176  G4HadronicProcess* capture = FindCaptureProcess();
177  capture->AddDataSet((G4NeutronCaptureXS*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronCaptureXS::Default_Name()));
178  hp->RegisterMe(theBERT);
179  hp->RegisterMe(theBIC);
180  capture->RegisterMe(new G4NeutronRadCapture());
181 
182  } else if(pname == "pi-" || pname == "pi+") {
183  G4HadronicProcess* hp = FindInelasticProcess(particle);
184  hp->AddDataSet(new G4BGGPionInelasticXS(particle));
185  //hp->AddDataSet(new G4CrossSectionPairGG((G4PiNuclearCrossSection*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4PiNuclearCrossSection::Default_Name()), 91*GeV));
186  //hp->RegisterMe(theQGSP);
187  hp->RegisterMe(theFTFP);
188  hp->RegisterMe(theBERT1);
189 
190  } else if(pname == "kaon-" ) {
191  G4HadronicProcess* hp = FindInelasticProcess(particle);
192  hp->RegisterMe(theFTFP1);
193  hp->RegisterMe(theBERT1);
194  hp->AddDataSet(kaonxs);
195 
196  } else if(pname == "kaon+" ) {
197  G4HadronicProcess* hp = FindInelasticProcess(particle);
198  hp->RegisterMe(theFTFP1);
199  hp->RegisterMe(theBERT1);
200  hp->AddDataSet(kaonxs);
201 
202  } else if(pname == "kaon0S" ||
203  pname == "kaon0L") {
204  G4HadronicProcess* hp = FindInelasticProcess(particle);
205  hp->RegisterMe(theFTFP1);
206  hp->RegisterMe(theBERT1);
207  hp->AddDataSet(kaonxs);
208 
209  } else if(pname == "lambda" ||
210  pname == "omega-" ||
211  pname == "sigma-" ||
212  pname == "sigma+" ||
213  pname == "sigma0" ||
214  pname == "xi-" ||
215  pname == "xi0") {
216  G4HadronicProcess* hp = FindInelasticProcess(particle);
217  hp->RegisterMe(theFTFP1);
218  hp->RegisterMe(theBERT1);
219  hp->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsHyperonInelasticXS::Default_Name()));
220 
221  } else if(pname == "anti_alpha" ||
222  pname == "anti_deuteron"||
223  pname == "anti_He3" ||
224  pname == "anti_proton" ||
225  pname == "anti_triton" ||
226  pname == "anti_lambda" ||
227  pname == "anti_neutron" ||
228  pname == "anti_omega-" ||
229  pname == "anti_sigma-" ||
230  pname == "anti_sigma+" ||
231  pname == "anti_xi-" ||
232  pname == "anti_xi0"
233  ) {
234 
235  G4HadronicProcess* hp = FindInelasticProcess(particle);
236  hp->RegisterMe(theFTFP2);
237  hp->AddDataSet(anucxs);
238 
239  }
240  }
241 }
G4VHadronPhysics::FindCaptureProcess
G4HadronicProcess * FindCaptureProcess()
Definition: G4VHadronPhysics.cc:287
p
const char * p
Definition: xmltok.h:285
G4NeutronInelasticXS::Default_Name
static const char * Default_Name()
Definition: G4NeutronInelasticXS.hh:70
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4HadronicProcess::RegisterMe
void RegisterMe(G4HadronicInteraction *a)
Definition: G4HadronicProcess.cc:153
TeV
static constexpr double TeV
Definition: G4SIunits.hh:218
G4VPhysicsConstructor::GetParticleIterator
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
Definition: G4VPhysicsConstructor.cc:78
G4HadronicProcess::AddDataSet
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
Definition: G4HadronicProcess.hh:111
G4cout
G4GLOB_DLL std::ostream G4cout
G4CrossSectionDataSetRegistry::Instance
static G4CrossSectionDataSetRegistry * Instance()
Definition: G4CrossSectionDataSetRegistry.cc:79
G4VHadronPhysics::FindInelasticProcess
G4HadronicProcess * FindInelasticProcess(const G4String &)
Definition: G4VHadronPhysics.cc:217
eplot.pname
string pname
Definition: eplot.py:33
emax
static const G4double emax
Definition: G4ChatterjeeCrossSection.hh:39
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition: G4HadronicInteractionRegistry.cc:117
G4VHadronPhysics::NewModel
G4HadronicInteraction * NewModel(G4HadronicInteraction *, G4double emin, G4double emax)
Definition: G4VHadronPhysics.cc:116
G4NeutronCaptureXS::Default_Name
static const char * Default_Name()
Definition: G4NeutronCaptureXS.hh:69
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition: G4HadronicInteractionRegistry.cc:36
GeV
static constexpr double GeV
Definition: G4SIunits.hh:217
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4VHadronPhysics::BuildModel
G4HadronicInteraction * BuildModel(G4VHadronModelBuilder *, G4double emin, G4double emax)
Definition: G4VHadronPhysics.cc:96

Here is the call graph for this function:

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