GammaPhysics.cxx

Go to the documentation of this file.

 
#include "SimCore/GammaPhysics.h"
 
#include "G4ProcessManager.hh"
#include "G4VProcess.hh"
#include "SimCore/PhotoNuclearModels/PhotoNuclearModel.h"
 
namespace simcore {
 
GammaPhysics::GammaPhysics(const G4String& name,
                           const framework::config::Parameters& parameters)
    : G4VPhysicsConstructor(name),
      model_parameters_{
          parameters.get<framework::config::Parameters>("photonuclear_model")} {
}
 
void GammaPhysics::ConstructParticle() {}
 
void GammaPhysics::ConstructProcess() {
  G4ProcessManager* process_manager = G4Gamma::Gamma()->GetProcessManager();
  if (process_manager == 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(
      model_parameters_.get<std::string>("class_name"),
      model_parameters_.get<std::string>("instance_name"), model_parameters_);
  if (not pn) {
    EXCEPTION_RAISE("UnableToCreate",
                    "Unable to create a PhotoNuclearModel of type " +
                        model_parameters_.get<std::string>("class_name"));
  }
  pn.value()->removeExistingModel(process_manager);
  pn.value()->constructGammaProcess(process_manager);
  const auto processes{process_manager->GetProcessList()};
  for (int i{0}; i < processes->size(); i++) {
    const auto process{(*processes)[i]};
    if (process->GetProcessName() == "photonNuclear") {
      process_manager->SetProcessOrderingToFirst(
          process, G4ProcessVectorDoItIndex::idxAll);
    }
  }
  // Add the gamma -> mumu to the physics list.
  process_manager->AddDiscreteProcess(&gamma_conv_process_);
}
}  // namespace simcore