LaserOpticalPhysicsList.cc

Go to the documentation of this file.

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

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

#include "G4Cerenkov.hh" 
#include "G4Scintillation.hh" 
#include "G4OpAbsorption.hh" 
#include "G4OpRayleigh.hh" 
#include "G4OpBoundaryProcess.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;
}