OpNovicePhysicsList Class Reference

#include <OpNovicePhysicsList.hh>

Inheritance diagram for OpNovicePhysicsList:

Collaboration diagram for OpNovicePhysicsList:

Public Member Functions
	OpNovicePhysicsList ()
virtual	~OpNovicePhysicsList ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
virtual void	SetCuts ()
void	ConstructDecay ()
void	ConstructEM ()
void	ConstructOp ()
void	SetVerbose (G4int)
void	SetNbOfPhotonsCerenkov (G4int)
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 ()

Private Attributes
OpNovicePhysicsListMessenger *	fMessenger

Static Private Attributes
static G4ThreadLocal G4int	fVerboseLevel = 1
static G4ThreadLocal G4int	fMaxNumPhotonStep = 20
static G4ThreadLocal G4Cerenkov *	fCerenkovProcess = 0
static G4ThreadLocal G4Scintillation *	fScintillationProcess = 0
static G4ThreadLocal G4OpAbsorption *	fAbsorptionProcess = 0
static G4ThreadLocal G4OpRayleigh *	fRayleighScatteringProcess = 0
static G4ThreadLocal G4OpMieHG *	fMieHGScatteringProcess = 0
static G4ThreadLocal G4OpBoundaryProcess *	fBoundaryProcess = 0

Additional Inherited Members
Static Public Member Functions inherited from G4VUserPhysicsList
static const G4VUPLManager &	GetSubInstanceManager ()
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 ()
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
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 G4VUserPhysicsList
static G4RUN_DLL G4VUPLManager	subInstanceManager

Detailed Description

Definition at line 48 of file OpNovicePhysicsList.hh.

Constructor & Destructor Documentation

OpNovicePhysicsList()

OpNovicePhysicsList::OpNovicePhysicsList

(

)

Definition at line 69 of file OpNovicePhysicsList.cc.

 : G4VUserPhysicsList()
{
  fMessenger = new OpNovicePhysicsListMessenger(this);
}

~OpNovicePhysicsList()

OpNovicePhysicsList::~OpNovicePhysicsList ( )

virtual

Definition at line 77 of file OpNovicePhysicsList.cc.

{ delete fMessenger; }

Member Function Documentation

ConstructDecay()

void OpNovicePhysicsList::ConstructDecay

(

)

Definition at line 120 of file OpNovicePhysicsList.cc.

{
  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    if (theDecayProcess->IsApplicable(*particle)) {
      pmanager ->AddProcess(theDecayProcess);
      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructEM()

void OpNovicePhysicsList::ConstructEM

(

)

Definition at line 160 of file OpNovicePhysicsList.cc.

{
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();

    if (particleName == "gamma") {
    // gamma
      // Construct processes for gamma
      pmanager->AddDiscreteProcess(new G4GammaConversion());
      pmanager->AddDiscreteProcess(new G4ComptonScattering());
      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());

    } else if (particleName == "e-") {
    //electron
      // Construct processes for electron
      pmanager->AddProcess(new G4eMultipleScattering(),-1, 1, 1);
      pmanager->AddProcess(new G4eIonisation(),       -1, 2, 2);
      pmanager->AddProcess(new G4eBremsstrahlung(),   -1, 3, 3);

    } else if (particleName == "e+") {
    //positron
      // Construct processes for positron
      pmanager->AddProcess(new G4eMultipleScattering(),-1, 1, 1);
      pmanager->AddProcess(new G4eIonisation(),       -1, 2, 2);
      pmanager->AddProcess(new G4eBremsstrahlung(),   -1, 3, 3);
      pmanager->AddProcess(new G4eplusAnnihilation(),  0,-1, 4);

    } else if( particleName == "mu+" ||
               particleName == "mu-"    ) {
    //muon
     // Construct processes for muon
     pmanager->AddProcess(new G4MuMultipleScattering(),-1, 1, 1);
     pmanager->AddProcess(new G4MuIonisation(),      -1, 2, 2);
     pmanager->AddProcess(new G4MuBremsstrahlung(),  -1, 3, 3);
     pmanager->AddProcess(new G4MuPairProduction(),  -1, 4, 4);

    } else {
      if ((particle->GetPDGCharge() != 0.0) &&
          (particle->GetParticleName() != "chargedgeantino") &&
          !particle->IsShortLived()) {
       // all others charged particles except geantino
       pmanager->AddProcess(new G4hMultipleScattering(),-1,1,1);
       pmanager->AddProcess(new G4hIonisation(),       -1,2,2);
     }
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructOp()

void OpNovicePhysicsList::ConstructOp

(

)

Definition at line 215 of file OpNovicePhysicsList.cc.

{
  fCerenkovProcess = new G4Cerenkov("Cerenkov");
  fCerenkovProcess->SetMaxNumPhotonsPerStep(fMaxNumPhotonStep);
  fCerenkovProcess->SetMaxBetaChangePerStep(10.0);
  fCerenkovProcess->SetTrackSecondariesFirst(true);
  fScintillationProcess = new G4Scintillation("Scintillation");
  fScintillationProcess->SetScintillationYieldFactor(1.);
  fScintillationProcess->SetTrackSecondariesFirst(true);
  fAbsorptionProcess = new G4OpAbsorption();
  fRayleighScatteringProcess = new G4OpRayleigh();
  fMieHGScatteringProcess = new G4OpMieHG();
  fBoundaryProcess = new G4OpBoundaryProcess();

  fCerenkovProcess->SetVerboseLevel(fVerboseLevel);
  fScintillationProcess->SetVerboseLevel(fVerboseLevel);
  fAbsorptionProcess->SetVerboseLevel(fVerboseLevel);
  fRayleighScatteringProcess->SetVerboseLevel(fVerboseLevel);
  fMieHGScatteringProcess->SetVerboseLevel(fVerboseLevel);
  fBoundaryProcess->SetVerboseLevel(fVerboseLevel);
  
  // Use Birks Correction in the Scintillation process
  if(G4Threading::IsMasterThread())
  {
    G4EmSaturation* emSaturation =
              G4LossTableManager::Instance()->EmSaturation();
      fScintillationProcess->AddSaturation(emSaturation);
  }

  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();
    if (fCerenkovProcess->IsApplicable(*particle)) {
      pmanager->AddProcess(fCerenkovProcess);
      pmanager->SetProcessOrdering(fCerenkovProcess,idxPostStep);
    }
    if (fScintillationProcess->IsApplicable(*particle)) {
      pmanager->AddProcess(fScintillationProcess);
      pmanager->SetProcessOrderingToLast(fScintillationProcess, idxAtRest);
      pmanager->SetProcessOrderingToLast(fScintillationProcess, idxPostStep);
    }
    if (particleName == "opticalphoton") {
      G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl;
      pmanager->AddDiscreteProcess(fAbsorptionProcess);
      pmanager->AddDiscreteProcess(fRayleighScatteringProcess);
      pmanager->AddDiscreteProcess(fMieHGScatteringProcess);
      pmanager->AddDiscreteProcess(fBoundaryProcess);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructParticle()

void OpNovicePhysicsList::ConstructParticle ( void  )

virtual

Implements G4VUserPhysicsList.

Definition at line 81 of file OpNovicePhysicsList.cc.

{
  // In this method, static member functions should be called
  // for all particles which you want to use.
  // This ensures that objects of these particle types will be
  // created in the program.

  G4BosonConstructor bConstructor;
  bConstructor.ConstructParticle();

  G4LeptonConstructor lConstructor;
  lConstructor.ConstructParticle();

  G4MesonConstructor mConstructor;
  mConstructor.ConstructParticle();

  G4BaryonConstructor rConstructor;
  rConstructor.ConstructParticle();

  G4IonConstructor iConstructor;
  iConstructor.ConstructParticle(); 
}

Here is the call graph for this function:

ConstructProcess()

void OpNovicePhysicsList::ConstructProcess ( void  )

virtual

Implements G4VUserPhysicsList.

Definition at line 106 of file OpNovicePhysicsList.cc.

{
  AddTransportation();
  ConstructDecay();
  ConstructEM();
  ConstructOp();
}

Here is the call graph for this function:

SetCuts()

void OpNovicePhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 294 of file OpNovicePhysicsList.cc.

{
  //  " G4VUserPhysicsList::SetCutsWithDefault" method sets
  //   the default cut value for all particle types
  //
  SetCutsWithDefault();

  if (verboseLevel>0) DumpCutValuesTable();
}

Here is the call graph for this function:

SetNbOfPhotonsCerenkov()

void OpNovicePhysicsList::SetNbOfPhotonsCerenkov

(

G4int

MaxNumber

)

Definition at line 285 of file OpNovicePhysicsList.cc.

{
  fMaxNumPhotonStep = MaxNumber;

  fCerenkovProcess->SetMaxNumPhotonsPerStep(fMaxNumPhotonStep);
}

Here is the call graph for this function:

SetVerbose()

void OpNovicePhysicsList::SetVerbose

(

G4int

verbose

)

Definition at line 271 of file OpNovicePhysicsList.cc.

{
  fVerboseLevel = verbose;

  fCerenkovProcess->SetVerboseLevel(fVerboseLevel);
  fScintillationProcess->SetVerboseLevel(fVerboseLevel);
  fAbsorptionProcess->SetVerboseLevel(fVerboseLevel);
  fRayleighScatteringProcess->SetVerboseLevel(fVerboseLevel);
  fMieHGScatteringProcess->SetVerboseLevel(fVerboseLevel);
  fBoundaryProcess->SetVerboseLevel(fVerboseLevel);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

fAbsorptionProcess

G4ThreadLocal G4OpAbsorption * OpNovicePhysicsList::fAbsorptionProcess = 0

staticprivate

Definition at line 80 of file OpNovicePhysicsList.hh.

fBoundaryProcess

G4ThreadLocal G4OpBoundaryProcess * OpNovicePhysicsList::fBoundaryProcess = 0

staticprivate

Definition at line 83 of file OpNovicePhysicsList.hh.

fCerenkovProcess

G4ThreadLocal G4Cerenkov * OpNovicePhysicsList::fCerenkovProcess = 0

staticprivate

Definition at line 78 of file OpNovicePhysicsList.hh.

fMaxNumPhotonStep

G4ThreadLocal G4int OpNovicePhysicsList::fMaxNumPhotonStep = 20

staticprivate

Definition at line 76 of file OpNovicePhysicsList.hh.

fMessenger

OpNovicePhysicsListMessenger* OpNovicePhysicsList::fMessenger

private

Definition at line 73 of file OpNovicePhysicsList.hh.

fMieHGScatteringProcess

G4ThreadLocal G4OpMieHG * OpNovicePhysicsList::fMieHGScatteringProcess = 0

staticprivate

Definition at line 82 of file OpNovicePhysicsList.hh.

fRayleighScatteringProcess

G4ThreadLocal G4OpRayleigh * OpNovicePhysicsList::fRayleighScatteringProcess = 0

staticprivate

Definition at line 81 of file OpNovicePhysicsList.hh.

fScintillationProcess

G4ThreadLocal G4Scintillation * OpNovicePhysicsList::fScintillationProcess = 0

staticprivate

Definition at line 79 of file OpNovicePhysicsList.hh.

fVerboseLevel

G4ThreadLocal G4int OpNovicePhysicsList::fVerboseLevel = 1

staticprivate

Definition at line 75 of file OpNovicePhysicsList.hh.

---

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

• OpNovicePhysicsList.hh
• OpNovicePhysicsList.cc