A2PrimaryGeneratorMessenger: added command /A2/generator/SetBeamEnergy

Author: werthm

README.md

@@ -156,6 +156,7 @@ Command                                | Meaning
 `/A2/generator/SetTMax 450 MeV`        | maximum kinetic energy for a particle in the phase-space generator
 `/A2/generator/SetThetaMin 0 deg`      | minimum polar angle for a particle in the phase-space generator
 `/A2/generator/SetThetaMax 120 deg`    | maximum polar angle for a particle in the phase-space generator
+`/A2/generator/SetBeamEnergy 0.2 GeV`  | energy of incoming beam (overwrites beam energy in event file)
 `/A2/generator/SetBeamXSigma 10 mm`    | x-sigma of incoming beam
 `/A2/generator/SetBeamYSigma 10 mm`    | y-sigma of incoming beam
 `/A2/generator/SetBeamDiameter 10 mm`  | beam-spot diameter on target

include/A2PrimaryGeneratorMessenger.hh

@@ -38,6 +38,7 @@ class A2PrimaryGeneratorMessenger: public G4UImessenger
   G4UIcmdWithADoubleAndUnit* SetTmaxCmd;
   G4UIcmdWithADoubleAndUnit* SetThetaminCmd;
   G4UIcmdWithADoubleAndUnit* SetThetamaxCmd;
+  G4UIcmdWithADoubleAndUnit* SetBeamEnergyCmd;
   G4UIcmdWithADoubleAndUnit* SetBeamXSigmaCmd;
   G4UIcmdWithADoubleAndUnit* SetBeamYSigmaCmd;
   G4UIcmdWithADoubleAndUnit* SetBeamDiameterCmd;

src/A2PrimaryGeneratorMessenger.cc

@@ -71,6 +71,12 @@ A2PrimaryGeneratorMessenger::A2PrimaryGeneratorMessenger(
   SetThetamaxCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
   SetThetamaxCmd->SetUnitCategory("Angle");
 
+  SetBeamEnergyCmd = new G4UIcmdWithADoubleAndUnit("/A2/generator/SetBeamEnergy",this);
+  SetBeamEnergyCmd->SetGuidance("Set the energy of the photon beam");
+  SetBeamEnergyCmd->SetParameterName("BeamEnergy",false);
+  SetBeamEnergyCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
+  SetBeamEnergyCmd->SetUnitCategory("Energy");
+
   SetBeamXSigmaCmd = new G4UIcmdWithADoubleAndUnit("/A2/generator/SetBeamXSigma",this);
   SetBeamXSigmaCmd->SetGuidance("Set the X width of the photon beam");
   SetBeamXSigmaCmd->SetParameterName("BeamXSigma",false);
@@ -130,6 +136,7 @@ A2PrimaryGeneratorMessenger::~A2PrimaryGeneratorMessenger()
   delete SetThetamaxCmd;
   delete SetModeCmd;
   delete SetSeedCmd;
+  delete SetBeamEnergyCmd;
   delete SetBeamXSigmaCmd;
   delete SetBeamYSigmaCmd;
   delete SetBeamDiameterCmd;
@@ -174,6 +181,9 @@ void A2PrimaryGeneratorMessenger::SetNewValue(
   if( command == SetThetamaxCmd )
      { A2Action->SetThetamax(SetThetamaxCmd->GetNewDoubleValue(newValue));}
 
+   if( command == SetBeamEnergyCmd )
+     { A2Action->SetBeamEnergy(SetBeamEnergyCmd->GetNewDoubleValue(newValue));}
+
    if( command == SetBeamXSigmaCmd )
      { A2Action->SetBeamXSigma(SetBeamXSigmaCmd->GetNewDoubleValue(newValue));}