GammaPhysics.cxx

Go to the documentation of this file.

 
#include "SimCore/GammaPhysics.h"
 
namespace simcore {
 
GammaPhysics::GammaPhysics(const G4String& name,
                           const framework::config::Parameters& parameters)
    : G4VPhysicsConstructor(name),
      modelParameters{parameters.getParameter<framework::config::Parameters>(
          "photonuclear_model")} {}
 
void GammaPhysics::ConstructParticle() {}
 
void GammaPhysics::ConstructProcess() {
  G4ProcessManager* processManager = G4Gamma::Gamma()->GetProcessManager();
  if (processManager == nullptr) {
    EXCEPTION_RAISE("GammaPhysics",
                    "Was unable to access the process manager for photons, "
                    "something is very wrong!");
  }
  // configure our PN model based on runtime parameters
  auto pn = PhotoNuclearModel::Factory::get().make(
      modelParameters.getParameter<std::string>("class_name"),
      modelParameters.getParameter<std::string>("instance_name"),
      modelParameters);
  pn->removeExistingModel(processManager);
  pn->ConstructGammaProcess(processManager);
  const auto processes{processManager->GetProcessList()};
  for (int i{0}; i < processes->size(); i++) {
    const auto process{(*processes)[i]};
    if (process->GetProcessName() == "photonNuclear") {
      processManager->SetProcessOrderingToFirst(
          process, G4ProcessVectorDoItIndex::idxAll);
    }
  }
  // Add the gamma -> mumu to the physics list.
  processManager->AddDiscreteProcess(&gammaConvProcess);
}
}  // namespace simcore