#include <GammaRayTelHadronPhysics.hh>

Inheritance diagram for GammaRayTelHadronPhysics:

Public Member Functions
	GammaRayTelHadronPhysics (const G4String &name="hadron")

virtual	~GammaRayTelHadronPhysics ()

void	ConstructParticle ()

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

Protected Attributes
G4HadronElasticProcess	theElasticProcess

G4HadronElastic *	theElasticModel

G4PionPlusInelasticProcess	thePionPlusInelastic

G4hMultipleScattering	thePionPlusMult

G4hIonisation	thePionPlusIonisation

G4PionMinusInelasticProcess	thePionMinusInelastic

G4hMultipleScattering	thePionMinusMult

G4hIonisation	thePionMinusIonisation

G4PiMinusAbsorptionBertini	thePionMinusAbsorption

G4TheoFSGenerator *	theModel

G4ExcitationHandler	theHandler

G4PreCompoundModel *	thePreEquilib

G4GeneratorPrecompoundInterface *	theCascade

G4FTFModel *	theStringModel

G4QGSMFragmentation	theFragmentation

G4ExcitedStringDecay *	theStringDecay

G4KaonPlusInelasticProcess	theKaonPlusInelastic

G4hMultipleScattering	theKaonPlusMult

G4hIonisation	theKaonPlusIonisation

G4KaonMinusInelasticProcess	theKaonMinusInelastic

G4hMultipleScattering	theKaonMinusMult

G4hIonisation	theKaonMinusIonisation

G4KaonMinusAbsorptionBertini	theKaonMinusAbsorption

G4KaonZeroLInelasticProcess	theKaonZeroLInelastic

G4KaonZeroSInelasticProcess	theKaonZeroSInelastic

G4ProtonInelasticProcess	theProtonInelastic

G4hMultipleScattering	theProtonMult

G4hIonisation	theProtonIonisation

G4AntiProtonInelasticProcess	theAntiProtonInelastic

G4hMultipleScattering	theAntiProtonMult

G4hIonisation	theAntiProtonIonisation

G4AntiProtonAbsorptionFritiof	theAntiProtonAnnihilation

G4NeutronInelasticProcess	theNeutronInelastic

G4HadronFissionProcess	theNeutronFission

G4LFission *	theNeutronFissionModel

G4HadronCaptureProcess *	theNeutronCapture

G4AntiNeutronInelasticProcess	theAntiNeutronInelastic

Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel

G4String	namePhysics

G4int	typePhysics

G4ParticleTable *	theParticleTable

G4int	g4vpcInstanceID

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)

Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 95 of file GammaRayTelHadronPhysics.hh.

Constructor & Destructor Documentation

GammaRayTelHadronPhysics::GammaRayTelHadronPhysics ( const G4String & name = "hadron" )

Definition at line 39 of file GammaRayTelHadronPhysics.cc.

40 : G4VPhysicsConstructor(name)

41 {;}

G4VPhysicsConstructor::G4VPhysicsConstructor

G4VPhysicsConstructor(const G4String &="")

Definition: G4VPhysicsConstructor.cc:60

GammaRayTelHadronPhysics::~GammaRayTelHadronPhysics ( )

virtual

Definition at line 43 of file GammaRayTelHadronPhysics.cc.

References theStringDecay.

 {
   delete theStringDecay;
 }

Member Function Documentation

void GammaRayTelHadronPhysics::ConstructParticle ( void )

inlinevirtual

Implements G4VPhysicsConstructor.

Definition at line 104 of file GammaRayTelHadronPhysics.hh.

104 {};

void GammaRayTelHadronPhysics::ConstructProcess ( void )

virtual

Implements G4VPhysicsConstructor.

Definition at line 52 of file GammaRayTelHadronPhysics.cc.

References G4HadronicProcess::AddDataSet(), G4ProcessManager::AddDiscreteProcess(), G4ProcessManager::AddProcess(), G4ProcessManager::AddRestProcess(), G4AntiNeutron::AntiNeutron(), G4AntiProton::AntiProton(), G4ChipsKaonMinusInelasticXS::Default_Name(), G4ChipsKaonPlusInelasticXS::Default_Name(), G4ChipsKaonZeroInelasticXS::Default_Name(), python.hepunit::eV, G4ParticleDefinition::GetProcessManager(), python.hepunit::GeV, idxAlongStep, idxPostStep, G4CrossSectionDataSetRegistry::Instance(), G4InuclParticleNames::kaonMinus, G4KaonMinus::KaonMinusDefinition(), G4InuclParticleNames::kaonPlus, G4KaonPlus::KaonPlusDefinition(), G4KaonZeroLong::KaonZeroLong(), G4KaonZeroShort::KaonZeroShort(), G4Neutron::Neutron(), ordDefault, ordInActive, G4InuclParticleNames::pion(), G4InuclParticleNames::pionMinus, G4PionMinus::PionMinusDefinition(), G4PionPlus::PionPlusDefinition(), G4Proton::Proton(), G4HadronicProcess::RegisterMe(), G4VPartonStringModel::SetFragmentationModel(), G4TheoFSGenerator::SetHighEnergyGenerator(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), G4ProcessManager::SetProcessOrdering(), G4TheoFSGenerator::SetTransport(), python.hepunit::TeV, theAntiNeutronInelastic, theAntiProtonAnnihilation, theAntiProtonInelastic, theAntiProtonIonisation, theAntiProtonMult, theCascade, theElasticModel, theElasticProcess, theKaonMinusAbsorption, theKaonMinusInelastic, theKaonMinusIonisation, theKaonMinusMult, theKaonPlusInelastic, theKaonPlusIonisation, theKaonPlusMult, theKaonZeroLInelastic, theKaonZeroSInelastic, theModel, theNeutronCapture, theNeutronFission, theNeutronFissionModel, theNeutronInelastic, thePionMinusAbsorption, thePionMinusInelastic, thePionMinusIonisation, thePionMinusMult, thePionPlusInelastic, thePionPlusIonisation, thePionPlusMult, thePreEquilib, theProtonInelastic, theProtonIonisation, theProtonMult, theStringDecay, and theStringModel.

 {
 
   G4ProcessManager * pManager = 0;
   /*
   G4cout << "" << G4endl;
   G4cout << "You are using the GammaRayTelHadronPhysics" << G4endl;
   G4cout << " - Note that this hadronic physics list is not optimized for any particular usage" << G4endl;
   G4cout << " - If you wish to have a starting point tailored for a particular area of work," << G4endl;
   G4cout << "   please use one of the available physics lists by use-case." << G4endl;
   G4cout << "" << G4endl;
   */
 
   // Elastic Process
   theElasticModel = new G4HadronElastic();
   theElasticProcess.RegisterMe(theElasticModel);
 
   const G4double theBERTMin =   0.0*GeV;
   const G4double theBERTMax =   5.0*GeV;
   const G4double theFTFMin =    4.0*GeV;
   const G4double theFTFMax =  100.0*TeV;
   
   theStringModel = new G4FTFModel;
   theStringDecay = new G4ExcitedStringDecay( new G4LundStringFragmentation );
   theStringModel->SetFragmentationModel( theStringDecay );
   thePreEquilib = new G4PreCompoundModel( new G4ExcitationHandler );
   theCascade = new G4GeneratorPrecompoundInterface( thePreEquilib );
 
   theModel = new G4TheoFSGenerator( "FTFP" );
   theModel->SetHighEnergyGenerator( theStringModel );
   theModel->SetTransport( theCascade );
   theModel->SetMinEnergy( theFTFMin );
   theModel->SetMaxEnergy( theFTFMax );
 
   G4TheoFSGenerator* theModelDownToZero =  new G4TheoFSGenerator( "FTFP" );
   theModelDownToZero->SetHighEnergyGenerator( theStringModel );
   theModelDownToZero->SetTransport( theCascade );
   theModelDownToZero->SetMinEnergy(0*eV);
   theModelDownToZero->SetMaxEnergy(theFTFMax );
 
   G4CascadeInterface * theBERTModel = new G4CascadeInterface;
   theBERTModel->SetMinEnergy( theBERTMin );
   theBERTModel->SetMaxEnergy( theBERTMax );
 
   // pi+ and pi-      
 
   G4VCrossSectionDataSet * thePiData = new G4CrossSectionPairGG( new G4PiNuclearCrossSection, 91*GeV );
   G4VCrossSectionDataSet * theAntiNucleonData = new G4CrossSectionInelastic( new G4ComponentAntiNuclNuclearXS );
 
 
   // PionPlus
   G4ParticleDefinition* pion = G4PionPlus::PionPlusDefinition();
   pManager = pion ->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   thePionPlusInelastic.AddDataSet(thePiData);
   thePionPlusInelastic.RegisterMe(theBERTModel);
   thePionPlusInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&thePionPlusInelastic);
 
   pManager->AddProcess(&thePionPlusIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&thePionPlusMult);
   pManager->SetProcessOrdering(&thePionPlusMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&thePionPlusMult, idxPostStep, 1);
 
   // PionMinus
  G4ParticleDefinition* pionMinus = G4PionMinus::PionMinusDefinition();
   pManager = pionMinus -> GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   thePionMinusInelastic.AddDataSet(thePiData);
   thePionMinusInelastic.RegisterMe(theBERTModel);
   thePionMinusInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&thePionMinusInelastic);
 
   pManager->AddProcess(&thePionMinusIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&thePionMinusMult);
   pManager->SetProcessOrdering(&thePionMinusMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&thePionMinusMult, idxPostStep, 1);
 
   pManager->AddRestProcess(&thePionMinusAbsorption, ordDefault);
 
   // KaonPlus
   G4ParticleDefinition* kaonPlus = G4KaonPlus::KaonPlusDefinition();
   pManager = kaonPlus->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theKaonPlusInelastic.AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
                    GetCrossSectionDataSet(G4ChipsKaonPlusInelasticXS::Default_Name()));
   theKaonPlusInelastic.RegisterMe(theBERTModel);
   theKaonPlusInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&theKaonPlusInelastic);
 
   pManager->AddProcess(&theKaonPlusIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&theKaonPlusMult);
   pManager->SetProcessOrdering(&theKaonPlusMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&theKaonPlusMult, idxPostStep, 1);
 
   // KaonMinus
   G4ParticleDefinition* kaonMinus = G4KaonMinus::KaonMinusDefinition();
   pManager = kaonMinus->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theKaonMinusInelastic.AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
                    GetCrossSectionDataSet(G4ChipsKaonMinusInelasticXS::Default_Name()));
   theKaonMinusInelastic.RegisterMe(theBERTModel);
   theKaonMinusInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&theKaonMinusInelastic);
 
   pManager->AddProcess(&theKaonMinusIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&theKaonMinusMult);
   pManager->SetProcessOrdering(&theKaonMinusMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&theKaonMinusMult, idxPostStep, 1);
 
   pManager->AddRestProcess(&theKaonMinusAbsorption, ordDefault);
 
   // KaonZeroL
   pManager = G4KaonZeroLong::KaonZeroLong()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theKaonZeroLInelastic.AddDataSet( G4CrossSectionDataSetRegistry::Instance()
                     ->GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()));
   theKaonZeroLInelastic.RegisterMe(theBERTModel);
   theKaonZeroLInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&theKaonZeroLInelastic);
  
   // KaonZeroS
   pManager = G4KaonZeroShort::KaonZeroShort()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theKaonZeroSInelastic.AddDataSet( G4CrossSectionDataSetRegistry::Instance()
                     ->GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()));
   theKaonZeroSInelastic.RegisterMe(theBERTModel);
   theKaonZeroSInelastic.RegisterMe(theModel); 
   pManager->AddDiscreteProcess(&theKaonZeroSInelastic);
 
   // Proton
   pManager = G4Proton::Proton()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theProtonInelastic.AddDataSet(new G4BGGNucleonInelasticXS( G4Proton::Proton() ) );
   theProtonInelastic.RegisterMe(theBERTModel);
   theProtonInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&theProtonInelastic);
 
   pManager->AddProcess(&theProtonIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&theProtonMult);
   pManager->SetProcessOrdering(&theProtonMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&theProtonMult, idxPostStep, 1);
 
   // anti-Proton
   pManager = G4AntiProton::AntiProton()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theAntiProtonInelastic.AddDataSet( theAntiNucleonData );
   theAntiProtonInelastic.RegisterMe(theModelDownToZero);
   pManager->AddDiscreteProcess(&theAntiProtonInelastic);
 
   pManager->AddProcess(&theAntiProtonIonisation, ordInActive,2, 2);
 
   pManager->AddProcess(&theAntiProtonMult);
   pManager->SetProcessOrdering(&theAntiProtonMult, idxAlongStep, 1);
   pManager->SetProcessOrdering(&theAntiProtonMult, idxPostStep, 1);
 
   pManager->AddRestProcess(&theAntiProtonAnnihilation);
 
   // Neutron
   pManager = G4Neutron::Neutron()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);
   theNeutronInelastic.AddDataSet( new G4BGGNucleonInelasticXS( G4Neutron::Neutron() ));
   theNeutronInelastic.RegisterMe(theBERTModel);
   theNeutronInelastic.RegisterMe(theModel);
   pManager->AddDiscreteProcess(&theNeutronInelastic);
   
   theNeutronFissionModel = new G4LFission();
   theNeutronFission.RegisterMe(theNeutronFissionModel);
   pManager->AddDiscreteProcess(&theNeutronFission);
 
   theNeutronCapture = new G4HadronCaptureProcess();
   theNeutronCapture->AddDataSet(new G4NeutronCaptureXS()); 
   pManager->AddDiscreteProcess(theNeutronCapture);
 
   // AntiNeutron
   pManager = G4AntiNeutron::AntiNeutron()->GetProcessManager();
   // add process
   pManager->AddDiscreteProcess(&theElasticProcess);  
   theAntiNeutronInelastic.AddDataSet( theAntiNucleonData );
   theAntiNeutronInelastic.RegisterMe(theModelDownToZero);
   pManager->AddDiscreteProcess(&theAntiNeutronInelastic);
 
 }