ExN03PhysicsList Class Reference

#include <ExN03PhysicsList.hh>

Inheritance diagram for ExN03PhysicsList:

Public Member Functions
	ExN03PhysicsList ()
	~ExN03PhysicsList ()

Protected Member Functions
void	ConstructParticle ()
void	ConstructProcess ()
void	SetCuts ()
void	ConstructBosons ()
void	ConstructLeptons ()
void	ConstructMesons ()
void	ConstructBaryons ()
void	ConstructGeneral ()
void	ConstructEM ()

Detailed Description

Definition at line 40 of file ExN03PhysicsList.hh.

Constructor & Destructor Documentation

ExN03PhysicsList()

ExN03PhysicsList::ExN03PhysicsList

(

)

Definition at line 45 of file ExN03PhysicsList.cxx.

                                   : G4VUserPhysicsList() {
  defaultCutValue = 1.0 * mm;
  SetVerboseLevel( 1 );
}

~ExN03PhysicsList()

ExN03PhysicsList::~ExN03PhysicsList

(

)

Definition at line 52 of file ExN03PhysicsList.cxx.

{}

Member Function Documentation

ConstructBaryons()

void ExN03PhysicsList::ConstructBaryons ( )

protected

Definition at line 116 of file ExN03PhysicsList.cxx.

                                        {
  //  barions
  G4Proton::ProtonDefinition();
  G4AntiProton::AntiProtonDefinition();
  G4Neutron::NeutronDefinition();
  G4AntiNeutron::AntiNeutronDefinition();
}

Referenced by ConstructParticle().

ConstructBosons()

void ExN03PhysicsList::ConstructBosons ( )

protected

Definition at line 70 of file ExN03PhysicsList.cxx.

                                       {
  // pseudo-particles
  G4Geantino::GeantinoDefinition();
  G4ChargedGeantino::ChargedGeantinoDefinition();
 
  // gamma
  G4Gamma::GammaDefinition();
 
  // optical photon
  G4OpticalPhoton::OpticalPhotonDefinition();
}

Referenced by ConstructParticle().

ConstructEM()

void ExN03PhysicsList::ConstructEM ( )

protected

Definition at line 152 of file ExN03PhysicsList.cxx.

                                   {
  theParticleIterator->reset();
  while ( ( *theParticleIterator )() )
  {
    G4ParticleDefinition* particle     = theParticleIterator->value();
    G4ProcessManager*     pmanager     = particle->GetProcessManager();
    G4String              particleName = particle->GetParticleName();
 
    if ( particleName == "gamma" )
    {
      // gamma
      pmanager->AddDiscreteProcess( new G4GammaConversion() );
      pmanager->AddDiscreteProcess( new G4ComptonScattering() );
      pmanager->AddDiscreteProcess( new G4PhotoElectricEffect() );
    }
    else if ( particleName == "e-" )
    {
      // electron
      G4VProcess* theeminusMultipleScattering = new G4MultipleScattering();
      G4VProcess* theeminusIonisation         = new G4eIonisation();
      G4VProcess* theeminusBremsstrahlung     = new G4eBremsstrahlung();
      //
      // add processes
      pmanager->AddProcess( theeminusMultipleScattering );
      pmanager->AddProcess( theeminusIonisation );
      pmanager->AddProcess( theeminusBremsstrahlung );
      //
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering( theeminusMultipleScattering, idxAlongStep, 1 );
      pmanager->SetProcessOrdering( theeminusIonisation, idxAlongStep, 2 );
      //
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering( theeminusMultipleScattering, idxPostStep, 1 );
      pmanager->SetProcessOrdering( theeminusIonisation, idxPostStep, 2 );
      pmanager->SetProcessOrdering( theeminusBremsstrahlung, idxPostStep, 3 );
    }
    else if ( particleName == "e+" )
    {
      // positron
      G4VProcess* theeplusMultipleScattering = new G4MultipleScattering();
      G4VProcess* theeplusIonisation         = new G4eIonisation();
      G4VProcess* theeplusBremsstrahlung     = new G4eBremsstrahlung();
      G4VProcess* theeplusAnnihilation       = new G4eplusAnnihilation();
      //
      // add processes
      pmanager->AddProcess( theeplusMultipleScattering );
      pmanager->AddProcess( theeplusIonisation );
      pmanager->AddProcess( theeplusBremsstrahlung );
      pmanager->AddProcess( theeplusAnnihilation );
      //
      // set ordering for AtRestDoIt
      pmanager->SetProcessOrderingToFirst( theeplusAnnihilation, idxAtRest );
      //
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering( theeplusMultipleScattering, idxAlongStep, 1 );
      pmanager->SetProcessOrdering( theeplusIonisation, idxAlongStep, 2 );
      //
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering( theeplusMultipleScattering, idxPostStep, 1 );
      pmanager->SetProcessOrdering( theeplusIonisation, idxPostStep, 2 );
      pmanager->SetProcessOrdering( theeplusBremsstrahlung, idxPostStep, 3 );
      pmanager->SetProcessOrdering( theeplusAnnihilation, idxPostStep, 4 );
    }
    else if ( particleName == "mu+" || particleName == "mu-" )
    {
      // muon
      G4VProcess* aMultipleScattering = new G4MultipleScattering();
      G4VProcess* aBremsstrahlung     = new G4MuBremsstrahlung();
      G4VProcess* aPairProduction     = new G4MuPairProduction();
      G4VProcess* anIonisation        = new G4MuIonisation();
      //
      // add processes
      pmanager->AddProcess( anIonisation );
      pmanager->AddProcess( aMultipleScattering );
      pmanager->AddProcess( aBremsstrahlung );
      pmanager->AddProcess( aPairProduction );
      //
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering( aMultipleScattering, idxAlongStep, 1 );
      pmanager->SetProcessOrdering( anIonisation, idxAlongStep, 2 );
      //
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering( aMultipleScattering, idxPostStep, 1 );
      pmanager->SetProcessOrdering( anIonisation, idxPostStep, 2 );
      pmanager->SetProcessOrdering( aBremsstrahlung, idxPostStep, 3 );
      pmanager->SetProcessOrdering( aPairProduction, idxPostStep, 4 );
    }
    else if ( ( !particle->IsShortLived() ) && ( particle->GetPDGCharge() != 0.0 ) &&
              ( particle->GetParticleName() != "chargedgeantino" ) )
    {
      // all others charged particles except geantino
      G4VProcess* aMultipleScattering = new G4MultipleScattering();
      G4VProcess* anIonisation        = new G4hIonisation();
      //
      // add processes
      pmanager->AddProcess( anIonisation );
      pmanager->AddProcess( aMultipleScattering );
      //
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering( aMultipleScattering, idxAlongStep, 1 );
      pmanager->SetProcessOrdering( anIonisation, idxAlongStep, 2 );
      //
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering( aMultipleScattering, idxPostStep, 1 );
      pmanager->SetProcessOrdering( anIonisation, idxPostStep, 2 );
    }
  }
}

Referenced by ConstructProcess().

ConstructGeneral()

void ExN03PhysicsList::ConstructGeneral ( )

protected

Definition at line 265 of file ExN03PhysicsList.cxx.

                                        {
  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while ( ( *theParticleIterator )() )
  {
    G4ParticleDefinition* particle = theParticleIterator->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 );
    }
  }
}

Referenced by ConstructProcess().

ConstructLeptons()

void ExN03PhysicsList::ConstructLeptons ( )

protected

Definition at line 84 of file ExN03PhysicsList.cxx.

                                        {
  // leptons
  G4Electron::ElectronDefinition();
  G4Positron::PositronDefinition();
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();
 
  G4NeutrinoE::NeutrinoEDefinition();
  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
  G4NeutrinoMu::NeutrinoMuDefinition();
  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
}

Referenced by ConstructParticle().

ConstructMesons()

void ExN03PhysicsList::ConstructMesons ( )

protected

Definition at line 99 of file ExN03PhysicsList.cxx.

                                       {
  //  mesons
  G4PionPlus::PionPlusDefinition();
  G4PionMinus::PionMinusDefinition();
  G4PionZero::PionZeroDefinition();
  G4Eta::EtaDefinition();
  G4EtaPrime::EtaPrimeDefinition();
  G4KaonPlus::KaonPlusDefinition();
  G4KaonMinus::KaonMinusDefinition();
  G4KaonZero::KaonZeroDefinition();
  G4AntiKaonZero::AntiKaonZeroDefinition();
  G4KaonZeroLong::KaonZeroLongDefinition();
  G4KaonZeroShort::KaonZeroShortDefinition();
}

Referenced by ConstructParticle().

ConstructParticle()

void ExN03PhysicsList::ConstructParticle ( )

protected

Definition at line 56 of file ExN03PhysicsList.cxx.

                                         {
  // 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.
 
  ConstructBosons();
  ConstructLeptons();
  ConstructMesons();
  ConstructBaryons();
}

ConstructProcess()

void ExN03PhysicsList::ConstructProcess ( )

protected

Definition at line 126 of file ExN03PhysicsList.cxx.

                                        {
  AddTransportation();
  ConstructEM();
  ConstructGeneral();
}

SetCuts()

void ExN03PhysicsList::SetCuts ( )

protected

Definition at line 285 of file ExN03PhysicsList.cxx.

                               {
  if ( verboseLevel > 0 )
  {
    G4cout << "ExN03PhysicsList::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( defaultCutValue, "gamma" );
  SetCutValue( defaultCutValue, "e-" );
  SetCutValue( defaultCutValue, "e+" );
 
  if ( verboseLevel > 0 ) DumpCutValuesTable();
}

---

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

• ExN03PhysicsList.hh
• ExN03PhysicsList.cxx