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 ()
virtual void	ConstructProcess ()
Public Member Functions inherited from G4VHadronPhysics
	G4VHadronPhysics (const G4String &name="hInelastic", G4int verbose=0)
virtual	~G4VHadronPhysics ()
virtual void	ConstructParticle ()
G4HadronicInteraction *	BuildModel (G4VHadronModelBuilder *, G4double emin, G4double emax)
G4HadronicInteraction *	NewModel (G4HadronicInteraction *, G4double emin, G4double emax)
void	AddInelasticCrossSection (const G4String &, G4VCrossSectionDataSet *)
void	AddInelasticCrossSection (const G4ParticleDefinition *, G4VCrossSectionDataSet *)
void	AddElasticCrossSection (const G4String &, G4VCrossSectionDataSet *)
void	AddElasticCrossSection (const G4ParticleDefinition *, G4VCrossSectionDataSet *)
void	AddCaptureCrossSection (G4VCrossSectionDataSet *)
void	AddFissionCrossSection (G4VCrossSectionDataSet *)
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
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const

Additional Inherited Members
Protected Member Functions inherited from G4VHadronPhysics
G4HadronicProcess *	FindInelasticProcess (const G4String &)
G4HadronicProcess *	FindInelasticProcess (const G4ParticleDefinition *)
G4HadronicProcess *	FindElasticProcess (const G4String &)
G4HadronicProcess *	FindElasticProcess (const G4ParticleDefinition *)
G4HadronicProcess *	FindCaptureProcess ()
G4HadronicProcess *	FindFissionProcess ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4ParticleTable::G4PTblDicIterator *	theParticleIterator
G4PhysicsListHelper *	thePLHelper

Detailed Description

Definition at line 47 of file G4HadronInelasticQBBC.hh.

Constructor & Destructor Documentation

G4HadronInelasticQBBC() [1/2]

G4HadronInelasticQBBC::G4HadronInelasticQBBC

(

G4int

ver = 1

)

Definition at line 90 of file G4HadronInelasticQBBC.cc.

  : G4VHadronPhysics("hInelastic"),verbose(ver),wasActivated(false)
{
  htype = "QBBC";
  theAntiNuclXS = 0;
}

G4HadronInelasticQBBC() [2/2]

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 97 of file G4HadronInelasticQBBC.cc.

  : G4VHadronPhysics("hInelastic"),verbose(ver),wasActivated(false)
{
  htype = name;
  theAntiNuclXS = 0;
}

~G4HadronInelasticQBBC()

G4HadronInelasticQBBC::~G4HadronInelasticQBBC ( )

virtual

Definition at line 105 of file G4HadronInelasticQBBC.cc.

{
  delete theAntiNuclXS;
}

Member Function Documentation

ConstructProcess()

void G4HadronInelasticQBBC::ConstructProcess ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 110 of file G4HadronInelasticQBBC.cc.

{
  if(wasActivated) return;
  wasActivated = true;
 
  if(verbose > 1) {
    G4cout << "### HadronInelasticQBBC Construct Process with type <"
       << htype << ">" << G4endl;
  }
 
  G4double emax = 100.*TeV;
 
 
  //G4cout << "G4HadronInelasticQBBC::ConstructProcess new PRECO"<< G4endl;
 
  // PreCompound and Evaporation models are instantiated here
  G4PreCompoundModel* thePreCompound = 0;
  G4HadronicInteraction* p =
    G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
  thePreCompound = static_cast<G4PreCompoundModel*>(p);
  if(!thePreCompound) { thePreCompound = new G4PreCompoundModel(); }
  //G4ExcitationHandler* handler = thePreCompound->GetExcitationHandler();
 
  // configure models
  //G4HadronicInteraction* theQGSP = 
  //  BuildModel(new G4QGSBuilder("QGSP",thePreCompound,true,false),12.5*GeV,emax);
  G4HadronicInteraction* theFTFP = 
    BuildModel(new G4FTFBuilder("FTFP",thePreCompound),3.0*GeV,emax);
  G4HadronicInteraction* theFTFP1 = 
    BuildModel(new G4FTFBuilder("FTFP",thePreCompound),3.0*GeV,emax);
  G4HadronicInteraction* theFTFP2 = 
    BuildModel(new G4FTFBuilder("FTFP",thePreCompound),0.0,emax);
 
  G4HadronicInteraction* theBERT = 
    NewModel(new G4CascadeInterface(),1.0*GeV,12.0*GeV);
  G4HadronicInteraction* theBERT1 = 
    NewModel(new G4CascadeInterface(),0.0*GeV,12.0*GeV);
 
  //G4cout << "G4HadronInelasticQBBC::ConstructProcess new Binary"<< G4endl;
  G4BinaryCascade* bic = new G4BinaryCascade(thePreCompound);
  G4HadronicInteraction* theBIC = NewModel(bic,0.0,1.5*GeV);
 
  G4QInelastic* theCHIPS = new G4QInelastic();
  G4HadronicProcessStore* store = G4HadronicProcessStore::Instance();
  store->RegisterExtraProcess(theCHIPS);
 
  // cross sections
  theAntiNuclXS = new G4ComponentAntiNuclNuclearXS();
  G4CrossSectionInelastic* anucxs = 
    new G4CrossSectionInelastic(theAntiNuclXS);
 
  // loop over particles
  theParticleIterator->reset();
  while( (*theParticleIterator)() ) {
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4String pname = particle->GetParticleName();
    //G4ProcessManager* pmanager = particle->GetProcessManager();
    if(verbose > 1) { 
      G4cout << "### HadronInelasticQBBC:  " << pname << G4endl; 
    }
 
    //
    // model and X-section configuration per particle type
    //
    if(pname == "proton") {
      G4HadronicProcess* hp = FindInelasticProcess(particle);
      hp->AddDataSet(new G4BGGNucleonInelasticXS(particle));
      
      //hp->RegisterMe(theQGSP);
      hp->RegisterMe(theFTFP);
      hp->RegisterMe(theBERT);
      hp->RegisterMe(theBIC);
 
    } else if(pname == "neutron") {
      G4HadronicProcess* hp = FindInelasticProcess(particle);
      hp->AddDataSet(new G4NeutronInelasticXS());
      //hp->RegisterMe(theQGSP);
      hp->RegisterMe(theFTFP);
       
      G4HadronicProcess* capture = FindCaptureProcess();
      capture->AddDataSet(new G4NeutronCaptureXS());
      hp->RegisterMe(theBERT);
      hp->RegisterMe(theBIC);
      capture->RegisterMe(new G4NeutronRadCapture());
 
    } else if(pname == "pi-" || pname == "pi+") {
      G4HadronicProcess* hp = FindInelasticProcess(particle);
//      hp->AddDataSet(new G4BGGPionInelasticXS(particle));
      hp->AddDataSet(new G4CrossSectionPairGG(new G4PiNuclearCrossSection(), 91*GeV));
      //hp->RegisterMe(theQGSP);
      hp->RegisterMe(theFTFP);
      hp->RegisterMe(theBERT1);
 
    } else if(pname == "kaon-" ) {
      G4HadronicProcess* hp = FindInelasticProcess(particle);
      hp->RegisterMe(theFTFP1);
      hp->RegisterMe(theBERT1);
      hp->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsKaonMinusInelasticXS::Default_Name()));
 
    } else if(pname == "kaon+" ) {
        G4HadronicProcess* hp = FindInelasticProcess(particle);
        hp->RegisterMe(theFTFP1);
        hp->RegisterMe(theBERT1);
        hp->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsKaonPlusInelasticXS::Default_Name()));
 
    } else if(pname == "kaon0S"    ||
              pname == "kaon0L") {
        G4HadronicProcess* hp = FindInelasticProcess(particle);
        hp->RegisterMe(theFTFP1);
        hp->RegisterMe(theBERT1);
        hp->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()));
        
    } else if(pname == "lambda"    ||
              pname == "omega-"    ||
          pname == "sigma-"    || 
          pname == "sigma+"    || 
          pname == "sigma0"    || 
          pname == "xi-"       || 
          pname == "xi0") {
      G4HadronicProcess* hp = FindInelasticProcess(particle);
      hp->RegisterMe(theFTFP1);
      hp->RegisterMe(theBERT1);
      hp->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsHyperonInelasticXS::Default_Name()));
 
    } else if(pname == "anti_alpha"   ||
          pname == "anti_deuteron"||
              pname == "anti_He3"     ||
          pname == "anti_proton"  || 
              pname == "anti_triton"  ||  
          pname == "anti_lambda"  ||
              pname == "anti_neutron" ||
              pname == "anti_omega-"  || 
              pname == "anti_sigma-"  || 
              pname == "anti_sigma+"  || 
              pname == "anti_xi-"     || 
              pname == "anti_xi0"     
          ) {
 
      G4HadronicProcess* hp = FindInelasticProcess(particle);
      hp->RegisterMe(theFTFP2);
      hp->AddDataSet(anucxs);
 
      //pmanager->AddDiscreteProcess(theCHIPS);
      //store->RegisterParticleForExtraProcess(theCHIPS,particle);
    } 
  }
}

---

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

• G4HadronInelasticQBBC.hh
• G4HadronInelasticQBBC.cc