Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4HadronInelasticQBBC Class Reference

#include <G4HadronInelasticQBBC.hh>

+ Inheritance diagram for G4HadronInelasticQBBC:

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)
 
 ~G4HadronInelasticQBBC () override
 
void ConstructProcess () override
 
 G4HadronInelasticQBBC (G4HadronInelasticQBBC &)=delete
 
G4HadronInelasticQBBCoperator= (const G4HadronInelasticQBBC &right)=delete
 
- Public Member Functions inherited from G4VHadronPhysics
 G4VHadronPhysics (const G4String &name="hInelastic", G4int verbose=0)
 
 ~G4VHadronPhysics () override
 
void ConstructParticle () override
 
G4HadronicInteractionBuildModel (G4VHadronModelBuilder *, G4double emin, G4double emax)
 
G4HadronicInteractionNewModel (G4HadronicInteraction *, G4double emin, G4double emax)
 
 G4VHadronPhysics (G4VHadronPhysics &)=delete
 
G4VHadronPhysicsoperator= (const G4VHadronPhysics &right)=delete
 
- Public Member Functions inherited from G4VPhysicsConstructor
 G4VPhysicsConstructor (const G4String &="")
 
 G4VPhysicsConstructor (const G4String &name, G4int physics_type)
 
virtual ~G4VPhysicsConstructor ()
 
virtual void ConstructParticle ()=0
 
virtual void ConstructProcess ()=0
 
void SetPhysicsName (const G4String &="")
 
const G4StringGetPhysicsName () const
 
void SetPhysicsType (G4int)
 
G4int GetPhysicsType () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
G4int GetInstanceID () const
 
virtual void TerminateWorker ()
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManagerGetSubInstanceManager ()
 
- Protected Types inherited from G4VPhysicsConstructor
using PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V
 
- Protected Member Functions inherited from G4VPhysicsConstructor
G4bool RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
 
G4ParticleTable::G4PTblDicIteratorGetParticleIterator () const
 
PhysicsBuilder_V GetBuilders () const
 
void AddBuilder (G4PhysicsBuilderInterface *bld)
 
- 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 44 of file G4HadronInelasticQBBC.hh.

Constructor & Destructor Documentation

◆ G4HadronInelasticQBBC() [1/3]

G4HadronInelasticQBBC::G4HadronInelasticQBBC ( G4int  ver = 1)

Definition at line 79 of file G4HadronInelasticQBBC.cc.

80 : G4VHadronPhysics("hInelasticQBBC"),verbose(ver)
81{
83}
static G4HadronicParameters * Instance()
void SetEnableBCParticles(G4bool val)

◆ G4HadronInelasticQBBC() [2/3]

G4HadronInelasticQBBC::G4HadronInelasticQBBC ( const G4String name,
G4int  ver = 1,
G4bool  ftf = false,
G4bool  bert = false,
G4bool  chips = false,
G4bool  hp = false,
G4bool  glauber = false 
)

◆ ~G4HadronInelasticQBBC()

G4HadronInelasticQBBC::~G4HadronInelasticQBBC ( )
override

Definition at line 89 of file G4HadronInelasticQBBC.cc.

90{}

◆ G4HadronInelasticQBBC() [3/3]

G4HadronInelasticQBBC::G4HadronInelasticQBBC ( G4HadronInelasticQBBC )
delete

Member Function Documentation

◆ ConstructProcess()

void G4HadronInelasticQBBC::ConstructProcess ( )
overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 92 of file G4HadronInelasticQBBC.cc.

93{
95 G4bool useFactorXS = param->ApplyFactorXS();
97
98 // configure models
99 const G4double eminFtf = param->GetMinEnergyTransitionFTF_Cascade();
100 const G4double eminBert = 1.0*CLHEP::GeV;
101 const G4double emaxBic = 1.5*CLHEP::GeV;
102 const G4double emaxBert = param->GetMaxEnergyTransitionFTF_Cascade();
103 const G4double emaxBertPions = 12.*CLHEP::GeV;
104 const G4double emax = param->GetMaxEnergy();
105
106 if(G4Threading::IsMasterThread() && verbose > 0) {
107 G4cout << "### HadronInelasticQBBC Construct Process:\n"
108 << " Emin(FTFP)= " << eminFtf/CLHEP::GeV
109 << " GeV; Emax(FTFP)= " << emax/CLHEP::GeV << " GeV\n"
110 << " Emin(BERT)= " << eminBert/CLHEP::GeV
111 << " GeV; Emax(BERT)= " << emaxBert/CLHEP::GeV
112 << " GeV; Emax(BERTpions)= " << emaxBertPions/CLHEP::GeV
113 << " GeV;\n" << " Emin(BIC) = 0 GeV; Emax(BIC)= "
114 << emaxBic/CLHEP::GeV << " GeV." << G4endl;
115 }
116
117 // PreCompound and Evaporation models are instantiated here
118 G4PreCompoundModel* thePreCompound = nullptr;
121 thePreCompound = static_cast<G4PreCompoundModel*>(p);
122 if(!thePreCompound) { thePreCompound = new G4PreCompoundModel(); }
123
124 auto theFTFP = new G4TheoFSGenerator("FTFP");
125 auto theStringModel = new G4FTFModel();
126 theStringModel->SetFragmentationModel(new G4ExcitedStringDecay());
127 theFTFP->SetHighEnergyGenerator( theStringModel );
128 theFTFP->SetTransport( new G4GeneratorPrecompoundInterface() );
129 theFTFP->SetMinEnergy( eminFtf );
130 theFTFP->SetMaxEnergy( emax );
131
132 auto theBERT = new G4CascadeInterface();
133 theBERT->SetMinEnergy( eminBert );
134 theBERT->SetMaxEnergy( emaxBert );
135 theBERT->usePreCompoundDeexcitation();
136
137 auto theBERT1 = new G4CascadeInterface();
138 theBERT1->SetMinEnergy( eminBert );
139 theBERT1->SetMaxEnergy( emaxBertPions );
140 theBERT1->usePreCompoundDeexcitation();
141
142 auto theBIC = new G4BinaryCascade(thePreCompound);
143 theBIC->SetMaxEnergy( emaxBic );
144
145 // p
147 G4HadronicProcess* hp =
148 new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
149 hp->AddDataSet(new G4ParticleInelasticXS(particle));
150 hp->RegisterMe(theFTFP);
151 hp->RegisterMe(theBERT);
152 hp->RegisterMe(theBIC);
153 ph->RegisterProcess(hp, particle);
154 if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
155
156 // n
157 particle = G4Neutron::Neutron();
158 hp = new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
160 hp->RegisterMe(theFTFP);
161 hp->RegisterMe(theBERT);
162 hp->RegisterMe(theBIC);
163 ph->RegisterProcess(hp, particle);
164 if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
165
166 hp = new G4HadronCaptureProcess("nCapture");
168 ph->RegisterProcess(hp, particle);
169
170 // pi+
171 particle = G4PionPlus::PionPlus();
172 hp = new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
173 hp->AddDataSet(new G4BGGPionInelasticXS(particle));
174 hp->RegisterMe(theFTFP);
175 hp->RegisterMe(theBERT1);
176 hp->RegisterMe(theBIC);
177 ph->RegisterProcess(hp, particle);
178 if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
179
180 // pi-
181 particle = G4PionMinus::PionMinus();
182 hp = new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
183 hp->AddDataSet(new G4BGGPionInelasticXS(particle));
184 hp->RegisterMe(theFTFP);
185 hp->RegisterMe(theBERT1);
186 hp->RegisterMe(theBIC);
187 ph->RegisterProcess(hp, particle);
188 if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
189
190 // kaons
192
193 // high energy particles
194 if( emax > param->EnergyThresholdForHeavyHadrons() ) {
195
196 // pbar, nbar, anti light ions
198
199 // hyperons
201
202 // b-, c- baryons and mesons
203 if( param->EnableBCParticles() ) {
205 }
206 }
207}
double G4double
Definition: G4Types.hh:83
bool G4bool
Definition: G4Types.hh:86
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
static void BuildBCHadronsFTFP_BERT()
static void BuildHyperonsFTFP_BERT()
static void BuildKaonsFTFP_BERT()
static void BuildAntiLightIonsFTFP()
G4HadronicInteraction * FindModel(const G4String &name)
static G4HadronicInteractionRegistry * Instance()
G4double GetMinEnergyTransitionFTF_Cascade() const
G4double GetMaxEnergyTransitionFTF_Cascade() const
G4bool EnableBCParticles() const
G4double EnergyThresholdForHeavyHadrons() const
G4double XSFactorPionInelastic() const
G4double GetMaxEnergy() const
G4double XSFactorNucleonInelastic() const
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
void MultiplyCrossSectionBy(G4double factor)
void RegisterMe(G4HadronicInteraction *a)
static G4Neutron * Neutron()
Definition: G4Neutron.cc:103
const G4String & GetParticleName() const
G4bool RegisterProcess(G4VProcess *process, G4ParticleDefinition *particle)
static G4PhysicsListHelper * GetPhysicsListHelper()
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:97
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:97
static G4Proton * Proton()
Definition: G4Proton.cc:92
G4bool IsMasterThread()
Definition: G4Threading.cc:124

◆ operator=()

G4HadronInelasticQBBC & G4HadronInelasticQBBC::operator= ( const G4HadronInelasticQBBC right)
delete

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