LaserOpticalPhysicsList.cc

Go to the documentation of this file.

#include "Alignment/LaserAlignmentSimulation/interface/LaserOpticalPhysicsList.h"
#include "SimG4Core/PhysicsLists/interface/EmParticleList.h"

#include "G4ParticleTable.hh"
#include "G4ProcessManager.hh"

#include "G4Cerenkov.hh"
#include "G4OpAbsorption.hh"
#include "G4OpBoundaryProcess.hh"
#include "G4OpRayleigh.hh"
#include "G4Scintillation.hh"

LaserOpticalPhysicsList::LaserOpticalPhysicsList(const G4String &name)
    : G4VPhysicsConstructor(name),
      wasActivated(false),
      theScintProcess(),
      theCerenkovProcess(),
      theAbsorptionProcess(),
      theRayleighScattering(),
      theBoundaryProcess(),
      theWLSProcess() {
  if (verboseLevel > 0)
    std::cout << "<LaserOpticalPhysicsList::LaserOpticalPhysicsList(...)> "
                 "entering constructor ..."
              << std::endl;
}

LaserOpticalPhysicsList::~LaserOpticalPhysicsList() {
  if (verboseLevel > 0) {
    std::cout << "<LaserOpticalPhysicsList::~LaserOpticalPhysicsList()> "
                 "entering destructor ... "
              << std::endl;
    std::cout << "  deleting the processes ... ";
  }
  if (theWLSProcess != nullptr) {
    delete theWLSProcess;
  }
  if (theBoundaryProcess != nullptr) {
    delete theBoundaryProcess;
  }
  if (theRayleighScattering != nullptr) {
    delete theRayleighScattering;
  }
  if (theAbsorptionProcess != nullptr) {
    delete theAbsorptionProcess;
  }
  if (theScintProcess != nullptr) {
    delete theScintProcess;
  }
  if (verboseLevel > 0)
    std::cout << " done " << std::endl;
}

void LaserOpticalPhysicsList::ConstructParticle() {
  if (verboseLevel > 0)
    std::cout << "<LaserOpticalPhysicsList::ConstructParticle()>: constructing "
                 "the optical photon ... "
              << std::endl;

  // optical photon
  G4OpticalPhoton::OpticalPhotonDefinition();
}

void LaserOpticalPhysicsList::ConstructProcess() {
  if (verboseLevel > 0)
    std::cout << "<LaserOpticalPhysicsList::ConstructProcess()>: constructing "
                 "the physics ... "
              << std::endl;

  theScintProcess = new G4Scintillation();
  //  theCerenkovProcess=new G4Cerenkov();
  theAbsorptionProcess = new G4OpAbsorption();
  theRayleighScattering = new G4OpRayleigh();
  theBoundaryProcess = new G4OpBoundaryProcess("OpBoundary");
  theWLSProcess = new G4OpWLS();

  // set the verbosity level
  theAbsorptionProcess->SetVerboseLevel(verboseLevel);
  theBoundaryProcess->SetVerboseLevel(verboseLevel);

  G4ProcessManager *pManager = nullptr;

  pManager = G4OpticalPhoton::OpticalPhoton()->GetProcessManager();
  pManager->AddDiscreteProcess(theAbsorptionProcess);
  pManager->AddDiscreteProcess(theRayleighScattering);
  // theBoundaryProcess->SetModel(unified);
  pManager->AddDiscreteProcess(theBoundaryProcess);
  pManager->AddDiscreteProcess(theWLSProcess);

  theScintProcess->SetScintillationYieldFactor(1.);
  theScintProcess->SetScintillationExcitationRatio(0.0);
  theScintProcess->SetTrackSecondariesFirst(true);

  G4ParticleTable *table = G4ParticleTable::GetParticleTable();
  EmParticleList emList;
  for (const auto &particleName : emList.PartNames()) {
    G4ParticleDefinition *particle = table->FindParticle(particleName);
    pManager = particle->GetProcessManager();
    if (theScintProcess->IsApplicable(*particle)) {
      pManager->AddProcess(theScintProcess);
      pManager->SetProcessOrderingToLast(theScintProcess, idxAtRest);
      pManager->SetProcessOrderingToLast(theScintProcess, idxPostStep);
    }
  }

  wasActivated = true;
}