larg4::OpticalPhysics Class Reference

#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.

    : 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 94 of file OpticalPhysics.cxx.

  {
    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 129 of file OpticalPhysics.cxx.

  {
    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);
      }
    }
  }

Member Data Documentation

bool larg4::OpticalPhysics::fSimpleBoundary

private

Definition at line 88 of file OpticalPhysics.hh.

G4OpAbsorption* larg4::OpticalPhysics::fTheAbsorptionProcess

private

Definition at line 91 of file OpticalPhysics.hh.

OpBoundaryProcessSimple* larg4::OpticalPhysics::fTheBoundaryProcess

private

Definition at line 93 of file OpticalPhysics.hh.

G4OpBoundaryProcess* larg4::OpticalPhysics::fTheBoundaryProcess_g4

private

Definition at line 94 of file OpticalPhysics.hh.

G4Cerenkov* larg4::OpticalPhysics::fTheCerenkovProcess

private

Definition at line 89 of file OpticalPhysics.hh.

G4OpRayleigh* larg4::OpticalPhysics::fTheRayleighScatteringProcess

private

Definition at line 92 of file OpticalPhysics.hh.

G4Scintillation* larg4::OpticalPhysics::fTheScintillationProcess

private

Definition at line 90 of file OpticalPhysics.hh.

G4OpWLS* larg4::OpticalPhysics::fTheWLSProcess

private

Definition at line 95 of file OpticalPhysics.hh.

---

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

• larsim/larsim/LegacyLArG4/OpticalPhysics.hh
• larsim/larsim/LegacyLArG4/OpticalPhysics.cxx