78 G4HadronPhysicsQGSP_BERT_HP::tpdata = 0;
126 void G4HadronPhysicsQGSP_BERT_HP::CreateModels()
128 G4bool quasiElasticFTF=
false;
129 G4bool quasiElasticQGS=
true;
138 tpdata->theNeutrons->RegisterMe(tpdata->theQGSPNeutron=
new G4QGSPNeutronBuilder(quasiElasticQGS));
139 tpdata->theQGSPNeutron->SetMinEnergy(minQGSP);
140 tpdata->theNeutrons->RegisterMe(tpdata->theFTFPNeutron=
new G4FTFPNeutronBuilder(quasiElasticFTF));
141 tpdata->theFTFPNeutron->SetMinEnergy(minFTFP);
142 tpdata->theFTFPNeutron->SetMaxEnergy(maxFTFP);
145 tpdata->theBertiniNeutron->SetMinEnergy(maxHP);
146 tpdata->theBertiniNeutron->SetMaxEnergy(maxBERT);
152 tpdata->theQGSPPro->SetMinEnergy(minQGSP);
154 tpdata->theFTFPPro->SetMinEnergy(minFTFP);
155 tpdata->theFTFPPro->SetMaxEnergy(maxFTFP);
157 tpdata->theBertiniPro->SetMaxEnergy(maxBERT);
160 tpdata->thePiK->RegisterMe(tpdata->theQGSPPiK=
new G4QGSPPiKBuilder(quasiElasticQGS));
161 tpdata->theQGSPPiK->SetMinEnergy(minQGSP);
162 tpdata->thePiK->RegisterMe(tpdata->theFTFPPiK=
new G4FTFPPiKBuilder(quasiElasticFTF));
163 tpdata->theFTFPPiK->SetMinEnergy(minFTFP);
164 tpdata->theFTFPPiK->SetMaxEnergy(maxFTFP);
166 tpdata->theBertiniPiK->SetMaxEnergy(maxBERT);
178 delete tpdata->theHPNeutron;
179 delete tpdata->theBertiniNeutron;
180 delete tpdata->theQGSPNeutron;
181 delete tpdata->theFTFPNeutron;
182 delete tpdata->theNeutrons;
183 delete tpdata->theBertiniPiK;
184 delete tpdata->theQGSPPiK;
185 delete tpdata->theFTFPPiK;
186 delete tpdata->thePiK;
187 delete tpdata->theBertiniPro;
188 delete tpdata->theQGSPPro;
189 delete tpdata->theFTFPPro;
190 delete tpdata->thePro;
191 delete tpdata->theFTFPAntiBaryon;
192 delete tpdata->theAntiBaryon;
193 delete tpdata->theHyperon;
195 delete tpdata; tpdata = 0;
216 if ( tpdata == 0 ) tpdata =
new ThreadPrivate;
218 tpdata->theNeutrons->Build();
219 tpdata->thePro->Build();
220 tpdata->thePiK->Build();
230 tpdata->theHyperon->Build();
231 tpdata->theAntiBaryon->Build();
238 for (
size_t i=0; i < static_cast<size_t>(pv->
size()); ++i ) {
239 if (
fCapture == ((*pv)[i])->GetProcessSubType() ) {
241 }
else if (
fFission == ((*pv)[i])->GetProcessSubType() ) {
250 capture->
AddDataSet(tpdata->xsNeutronCaptureXS);
G4VCrossSectionDataSet * GetCrossSectionDataSet(const G4String &name, G4bool warning=true)
virtual ~G4HadronPhysicsQGSP_BERT_HP()
virtual void ConstructParticle()
G4HadronPhysicsQGSP_BERT_HP(G4int verbose=1)
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
static G4KaonZeroLong * KaonZeroLong()
static void ConstructParticle()
static void ConstructParticle()
static void ConstructParticle()
static G4KaonMinus * KaonMinus()
static void ConstructParticle()
void RegisterMe(G4HadronicInteraction *a)
void SetMinEnergy(G4double anEnergy)
static G4HadronicProcess * FindInelasticProcess(const G4ParticleDefinition *)
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
static G4KaonZeroShort * KaonZeroShort()
static G4CrossSectionDataSetRegistry * Instance()
static G4Neutron * Neutron()
static const char * Default_Name()
G4ProcessManager * GetProcessManager() const
static constexpr double GeV
virtual void ConstructProcess()
static constexpr double MeV
static G4KaonPlus * KaonPlus()
G4ProcessVector * GetProcessList() const
#define G4_DECLARE_PHYSCONSTR_FACTORY(physics_constructor)