File indexing completed on 2025-02-23 09:22:11
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028 #include "PrimaryGeneratorMessenger.hh"
0029
0030 #include "PrimaryGeneratorAction.hh"
0031
0032 #include "G4ParticleTable.hh"
0033 #include "G4SingleParticleSource.hh"
0034 #include "G4UIcmdWithADoubleAndUnit.hh"
0035 #include "G4UIcmdWithAString.hh"
0036 #include "G4UIdirectory.hh"
0037
0038
0039 PrimaryGeneratorMessenger::PrimaryGeneratorMessenger(PrimaryGeneratorAction* pPrimaryGenerator)
0040 : G4UImessenger(), fpPrimaryGenerator(pPrimaryGenerator)
0041 {
0042 fparticle = std::make_unique<G4UIcmdWithAString>("/UHDR/source/particle", this);
0043 fparticle->SetGuidance("Add time structure.");
0044 fparticle->SetGuidance("e-, proton, alpha");
0045
0046 fenergy = std::make_unique<G4UIcmdWithADoubleAndUnit>("/UHDR/source/energy", this);
0047 fenergy->SetGuidance("Sets a monocromatic energy (same as gps/energy)");
0048 fenergy->SetParameterName("monoenergy", false, false);
0049 fenergy->SetDefaultUnit("keV");
0050 }
0051
0052
0053
0054 PrimaryGeneratorMessenger::~PrimaryGeneratorMessenger() = default;
0055
0056
0057
0058 void PrimaryGeneratorMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
0059 {
0060 auto particleGun = fpPrimaryGenerator->GetSPGun();
0061 if (command == fparticle.get()) {
0062 G4ExceptionDescription ed;
0063 G4ParticleDefinition* pd = G4ParticleTable::GetParticleTable()->FindParticle(newValue);
0064 if (pd != nullptr) {
0065 particleGun->SetParticleDefinition(pd);
0066 }
0067 else {
0068 ed << "Particle [" << newValue << "] is not found.";
0069 command->CommandFailed(ed);
0070 }
0071 }
0072 else if (command == fenergy.get()) {
0073 auto pEnergyDis = particleGun->GetEneDist();
0074 pEnergyDis->SetMonoEnergy(fenergy->GetNewDoubleValue(newValue));
0075 }
0076 }
0077
0078