Geant4CerenkovPhysics.cpp

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
// Author     : M.Frank
//
//==========================================================================
 
#ifndef DDG4_GEANT4CERENKOVPHYSICS_H
#define DDG4_GEANT4CERENKOVPHYSICS_H 1
 
#include <DDG4/Geant4PhysicsList.h>
 
#include <G4ParticleTableIterator.hh>
#include <G4ParticleDefinition.hh>
#include <G4ParticleTypes.hh>
#include <G4ParticleTable.hh>
#include <G4ProcessManager.hh>
#include <G4Version.hh>
 
#if G4VERSION_NUMBER >= 1070
#include <G4OpticalParameters.hh>
#endif
 
#include <G4Cerenkov.hh>
 
namespace dd4hep {
 
  namespace sim {
 
 
    class Geant4CerenkovPhysics : public Geant4PhysicsList    {
    public:
      Geant4CerenkovPhysics() = delete;
      Geant4CerenkovPhysics(const Geant4CerenkovPhysics&) = delete;
      Geant4CerenkovPhysics(Geant4Context* ctxt, const std::string& nam)
        : Geant4PhysicsList(ctxt, nam)
      {
        declareProperty("MaxNumPhotonsPerStep",  m_maxNumPhotonsPerStep = 0);
        declareProperty("MaxBetaChangePerStep",  m_maxBetaChangePerStep = 0.0);
        declareProperty("TrackSecondariesFirst", m_trackSecondariesFirst = false);
        declareProperty("StackPhotons",          m_stackPhotons = true);
        declareProperty("VerboseLevel",          m_verbosity = 0);
      }
      virtual ~Geant4CerenkovPhysics() = default;
      virtual void constructProcesses(G4VUserPhysicsList* physics_list)   { 
        this->Geant4PhysicsList::constructProcesses(physics_list);
        info("+++ Constructing: maxNumPhotonsPerStep:%d maxBeta:%f "
             "track secondaries:%s stack photons:%s track secondaries:%s",
             m_maxNumPhotonsPerStep, m_maxBetaChangePerStep,
             yes_no(m_trackSecondariesFirst), yes_no(m_stackPhotons),
             yes_no(m_trackSecondariesFirst));
        G4Cerenkov* process = new G4Cerenkov(name());
#if G4VERSION_NUMBER >= 1070
        G4OpticalParameters* params = G4OpticalParameters::Instance();
        params->SetCerenkovVerboseLevel(m_verbosity);
        params->SetCerenkovMaxPhotonsPerStep(m_maxNumPhotonsPerStep);
        params->SetCerenkovMaxBetaChange(m_maxBetaChangePerStep);
        params->SetCerenkovTrackSecondariesFirst(m_trackSecondariesFirst);
        params->SetCerenkovStackPhotons(m_stackPhotons);
#else
        process->SetVerboseLevel(m_verbosity);
        process->SetMaxNumPhotonsPerStep(m_maxNumPhotonsPerStep);
        process->SetMaxBetaChangePerStep(m_maxBetaChangePerStep);
        process->SetTrackSecondariesFirst(m_trackSecondariesFirst);
#if G4VERSION_NUMBER > 1030
        process->SetStackPhotons(m_stackPhotons);
#endif
#endif
        auto pit = G4ParticleTable::GetParticleTable()->GetIterator();
        pit->reset();
        while( (*pit)() ){
          G4ParticleDefinition* particle = pit->value();
          if (process->IsApplicable(*particle)) {
            G4ProcessManager* pmanager = particle->GetProcessManager();
            pmanager->AddProcess(process);
            pmanager->SetProcessOrdering(process,idxPostStep);
          }
        }
      }
 
    private:
      double      m_maxBetaChangePerStep;
      int         m_maxNumPhotonsPerStep;
      int         m_verbosity;
      bool        m_trackSecondariesFirst;
      bool        m_stackPhotons;
    };
  }
}
#endif   // DDG4_GEANT4CERENKOVPHYSICS_H
 
#include <DDG4/Factories.h>
using namespace dd4hep::sim;
DECLARE_GEANT4ACTION(Geant4CerenkovPhysics)