HadrontherapyPhysicsList Class Reference

#include <HadrontherapyPhysicsList.hh>

Inheritance diagram for HadrontherapyPhysicsList:

Public Member Functions
	HadrontherapyPhysicsList ()
virtual	~HadrontherapyPhysicsList ()
void	ConstructParticle ()
void	SetCuts ()
void	SetCutForGamma (G4double)
void	SetCutForElectron (G4double)
void	SetCutForPositron (G4double)
void	SetDetectorCut (G4double cut)
void	AddPhysicsList (const G4String &name)
void	ConstructProcess ()
void	AddStepMax ()
HadrontherapyStepMax *	GetStepMaxProcess ()
void	AddPackage (const G4String &name)
Public Member Functions inherited from G4VModularPhysicsList
	G4VModularPhysicsList ()
virtual	~G4VModularPhysicsList ()
void	RegisterPhysics (G4VPhysicsConstructor *)
const G4VPhysicsConstructor *	GetPhysics (G4int index) const
const G4VPhysicsConstructor *	GetPhysics (const G4String &name) const
const G4VPhysicsConstructor *	GetPhysicsWithType (G4int physics_type) const
void	ReplacePhysics (G4VPhysicsConstructor *)
void	RemovePhysics (G4VPhysicsConstructor *)
void	RemovePhysics (G4int type)
void	RemovePhysics (const G4String &name)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const
Public Member Functions inherited from G4VUserPhysicsList
	G4VUserPhysicsList ()
virtual	~G4VUserPhysicsList ()
	G4VUserPhysicsList (const G4VUserPhysicsList &)
G4VUserPhysicsList &	operator= (const G4VUserPhysicsList &)
void	Construct ()
void	UseCoupledTransportation (G4bool vl=true)
void	SetDefaultCutValue (G4double newCutValue)
G4double	GetDefaultCutValue () const
void	BuildPhysicsTable ()
void	PreparePhysicsTable (G4ParticleDefinition *)
void	BuildPhysicsTable (G4ParticleDefinition *)
G4bool	StorePhysicsTable (const G4String &directory=".")
G4bool	IsPhysicsTableRetrieved () const
G4bool	IsStoredInAscii () const
const G4String &	GetPhysicsTableDirectory () const
void	SetPhysicsTableRetrieved (const G4String &directory="")
void	SetStoredInAscii ()
void	ResetPhysicsTableRetrieved ()
void	ResetStoredInAscii ()
void	DumpList () const
void	DumpCutValuesTable (G4int flag=1)
void	DumpCutValuesTableIfRequested ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
void	SetCutsWithDefault ()
void	SetCutValue (G4double aCut, const G4String &pname)
G4double	GetCutValue (const G4String &pname) const
void	SetCutValue (G4double aCut, const G4String &pname, const G4String &rname)
void	SetParticleCuts (G4double cut, G4ParticleDefinition *particle, G4Region *region=0)
void	SetParticleCuts (G4double cut, const G4String &particleName, G4Region *region=0)
void	SetCutsForRegion (G4double aCut, const G4String &rname)
void	ResetCuts ()
	obsolete methods More...
void	SetApplyCuts (G4bool value, const G4String &name)
G4bool	GetApplyCuts (const G4String &name) const
void	RemoveProcessManager ()
void	AddProcessManager (G4ParticleDefinition *newParticle, G4ProcessManager *newManager=0)
void	CheckParticleList ()
void	DisableCheckParticleList ()
G4int	GetInstanceID () const
void	InitializeWorker ()

Additional Inherited Members
Static Public Member Functions inherited from G4VModularPhysicsList
static const G4VMPLManager &	GetSubInstanceManager ()
Static Public Member Functions inherited from G4VUserPhysicsList
static const G4VUPLManager &	GetSubInstanceManager ()
Protected Types inherited from G4VModularPhysicsList
typedef G4VMPLData::G4PhysConstVectorData	G4PhysConstVector
Protected Member Functions inherited from G4VModularPhysicsList
	G4VModularPhysicsList (const G4VModularPhysicsList &)
G4VModularPhysicsList &	operator= (const G4VModularPhysicsList &)
Protected Member Functions inherited from G4VUserPhysicsList
void	AddTransportation ()
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
void	BuildIntegralPhysicsTable (G4VProcess *, G4ParticleDefinition *)
virtual void	RetrievePhysicsTable (G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)
void	InitializeProcessManager ()
Protected Attributes inherited from G4VModularPhysicsList
G4int	verboseLevel
G4int	g4vmplInstanceID
Protected Attributes inherited from G4VUserPhysicsList
G4ParticleTable *	theParticleTable
G4int	verboseLevel
G4double	defaultCutValue
G4bool	isSetDefaultCutValue
G4ProductionCutsTable *	fCutsTable
G4bool	fRetrievePhysicsTable
G4bool	fStoredInAscii
G4bool	fIsCheckedForRetrievePhysicsTable
G4bool	fIsRestoredCutValues
G4String	directoryPhysicsTable
G4bool	fDisableCheckParticleList
G4int	g4vuplInstanceID
Static Protected Attributes inherited from G4VModularPhysicsList
static G4RUN_DLL G4VMPLManager	G4VMPLsubInstanceManager
Static Protected Attributes inherited from G4VUserPhysicsList
static G4RUN_DLL G4VUPLManager	subInstanceManager

Detailed Description

Definition at line 45 of file HadrontherapyPhysicsList.hh.

Constructor & Destructor Documentation

HadrontherapyPhysicsList::HadrontherapyPhysicsList

(

)

Definition at line 85 of file HadrontherapyPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4LossTableManager::Instance(), python.hepunit::mm, and G4VModularPhysicsList::SetVerboseLevel().

                                                   : G4VModularPhysicsList()
{
  G4LossTableManager::Instance();
  defaultCutValue = 1.*mm;
  cutForGamma     = defaultCutValue;
  cutForElectron  = defaultCutValue;
  cutForPositron  = defaultCutValue;

  helIsRegistered  = false;
  bicIsRegistered  = false;
  biciIsRegistered = false;
  locIonIonInelasticIsRegistered = false;
  radioactiveDecayIsRegistered = false;

  stepMaxProcess  = 0;

  pMessenger = new HadrontherapyPhysicsListMessenger(this);

  SetVerboseLevel(1);

  // ******     Definition of some defaults for the physics     *****
  // ******     in case no physics is called by the macro file  *****
  // EM physics
  emPhysicsList = new G4EmStandardPhysics_option3(1);
  emName = G4String("emstandard_opt3");

  // Decay physics and all particles
  decPhysicsList = new G4DecayPhysics();
  raddecayList = new G4RadioactiveDecayPhysics();
}

HadrontherapyPhysicsList::~HadrontherapyPhysicsList ( )

virtual

Definition at line 117 of file HadrontherapyPhysicsList.cc.

{
  delete pMessenger;
  delete emPhysicsList;
  delete decPhysicsList;
  delete raddecayList;

  for(size_t i=0; i<hadronPhys.size(); i++) {delete hadronPhys[i];}
}

Member Function Documentation

void HadrontherapyPhysicsList::AddPackage

(

const G4String &

name

)

void HadrontherapyPhysicsList::AddPhysicsList

(

const G4String &

name

)

Definition at line 164 of file HadrontherapyPhysicsList.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), and G4VModularPhysicsList::verboseLevel.

Referenced by HadrontherapyPhysicsListMessenger::SetNewValue().

{

  if (verboseLevel>1) {
    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
  }
  if (name == emName) return;

  /////////////////////////////////////////////////////////////////////////////
  //   ELECTROMAGNETIC MODELS
  /////////////////////////////////////////////////////////////////////////////
  if (name == "standard_opt3") {
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmStandardPhysics_option3();
    G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmStandardPhysics_option3" << G4endl;

  } else if (name == "LowE_Livermore") {
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmLivermorePhysics();
    G4RunManager::GetRunManager()-> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmLivermorePhysics" << G4endl;

  } else if (name == "LowE_Penelope") {
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmPenelopePhysics();
    G4RunManager::GetRunManager()-> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmPenelopePhysics" << G4endl;
    
  } else if (name == "local_ion_ion_inelastic") {
    hadronPhys.push_back(new LocalIonIonInelasticPhysic());
    locIonIonInelasticIsRegistered = true;

  } else if (name == "QGSP_BIC_EMY") {
    AddPhysicsList("emstandard_opt3");
    hadronPhys.push_back( new G4HadronPhysicsQGSP_BIC());
    hadronPhys.push_back( new G4EmExtraPhysics());
    hadronPhys.push_back( new G4HadronElasticPhysics());
    hadronPhys.push_back( new G4StoppingPhysics());
    hadronPhys.push_back( new G4IonBinaryCascadePhysics());
    hadronPhys.push_back( new G4NeutronTrackingCut());
    
  } else { 
    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
       << " is not defined"
       << G4endl;
  }
}

void HadrontherapyPhysicsList::AddStepMax

(

)

Definition at line 217 of file HadrontherapyPhysicsList.cc.

References G4ProcessManager::AddDiscreteProcess(), G4ParticleDefinition::GetProcessManager(), HadrontherapyStepMax::IsApplicable(), and theParticleIterator.

Referenced by ConstructProcess().

{
  // Step limitation seen as a process
  stepMaxProcess = new HadrontherapyStepMax();

  theParticleIterator->reset();
  while ((*theParticleIterator)()){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();

    if (stepMaxProcess->IsApplicable(*particle) && pmanager)
      {
    pmanager ->AddDiscreteProcess(stepMaxProcess);
      }
  }
}

void HadrontherapyPhysicsList::ConstructParticle ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 128 of file HadrontherapyPhysicsList.cc.

References G4VPhysicsConstructor::ConstructParticle().

{
  decPhysicsList->ConstructParticle();
}

void HadrontherapyPhysicsList::ConstructProcess ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 134 of file HadrontherapyPhysicsList.cc.

References G4EmConfigurator::AddModels(), AddStepMax(), G4VUserPhysicsList::AddTransportation(), G4ParallelWorldPhysics::ConstructProcess(), and G4VPhysicsConstructor::ConstructProcess().

{
  // transportation
  AddTransportation();

  // electromagnetic physics list
  emPhysicsList->ConstructProcess();
  em_config.AddModels();

  // decay physics list
  decPhysicsList->ConstructProcess();
  raddecayList->ConstructProcess();

  // hadronic physics lists
  for(size_t i=0; i<hadronPhys.size(); i++) {
    hadronPhys[i] -> ConstructProcess();
  }

  // step limitation (as a full process)
  //
  AddStepMax();

  //Parallel world sensitivity
  G4ParallelWorldPhysics* pWorld = new G4ParallelWorldPhysics("DetectorROGeometry");
  pWorld->ConstructProcess();

  return;
}

HadrontherapyStepMax* HadrontherapyPhysicsList::GetStepMaxProcess ( )

inline

Definition at line 63 of file HadrontherapyPhysicsList.hh.

{return stepMaxProcess;};

void HadrontherapyPhysicsList::SetCutForElectron

(

G4double

cut

)

Definition at line 262 of file HadrontherapyPhysicsList.cc.

References G4Electron::Electron(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by HadrontherapyPhysicsListMessenger::SetNewValue().

{
  cutForElectron = cut;
  SetParticleCuts(cutForElectron, G4Electron::Electron());
}

void HadrontherapyPhysicsList::SetCutForGamma

(

G4double

cut

)

Definition at line 255 of file HadrontherapyPhysicsList.cc.

References G4Gamma::Gamma(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by HadrontherapyPhysicsListMessenger::SetNewValue().

{
  cutForGamma = cut;
  SetParticleCuts(cutForGamma, G4Gamma::Gamma());
}

void HadrontherapyPhysicsList::SetCutForPositron

(

G4double

cut

)

Definition at line 269 of file HadrontherapyPhysicsList.cc.

References G4Positron::Positron(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by HadrontherapyPhysicsListMessenger::SetNewValue().

{
  cutForPositron = cut;
  SetParticleCuts(cutForPositron, G4Positron::Positron());
}

void HadrontherapyPhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 235 of file HadrontherapyPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4VUserPhysicsList::DumpCutValuesTable(), G4BestUnit, G4cout, G4endl, G4VUserPhysicsList::SetCutValue(), SetDetectorCut(), and G4VModularPhysicsList::verboseLevel.

{

  if (verboseLevel >0){
    G4cout << "PhysicsList::SetCuts:";
    G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
  }

  // set cut values for gamma at first and for e- second and next for e+,
  // because some processes for e+/e- need cut values for gamma
  SetCutValue(cutForGamma, "gamma");
  SetCutValue(cutForElectron, "e-");
  SetCutValue(cutForPositron, "e+");

  // Set cuts for detector
  SetDetectorCut(defaultCutValue); 
  if (verboseLevel>0) DumpCutValuesTable();
}

void HadrontherapyPhysicsList::SetDetectorCut

(

G4double

cut

)

Definition at line 275 of file HadrontherapyPhysicsList.cc.

References G4ProductionCuts::GetIndex(), G4RegionStore::GetInstance(), and G4RegionStore::GetRegion().

Referenced by SetCuts().

{
  G4String regionName = "DetectorLog";
  G4Region* region = G4RegionStore::GetInstance()->GetRegion(regionName);

  G4ProductionCuts* cuts = new G4ProductionCuts ;
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("gamma"));
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e-"));
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e+"));
  region -> SetProductionCuts(cuts);
}

---

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

• HadrontherapyPhysicsList.hh
• HadrontherapyPhysicsList.cc