#include "FastOpticalPhysics.h"

Inheritance diagram for larg4::FastOpticalPhysics:

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

virtual	~FastOpticalPhysics ()

virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

Private Attributes
G4int	verbose

G4Cerenkov *	fTheCerenkovProcess

OpFastScintillation *	fTheScintillationProcess

G4OpAbsorption *	fTheAbsorptionProcess

G4OpRayleigh *	fTheRayleighScatteringProcess

OpBoundaryProcessSimple *	fTheBoundaryProcess

G4OpWLS *	fTheWLSProcess

Detailed Description

Definition at line 82 of file FastOpticalPhysics.h.

Constructor & Destructor Documentation

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

Definition at line 121 of file FastOpticalPhysics.cxx.

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

larg4::FastOpticalPhysics::~FastOpticalPhysics ( )

virtual

Definition at line 130 of file FastOpticalPhysics.cxx.

131 {}

Member Function Documentation

void larg4::FastOpticalPhysics::ConstructParticle ( )

virtual

Definition at line 134 of file FastOpticalPhysics.cxx.

References LOG_DEBUG.

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

virtual

Definition at line 168 of file FastOpticalPhysics.cxx.

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

     {
     // Add standard EM Processes
     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       = new 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)) {
         pmanager->AddProcess(fTheScintillationProcess);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxAtRest);
         pmanager->SetProcessOrderingToLast(fTheScintillationProcess, idxPostStep);
 
       }
   
      if (particleName == "opticalphoton") {
        mf::LogInfo("FastOptical") << " AddDiscreteProcess to OpticalPhoton ";
        pmanager->AddDiscreteProcess(fTheAbsorptionProcess);
        pmanager->AddDiscreteProcess(fTheRayleighScatteringProcess);
        pmanager->AddDiscreteProcess(fTheBoundaryProcess);
        pmanager->AddDiscreteProcess(fTheWLSProcess);
      }
     }
     
   }