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)
virtual	~G4HadronInelasticQBBC ()=default
void	ConstructProcess () override
	G4HadronInelasticQBBC (G4HadronInelasticQBBC &)=delete
G4HadronInelasticQBBC &	operator= (const G4HadronInelasticQBBC &right)=delete
Public Member Functions inherited from G4VHadronPhysics
	G4VHadronPhysics (const G4String &name="hInelastic", G4int verbose=0)
virtual	~G4VHadronPhysics ()
void	ConstructParticle () override
G4HadronicInteraction *	BuildModel (G4VHadronModelBuilder *, G4double emin, G4double emax)
G4HadronicInteraction *	NewModel (G4HadronicInteraction *, G4double emin, G4double emax)
	G4VHadronPhysics (G4VHadronPhysics &)=delete
G4VHadronPhysics &	operator= (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 G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Types inherited from G4VPhysicsConstructor
using	PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
PhysicsBuilder_V	GetBuilders () const
void	AddBuilder (G4PhysicsBuilderInterface *bld)
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel = 0
G4String	namePhysics = ""
G4int	typePhysics = 0
G4ParticleTable *	theParticleTable = nullptr
G4int	g4vpcInstanceID = 0
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 82 of file G4HadronInelasticQBBC.cc.

  : G4VHadronPhysics("hInelasticQBBC")
{
  SetPhysicsType(bHadronInelastic);
  auto param = G4HadronicParameters::Instance();
  param->SetEnableBCParticles(true);
  param->SetEnableNeutronGeneralProcess(true);
  param->SetVerboseLevel(ver);
}

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
	)

Definition at line 92 of file G4HadronInelasticQBBC.cc.

                                           : G4HadronInelasticQBBC(ver)
{}

~G4HadronInelasticQBBC()

virtual G4HadronInelasticQBBC::~G4HadronInelasticQBBC ( )

virtualdefault

G4HadronInelasticQBBC() [3/3]

G4HadronInelasticQBBC::G4HadronInelasticQBBC ( G4HadronInelasticQBBC &  )

delete

Member Function Documentation

ConstructProcess()

void G4HadronInelasticQBBC::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 96 of file G4HadronInelasticQBBC.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  // configure models
  const G4double eminFtf = param->GetMinEnergyTransitionFTF_Cascade();
  const G4double eminBert = 1.0*CLHEP::GeV;
  const G4double emaxBic  = 1.5*CLHEP::GeV;
  const G4double emaxBert = param->GetMaxEnergyTransitionFTF_Cascade();
  const G4double emaxBertPions = 12.*CLHEP::GeV;
  const G4double emax = param->GetMaxEnergy();
 
  if(G4Threading::IsMasterThread() && param->GetVerboseLevel() > 0) {
    G4cout << "### HadronInelasticQBBC Construct Process:\n"
           << "    Emin(FTFP)= " << eminFtf/CLHEP::GeV 
           << " GeV; Emax(FTFP)= " << emax/CLHEP::GeV << " GeV\n"
       << "    Emin(BERT)= " << eminBert/CLHEP::GeV
       << " GeV; Emax(BERT)= " << emaxBert/CLHEP::GeV
           << " GeV; Emax(BERTpions)= " << emaxBertPions/CLHEP::GeV  
           << " GeV;\n" << "    Emin(BIC) = 0 GeV; Emax(BIC)= " 
           << emaxBic/CLHEP::GeV << " GeV." << G4endl;
  }
 
  // PreCompound and Evaporation models are instantiated here
  G4PreCompoundModel* thePreCompound = nullptr;
  G4HadronicInteraction* p =
    G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
  thePreCompound = static_cast<G4PreCompoundModel*>(p);
  if(!thePreCompound) { thePreCompound = new G4PreCompoundModel(); }
 
  auto theFTFP = new G4TheoFSGenerator("FTFP");
  auto theStringModel = new G4FTFModel();
  theStringModel->SetFragmentationModel(new G4ExcitedStringDecay());
  theFTFP->SetHighEnergyGenerator( theStringModel );
  theFTFP->SetTransport( new G4GeneratorPrecompoundInterface() );
  theFTFP->SetMinEnergy( eminFtf );
  theFTFP->SetMaxEnergy( emax );
 
  auto theBERT = new G4CascadeInterface();
  theBERT->SetMinEnergy( eminBert );
  theBERT->SetMaxEnergy( emaxBert );
  theBERT->usePreCompoundDeexcitation();
 
  auto theBERT1 = new G4CascadeInterface();
  theBERT1->SetMinEnergy( eminBert );
  theBERT1->SetMaxEnergy( emaxBertPions );
  theBERT1->usePreCompoundDeexcitation();
 
  auto theBIC = new G4BinaryCascade(thePreCompound);
  theBIC->SetMaxEnergy( emaxBic );
 
  // p
  G4ParticleDefinition* particle = G4Proton::Proton();
  G4HadronicProcess* hp = 
    new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
  hp->AddDataSet(new G4ParticleInelasticXS(particle));
  hp->RegisterMe(theFTFP);
  hp->RegisterMe(theBERT);
  hp->RegisterMe(theBIC);
  ph->RegisterProcess(hp, particle);
  if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
 
  // n
  particle = G4Neutron::Neutron();
  G4HadronicProcess* ni = 
    new G4HadronInelasticProcess( "neutronInelastic", particle );
  ni->RegisterMe(theFTFP);
  ni->RegisterMe(theBERT);
  ni->RegisterMe(theBIC);
  G4HadProcesses::BuildNeutronInelasticAndCapture(ni);
 
  // pi+
  particle = G4PionPlus::PionPlus();
  hp = new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
  hp->AddDataSet(new G4BGGPionInelasticXS(particle));
  hp->RegisterMe(theFTFP);
  hp->RegisterMe(theBERT1);
  hp->RegisterMe(theBIC);
  ph->RegisterProcess(hp, particle);
  if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
 
  // pi-
  particle = G4PionMinus::PionMinus();
  hp = new G4HadronInelasticProcess( particle->GetParticleName()+"Inelastic", particle );
  hp->AddDataSet(new G4BGGPionInelasticXS(particle));
  hp->RegisterMe(theFTFP);
  hp->RegisterMe(theBERT1);
  hp->RegisterMe(theBIC);
  ph->RegisterProcess(hp, particle);
  if( useFactorXS ) hp->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
 
  // kaons
  G4HadronicBuilder::BuildKaonsFTFP_BERT();
 
  // high energy particles
  if( emax > param->EnergyThresholdForHeavyHadrons() ) {
 
    // pbar, nbar, anti light ions
    G4HadronicBuilder::BuildAntiLightIonsFTFP();
 
    // hyperons
    G4HadronicBuilder::BuildHyperonsFTFP_BERT();
 
    // b-, c- baryons and mesons
    if( param->EnableBCParticles() ) {
      G4HadronicBuilder::BuildBCHadronsFTFP_BERT();
    }
  }
}

operator=()

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

delete

---

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

• G4HadronInelasticQBBC.hh
• G4HadronInelasticQBBC.cc