IonPhysics Class Reference

Ion hadronic physics constructor for Belle II physics list. More...

#include <IonPhysics.h>

Inheritance diagram for IonPhysics:

Collaboration diagram for IonPhysics:

Public Member Functions
virtual void	ConstructParticle () override
	Build processes, models, cross sections used in physics list.
virtual void	ConstructProcess () override
	Build all particle types used in physics list (empty in this class)

Private Attributes
G4TheoFSGenerator *	m_ftfp
	Final state generator for QCD string models.
G4FTFModel *	m_stringModel
	Fritiof string model.
G4ExcitedStringDecay *	m_stringDecay
	Model to decay strings into hadrons.
G4LundStringFragmentation *	m_fragModel
	Lund string fragmentation model.
G4GeneratorPrecompoundInterface *	m_preCompoundModel
	Precompound model to deexcite post-collision nucleus.
G4VComponentCrossSection *	m_theGGNuclNuclXS
	Cross section set for inelastic nucleus-nucleus collisions.
G4ComponentGGNuclNuclXsc *	m_ionGGXS
	Cross section set for elastic nucleus-nucleus collisions.

Detailed Description

Ion hadronic physics constructor for Belle II physics list.

Definition at line 33 of file IonPhysics.h.

Constructor & Destructor Documentation

IonPhysics()

IonPhysics

(

)

Definition at line 35 of file IonPhysics.cc.

  : m_ftfp(nullptr), m_stringModel(nullptr), m_stringDecay(nullptr),
    m_fragModel(nullptr), m_preCompoundModel(nullptr)
{}

~IonPhysics()

~IonPhysics

(

)

Definition at line 41 of file IonPhysics.cc.

{
  delete m_stringDecay;
  delete m_stringModel;
  delete m_fragModel;
  delete m_preCompoundModel;
 
  delete m_theGGNuclNuclXS;
  delete m_ionGGXS;
}

Member Function Documentation

ConstructParticle()

void ConstructParticle ( )

overridevirtual

Build processes, models, cross sections used in physics list.

Definition at line 53 of file IonPhysics.cc.

{}

ConstructProcess()

void ConstructProcess ( )

overridevirtual

Build all particle types used in physics list (empty in this class)

Definition at line 57 of file IonPhysics.cc.

{
  G4ProcessManager* procMan = 0;
 
  // Elastic model for generic ions (z > 2)
  G4NuclNuclDiffuseElastic* ionElastic = new G4NuclNuclDiffuseElastic;
  ionElastic->SetMinEnergy(0.0);
 
  // FTFP    ==>>   eventually replace this with new class FTFPInterface
  m_ftfp = new G4TheoFSGenerator("FTFP");
  m_stringModel = new G4FTFModel;
  m_stringDecay =
    new G4ExcitedStringDecay(m_fragModel = new G4LundStringFragmentation);
  m_stringModel->SetFragmentationModel(m_stringDecay);
  m_preCompoundModel = new G4GeneratorPrecompoundInterface();
 
  m_ftfp->SetHighEnergyGenerator(m_stringModel);
  m_ftfp->SetTransport(m_preCompoundModel);
  m_ftfp->SetMinEnergy(10.01 * GeV);
  m_ftfp->SetMaxEnergy(1.0 * TeV);
 
  // QMD model
  G4QMDReaction* qmd = new G4QMDReaction;
  qmd->SetMinEnergy(100.0 * MeV);
  qmd->SetMaxEnergy(10.0 * GeV);
 
  // BIC ion model
  G4HadronicInteraction* p =
    G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
  G4PreCompoundModel* thePreCompound = static_cast<G4PreCompoundModel*>(p);
  if (!thePreCompound) { thePreCompound = new G4PreCompoundModel; }
 
  G4BinaryLightIonReaction* ionBC = new G4BinaryLightIonReaction(thePreCompound);
  ionBC->SetMinEnergy(0.0 * MeV);
  ionBC->SetMaxEnergy(110.0 * MeV);
 
  // Elastic cross section set
  m_ionGGXS = new G4ComponentGGNuclNuclXsc;
  G4VCrossSectionDataSet* ionElasticXS = new G4CrossSectionElastic(m_ionGGXS);
  ionElasticXS->SetMinKinEnergy(0.0);
 
  // Inelastic cross section set
  m_theGGNuclNuclXS = new G4ComponentGGNuclNuclXsc();
  G4VCrossSectionDataSet* nuclNuclXS =
    new G4CrossSectionInelastic(m_theGGNuclNuclXS);
 
  //   Deuteron                                                               //
 
  procMan = G4Deuteron::Deuteron()->GetProcessManager();
 
  // elastic
  // no model available
 
  // inelastic
  G4HadronInelasticProcess* deutProcInel =
    new G4HadronInelasticProcess("DeuteronInelProcess", G4Deuteron::Deuteron());
  deutProcInel->RegisterMe(ionBC);
  deutProcInel->RegisterMe(qmd);
  deutProcInel->RegisterMe(m_ftfp);
  deutProcInel->AddDataSet(nuclNuclXS);
  procMan->AddDiscreteProcess(deutProcInel);
 
  //   Triton                                                                 //
 
  procMan = G4Triton::Triton()->GetProcessManager();
 
  // elastic
  // no model available
 
  // inelastic
  G4HadronInelasticProcess* tritProcInel =
    new G4HadronInelasticProcess("TritonInelProcess", G4Triton::Triton());
  tritProcInel->RegisterMe(ionBC);
  tritProcInel->RegisterMe(qmd);
  tritProcInel->RegisterMe(m_ftfp);
  tritProcInel->AddDataSet(nuclNuclXS);
  procMan->AddDiscreteProcess(tritProcInel);
 
  //   He3                                                                    //
 
  procMan = G4He3::He3()->GetProcessManager();
 
  // elastic
  // no model available
 
  // inelastic
  G4HadronInelasticProcess* he3ProcInel =
    new G4HadronInelasticProcess("He3InelProcess", G4He3::He3());
  he3ProcInel->RegisterMe(ionBC);
  he3ProcInel->RegisterMe(qmd);
  he3ProcInel->RegisterMe(m_ftfp);
  he3ProcInel->AddDataSet(nuclNuclXS);
  procMan->AddDiscreteProcess(he3ProcInel);
 
  //   Alpha                                                                  //
 
  procMan = G4Alpha::Alpha()->GetProcessManager();
 
  // elastic
  // no model available
 
  // inelastic
  G4HadronInelasticProcess* alphProcInel =
    new G4HadronInelasticProcess("AlphaInelProcess", G4Alpha::Alpha());
  alphProcInel->RegisterMe(ionBC);
  alphProcInel->RegisterMe(qmd);
  alphProcInel->RegisterMe(m_ftfp);
  alphProcInel->AddDataSet(nuclNuclXS);
  procMan->AddDiscreteProcess(alphProcInel);
 
  //   Generic ion                                                            //
 
  procMan = G4GenericIon::GenericIon()->GetProcessManager();
 
  // elastic
  G4HadronElasticProcess* ionProcEl = new G4HadronElasticProcess;
  ionProcEl->RegisterMe(ionElastic);
  ionProcEl->AddDataSet(ionElasticXS);
  procMan->AddDiscreteProcess(ionProcEl);
 
  // inelastic
  G4HadronInelasticProcess* genIonProcInel =
    new G4HadronInelasticProcess("IonInelProcess", G4GenericIon::GenericIon());
  genIonProcInel->RegisterMe(ionBC);
  genIonProcInel->RegisterMe(qmd);
  genIonProcInel->RegisterMe(m_ftfp);
  genIonProcInel->AddDataSet(nuclNuclXS);
  procMan->AddDiscreteProcess(genIonProcInel);
 
}

Member Data Documentation

m_fragModel

G4LundStringFragmentation* m_fragModel

private

Lund string fragmentation model.

Definition at line 55 of file IonPhysics.h.

m_ftfp

G4TheoFSGenerator* m_ftfp

private

Final state generator for QCD string models.

Definition at line 46 of file IonPhysics.h.

m_ionGGXS

G4ComponentGGNuclNuclXsc* m_ionGGXS

private

Cross section set for elastic nucleus-nucleus collisions.

Definition at line 64 of file IonPhysics.h.

m_preCompoundModel

G4GeneratorPrecompoundInterface* m_preCompoundModel

private

Precompound model to deexcite post-collision nucleus.

Definition at line 58 of file IonPhysics.h.

m_stringDecay

G4ExcitedStringDecay* m_stringDecay

private

Model to decay strings into hadrons.

Definition at line 52 of file IonPhysics.h.

m_stringModel

G4FTFModel* m_stringModel

private

Fritiof string model.

Definition at line 49 of file IonPhysics.h.

m_theGGNuclNuclXS

G4VComponentCrossSection* m_theGGNuclNuclXS

private

Cross section set for inelastic nucleus-nucleus collisions.

Definition at line 61 of file IonPhysics.h.

---

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

• simulation/physicslist/include/IonPhysics.h
• simulation/physicslist/src/IonPhysics.cc