19#include "G4Version.hh"
22#include "BDSException.hh"
23#include "BDSGlobalConstants.hh"
24#include "BDSModularPhysicsList.hh"
25#include "BDSIonDefinition.hh"
26#include "BDSParticleDefinition.hh"
27#if G4VERSION_NUMBER > 1039
28#include "BDSPhysicsChannelling.hh"
30#include "BDSPhysicsCutsAndLimits.hh"
31#include "BDSPhysicsEMDissociation.hh"
32#include "BDSPhysicsMuonSplitting.hh"
33#include "BDSPhysicsUtilities.hh"
34#include "BDSUtilities.hh"
35#include "BDSWarning.hh"
36#include "BDSEmStandardPhysicsOp4Channelling.hh"
38#include "FTFP_BERT.hh"
40#include "G4AntiNeutrinoE.hh"
41#include "G4AntiNeutrinoMu.hh"
42#include "G4AntiNeutrinoTau.hh"
43#include "G4AntiNeutron.hh"
44#include "G4AntiProton.hh"
45#include "G4Electron.hh"
46#include "G4EmParameters.hh"
47#include "G4EmStandardPhysics_option4.hh"
48#include "G4EmStandardPhysicsSS.hh"
49#include "G4DynamicParticle.hh"
51#include "G4GenericBiasingPhysics.hh"
52#include "G4GenericIon.hh"
53#include "G4IonElasticPhysics.hh"
54#include "G4IonTable.hh"
55#include "G4KaonMinus.hh"
56#include "G4KaonPlus.hh"
57#include "G4KaonZero.hh"
58#include "G4KaonZeroLong.hh"
59#include "G4KaonZeroShort.hh"
60#include "G4LeptonConstructor.hh"
61#include "G4MuonMinus.hh"
62#include "G4MuonPlus.hh"
63#include "G4NeutrinoE.hh"
64#include "G4NeutrinoMu.hh"
65#include "G4NeutrinoTau.hh"
66#include "G4Neutron.hh"
67#include "G4ParticleTable.hh"
68#include "G4ParticleTableIterator.hh"
69#include "G4PionMinus.hh"
70#include "G4PionPlus.hh"
71#include "G4PionZero.hh"
72#include "G4Positron.hh"
73#include "G4ProductionCutsTable.hh"
75#include "G4PhysListFactory.hh"
76#include "G4ProcessManager.hh"
77#include "G4ProcessVector.hh"
80#include "G4UImanager.hh"
81#if G4VERSION_NUMBER > 1049
82#include "G4ParticleDefinition.hh"
83#include "G4CoupledTransportation.hh"
84#include "G4Transportation.hh"
88#if G4VERSION_NUMBER > 1069
89#include "G4HadronicParameters.hh"
92#include "parser/beam.h"
93#include "parser/fastlist.h"
94#include "parser/physicsbiasing.h"
105 return G4IonTable::IsIon(particle) && particle!=G4Proton::Definition();
110 return BDS::IsIon(particle->GetDefinition()) || particle->GetTotalOccupancy()>0;
119 G4VModularPhysicsList* result =
nullptr;
126 G4bool completePhysics =
BDS::StrContains(physicsListNameLower,
"complete");
127 if (useGeant4Physics)
130 G4String geant4PhysicsList = physicsList.substr(2);
131 G4PhysListFactory factory;
132 if (!factory.IsReferencePhysList(geant4PhysicsList))
134 G4cerr <<
"Unknown Geant4 physics list \"" << geant4PhysicsList <<
"\"" << G4endl;
135 G4cout <<
"Available Geant4 hadronic physics lists:" << G4endl;
136 for (
const auto &name : factory.AvailablePhysLists())
137 {G4cout <<
"\"" << name <<
"\"" << G4endl;}
138 G4cout <<
"Available Geant4 EM physics lists:" << G4endl;
139 for (
const auto &name : factory.AvailablePhysListsEM())
140 {G4cout <<
"\"" << name <<
"\"" << G4endl;}
141 throw BDSException(__METHOD_NAME__,
"Unknown Geant4 physics list \"" + geant4PhysicsList +
"\"");
145 result = factory.GetReferencePhysList(geant4PhysicsList);
146 if (g->G4PhysicsUseBDSIMRangeCuts())
148 if (g->MinimumKineticEnergy() > 0 || g->G4PhysicsUseBDSIMCutsAndLimits())
150 G4cout <<
"\nAdding cuts and limits physics process to Geant4 reference physics list" << G4endl;
151 G4cout <<
"This is to enforce BDSIM range cuts and the minimumKinetic energy option.\n";
152 G4cout <<
"This is done by default for the functionality of BDSIM tracking and should not affect the physics greatly.\n";
153 G4cout <<
"See the BDSIM manual about Geant4 reference physics lists for details." << G4endl;
156 else if (!g->G4PhysicsUseBDSIMCutsAndLimits() && g->Circular())
158 G4String message =
"g4PhysicsUseBDSIMCutsAndLimits turned off but using a circular machine - circular mechanics will be broken";
159 BDS::Warning(__METHOD_NAME__, message);
163 else if (completePhysics)
167 G4cout <<
"Constructing \"" << physicsListNameLower <<
"\" complete physics list" << G4endl;
168#if G4VERSION_NUMBER > 1039
176 r->SetVerboseLevel(verbosity);
179 throw BDSException(__METHOD_NAME__,
"Channel physics is not supported with Geant4 versions less than 10.4");
183 {
throw BDSException(__METHOD_NAME__,
"Unknown 'complete' physics list \"" + physicsList +
"\"");}
193 result->ConstructParticle();
194 result->SetVerboseLevel(verbosity);
197 G4UImanager* UIManager = G4UImanager::GetUIpointer();
199 if (!physicsMacro.empty())
202 G4String physicsMacroFull =
BDS::GetFullPath(physicsMacro,
false, setInExecOptions);
203 G4cout <<
"Applying geant4 physics macro file: " << physicsMacroFull << G4endl;
204 UIManager->ApplyCommand(
"/control/execute " + physicsMacroFull);
207 G4VUserPhysicsList* resultAsUserPhysicsList =
dynamic_cast<G4VUserPhysicsList*
>(result);
208 if (resultAsUserPhysicsList)
210 resultAsUserPhysicsList->DumpCutValuesTable(verbosity);
211 resultAsUserPhysicsList->SetVerboseLevel(verbosity);
217 const std::set<std::string>& keys)
220 for (
const auto& k : keys)
221 {nSet += beamDefinition.
HasBeenSet(k) ? 1 : 0;}
226 const std::set<std::string>& keys,
228 G4bool warnZeroParamsSet,
229 const G4String& unitString)
233 G4cerr <<
"Beam> More than one parameter set - there should only be one" << G4endl;
234 for (
const auto& k : keys)
235 {G4cerr << std::setw(14) << std::left << k <<
": " << std::setw(7) << std::right << beamDefinition.
get_value(k) <<
" " << unitString << G4endl;}
236 throw BDSException(__METHOD_NAME__,
"conflicting parameters set");
238 else if (nSet == 0 && warnZeroParamsSet)
240 G4cerr <<
"Beam> One of the following required to be set" << G4endl;
241 for (
const auto &k : keys)
242 {G4cerr << std::setw(14) << std::left << k <<
": " << std::setw(7) << std::right << beamDefinition.
get_value(k) <<
" " << unitString << G4endl;}
243 throw BDSException(__METHOD_NAME__,
"insufficient parameters set");
251 G4bool& beamDifferentFromDesignParticle)
253 if (beamDefinition.
particle.empty())
254 {
throw BDSException(
"Beam> no \"particle\" specified (required).");}
257 std::set<std::string> keysDesign = {
"energy",
"momentum",
"kineticEnergy"};
260 std::set<std::string> keysParticle = {
"E0",
"P0",
"Ek0"};
271 beamDifferentFromDesignParticle = nSetParticle > 0 || beamParticleName != beamDefinition.
particle;
272 if (nSetParticle > 0)
276 beamDefinition.
E0 * CLHEP::GeV,
277 beamDefinition.
Ek0 * CLHEP::GeV,
278 beamDefinition.
P0 * CLHEP::GeV,
296 G4double totalEnergyIn,
297 G4double kineticEnergyIn,
304 std::map<G4String, G4String> commonSubstitutions = { {
"photon",
"gamma"},
310 {
"antiproton",
"anti_proton"},
311 {
"anti-proton",
"anti_proton"} };
313 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
314 std::regex ionParticle(
"(ion\\s).*");
315 if (std::regex_match(particleName, ionParticle))
317 G4GenericIon::GenericIonDefinition();
320 G4IonTable* ionTable = particleTable->GetIonTable();
322 G4int ionPDGID = G4IonTable::GetNucleusEncoding(ionDef->Z(), ionDef->A());
323 G4double mass = ionTable->GetIonMass(ionDef->Z(), ionDef->A());
324 mass += ionDef->NElectrons()*G4Electron::Definition()->GetPDGMass();
325 G4double charge = ionDef->Charge();
327 totalEnergyIn, kineticEnergyIn, momentumIn, ffact, ionDef, ionPDGID);
333 auto searchName = commonSubstitutions.find(particleName);
334 if (searchName != commonSubstitutions.end())
336 G4cout <<
"Substituting particle name \"" << particleName <<
"\" for the Geant4 name: \"" << searchName->second <<
"\"" << G4endl;
337 particleName = searchName->second;
341 G4ParticleDefinition* particleDef =
nullptr;
348 int particleID = std::stoi(particleName);
352 G4ParticleTable::G4PTblEncodingDictionary* encoding = G4ParticleTable::fEncodingDictionary;
353 auto search = encoding->find(particleID);
354 if (search != encoding->end())
355 {particleDef = search->second;}
357 {
throw BDSException(__METHOD_NAME__,
"PDG ID \"" + particleName +
"not found in particle table");}
359 catch (
const std::logic_error&)
361 particleDef = particleTable->FindParticle(particleName);
363 {particleDef = particleTable->FindParticle(particleNameIn);}
368 throw BDSException(__METHOD_NAME__,
"Particle \"" + particleName +
"\" not found.");
370 particleDefB =
new BDSParticleDefinition(particleDef, totalEnergyIn, kineticEnergyIn, momentumIn, ffact);
379 if (name ==
"proton")
380 {G4Proton::ProtonDefinition();}
381 else if (name ==
"anti_proton")
382 {G4AntiProton::AntiProtonDefinition();}
383 else if (name ==
"e-")
384 {G4Electron::ElectronDefinition();}
385 else if (name ==
"e+")
386 {G4Positron::PositronDefinition();}
387 else if (name ==
"pi-")
388 {G4PionMinus::PionMinusDefinition();}
389 else if (name ==
"pi+")
390 {G4PionPlus::PionPlusDefinition();}
391 else if (name ==
"pi0")
392 {G4PionZero::PionZeroDefinition();}
393 else if (name ==
"neutron")
394 {G4Neutron::NeutronDefinition();}
395 else if (name ==
"photon" || name ==
"gamma")
397 else if (name ==
"mu-")
398 {G4MuonMinus::MuonMinusDefinition();}
399 else if (name ==
"mu+")
400 {G4MuonPlus::MuonPlusDefinition();}
401 else if (name ==
"kaon-")
402 {G4KaonMinus::KaonMinusDefinition();}
403 else if (name ==
"kaon+")
404 {G4KaonPlus::KaonPlusDefinition();}
405 else if (name ==
"kaon0")
406 {G4KaonZero::KaonZeroDefinition();}
407 else if (name ==
"kaon0l")
408 {G4KaonZeroLong::KaonZeroLongDefinition();}
409 else if (name ==
"kaon0s")
410 {G4KaonZeroShort::KaonZeroShortDefinition();}
411 else if (name ==
"nu_e")
412 {G4NeutrinoE::NeutrinoEDefinition();}
413 else if (name ==
"anti_nu_e")
414 {G4AntiNeutrinoE::AntiNeutrinoEDefinition();}
415 else if (name ==
"nu_mu")
416 {G4NeutrinoMu::NeutrinoMuDefinition();}
417 else if (name ==
"anti_nu_mu")
418 {G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();}
419 else if (name ==
"nu_tau")
420 {G4NeutrinoTau::NeutrinoTauDefinition();}
421 else if (name ==
"anti_nu_tau")
422 {G4AntiNeutrinoTau::AntiNeutrinoTauDefinition();}
425 G4String msg =
"Unknown common beam particle type \"" + name;
426 msg +=
"\" - if it doesn't work, include all \"all_particles\" in the physicsList option.";
427 BDS::Warning(__METHOD_NAME__, msg);
434 G4Electron::ElectronDefinition();
435 G4Positron::PositronDefinition();
436 G4NeutrinoE::NeutrinoEDefinition();
437 G4AntiNeutrinoE::AntiNeutrinoEDefinition();
440 G4Proton::ProtonDefinition();
441 G4AntiProton::AntiProtonDefinition();
442 G4Neutron::NeutronDefinition();
443 G4AntiNeutron::AntiNeutronDefinition();
451 G4LeptonConstructor::ConstructParticle();
452 G4PionPlus::PionPlusDefinition();
453 G4PionMinus::PionMinusDefinition();
454 G4KaonPlus::KaonPlusDefinition();
455 G4KaonMinus::KaonMinusDefinition();
456 G4GenericIon::GenericIonDefinition();
461 G4cout <<
"Register physics processes by name for the primary particle \""
462 << primaryParticleName <<
"\":" << G4endl;
464 auto particle = G4ParticleTable::GetParticleTable()->FindParticle(primaryParticleName);
468 G4cout << __METHOD_NAME__ <<
"primary particle not defined yet - could be ion" << G4endl;
472 auto pl = particle->GetProcessManager()->GetProcessList();
473 for (G4int i = 0; i < (G4int)pl->length(); i++)
474 {G4cout <<
"\"" << (*pl)[i]->GetProcessName() <<
"\"" << G4endl;}
480 std::set<const G4ParticleDefinition*> particlesToBias;
481 for (
const auto& b : biases)
483 const G4ParticleDefinition* particle =
nullptr;
484 G4String particleName = G4String(b.particle);
486 particle = G4ParticleTable::GetParticleTable()->FindParticle(particleName);
489 {particlesToBias.insert(particle);}
491 {
throw BDSException(__METHOD_NAME__,
"Unknown particle type for biasing: \"" + particleName +
"\"");}
494 if (particlesToBias.empty())
497 G4GenericBiasingPhysics* physBias =
new G4GenericBiasingPhysics();
498 for (
auto part : particlesToBias)
500 const G4String& particleName = part->GetParticleName();
501 G4cout << __METHOD_NAME__ <<
"wrapping \"" << particleName <<
"\" for biasing" << G4endl;
502 physBias->Bias(particleName);
510 G4int muonSplittingFactor = globals->MuonSplittingFactor();
511 if (muonSplittingFactor > 1)
513 G4int muonSplittingFactor2 = globals->MuonSplittingFactor2();
514 G4double muonSplittingThresholdParentEk = globals->MuonSplittingThresholdParentEk();
515 G4double muonSplittingThresholdParentEk2 = globals->MuonSplittingThresholdParentEk2();
516 G4cout <<
"BDSPhysicsMuonSplitting -> using muon splitting wrapper -> factor of: " << muonSplittingFactor << G4endl;
517 if (muonSplittingThresholdParentEk > 0)
518 {G4cout <<
"BDSPhysicsMuonSplitting -> minimum parent kinetic energy: " << muonSplittingThresholdParentEk / CLHEP::GeV <<
" GeV" << G4endl;}
519 if (muonSplittingFactor2 > 1)
521 G4cout <<
"BDSPhysicsMuonSplitting -> factor #2: " << muonSplittingFactor2 <<
" for muons above "
522 << muonSplittingThresholdParentEk / CLHEP::GeV <<
" GeV" << G4endl;
524 G4bool excludeW1P = globals->MuonSplittingExcludeWeight1Particles();
526 muonSplittingFactor2, muonSplittingThresholdParentEk2,
527 excludeW1P, globals->MuonSplittingExclusionWeight()));
533 G4cout << __METHOD_NAME__ <<
"Defined particles: " << G4endl;
534 auto it = G4ParticleTable::GetParticleTable()->GetIterator();
537 {G4cout << it->value()->GetParticleName() <<
" ";}
541#if G4VERSION_NUMBER > 1039
547 G4VModularPhysicsList* physlist =
new FTFP_BERT();
548 physlist->RegisterPhysics(
new G4IonElasticPhysics());
549 G4GenericBiasingPhysics* biasingPhysics =
new G4GenericBiasingPhysics();
558 physlist->ReplacePhysics(
new G4EmStandardPhysics_option4());
562 physlist->ReplacePhysics(
new G4EmStandardPhysicsSS());
568 G4cout <<
"Adding EM Dissocation to crystal channelling physics list" << G4endl;
572 biasingPhysics->PhysicsBiasAllCharged();
573 physlist->RegisterPhysics(biasingPhysics);
577 G4cout <<
"\nWARNING - adding cuts and limits physics process to \"COMPLETE\" physics list" << G4endl;
578 G4cout <<
"This is to enforce BDSIM range cuts and the minimumKinetic energy option.\n";
579 G4cout <<
"This can be turned off by setting option, g4PhysicsUseBDSIMCutsAndLimits=0;\n" << G4endl;
591 if (globals->DefaultRangeCutsSet())
593 G4cout << __METHOD_NAME__ <<
"Default production range cut " << physicsList->GetDefaultCutValue() <<
" mm" << G4endl;
594 physicsList->SetDefaultCutValue(globals->DefaultRangeCut());
596 if (globals->ProdCutPhotonsSet())
598 G4cout << __METHOD_NAME__ <<
"Photon production range cut " << physicsList->GetCutValue(
"gamma") <<
" mm" << G4endl;
599 physicsList->SetCutValue(globals->ProdCutPhotons(),
"gamma");
601 if (globals->ProdCutElectronsSet())
603 G4cout << __METHOD_NAME__ <<
"Electron production range cut " << physicsList->GetCutValue(
"e-") <<
" mm" << G4endl;
604 physicsList->SetCutValue(globals->ProdCutElectrons(),
"e-");
606 if (globals->ProdCutPositronsSet())
608 G4cout << __METHOD_NAME__ <<
"Positron production range cut " << physicsList->GetCutValue(
"e+") <<
" mm" << G4endl;
609 physicsList->SetCutValue(globals->ProdCutPositrons(),
"e+");
611 if (globals->ProdCutProtonsSet())
613 G4cout << __METHOD_NAME__ <<
"Proton production range cut " << physicsList->GetCutValue(
"proton") <<
" mm" << G4endl;
614 physicsList->SetCutValue(globals->ProdCutProtons(),
"proton");
618 G4cout << __METHOD_NAME__ <<
"Range cuts from inspection of the physics list" << G4endl;
619 G4cout << __METHOD_NAME__ <<
"Default production range cut " << physicsList->GetDefaultCutValue() <<
" mm" << G4endl;
622 G4cout << __METHOD_NAME__ <<
"List of all constructed particles by physics lists" << G4endl;
623 for (
auto particle : *G4ParticleTable::fDictionary)
624 {G4cout << particle.second->GetParticleName() <<
", ";}
628 physicsList->DumpCutValuesTable(verbosity);
634 G4double energyLimitLow = globals->PhysicsEnergyLimitLow();
635 G4double energyLimitHigh = globals->PhysicsEnergyLimitHigh();
639 if (setEnergyLimitLow || setEnergyLimitHigh)
641 auto table = G4ProductionCutsTable::GetProductionCutsTable();
642 G4double defaultEnergyLimitLow = table->GetLowEdgeEnergy();
643 G4double defaultEnergyLimitHigh = table->GetHighEdgeEnergy();
644 G4double elLow = setEnergyLimitLow ? energyLimitLow : defaultEnergyLimitLow;
645 G4double elHigh = setEnergyLimitHigh ? energyLimitHigh : defaultEnergyLimitHigh;
646 table->SetEnergyRange(elLow, elHigh);
647 if (setEnergyLimitLow)
649 G4cout << __METHOD_NAME__ <<
"set EM physics low energy limit: "
650 << elLow/CLHEP::MeV <<
" MeV" << G4endl;
652 if (setEnergyLimitHigh)
654 G4cout << __METHOD_NAME__ <<
"set high energy limit: "
655 << elHigh/CLHEP::TeV <<
" TeV" << G4endl;
656 if (elHigh > G4EmParameters::Instance()->MaxKinEnergy())
658 G4cout << __METHOD_NAME__ <<
"set EM physics Ek limit to " << elHigh/CLHEP::TeV <<
" TeV" << G4endl;
659 G4EmParameters::Instance()->SetMaxEnergy(elHigh);
661#if G4VERSION_NUMBER > 1069
663 if (elHigh > G4HadronicParameters::Instance()->GetMaxEnergy())
665 G4cout << __METHOD_NAME__ <<
"set hadronic physics Ek limit to " << elHigh/CLHEP::TeV <<
" TeV" << G4endl;
666 G4HadronicParameters::Instance()->SetMaxEnergy(elHigh);
673#if G4VERSION_NUMBER > 1049
688 G4double warningEnergy = 1.0 * CLHEP::kiloelectronvolt;
689 G4double importantEnergy = 10.0 * CLHEP::kiloelectronvolt;
690 G4double numberOfTrials = 1500;
692 auto transport = transportPair.first;
693 auto coupledTransport = transportPair.second;
698 transport->SetThresholdWarningEnergy(warningEnergy);
699 transport->SetThresholdImportantEnergy(importantEnergy);
700 transport->SetThresholdTrials(numberOfTrials);
702 else if(coupledTransport)
705 coupledTransport->SetThresholdWarningEnergy(warningEnergy);
706 coupledTransport->SetThresholdImportantEnergy(importantEnergy);
707 coupledTransport->SetThresholdTrials(numberOfTrials);
713 const auto* partPM = particleDef->GetProcessManager();
715 G4VProcess* partTransport = partPM->GetProcess(
"Transportation");
716 auto transport =
dynamic_cast<G4Transportation*
>(partTransport);
718 partTransport = partPM->GetProcess(
"CoupledTransportation");
719 auto coupledTransport =
dynamic_cast<G4CoupledTransportation*
>(partTransport);
721 return std::make_pair(transport, coupledTransport);
Copy of crystal channelling em physics from Geant4.
General exception with possible name of object and message.
A class that holds global options and constants.
static BDSGlobalConstants * Instance()
Access method.
Class to parse an ion particle definition.
Modular physics list based on Geant4 constructors.
Wrapper for particle definition.
G4ParticleDefinition * ParticleDefinition() const
Accessor.
Channelling physics process.
Physics processes required for user tracking limits.
Electromagnetic dissociation for high energy ions.
High energy muon processes.
std::string particle
beam parameters
double beamKineticEnergy
beam parameters
double beamEnergy
beam parameters
double Ek0
initial beam centroid
std::string beamParticleName
beam parameters
double P0
initial beam centroid
double beamMomentum
beam parameters
double E0
initial beam centroid
double get_value(std::string name) const
get method (only for doubles)
bool HasBeenSet(const std::string &name) const
Whether a parameter has been set using the set_value method or not.
List with Efficient Lookup.
void SetRangeCuts(G4VModularPhysicsList *physicsList, G4int verbosity=1)
void ConstructBeamParticleG4(const G4String &name)
void PrintDefinedParticles()
G4bool StrContains(const G4String &str, const G4String &test)
Utility function to simplify lots of syntax changes for pedantic g4 changes.
std::pair< G4Transportation *, G4CoupledTransportation * > FindTransportation(const G4ParticleDefinition *particleDef)
Taken from Geant4 field01 example. Get the two possible transportation processes.
G4String GetFullPath(G4String filename, bool excludeNameFromPath=false, bool useCWDForPrefix=false)
BDSParticleDefinition * ConstructParticleDefinition(const G4String &particleNameIn, G4double totalEnergyIn, G4double kineticEnergyIn, G4double momentumIn, G4double ffact=1)
G4String LowerCase(const G4String &str)
Utility function to simplify lots of syntax changes for pedantic g4 changes.
void ConflictingParametersSet(const GMAD::Beam &beamDefinition, const std::set< std::string > &keys, G4int nSet, G4bool warnZeroParamsSet=true, const G4String &unitString="")
Throw an exception if too few or too many parameters are set for the supplied keys.
void PrintPrimaryParticleProcesses(const G4String &primaryParticleName)
void BuildMuonBiasing(G4VModularPhysicsList *physicsList)
Build muon splitting biasing and wrap the various processes in the physics list.
void CheckAndSetEnergyValidityRange()
void FixGeant105ThreshholdsForBeamParticle(const BDSParticleDefinition *particleDefinition)
Apply FixGeant105ThreshholdsForParticle to the beam particle definition.
void FixGeant105ThreshholdsForParticle(const G4ParticleDefinition *particleDefinition)
void ConstructExtendedParticleSet()
void ConstructDesignAndBeamParticle(const GMAD::Beam &beamDefinition, G4double ffact, BDSParticleDefinition *&designParticle, BDSParticleDefinition *&beamParticle, G4bool &beamDifferentFromDesignParticle)
G4bool IsFinite(G4double value, G4double tolerance=std::numeric_limits< double >::epsilon())
G4int NBeamParametersSet(const GMAD::Beam &beamDefinition, const std::set< std::string > &keys)
Count how many out of the set of keys in a beam instance are set.
G4VModularPhysicsList * BuildPhysics(const G4String &physicsList, G4int verbosity=1)
G4bool IsIon(const G4ParticleDefinition *particle)
Whether a particle is an ion. A proton is counted NOT as an ion.
void ConstructMinimumParticleSet()
G4VModularPhysicsList * ChannellingPhysicsComplete(G4bool useEMD=false, G4bool regular=false, G4bool em4=false, G4bool emss=false)
Build the physics required for channelling to work correctly.