#include "OpticalPhysics.hh"

Inheritance diagram for larg4::OpticalPhysics:

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

virtual	~OpticalPhysics ()

virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

Private Attributes
G4int	verbose

G4Cerenkov *	fTheCerenkovProcess

G4Scintillation *	fTheScintillationProcess

G4OpAbsorption *	fTheAbsorptionProcess

G4OpRayleigh *	fTheRayleighScatteringProcess

OpBoundaryProcessSimple *	fTheBoundaryProcess

G4OpBoundaryProcess *	fTheBoundaryProcess_g4

G4OpWLS *	fTheWLSProcess

Detailed Description

Definition at line 81 of file OpticalPhysics.hh.

Constructor & Destructor Documentation

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

Definition at line 116 of file OpticalPhysics.cxx.

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

larg4::OpticalPhysics::~OpticalPhysics ( )

virtual

Definition at line 125 of file OpticalPhysics.cxx.

126 {}

Member Function Documentation

void larg4::OpticalPhysics::ConstructParticle ( )

virtual

Definition at line 129 of file OpticalPhysics.cxx.

References LOG_DEBUG.

   {
     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 ( )

virtual

Definition at line 163 of file OpticalPhysics.cxx.

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

   {
 
     const detinfo::LArProperties* larp = lar::providerFrom<detinfo::LArPropertiesService>();
     const detinfo::DetectorProperties* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
     // Add standard EM Processes
     LOG_DEBUG("OpticalPhysics") << "PROCESSES BEING CONSTRUCTED IN OPTICAL PHYSICS";
     
     fTheCerenkovProcess            = new G4Cerenkov("Cerenkov");
     fTheScintillationProcess       = new G4Scintillation("Scintillation");
     fTheAbsorptionProcess          = new G4OpAbsorption();
     fTheRayleighScatteringProcess  = new G4OpRayleigh();
     if(detp->SimpleBoundary())
       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);
         //      mf::LogInfo("OpticalPhysics")<<"OpticalPhysics : Cerenkov applicable : " << particleName;
       }
       if (fTheScintillationProcess->IsApplicable(*particle)) {
         pmanager->AddProcess(fTheScintillationProcess);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxAtRest);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxPostStep);
         //      mf::LogInfo("OpticalPhysics")<<"OpticalPhysics : Scintillation applicable : " << particleName;
       }
   
      if (particleName == "opticalphoton") {
        mf::LogInfo("OpticalPhysics") << " AddDiscreteProcess to OpticalPhoton ";
         pmanager->AddDiscreteProcess(fTheAbsorptionProcess);
         pmanager->AddDiscreteProcess(fTheRayleighScatteringProcess);
         if(detp->SimpleBoundary())
           pmanager->AddDiscreteProcess(fTheBoundaryProcess);
         else
           pmanager->AddDiscreteProcess(fTheBoundaryProcess_g4);
         pmanager->AddDiscreteProcess(fTheWLSProcess);
       }
     }
     
   }