#include <G4HadronInelasticQBBC_ABLA.hh>

Inheritance diagram for G4HadronInelasticQBBC_ABLA:

Public Member Functions
	G4HadronInelasticQBBC_ABLA (G4int ver=1)

virtual	~G4HadronInelasticQBBC_ABLA ()=default

void	ConstructProcess () override

	G4HadronInelasticQBBC_ABLA (G4HadronInelasticQBBC_ABLA &)=delete

G4HadronInelasticQBBC_ABLA &	operator= (const G4HadronInelasticQBBC_ABLA &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 ()

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 54 of file G4HadronInelasticQBBC_ABLA.hh.

Constructor & Destructor Documentation

◆ G4HadronInelasticQBBC_ABLA() [1/2]

G4HadronInelasticQBBC_ABLA::G4HadronInelasticQBBC_ABLA ( G4int ver = 1 )

Definition at line 79 of file G4HadronInelasticQBBC_ABLA.cc.

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

◆ ~G4HadronInelasticQBBC_ABLA()

virtual G4HadronInelasticQBBC_ABLA::~G4HadronInelasticQBBC_ABLA ( )

virtualdefault

◆ G4HadronInelasticQBBC_ABLA() [2/2]

G4HadronInelasticQBBC_ABLA::G4HadronInelasticQBBC_ABLA ( G4HadronInelasticQBBC_ABLA & )

delete

Member Function Documentation

◆ ConstructProcess()

void G4HadronInelasticQBBC_ABLA::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 90 of file G4HadronInelasticQBBC_ABLA.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.0*CLHEP::GeV;
  const G4double emax = param->GetMaxEnergy();
  if ( G4Threading::IsMasterThread()  &&  param->GetVerboseLevel() > 0 ) {
    G4cout << "### HadronInelasticQBBC_ABLA 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;
  }
 
  G4PreCompoundModel* thePreCompound = nullptr;
  G4HadronicInteraction* p = G4HadronicInteractionRegistry::Instance()->FindModel( "PRECO" );
  thePreCompound = static_cast< G4PreCompoundModel* >( p );
  if ( thePreCompound == nullptr ) thePreCompound = new G4PreCompoundModel;
 
  G4AblaInterface* theAblaInterface = nullptr;
  G4HadronicInteraction* pAbla = G4HadronicInteractionRegistry::Instance()->FindModel( "ABLAXX" );
  theAblaInterface = static_cast< G4AblaInterface* >( pAbla );
  if ( theAblaInterface == nullptr ) theAblaInterface = new G4AblaInterface;
 
  auto theFTFP = new G4TheoFSGenerator( "FTFP" );
  auto theStringModel = new G4FTFModel;
  theStringModel->SetFragmentationModel( new G4ExcitedStringDecay );
  theFTFP->SetHighEnergyGenerator( theStringModel );
  // theFTFP->SetTransport( new G4GeneratorPrecompoundInterface );
  theFTFP->SetTransport( new G4GeneratorPrecompoundInterface( theAblaInterface ) );
  theFTFP->SetMinEnergy( eminFtf );
  theFTFP->SetMaxEnergy( emax );
 
  auto theBERT = new G4CascadeInterface;
  theBERT->SetMinEnergy( eminBert );
  theBERT->SetMaxEnergy( emaxBert );
  // theBERT->usePreCompoundDeexcitation();
  theBERT->useAblaDeexcitation();
 
  auto theBERT1 = new G4CascadeInterface;
  theBERT1->SetMinEnergy( eminBert );
  theBERT1->SetMaxEnergy( emaxBertPions );
  // theBERT1->usePreCompoundDeexcitation();
  theBERT1->useAblaDeexcitation();
  
  // auto theBIC = new G4BinaryCascade( thePreCompound );
  auto theBIC = new G4BinaryCascade( theAblaInterface );
  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_ABLA & G4HadronInelasticQBBC_ABLA::operator= ( const G4HadronInelasticQBBC_ABLA & right )