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 ()
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
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
G4int	g4vpcInstanceID
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 47 of file G4HadronInelasticQBBC.hh.

Constructor & Destructor Documentation

G4HadronInelasticQBBC::G4HadronInelasticQBBC

(

G4int

ver = 1

)

Definition at line 89 of file G4HadronInelasticQBBC.cc.

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

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

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

G4HadronInelasticQBBC::~G4HadronInelasticQBBC ( )

virtual

Definition at line 104 of file G4HadronInelasticQBBC.cc.

{
  delete theAntiNuclXS;
}

Member Function Documentation

void G4HadronInelasticQBBC::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 109 of file G4HadronInelasticQBBC.cc.

References G4HadronicProcess::AddDataSet(), aParticleIterator, G4VHadronPhysics::BuildModel(), G4ChipsHyperonInelasticXS::Default_Name(), G4ChipsKaonMinusInelasticXS::Default_Name(), G4ChipsKaonPlusInelasticXS::Default_Name(), G4ChipsKaonZeroInelasticXS::Default_Name(), G4VHadronPhysics::FindCaptureProcess(), G4VHadronPhysics::FindInelasticProcess(), G4HadronicInteractionRegistry::FindModel(), G4cout, G4endl, G4ParticleDefinition::GetParticleName(), python.hepunit::GeV, G4HadronicInteractionRegistry::Instance(), G4CrossSectionDataSetRegistry::Instance(), G4VHadronPhysics::NewModel(), eplot::pname, G4HadronicProcess::RegisterMe(), and python.hepunit::TeV.

{
  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);

  // cross sections
  theAntiNuclXS = new G4ComponentAntiNuclNuclearXS();
  G4CrossSectionInelastic* anucxs = 
    new G4CrossSectionInelastic(theAntiNuclXS);

  // loop over particles
  aParticleIterator->reset();
  while( (*aParticleIterator)() ) {
    G4ParticleDefinition* particle = aParticleIterator->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);

    } 
  }
}

---

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

• G4HadronInelasticQBBC.hh
• G4HadronInelasticQBBC.cc