extra gamma particle physics for simulation and sets up the photonuclear model to use from the configuration More...

Public Member Functions
	GammaPhysics (const G4String &name, const framework::config::Parameters &parameters)
	Class constructor.

virtual	~GammaPhysics ()=default
	Class destructor.

void	ConstructParticle () final
	Construct particles.

void	ConstructProcess () final
	We do two things for this call back during initialization.

Private Attributes
G4GammaConversionToMuons	gammaConvProcess
	The gamma to muons process.

framework::config::Parameters	modelParameters
	Parameters from the configuration to pass along to the photonuclear model.

Detailed Description

extra gamma particle physics for simulation and sets up the photonuclear model to use from the configuration

Note

Is responsible for selecting the photonuclear model from the python configuration. Currently adds gamma -> mumu reaction using the G4GammaConversionToMuons process. Also changes ordering of gamma processes such that photonNuclear and GammaToMuMu are called first.

Definition at line 36 of file GammaPhysics.h.

Constructor & Destructor Documentation

◆ GammaPhysics()

simcore::GammaPhysics::GammaPhysics	(	const G4String &	name,
		const framework::config::Parameters &	parameters
	)

Class constructor.

Parameters

name	The name of the physics.
parameters	The python configuration

Definition at line 12 of file GammaPhysics.cxx.

    : G4VPhysicsConstructor(name),
      modelParameters{parameters.getParameter<framework::config::Parameters>(
          "photonuclear_model")} {}

Member Function Documentation

◆ ConstructParticle()

void simcore::GammaPhysics::ConstructParticle ( )

final

Construct particles.

We don't do anything here since we are just attaching/updating the photon physics.

Definition at line 18 of file GammaPhysics.cxx.

18{}

◆ ConstructProcess()

void simcore::GammaPhysics::ConstructProcess ( )

final

We do two things for this call back during initialization.

We add the muon-conversion process to the photon's process table, enabling it to be one of the processes that can happen.
We move the photonNuclear process to be the first process in the photon's ordering so that the bias operator has a chance to lower its cross-section before any EM process is called (if need be).

Put the PN process first in the ordering in case PN biasing is happening.

Process ordering is a complicated concept and unfortunately its affect on biasing is poorly documented. What has been said is that some processes need to be called last[1]. In addition, the practical experience of working on defining a custom process for G4DarkBreM showed that sometimes Geant4 does not get through the full list of processes but it always starts at the beginning. For these reasons, we put the PN process first in the ordering so that we can insure it is always check by Geant4 before continuing.

[1] https://indico.cern.ch/event/58317/contributions/2047449/attachments/992300/1411062/PartDef-ProcMan-AS.pdf

Definition at line 20 of file GammaPhysics.cxx.

                                    {
  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);
}

References gammaConvProcess, simcore::Factory< Prototype, PrototypePtr, PrototypeConstructorArgs >::get(), framework::config::Parameters::getParameter(), simcore::Factory< Prototype, PrototypePtr, PrototypeConstructorArgs >::make(), and modelParameters.

Member Data Documentation

◆ gammaConvProcess

G4GammaConversionToMuons simcore::GammaPhysics::gammaConvProcess

private

The gamma to muons process.

Definition at line 75 of file GammaPhysics.h.

Referenced by ConstructProcess().

◆ modelParameters

framework::config::Parameters simcore::GammaPhysics::modelParameters

private

Parameters from the configuration to pass along to the photonuclear model.

Definition at line 80 of file GammaPhysics.h.