GammaLeptoNuclearPhysics Class Reference

Gamma-nuclear, electro-nuclear and muon-nuclear physics constructor for Belle II physics list. More...

#include <GammaLeptoNuclearPhysics.h>

Inheritance diagram for GammaLeptoNuclearPhysics:

Public Member Functions
	GammaLeptoNuclearPhysics (const G4int verbosityLevel=0)
	Constructor.
	~GammaLeptoNuclearPhysics ()
	Destructor.
virtual void	ConstructProcess () override
	Build processes, models, cross sections used in physics list.
virtual void	ConstructParticle () override
	Build all particle types used in physics list (empty in this class)

Private Attributes
G4TheoFSGenerator *	m_qgsp
	Final state generator for QCD string models.
G4QGSModel< G4GammaParticipants > *	m_stringModel
	Quark Gluon String model.
G4ExcitedStringDecay *	m_stringDecay
	Model to decay strings into hadrons.
G4QGSMFragmentation *	m_fragModel
	Quark Gluon String fragmentation model.
G4GeneratorPrecompoundInterface *	m_preCompoundModel
	Precompound model to deexcite post-collision nucleus.
G4bool	m_useGammaNuclearXS
	Switch between GammaNuclearXS and PhotoNuclearXS.

Detailed Description

Gamma-nuclear, electro-nuclear and muon-nuclear physics constructor for Belle II physics list.

Definition at line 33 of file GammaLeptoNuclearPhysics.h.

Constructor & Destructor Documentation

GammaLeptoNuclearPhysics()

GammaLeptoNuclearPhysics

(

const G4int

verbosityLevel = 0

)

Constructor.

Definition at line 39 of file GammaLeptoNuclearPhysics.cc.

  : m_qgsp(nullptr), m_stringModel(nullptr), m_stringDecay(nullptr),
    m_fragModel(nullptr), m_preCompoundModel(nullptr), m_useGammaNuclearXS(false)
{
  G4HadronicParameters::Instance()->SetVerboseLevel(verbosityLevel);
}

~GammaLeptoNuclearPhysics()

~GammaLeptoNuclearPhysics

(

)

Destructor.

Definition at line 47 of file GammaLeptoNuclearPhysics.cc.

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

Member Function Documentation

ConstructParticle()

void ConstructParticle ( )

overridevirtual

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

Definition at line 125 of file GammaLeptoNuclearPhysics.cc.

{}

ConstructProcess()

void ConstructProcess ( )

overridevirtual

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

Definition at line 56 of file GammaLeptoNuclearPhysics.cc.

{
  // Use Bertini cascade for low energies
  G4CascadeInterface* theGammaReaction = new G4CascadeInterface;
  theGammaReaction->SetMinEnergy(0.0);
  theGammaReaction->SetMaxEnergy(3.5 * GeV);
 
  // Use QGSP for high energies
  m_qgsp = new G4TheoFSGenerator("QGSP");
  m_stringModel = new G4QGSModel<G4GammaParticipants>;
  m_stringDecay =
    new G4ExcitedStringDecay(m_fragModel = new G4QGSMFragmentation);
  m_stringModel->SetFragmentationModel(m_stringDecay);
  m_preCompoundModel = new G4GeneratorPrecompoundInterface();
 
  m_qgsp->SetHighEnergyGenerator(m_stringModel);
  m_qgsp->SetTransport(m_preCompoundModel);
  m_qgsp->SetMinEnergy(3 * GeV);
  m_qgsp->SetMaxEnergy(100 * TeV);
 
  // Lepto-nuclear models
  G4ElectroVDNuclearModel* evdn = new G4ElectroVDNuclearModel;
  G4MuonVDNuclearModel* mvdn = new G4MuonVDNuclearModel;
 
 
  G4ProcessManager* procMan = 0;
 
  // Gamma
  procMan = G4Gamma::Gamma()->GetProcessManager();
  G4HadronInelasticProcess* pnProc = new G4HadronInelasticProcess("photonNuclear", G4Gamma::Definition());
  auto xsreg = G4CrossSectionDataSetRegistry::Instance();
  G4VCrossSectionDataSet* xs = nullptr;
  if (m_useGammaNuclearXS) {
    xs = xsreg->GetCrossSectionDataSet("GammaNuclearXS");
    if (nullptr == xs) xs = new G4GammaNuclearXS();
  } else {
    xs = xsreg->GetCrossSectionDataSet("PhotoNuclearXS");
    if (nullptr == xs) xs = new G4PhotoNuclearCrossSection();
  }
  pnProc->AddDataSet(xs);
  pnProc->RegisterMe(theGammaReaction);
  pnProc->RegisterMe(m_qgsp);
  procMan->AddDiscreteProcess(pnProc);
 
  // Electron
  procMan = G4Electron::Electron()->GetProcessManager();
  G4ElectronNuclearProcess* emn = new G4ElectronNuclearProcess;
  emn->RegisterMe(evdn);
  procMan->AddDiscreteProcess(emn);
 
  // Positron
  procMan = G4Positron::Positron()->GetProcessManager();
  G4PositronNuclearProcess* epn = new G4PositronNuclearProcess;
  epn->RegisterMe(evdn);
  procMan->AddDiscreteProcess(epn);
 
  // Muon-
  procMan = G4MuonMinus::MuonMinus()->GetProcessManager();
  G4MuonNuclearProcess* mun = new G4MuonNuclearProcess;
  mun->RegisterMe(mvdn);
  procMan->AddDiscreteProcess(mun);
 
  // Muon+
  procMan = G4MuonPlus::MuonPlus()->GetProcessManager();
  procMan->AddDiscreteProcess(mun);
 
}

Member Data Documentation

m_fragModel

G4QGSMFragmentation* m_fragModel

private

Quark Gluon String fragmentation model.

Definition at line 57 of file GammaLeptoNuclearPhysics.h.

m_preCompoundModel

G4GeneratorPrecompoundInterface* m_preCompoundModel

private

Precompound model to deexcite post-collision nucleus.

Definition at line 60 of file GammaLeptoNuclearPhysics.h.

m_qgsp

G4TheoFSGenerator* m_qgsp

private

Final state generator for QCD string models.

Definition at line 48 of file GammaLeptoNuclearPhysics.h.

m_stringDecay

G4ExcitedStringDecay* m_stringDecay

private

Model to decay strings into hadrons.

Definition at line 54 of file GammaLeptoNuclearPhysics.h.

m_stringModel

G4QGSModel<G4GammaParticipants>* m_stringModel

private

Quark Gluon String model.

Definition at line 51 of file GammaLeptoNuclearPhysics.h.

m_useGammaNuclearXS

G4bool m_useGammaNuclearXS

private

Switch between GammaNuclearXS and PhotoNuclearXS.

Definition at line 63 of file GammaLeptoNuclearPhysics.h.

---

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

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