Tue 21 Aug 15:36:34 BST 2018

1 files changed, 327 insertions, 0 deletions

diff --git a/include/G4BeamTestHadronPhysics.h.backup b/include/G4BeamTestHadronPhysics.h.backup
new file mode 100644
index 0000000..dc2bf67
--- /dev/null
+++ b/include/G4BeamTestHadronPhysics.h.backup
@@ -0,0 +1,327 @@
+#ifndef G4TANKRESPONSE_G4BEAMTESTHADRONPHYSICS_H_INCLUDED
+#define G4TANKRESPONSE_G4BEAMTESTHADRONPHYSICS_H_INCLUDED
+
+#include <vector>
+
+#include <globals.hh>
+
+#include <G4VPhysicsConstructor.hh>
+
+#include <G4hMultipleScattering.hh>
+#include <G4hIonisation.hh>
+
+#include <G4HadronElasticProcess.hh>
+#include <G4HadronFissionProcess.hh>
+#include <G4HadronCaptureProcess.hh>
+
+#include <G4PionPlusInelasticProcess.hh>
+#include <G4PionMinusInelasticProcess.hh>
+#include <G4KaonPlusInelasticProcess.hh>
+#include <G4KaonZeroSInelasticProcess.hh>
+#include <G4KaonZeroLInelasticProcess.hh>
+#include <G4KaonMinusInelasticProcess.hh>
+#include <G4ProtonInelasticProcess.hh>
+#include <G4AntiProtonInelasticProcess.hh>
+#include <G4NeutronInelasticProcess.hh>
+#include <G4AntiNeutronInelasticProcess.hh>
+#include <G4LambdaInelasticProcess.hh>
+#include <G4AntiLambdaInelasticProcess.hh>
+#include <G4SigmaPlusInelasticProcess.hh>
+#include <G4SigmaMinusInelasticProcess.hh>
+#include <G4AntiSigmaPlusInelasticProcess.hh>
+#include <G4AntiSigmaMinusInelasticProcess.hh>
+#include <G4XiZeroInelasticProcess.hh>
+#include <G4XiMinusInelasticProcess.hh>
+#include <G4AntiXiZeroInelasticProcess.hh>
+#include <G4AntiXiMinusInelasticProcess.hh>
+#include <G4DeuteronInelasticProcess.hh>
+#include <G4TritonInelasticProcess.hh>
+#include <G4AlphaInelasticProcess.hh>
+#include <G4OmegaMinusInelasticProcess.hh>
+#include <G4AntiOmegaMinusInelasticProcess.hh>
+
+// Low-energy Models
+#include <G4LElastic.hh>   
+#include <G4LFission.hh>
+#include <G4LCapture.hh>
+
+#include <G4LEPionPlusInelastic.hh>
+#include <G4LEPionMinusInelastic.hh>
+#include <G4LEKaonPlusInelastic.hh>
+#include <G4LEKaonZeroSInelastic.hh>
+#include <G4LEKaonZeroLInelastic.hh>
+#include <G4LEKaonMinusInelastic.hh>
+#include <G4LEProtonInelastic.hh>
+#include <G4LEAntiProtonInelastic.hh>
+#include <G4LENeutronInelastic.hh>
+#include <G4LEAntiNeutronInelastic.hh>
+#include <G4LELambdaInelastic.hh>
+#include <G4LEAntiLambdaInelastic.hh>
+#include <G4LESigmaPlusInelastic.hh>
+#include <G4LESigmaMinusInelastic.hh>
+#include <G4LEAntiSigmaPlusInelastic.hh>
+#include <G4LEAntiSigmaMinusInelastic.hh>
+#include <G4LEXiZeroInelastic.hh>
+#include <G4LEXiMinusInelastic.hh>
+#include <G4LEAntiXiZeroInelastic.hh>
+#include <G4LEAntiXiMinusInelastic.hh>
+#include <G4LEDeuteronInelastic.hh>
+#include <G4LETritonInelastic.hh>
+#include <G4LEAlphaInelastic.hh>
+#include <G4LEOmegaMinusInelastic.hh>
+#include <G4LEAntiOmegaMinusInelastic.hh>
+
+// High-energy Models
+
+#include <G4HEPionPlusInelastic.hh>
+#include <G4HEPionMinusInelastic.hh>
+#include <G4HEKaonPlusInelastic.hh>
+#include <G4HEKaonZeroInelastic.hh>
+#include <G4HEKaonZeroInelastic.hh>
+#include <G4HEKaonMinusInelastic.hh>
+#include <G4HEProtonInelastic.hh>
+#include <G4HEAntiProtonInelastic.hh>
+#include <G4HENeutronInelastic.hh>
+#include <G4HEAntiNeutronInelastic.hh>
+#include <G4HELambdaInelastic.hh>
+#include <G4HEAntiLambdaInelastic.hh>
+#include <G4HESigmaPlusInelastic.hh>
+#include <G4HESigmaMinusInelastic.hh>
+#include <G4HEAntiSigmaPlusInelastic.hh>
+#include <G4HEAntiSigmaMinusInelastic.hh>
+#include <G4HEXiZeroInelastic.hh>
+#include <G4HEXiMinusInelastic.hh>
+#include <G4HEAntiXiZeroInelastic.hh>
+#include <G4HEAntiXiMinusInelastic.hh>
+#include <G4HEOmegaMinusInelastic.hh>
+#include <G4HEAntiOmegaMinusInelastic.hh>
+
+// Stopping processes
+#include <G4AntiProtonAnnihilationAtRest.hh>
+#include <G4AntiNeutronAnnihilationAtRest.hh>
+
+#ifdef TRIUMF_STOP_PIMINUS
+#include <G4PionMinusAbsorptionAtRest.hh>
+#else
+#include <G4PiMinusAbsorptionAtRest.hh>
+#endif
+#ifdef TRIUMF_STOP_KMINUS
+#include <G4KaonMinusAbsorption.hh>
+#else
+#include <G4KaonMinusAbsorptionAtRest.hh>
+#endif
+
+// quark gluon string model with chips afterburner.
+#include <G4TheoFSGenerator.hh>
+#include <G4ExcitationHandler.hh>
+#include <G4PreCompoundModel.hh>
+#include <G4GeneratorPrecompoundInterface.hh>
+#include <G4QGSModel.hh>
+#include <G4QGSParticipants.hh>
+#include <G4QGSMFragmentation.hh>
+#include <G4ExcitedStringDecay.hh>
+
+/**
+   @brief Hadron physics. Used only if Geant4 version is earlier than 4.10.
+   @author GEANT4/Peter Niessen
+   @date Sat Jul 24 23:53:47 EDT 2004
+
+   Many processes. You're encouraged to check the source.
+*/
+class G4BeamTestHadronPhysics : public G4VPhysicsConstructor {
+public: 
+
+  /**
+   * The constructor
+   */
+  G4BeamTestHadronPhysics();
+
+  /**
+   * The virtual destructor
+   */
+  virtual ~G4BeamTestHadronPhysics();
+
+  /**
+   * This method will be invoked in the Construct() method. 
+   * each particle type will be instantiated
+   */
+  virtual void ConstructParticle();
+ 
+  /**
+   * This method will be invoked in the Construct() method.
+   * each physics process will be instantiated and
+   * registered to the process manager of each particle type 
+   */
+  virtual void ConstructProcess();
+
+protected:
+
+  // Elastic Process
+  G4HadronElasticProcess elasticProcess_;
+  G4LElastic *elasticModel_;
+  
+  // Pi + 
+  G4PionPlusInelasticProcess pionPlusInelastic_;
+  G4LEPionPlusInelastic *lePionPlusModel_;
+  G4hMultipleScattering pionPlusMult_;
+  G4hIonisation pionPlusIonisation_;
+
+  // Pi -
+  G4PionMinusInelasticProcess pionMinusInelastic_;
+  G4LEPionMinusInelastic *lePionMinusModel_;
+  G4hMultipleScattering pionMinusMult_;
+  G4hIonisation pionMinusIonisation_;
+#ifdef TRIUMF_STOP_PIMINUS
+  G4PionMinusAbsorptionAtRest pionMinusAbsorption_;
+#else
+  G4PiMinusAbsorptionAtRest pionMinusAbsorption_;
+#endif
+
+  // pi+ and pi-
+  G4TheoFSGenerator theoModel_;
+  G4ExcitationHandler excitationHandler_;
+  G4PreCompoundModel *preEquilib_;
+  G4GeneratorPrecompoundInterface cascade_;
+  G4QGSModel< G4QGSParticipants > stringModel_;
+  G4QGSMFragmentation fragmentation_;
+  G4ExcitedStringDecay *stringDecay_;
+
+  // K + 
+  G4KaonPlusInelasticProcess kaonPlusInelastic_;
+  G4LEKaonPlusInelastic *leKaonPlusModel_;
+  G4HEKaonPlusInelastic *heKaonPlusModel_;
+  G4hMultipleScattering kaonPlusMult_;
+  G4hIonisation kaonPlusIonisation_;
+	
+  // K -
+  G4KaonMinusInelasticProcess kaonMinusInelastic_;
+  G4LEKaonMinusInelastic *leKaonMinusModel_;
+  G4HEKaonMinusInelastic *heKaonMinusModel_;
+  G4hMultipleScattering kaonMinusMult_;
+  G4hIonisation kaonMinusIonisation_;
+#ifdef TRIUMF_STOP_KMINUS
+  G4KaonMinusAbsorption kaonMinusAbsorption_;
+#else
+  G4PiMinusAbsorptionAtRest kaonMinusAbsorption_;
+#endif
+
+  // K0L
+  G4KaonZeroLInelasticProcess kaonZeroLInelastic_;
+  G4LEKaonZeroLInelastic *leKaonZeroLModel_;
+  G4HEKaonZeroInelastic *heKaonZeroLModel_;
+
+  // K0S
+  G4KaonZeroSInelasticProcess kaonZeroSInelastic_;
+  G4LEKaonZeroSInelastic *leKaonZeroSModel_;
+  G4HEKaonZeroInelastic *heKaonZeroSModel_;
+
+  // Proton
+  G4ProtonInelasticProcess protonInelastic_;
+  G4LEProtonInelastic *leProtonModel_;
+  G4HEProtonInelastic *heProtonModel_;
+  G4hMultipleScattering protonMult_;
+  G4hIonisation protonIonisation_;
+ 
+  // anti-proton
+  G4AntiProtonInelasticProcess antiProtonInelastic_;
+  G4LEAntiProtonInelastic *leAntiProtonModel_;
+  G4HEAntiProtonInelastic *heAntiProtonModel_;
+  G4hMultipleScattering antiProtonMult_;
+  G4hIonisation antiProtonIonisation_;
+  G4AntiProtonAnnihilationAtRest antiProtonAnnihilation_;
+    
+  // neutron
+  G4NeutronInelasticProcess neutronInelastic_;
+  G4LENeutronInelastic *leNeutronModel_;
+  G4HENeutronInelastic *heNeutronModel_;
+  G4HadronFissionProcess neutronFission_;
+  G4LFission *neutronFissionModel_;
+  G4HadronCaptureProcess neutronCapture_;
+  G4LCapture *neutronCaptureModel_;
+
+
+  // anti-neutron
+  G4AntiNeutronInelasticProcess antiNeutronInelastic_;
+  G4LEAntiNeutronInelastic *leAntiNeutronModel_;
+  G4HEAntiNeutronInelastic *heAntiNeutronModel_;
+  G4AntiNeutronAnnihilationAtRest antiNeutronAnnihilation_;
+   
+  // Lambda
+  G4LambdaInelasticProcess lambdaInelastic_;
+  G4LELambdaInelastic *leLambdaModel_;
+  G4HELambdaInelastic *heLambdaModel_;
+  
+  // AntiLambda
+  G4AntiLambdaInelasticProcess antiLambdaInelastic_;
+  G4LEAntiLambdaInelastic *leAntiLambdaModel_;
+  G4HEAntiLambdaInelastic *heAntiLambdaModel_;
+  
+  // SigmaMinus
+  G4SigmaMinusInelasticProcess sigmaMinusInelastic_;
+  G4LESigmaMinusInelastic *leSigmaMinusModel_;
+  G4HESigmaMinusInelastic *heSigmaMinusModel_;
+  G4hMultipleScattering sigmaMinusMult_;
+  G4hIonisation sigmaMinusIonisation_;
+  
+  // AntiSigmaMinus
+  G4AntiSigmaMinusInelasticProcess antiSigmaMinusInelastic_;
+  G4LEAntiSigmaMinusInelastic *leAntiSigmaMinusModel_;
+  G4HEAntiSigmaMinusInelastic *heAntiSigmaMinusModel_;
+  G4hMultipleScattering antiSigmaMinusMult_;
+  G4hIonisation antiSigmaMinusIonisation_;
+   
+  // SigmaPlus
+  G4SigmaPlusInelasticProcess sigmaPlusInelastic_;
+  G4LESigmaPlusInelastic *leSigmaPlusModel_;
+  G4HESigmaPlusInelastic *heSigmaPlusModel_;
+  G4hMultipleScattering sigmaPlusMult_;
+  G4hIonisation sigmaPlusIonisation_;
+  
+  // AntiSigmaPlus
+  G4AntiSigmaPlusInelasticProcess antiSigmaPlusInelastic_;
+  G4LEAntiSigmaPlusInelastic *leAntiSigmaPlusModel_;
+  G4HEAntiSigmaPlusInelastic *heAntiSigmaPlusModel_;
+  G4hMultipleScattering antiSigmaPlusMult_;
+  G4hIonisation antiSigmaPlusIonisation_;
+  
+  // XiZero
+  G4XiZeroInelasticProcess xiZeroInelastic_;
+  G4LEXiZeroInelastic *leXiZeroModel_;
+  G4HEXiZeroInelastic *heXiZeroModel_;
+  
+  // AntiXiZero
+  G4AntiXiZeroInelasticProcess  antiXiZeroInelastic_;
+  G4LEAntiXiZeroInelastic*  leAntiXiZeroModel_;
+  G4HEAntiXiZeroInelastic*  heAntiXiZeroModel_;
+  
+  // XiMinus
+  G4XiMinusInelasticProcess xiMinusInelastic_;
+  G4LEXiMinusInelastic *leXiMinusModel_;
+  G4HEXiMinusInelastic *heXiMinusModel_;
+  G4hMultipleScattering xiMinusMult_;
+  G4hIonisation xiMinusIonisation_;
+
+  // AntiXiMinus
+  G4AntiXiMinusInelasticProcess antiXiMinusInelastic_;
+  G4LEAntiXiMinusInelastic *leAntiXiMinusModel_;
+  G4HEAntiXiMinusInelastic *heAntiXiMinusModel_;
+  G4hMultipleScattering antiXiMinusMult_;
+  G4hIonisation antiXiMinusIonisation_;
+  
+  // OmegaMinus
+  G4OmegaMinusInelasticProcess omegaMinusInelastic_;
+  G4LEOmegaMinusInelastic *leOmegaMinusModel_;
+  G4HEOmegaMinusInelastic *heOmegaMinusModel_;
+  G4hMultipleScattering omegaMinusMult_;
+  G4hIonisation omegaMinusIonisation_;
+   
+  // AntiOmegaMinus
+  G4AntiOmegaMinusInelasticProcess antiOmegaMinusInelastic_;
+  G4LEAntiOmegaMinusInelastic *leAntiOmegaMinusModel_;
+  G4HEAntiOmegaMinusInelastic *heAntiOmegaMinusModel_;
+  G4hMultipleScattering antiOmegaMinusMult_;
+  G4hIonisation antiOmegaMinusIonisation_;
+};
+
+
+#endif // G4TANKRESPONSE_G4BEAMTESTHADRONPHYSICS_H_INCLUDED