#include "OpticalPhysics.hh"

Inheritance diagram for larg4::OpticalPhysics:

Public Member Functions
	OpticalPhysics (detinfo::DetectorPropertiesData const &detProp, G4int ver=0, G4String const &name="G4Optical")

Private Member Functions
void	ConstructParticle () override

void	ConstructProcess () override

Private Attributes
bool	fSimpleBoundary

G4Cerenkov *	fTheCerenkovProcess

G4Scintillation *	fTheScintillationProcess

G4OpAbsorption *	fTheAbsorptionProcess

G4OpRayleigh *	fTheRayleighScatteringProcess

OpBoundaryProcessSimple *	fTheBoundaryProcess

G4OpBoundaryProcess *	fTheBoundaryProcess_g4

G4OpWLS *	fTheWLSProcess

Detailed Description

Definition at line 78 of file OpticalPhysics.hh.

Constructor & Destructor Documentation

larg4::OpticalPhysics::OpticalPhysics	(	detinfo::DetectorPropertiesData const &	detProp,
		G4int	ver = `0`,
		G4String const &	name = `"G4Optical"`
	)

Definition at line 83 of file OpticalPhysics.cxx.

References detinfo::DetectorPropertiesData::SimpleBoundary().

     : G4VPhysicsConstructor(name), fSimpleBoundary{detProp.SimpleBoundary()}
   {
     G4LossTableManager::Instance();
     mf::LogInfo("OpticalPhysics") << "OBJECT BEING CONSTRUCTED IN OPTICAL PHYSICS";
   }

Member Function Documentation

void larg4::OpticalPhysics::ConstructParticle ( )

overrideprivate

Definition at line 93 of file OpticalPhysics.cxx.

References MF_LOG_DEBUG.

   {
     MF_LOG_DEBUG("OpticalPhysics") << "PARTICLES BEING CONSTRUCTED IN OPTICAL PHYSICS";
     // optical photon
     G4OpticalPhoton::OpticalPhotonDefinition();
 
     // gamma
     G4Gamma::Gamma();
 
     // leptons
     G4Electron::Electron();
     G4Positron::Positron();
     G4MuonPlus::MuonPlus();
     G4MuonMinus::MuonMinus();
 
     // mesons
     G4PionPlus::PionPlusDefinition();
     G4PionMinus::PionMinusDefinition();
     G4KaonPlus::KaonPlusDefinition();
     G4KaonMinus::KaonMinusDefinition();
 
     // barions
     G4Proton::Proton();
     G4AntiProton::AntiProton();
 
     // ions
     G4Deuteron::Deuteron();
     G4Triton::Triton();
     G4He3::He3();
     G4Alpha::Alpha();
     G4GenericIon::GenericIonDefinition();
   }

void larg4::OpticalPhysics::ConstructProcess ( )

overrideprivate

Definition at line 127 of file OpticalPhysics.cxx.

References detinfo::LArProperties::CerenkovLightEnabled(), fSimpleBoundary, fTheAbsorptionProcess, fTheBoundaryProcess, fTheBoundaryProcess_g4, fTheCerenkovProcess, fTheRayleighScatteringProcess, fTheScintillationProcess, fTheWLSProcess, and MF_LOG_DEBUG.

   {
     const detinfo::LArProperties* larp = lar::providerFrom<detinfo::LArPropertiesService>();
 
     // Add standard EM Processes
     MF_LOG_DEBUG("OpticalPhysics") << "PROCESSES BEING CONSTRUCTED IN OPTICAL PHYSICS";
 
     fTheCerenkovProcess = new G4Cerenkov("Cerenkov");
     fTheScintillationProcess = new G4Scintillation("Scintillation");
     fTheAbsorptionProcess = new G4OpAbsorption();
     fTheRayleighScatteringProcess = new G4OpRayleigh();
     if (fSimpleBoundary)
       fTheBoundaryProcess = new OpBoundaryProcessSimple();
     else
       fTheBoundaryProcess_g4 = new G4OpBoundaryProcess();
     fTheWLSProcess = new G4OpWLS();
 
     fTheCerenkovProcess->SetMaxNumPhotonsPerStep(700);
     fTheCerenkovProcess->SetMaxBetaChangePerStep(10.0);
     fTheCerenkovProcess->SetTrackSecondariesFirst(false);
 
     fTheScintillationProcess->SetScintillationYieldFactor(1.);
     fTheScintillationProcess->SetTrackSecondariesFirst(false);
 
     // Use Birks Correction in the Scintillation process
 
     G4EmSaturation* emSaturation = G4LossTableManager::Instance()->EmSaturation();
     fTheScintillationProcess->AddSaturation(emSaturation);
 
     bool CerenkovLightEnabled = larp->CerenkovLightEnabled();
 
     mf::LogInfo("OpticalPhysics") << "Cerenkov light enabled : " << CerenkovLightEnabled;
     static G4ParticleTable* fParticleTable = G4ParticleTable::GetParticleTable();
     G4ParticleTable::G4PTblDicIterator* aParticleIterator;
     aParticleIterator = fParticleTable->GetIterator();
 
     aParticleIterator->reset();
     while ((*aParticleIterator)()) {
       G4ParticleDefinition* particle = aParticleIterator->value();
       G4ProcessManager* pmanager = particle->GetProcessManager();
       G4String particleName = particle->GetParticleName();
       if (fTheCerenkovProcess->IsApplicable(*particle) && CerenkovLightEnabled) {
         pmanager->AddProcess(fTheCerenkovProcess);
         pmanager->SetProcessOrdering(fTheCerenkovProcess, idxPostStep);
       }
       if (fTheScintillationProcess->IsApplicable(*particle)) {
         pmanager->AddProcess(fTheScintillationProcess);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxAtRest);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxPostStep);
       }
 
       if (particleName == "opticalphoton") {
         mf::LogInfo("OpticalPhysics") << " AddDiscreteProcess to OpticalPhoton ";
         pmanager->AddDiscreteProcess(fTheAbsorptionProcess);
         pmanager->AddDiscreteProcess(fTheRayleighScatteringProcess);
         if (fSimpleBoundary)
           pmanager->AddDiscreteProcess(fTheBoundaryProcess);
         else
           pmanager->AddDiscreteProcess(fTheBoundaryProcess_g4);
         pmanager->AddDiscreteProcess(fTheWLSProcess);
       }
     }
   }