G4IonQMDPhysics Class Reference

#include <G4IonQMDPhysics.hh>

Inheritance diagram for G4IonQMDPhysics:

Collaboration diagram for G4IonQMDPhysics:

Public Member Functions
	G4IonQMDPhysics (G4int verb=0)
	G4IonQMDPhysics (const G4String &name, G4int ver=0)
virtual	~G4IonQMDPhysics ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
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
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const

Private Member Functions
void	AddProcess (const G4String &, G4ParticleDefinition *, G4BinaryLightIonReaction *, G4QMDReaction *, G4HadronicInteraction *)

Private Attributes
G4double	eminBIC
G4double	eminQMD
G4double	emaxQMD
G4double	emaxFTFP
G4double	overlap
G4int	verbose

Static Private Attributes
static G4ThreadLocal std::vector< G4HadronInelasticProcess * > *	p_list = 0
static G4ThreadLocal std::vector< G4HadronicInteraction * > *	model_list = 0
static G4ThreadLocal G4VCrossSectionDataSet *	theNuclNuclData = 0
static G4ThreadLocal G4VComponentCrossSection *	theGGNuclNuclXS = 0
static G4ThreadLocal G4BinaryLightIonReaction *	theIonBC = 0
static G4ThreadLocal G4HadronicInteraction *	theFTFP = 0
static G4ThreadLocal G4FTFBuilder *	theBuilder = 0
static G4ThreadLocal G4QMDReaction *	theQMD = 0
static G4ThreadLocal G4bool	wasActivated = false

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4int	g4vpcInstanceID
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 56 of file G4IonQMDPhysics.hh.

Constructor & Destructor Documentation

G4IonQMDPhysics() [1/2]

G4IonQMDPhysics::G4IonQMDPhysics

(

G4int

verb = 0

)

Definition at line 90 of file G4IonQMDPhysics.cc.

  :  G4VPhysicsConstructor("IonQMD"), verbose(ver)
{
  eminBIC  = 0.*MeV;
  eminQMD  = 100.*MeV;
  emaxQMD  = 10.*GeV;
  emaxFTFP = 1.*TeV;
  overlap  = 10*MeV;
  SetPhysicsType(bIons);
  if(verbose > 1) G4cout << "### G4IonQMDPhysics" << G4endl;
}

Here is the call graph for this function:

G4IonQMDPhysics() [2/2]

G4IonQMDPhysics::G4IonQMDPhysics	(	const G4String &	name,
		G4int	ver = 0
	)

Definition at line 102 of file G4IonQMDPhysics.cc.

  :  G4VPhysicsConstructor(name), verbose(ver)
{
  eminBIC  = 0.*MeV;
  eminQMD  = 100.*MeV;
  emaxQMD  = 10.*GeV;
  emaxFTFP = 1.*TeV;
  overlap  = 10*MeV;
  SetPhysicsType(bIons);
  if(verbose > 1) G4cout << "### G4IonQMDPhysics" << G4endl;
}

Here is the call graph for this function:

~G4IonQMDPhysics()

G4IonQMDPhysics::~G4IonQMDPhysics ( )

virtual

Definition at line 115 of file G4IonQMDPhysics.cc.

{
  if(wasActivated) {
    delete theBuilder; theBuilder = 0;
    delete theGGNuclNuclXS; theGGNuclNuclXS = 0; 
    delete theNuclNuclData; theNuclNuclData = 0;
    G4int i;
    if ( p_list ) {
      G4int n = p_list->size();
      for(i=0; i<n; i++) {delete (*p_list)[i];}
      delete p_list; p_list = 0;
    }
    if ( model_list ) {
      G4int n = model_list->size();
      for(i=0; i<n; i++) {delete (*model_list)[i];}
      delete model_list; model_list = 0;
    }
  }
}

Member Function Documentation

AddProcess()

void G4IonQMDPhysics::AddProcess ( const G4String &  name, G4ParticleDefinition *  p, G4BinaryLightIonReaction *  BIC, G4QMDReaction *  QMD, G4HadronicInteraction *  FTFP  )

private

Definition at line 166 of file G4IonQMDPhysics.cc.

{
  G4HadronInelasticProcess* hadi = new G4HadronInelasticProcess(name, p);
  if ( p_list == 0 ) p_list = new  std::vector<G4HadronInelasticProcess*>;
  p_list->push_back(hadi);
  G4ProcessManager* pManager = p->GetProcessManager();
  pManager->AddDiscreteProcess(hadi);

  hadi->AddDataSet(theNuclNuclData);

  BIC->SetMinEnergy(eminBIC);
  BIC->SetMaxEnergy(emaxQMD);  //reset when QMD is present 
  hadi->RegisterMe(BIC);

  if(QMD) {
    QMD->SetMinEnergy(eminQMD);
    BIC->SetMaxEnergy(eminQMD+overlap);
    QMD->SetMaxEnergy(emaxQMD);
    hadi->RegisterMe(QMD);
  }  

  if(FTFP) {
    FTFP->SetMinEnergy(emaxQMD - overlap);
    FTFP->SetMaxEnergy(emaxFTFP);
    hadi->RegisterMe(FTFP);
  }  

  if(verbose > 1) {
    G4cout << "Register " << hadi->GetProcessName()
       << " for " << p->GetParticleName() << G4endl
       << "       Binary Cascade for E(MeV)= " << eminBIC << " - " 
       << (QMD==0 ? emaxQMD : (eminQMD-overlap)) ;
    if(QMD) {
      G4cout  << G4endl <<"       QMD for E(MeV)= " << eminQMD << " - " << emaxQMD;
    }
    if(FTFP) {
      G4cout << G4endl<< "       FTFP for E(MeV)= " << emaxQMD-overlap << " - " << emaxFTFP;
    }
    G4cout << G4endl;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructParticle()

void G4IonQMDPhysics::ConstructParticle ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 212 of file G4IonQMDPhysics.cc.

{
  //  Construct light ions
  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

Here is the call graph for this function:

ConstructProcess()

void G4IonQMDPhysics::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 135 of file G4IonQMDPhysics.cc.

{
  if(wasActivated) return;
  wasActivated = true;

  G4HadronicInteraction* p =
    G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
  G4PreCompoundModel* thePreCompound = static_cast<G4PreCompoundModel*>(p); 
  if(!thePreCompound) { thePreCompound = new G4PreCompoundModel; }

  theIonBC = new G4BinaryLightIonReaction(thePreCompound);
  if ( model_list == 0 ) model_list = new std::vector<G4HadronicInteraction*>;
  model_list->push_back(theIonBC);

  theBuilder = new G4FTFBuilder("FTFP",thePreCompound);
  theFTFP = theBuilder->GetModel();
  model_list->push_back(theFTFP);

  theQMD= new G4QMDReaction();
  model_list->push_back(theQMD);

  theNuclNuclData = new G4CrossSectionInelastic( theGGNuclNuclXS = new G4ComponentGGNuclNuclXsc() );

  AddProcess("dInelastic", G4Deuteron::Deuteron(), theIonBC, theQMD, theFTFP);
  AddProcess("tInelastic", G4Triton::Triton(), theIonBC, theQMD, theFTFP);
  AddProcess("He3Inelastic", G4He3::He3(), theIonBC, theQMD, theFTFP);
  AddProcess("alphaInelastic", G4Alpha::Alpha(), theIonBC, theQMD, theFTFP);
  AddProcess("ionInelastic", G4GenericIon::GenericIon(), theIonBC, theQMD, theFTFP);

}

Here is the call graph for this function:

Member Data Documentation

emaxFTFP

G4double G4IonQMDPhysics::emaxFTFP

private

Definition at line 96 of file G4IonQMDPhysics.hh.

emaxQMD

G4double G4IonQMDPhysics::emaxQMD

private

Definition at line 94 of file G4IonQMDPhysics.hh.

eminBIC

G4double G4IonQMDPhysics::eminBIC

private

Definition at line 91 of file G4IonQMDPhysics.hh.

eminQMD

G4double G4IonQMDPhysics::eminQMD

private

Definition at line 93 of file G4IonQMDPhysics.hh.

model_list

G4ThreadLocal std::vector< G4HadronicInteraction * > * G4IonQMDPhysics::model_list = 0

staticprivate

Definition at line 81 of file G4IonQMDPhysics.hh.

overlap

G4double G4IonQMDPhysics::overlap

private

Definition at line 98 of file G4IonQMDPhysics.hh.

p_list

G4ThreadLocal std::vector< G4HadronInelasticProcess * > * G4IonQMDPhysics::p_list = 0

staticprivate

Definition at line 80 of file G4IonQMDPhysics.hh.

theBuilder

G4ThreadLocal G4FTFBuilder * G4IonQMDPhysics::theBuilder = 0

staticprivate

Definition at line 88 of file G4IonQMDPhysics.hh.

theFTFP

G4ThreadLocal G4HadronicInteraction * G4IonQMDPhysics::theFTFP = 0

staticprivate

Definition at line 87 of file G4IonQMDPhysics.hh.

theGGNuclNuclXS

G4ThreadLocal G4VComponentCrossSection * G4IonQMDPhysics::theGGNuclNuclXS = 0

staticprivate

Definition at line 84 of file G4IonQMDPhysics.hh.

theIonBC

G4ThreadLocal G4BinaryLightIonReaction * G4IonQMDPhysics::theIonBC = 0

staticprivate

Definition at line 86 of file G4IonQMDPhysics.hh.

theNuclNuclData

G4ThreadLocal G4VCrossSectionDataSet * G4IonQMDPhysics::theNuclNuclData = 0

staticprivate

Definition at line 83 of file G4IonQMDPhysics.hh.

theQMD

G4ThreadLocal G4QMDReaction * G4IonQMDPhysics::theQMD = 0

staticprivate

Definition at line 89 of file G4IonQMDPhysics.hh.

verbose

G4int G4IonQMDPhysics::verbose

private

Definition at line 100 of file G4IonQMDPhysics.hh.

wasActivated

G4ThreadLocal G4bool G4IonQMDPhysics::wasActivated = false

staticprivate

Definition at line 101 of file G4IonQMDPhysics.hh.

---

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

• G4IonQMDPhysics.hh
• G4IonQMDPhysics.cc