BDSModularPhysicsList.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2024.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSMODULARPHYSICSLIST_H
#define BDSMODULARPHYSICSLIST_H
 
#include "G4VModularPhysicsList.hh"
 
#include <map>
#include <string>
#include <vector>
 
class BDSGlobalConstants;
class BDSModularPhysicsList;
class BDSParticleDefinition;
class G4OpticalPhysics;
class G4VPhysicsConstructor;
 
namespace GMAD {
  class PhysicsBiasing;
  template<typename T> class FastList;
}
 
class BDSModularPhysicsList: public G4VModularPhysicsList
{
public:
  BDSModularPhysicsList() = delete;
  explicit BDSModularPhysicsList(const G4String& physicsList);
  virtual ~BDSModularPhysicsList();
 
  virtual void ConstructParticle();
 
  virtual void ConstructProcess();
  
  void Print();
 
  inline G4bool UsingIons() const {return usingIons;}
 
private:
  void ConstructAllLeptons();
 
  void ConstructAllShortLived();
 
  void ConstructAllMesons();
 
  void ConstructAllBaryons();
 
  void ConstructAllIons();
 
  G4bool constructedAllLeptons;
  G4bool constructedAllShortLived;
  G4bool constructedAllMesons;
  G4bool constructedAllBaryons;
  G4bool constructedAllIons;
 
  void ParsePhysicsList(const G4String& physListName);
 
  void ConfigurePhysics();
 
  void ConfigureOptical();
 
  void CheckIncompatiblePhysics(const G4String& singlePhysicsIn) const;
  
  typedef void(BDSModularPhysicsList::*Constructor)();
 
  G4String temporaryName;
 
  std::map<std::string, Constructor> physicsConstructors;
 
  std::vector<G4String> physicsLists;
 
  std::map<G4String, G4bool> physicsActivated;
 
  std::map<G4String, std::vector<G4String> > incompatible;
 
  std::map<G4String, G4String> aliasToOriginal;
  
  G4OpticalPhysics* opticalPhysics;
  std::vector<G4VPhysicsConstructor*> constructors;
 
  BDSGlobalConstants* globals;
 
  G4bool emWillBeUsed;
 
  mutable G4bool usingIons;
 
  G4bool particlesConstructed;
 
  void AllParticles();
  
  void AnnihiToMuMu();
  void ChargeExchange();
  void Cherenkov();
  void CutsAndLimits();
  void Decay();
  void DecayRadioactive();
  void Em();
  void EmExtra();
  void EmLivermore();
  void EmLivermorePolarised();
  void EmLowEP();
  void EmPenelope();
  void EmSS();
  void EmWVI();
  void Em1();
  void Em2();
  void Em3();
  void Em4();
  void FTFPBERT();
  void FTFPBERTHP();
  void FTFBIC();
  void GammaToMuMu();
  void HadronicElastic();
  void HadronicElasticD();
  void HadronicElasticH();
  void HadronicElasticHP();
  void HadronicElasticLEND();
  void HadronicElasticXS();
  void Ion();
  void IonBinary();
  void IonElastic();
  void IonElasticQMD();
  void IonEMDissociation();
  void IonINCLXX();
  void Ionisation();
  void LaserWire();
  void LaserPhotoDetachment();
  void LaserComptonScattering();
  void LaserCumulativePhotoDetachment();
  void LaserCumulativeCompton();
  void LaserIonExcitation();
  void Muon();
  void MuonInelastic();
  void NeutronTrackingCut();
  void Optical();
  void QGSPBERT();
  void QGSPBERTHP();
  void QGSPBIC();
  void QGSPBICHP();
  void Shielding();
  void Stopping();
  void SynchRad();
#if G4VERSION_NUMBER > 1019
  void EmGS();
  void DNAChemistry();
#endif
#if G4VERSION_NUMBER > 1020
  void DecaySpin();
#endif
#if G4VERSION_NUMBER > 1022
  void IonPHP();
#endif
#if G4VERSION_NUMBER > 1029
  void DecayMuonicAtom();
#endif
#if G4VERSION_NUMBER > 1039
  void Channelling();
  void DNA();
  void Radioactivation();
  void ShieldingLEND();
#endif
#if G4VERSION_NUMBER > 1119
  void XrayReflection();
#endif
};
 
#endif