#include "FastOpticalPhysics.h"

Inheritance diagram for larg4::FastOpticalPhysics:

Public Member Functions
	FastOpticalPhysics (G4int ver=0, const G4String &name="FastOptical")

void	ConstructParticle () override

void	ConstructProcess () override

Private Attributes
G4Cerenkov *	fTheCerenkovProcess

std::unique_ptr< OpFastScintillation >	fTheScintillationProcess

G4OpAbsorption *	fTheAbsorptionProcess

G4OpRayleigh *	fTheRayleighScatteringProcess

OpBoundaryProcessSimple *	fTheBoundaryProcess

G4OpWLS *	fTheWLSProcess

Detailed Description

Definition at line 79 of file FastOpticalPhysics.h.

Constructor & Destructor Documentation

larg4::FastOpticalPhysics::FastOpticalPhysics	(	G4int	ver = `0`,
		const G4String &	name = `"FastOptical"`
	)

Definition at line 83 of file FastOpticalPhysics.cxx.

     : G4VPhysicsConstructor(name)
   {
     G4LossTableManager::Instance();
     mf::LogInfo("FastOpticalPhysics") << "OBJECT BEING CONSTRUCTED IN OPTICAL PHYSICS";
   }

Member Function Documentation

void larg4::FastOpticalPhysics::ConstructParticle ( )

override

Definition at line 91 of file FastOpticalPhysics.cxx.

References MF_LOG_DEBUG.

   {
     MF_LOG_DEBUG("FastOpticalPhysics") << "PARTICLES BEING CONSTRUCTED IN FAST 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::FastOpticalPhysics::ConstructProcess ( )

override

Definition at line 125 of file FastOpticalPhysics.cxx.

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

   {
     // Add standard EM Processes
     MF_LOG_DEBUG("FastOpticalPhysics") << "PROCESSES BEING CONSTRUCTED IN OPTICAL PHYSICS";
 
     fTheCerenkovProcess = new G4Cerenkov("Cerenkov");
     fTheAbsorptionProcess = new G4OpAbsorption();
     fTheRayleighScatteringProcess = new G4OpRayleigh();
     fTheBoundaryProcess = new OpBoundaryProcessSimple();
     fTheWLSProcess = new G4OpWLS();
     fTheScintillationProcess = std::make_unique<OpFastScintillation>("FastScintillation");
 
     fTheCerenkovProcess->SetMaxNumPhotonsPerStep(700);
     fTheCerenkovProcess->SetMaxBetaChangePerStep(10.0);
     fTheCerenkovProcess->SetTrackSecondariesFirst(false);
 
     const detinfo::LArProperties* larp = lar::providerFrom<detinfo::LArPropertiesService>();
     bool CerenkovEnabled = larp->CerenkovLightEnabled();
 
     mf::LogInfo("FastOpticalPhysics") << "Cerenkov enabled : " << CerenkovEnabled;
     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) && CerenkovEnabled) {
         pmanager->AddProcess(fTheCerenkovProcess);
         pmanager->SetProcessOrdering(fTheCerenkovProcess, idxPostStep);
       }
       if (fTheScintillationProcess->IsApplicable(*particle)) {
         auto ptr = fTheScintillationProcess.get();
         pmanager->AddProcess(ptr);
         pmanager->SetProcessOrderingToLast(ptr, idxAtRest);
         pmanager->SetProcessOrderingToLast(ptr, idxPostStep);
       }
 
       if (particleName == "opticalphoton") {
         mf::LogInfo("FastOptical") << " AddDiscreteProcess to OpticalPhoton ";
         pmanager->AddDiscreteProcess(fTheAbsorptionProcess);
         pmanager->AddDiscreteProcess(fTheRayleighScatteringProcess);
         pmanager->AddDiscreteProcess(fTheBoundaryProcess);
         pmanager->AddDiscreteProcess(fTheWLSProcess);
       }
     }
   }