[u/mrichter/AliRoot.git] / TGeant4 / TG4PhysicsConstructorHadron.h

// $Id$
// Category: physics
//
// Author: I. Hrivnacova
//
// Class TG4PhysicsConstructorHadron
// ---------------------------------
// Constructor of hadron physics.
// According to ExN04HadronPhysics.hh, GEANT4 tag Name: geant4-03-02

#ifndef TG4_PHYSICS_CONSTRUCTOR_HADRON_H
#define TG4_PHYSICS_CONSTRUCTOR_HADRON_H

#include "TG4VPhysicsConstructor.h"

#include <g4std/vector>
#include <globals.hh>

#include <G4MultipleScattering.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

class TG4PhysicsConstructorHadron: public TG4VPhysicsConstructor
{
  typedef G4std::vector<G4VProcess*>  ProcessVector;

  public:
    TG4PhysicsConstructorHadron(const G4String& name = "Hadron");
    TG4PhysicsConstructorHadron(G4int verboseLevel, 
                                G4bool setEM, G4bool setHadron,
                                const G4String& name = "Hadron");
    // --> protected
    // TG4PhysicsConstructorHadron(const TG4PhysicsConstructorHadron& right);
    virtual ~TG4PhysicsConstructorHadron();

  protected:
    TG4PhysicsConstructorHadron(const TG4PhysicsConstructorHadron& right);
    
    // operators
    TG4PhysicsConstructorHadron& operator=(
                                const TG4PhysicsConstructorHadron& right);
  
    // methods
          // construct particle and physics
    virtual void ConstructParticle();
    virtual void ConstructProcess();

    // data members
         // Elastic Process
    G4HadronElasticProcess  fElasticProcess;            //hadron elastic process
    G4LElastic*             fElasticModel;              //elastic model
  
         // Pi + 
    G4PionPlusInelasticProcess fPionPlusInelastic;      //pi+ inel process
    G4LEPionPlusInelastic*     fLEPionPlusModel;        //pi+ LE inel model
    G4HEPionPlusInelastic*     fHEPionPlusModel;        //pi+ HE inel model
    G4MultipleScattering       fPionPlusMult;           //pi+ msc
    G4hIonisation              fPionPlusIonisation;     //pi+ ionisation

         // Pi -
    G4PionMinusInelasticProcess  fPionMinusInelastic;   //pi- inel process
    G4LEPionMinusInelastic*      fLEPionMinusModel;     //pi- LE inel model
    G4HEPionMinusInelastic*      fHEPionMinusModel;     //pi- HE inel model
    G4MultipleScattering         fPionMinusMult;        //pi- msc
    G4hIonisation                fPionMinusIonisation;  //pi- ionisation
#ifdef TRIUMF_STOP_PIMINUS
    G4PionMinusAbsorptionAtRest  fPionMinusAbsorption;  //pi- absorption
#else
    G4PiMinusAbsorptionAtRest    fPionMinusAbsorption;  //pi- absorption
#endif

         // K + 
    G4KaonPlusInelasticProcess  fKaonPlusInelastic;     //kaon+ inel process
    G4LEKaonPlusInelastic*      fLEKaonPlusModel;       //kaon+ LE inel model
    G4HEKaonPlusInelastic*      fHEKaonPlusModel;       //kaon+ HE inel model
    G4MultipleScattering        fKaonPlusMult;          //kaon+ msc
    G4hIonisation               fKaonPlusIonisation;    //kaon+ ionisation

         // K -
    G4KaonMinusInelasticProcess  fKaonMinusInelastic;   //kaon- inel process
    G4LEKaonMinusInelastic*      fLEKaonMinusModel;     //kaon- LE inel model
    G4HEKaonMinusInelastic*      fHEKaonMinusModel;     //kaon- HE inel model
    G4MultipleScattering         fKaonMinusMult;        //kaon- msc
    G4hIonisation                fKaonMinusIonisation;  //kaon- ionisation
#ifdef TRIUMF_STOP_KMINUS
    G4KaonMinusAbsorption        fKaonMinusAbsorption;  //kaon- absorption
#else
    G4PiMinusAbsorptionAtRest    fKaonMinusAbsorption;  //kaon- absorption
#endif

        // K0L
    G4KaonZeroLInelasticProcess  fKaonZeroLInelastic;   //kaon0 inel process
    G4LEKaonZeroLInelastic*      fLEKaonZeroLModel;     //kaon0 LE inel model
    G4HEKaonZeroInelastic*       fHEKaonZeroLModel;     //kaon0 HE inel model

        // K0S
    G4KaonZeroSInelasticProcess  fKaonZeroSInelastic;   //kaon0S inel process
    G4LEKaonZeroSInelastic*      fLEKaonZeroSModel;     //kaon0S LE inel model
    G4HEKaonZeroInelastic*       fHEKaonZeroSModel;     //kaon0S HE inel mode

        // Proton
    G4ProtonInelasticProcess  fProtonInelastic;         //p inel process
    G4LEProtonInelastic*      fLEProtonModel;           //p LE inel model
    G4HEProtonInelastic*      fHEProtonModel;           //p HE inel model
    G4MultipleScattering      fProtonMult;              //p msc
    G4hIonisation             fProtonIonisation;        //p ionisation
 
        // anti-proton
    G4AntiProtonInelasticProcess    fAntiProtonInelastic; //p_bar inel process
    G4LEAntiProtonInelastic*        fLEAntiProtonModel;   //p_bar LE inel model
    G4HEAntiProtonInelastic*        fHEAntiProtonModel;   //p_bar HE inel model
    G4MultipleScattering            fAntiProtonMult;      //p_bar msc
    G4hIonisation                   fAntiProtonIonisation;//p_bar ionisation
    G4AntiProtonAnnihilationAtRest  fAntiProtonAnnihilation;//p_bar annihilation
    
       // neutron
    G4NeutronInelasticProcess  fNeutronInelastic;       //n inel process
    G4LENeutronInelastic*      fLENeutronModel;         //n LE inel model
    G4HENeutronInelastic*      fHENeutronModel;         //n HE inel model
    G4HadronFissionProcess     fNeutronFission;         //n fission
    G4LFission*                fNeutronFissionModel;    //n fission model
    G4HadronCaptureProcess     fNeutronCapture;         //n capture
    G4LCapture*                fNeutronCaptureModel;    //n capture model

       // anti-neutron
    G4AntiNeutronInelasticProcess    fAntiNeutronInelastic;//n_bar inel process
    G4LEAntiNeutronInelastic*        fLEAntiNeutronModel;  //n_bar LE inel model
    G4HEAntiNeutronInelastic*        fHEAntiNeutronModel;  //n_bar HE inel model
    G4AntiNeutronAnnihilationAtRest  fAntiNeutronAnnihilation;//n_bar ionisation
     
       // Lambda
    G4LambdaInelasticProcess  fLambdaInelastic;         //lambda inel process
    G4LELambdaInelastic*      fLELambdaModel;           //lambda LE inel model
    G4HELambdaInelastic*      fHELambdaModel;           //lambda HE inel model
  
       // AntiLambda
    G4AntiLambdaInelasticProcess  fAntiLambdaInelastic; //lambda_bar inel process
    G4LEAntiLambdaInelastic*      fLEAntiLambdaModel;   //lambda_bar LE inel model
    G4HEAntiLambdaInelastic*      fHEAntiLambdaModel;   //lambda_bar HE inel model
  
       // SigmaMinus
    G4SigmaMinusInelasticProcess  fSigmaMinusInelastic; //sigma- inel process
    G4LESigmaMinusInelastic*      fLESigmaMinusModel;   //sigma- LE inel model
    G4HESigmaMinusInelastic*      fHESigmaMinusModel;   //sigma- HE inel model
    G4MultipleScattering          fSigmaMinusMult;      //sigma- msc
    G4hIonisation                 fSigmaMinusIonisation;//sigma- ionisation
  
       // AntiSigmaMinus
    G4AntiSigmaMinusInelasticProcess  fAntiSigmaMinusInelastic; //sigma-_bar inel process
    G4LEAntiSigmaMinusInelastic*      fLEAntiSigmaMinusModel;   //sigma-_bar LE inel model
    G4HEAntiSigmaMinusInelastic*      fHEAntiSigmaMinusModel;   //sigma-_bar HE inel model
    G4MultipleScattering              fAntiSigmaMinusMult;      //sigma-_bar msc
    G4hIonisation                     fAntiSigmaMinusIonisation;//sigma-_bar ionisation
   
       // SigmaPlus
    G4SigmaPlusInelasticProcess  fSigmaPlusInelastic;   //sigma+ inel process
    G4LESigmaPlusInelastic*      fLESigmaPlusModel;     //sigma+ LE inel model
    G4HESigmaPlusInelastic*      fHESigmaPlusModel;     //sigma+ HE inel model
    G4MultipleScattering         fSigmaPlusMult;        //sigma+ msc
    G4hIonisation                fSigmaPlusIonisation;  //sigma+ ionisation
  
       // AntiSigmaPlus
    G4AntiSigmaPlusInelasticProcess  fAntiSigmaPlusInelastic;  //sigma+_bar inel process
    G4LEAntiSigmaPlusInelastic*      fLEAntiSigmaPlusModel;    //sigma+_bar LE inel model
    G4HEAntiSigmaPlusInelastic*      fHEAntiSigmaPlusModel;    //sigma+_bar HE inel model
    G4MultipleScattering             fAntiSigmaPlusMult;       //sigma+_bar msc
    G4hIonisation                    fAntiSigmaPlusIonisation; //sigma+_bar ionisation
  
      // XiZero
    G4XiZeroInelasticProcess  fXiZeroInelastic;        //xi0 inel process
    G4LEXiZeroInelastic*      fLEXiZeroModel;          //xi0 LE inel model
    G4HEXiZeroInelastic*      fHEXiZeroModel;          //xi0 HE inel model
  
      // AntiXiZero
    G4AntiXiZeroInelasticProcess  fAntiXiZeroInelastic;//xi0_bar inel process
    G4LEAntiXiZeroInelastic*      fLEAntiXiZeroModel;  //xi0_bar LE inel model
    G4HEAntiXiZeroInelastic*      fHEAntiXiZeroModel;  //xi0_bar HE inel model
  
      // XiMinus
    G4XiMinusInelasticProcess  fXiMinusInelastic;      //xi- inel process
    G4LEXiMinusInelastic*      fLEXiMinusModel;        //xi- LE inel model
    G4HEXiMinusInelastic*      fHEXiMinusModel;        //xi- HE inel model
    G4MultipleScattering       fXiMinusMult;           //xi- msc
    G4hIonisation              fXiMinusIonisation;     //xi- ionisation

      // AntiXiMinus
    G4AntiXiMinusInelasticProcess  fAntiXiMinusInelastic; //xi-_bar inel process
    G4LEAntiXiMinusInelastic*      fLEAntiXiMinusModel;   //xi-_bar LE inel model
    G4HEAntiXiMinusInelastic*      fHEAntiXiMinusModel;   //xi-_bar HE inel model
    G4MultipleScattering           fAntiXiMinusMult;      //xi-_bar msc
    G4hIonisation                  fAntiXiMinusIonisation;//xi-_bar ionisation
  
      // OmegaMinus
    G4OmegaMinusInelasticProcess  fOmegaMinusInelastic;   //omega- inel process
    G4LEOmegaMinusInelastic*      fLEOmegaMinusModel;     //omega- LE inel model
    G4HEOmegaMinusInelastic*      fHEOmegaMinusModel;     //omega- HE inel model
    G4MultipleScattering          fOmegaMinusMult;        //omega- msc
    G4hIonisation                 fOmegaMinusIonisation;  //omega- ionisation
   
      // AntiOmegaMinus
    G4AntiOmegaMinusInelasticProcess  fAntiOmegaMinusInelastic; //omega-_bar inel process
    G4LEAntiOmegaMinusInelastic*      fLEAntiOmegaMinusModel;   //omega-_bar LE inel model
    G4HEAntiOmegaMinusInelastic*      fHEAntiOmegaMinusModel;   //omega-_bar HE inel model
    G4MultipleScattering              fAntiOmegaMinusMult;      //omega-_bar msc
    G4hIonisation                     fAntiOmegaMinusIonisation;//omega-_bar ionisation
    
      // Other
    ProcessVector  fOtherProcesses; //other process
    

  private:
    // methods

    // EM processes
    void ConstructEMProcessForPionPlus();
    void ConstructEMProcessForPionMinus();
    void ConstructEMProcessForKaonPlus();
    void ConstructEMProcessForKaonMinus();
    void ConstructEMProcessForProton();
    void ConstructEMProcessForAntiProton();
    void ConstructEMProcessForSigmaMinus();
    void ConstructEMProcessForAntiSigmaMinus();
    void ConstructEMProcessForSigmaPlus();
    void ConstructEMProcessForAntiSigmaPlus();
    void ConstructEMProcessForXiMinus();
    void ConstructEMProcessForAntiXiMinus();
    void ConstructEMProcessForOmegaMinus();
    void ConstructEMProcessForAntiOmegaMinus();
    void ConstructEMProcessForOther();

    // hadron processes
    void ConstructHadProcessForPionPlus();
    void ConstructHadProcessForPionMinus();
    void ConstructHadProcessForKaonPlus();
    void ConstructHadProcessForKaonMinus();
    void ConstructHadProcessForKaonZeroLong();
    void ConstructHadProcessForKaonZeroShort();
    void ConstructHadProcessForProton();
    void ConstructHadProcessForAntiProton();
    void ConstructHadProcessForNeutron();
    void ConstructHadProcessForAntiNeutron();
    void ConstructHadProcessForLambda();
    void ConstructHadProcessForAntiLambda();
    void ConstructHadProcessForSigmaMinus();
    void ConstructHadProcessForAntiSigmaMinus();
    void ConstructHadProcessForSigmaPlus();
    void ConstructHadProcessForAntiSigmaPlus();
    void ConstructHadProcessForXiMinus();
    void ConstructHadProcessForAntiXiMinus();
    void ConstructHadProcessForXiZero();
    void ConstructHadProcessForAntiXiZero();
    void ConstructHadProcessForOmegaMinus();
    void ConstructHadProcessForAntiOmegaMinus();
    void ConstructHadProcessForOther();

    // data members
    G4bool  fSetEM;    //if true - EM processes are constructed
    G4bool  fSetHadron;//if true - hadron processes are constructed
};

#endif //TG4_PHYSICS_CONSTRUCTOR_HADRON_H
Commit	Line	Data
b94a914b	1	// $Id$
	2	// Category: physics
	3	//
499b353a	4	// Author: I. Hrivnacova
	5	//
	6	// Class TG4PhysicsConstructorHadron
	7	// ---------------------------------
b94a914b	8	// Constructor of hadron physics.
696e37fa	9	// According to ExN04HadronPhysics.hh, GEANT4 tag Name: geant4-03-02
b94a914b	10
	11	#ifndef TG4_PHYSICS_CONSTRUCTOR_HADRON_H
	12	#define TG4_PHYSICS_CONSTRUCTOR_HADRON_H
	13
11f7ee82	14	#include "TG4VPhysicsConstructor.h"
11f7ee82	15
696e37fa	16	#include <g4std/vector>
b94a914b	17	#include <globals.hh>
b94a914b	18
696e37fa	19	#include <G4MultipleScattering.hh>
	20	#include <G4hIonisation.hh>
	21
	22	#include <G4HadronElasticProcess.hh>
	23	#include <G4HadronFissionProcess.hh>
	24	#include <G4HadronCaptureProcess.hh>
	25
	26	#include <G4PionPlusInelasticProcess.hh>
	27	#include <G4PionMinusInelasticProcess.hh>
	28	#include <G4KaonPlusInelasticProcess.hh>
	29	#include <G4KaonZeroSInelasticProcess.hh>
	30	#include <G4KaonZeroLInelasticProcess.hh>
	31	#include <G4KaonMinusInelasticProcess.hh>
	32	#include <G4ProtonInelasticProcess.hh>
	33	#include <G4AntiProtonInelasticProcess.hh>
	34	#include <G4NeutronInelasticProcess.hh>
	35	#include <G4AntiNeutronInelasticProcess.hh>
	36	#include <G4LambdaInelasticProcess.hh>
	37	#include <G4AntiLambdaInelasticProcess.hh>
	38	#include <G4SigmaPlusInelasticProcess.hh>
	39	#include <G4SigmaMinusInelasticProcess.hh>
	40	#include <G4AntiSigmaPlusInelasticProcess.hh>
	41	#include <G4AntiSigmaMinusInelasticProcess.hh>
	42	#include <G4XiZeroInelasticProcess.hh>
	43	#include <G4XiMinusInelasticProcess.hh>
	44	#include <G4AntiXiZeroInelasticProcess.hh>
	45	#include <G4AntiXiMinusInelasticProcess.hh>
	46	#include <G4DeuteronInelasticProcess.hh>
	47	#include <G4TritonInelasticProcess.hh>
	48	#include <G4AlphaInelasticProcess.hh>
	49	#include <G4OmegaMinusInelasticProcess.hh>
	50	#include <G4AntiOmegaMinusInelasticProcess.hh>
	51
	52	// Low-energy Models
	53	#include <G4LElastic.hh>
	54	#include <G4LFission.hh>
	55	#include <G4LCapture.hh>
	56
	57	#include <G4LEPionPlusInelastic.hh>
	58	#include <G4LEPionMinusInelastic.hh>
	59	#include <G4LEKaonPlusInelastic.hh>
	60	#include <G4LEKaonZeroSInelastic.hh>
	61	#include <G4LEKaonZeroLInelastic.hh>
	62	#include <G4LEKaonMinusInelastic.hh>
	63	#include <G4LEProtonInelastic.hh>
	64	#include <G4LEAntiProtonInelastic.hh>
	65	#include <G4LENeutronInelastic.hh>
	66	#include <G4LEAntiNeutronInelastic.hh>
	67	#include <G4LELambdaInelastic.hh>
	68	#include <G4LEAntiLambdaInelastic.hh>
	69	#include <G4LESigmaPlusInelastic.hh>
	70	#include <G4LESigmaMinusInelastic.hh>
	71	#include <G4LEAntiSigmaPlusInelastic.hh>
	72	#include <G4LEAntiSigmaMinusInelastic.hh>
	73	#include <G4LEXiZeroInelastic.hh>
	74	#include <G4LEXiMinusInelastic.hh>
	75	#include <G4LEAntiXiZeroInelastic.hh>
	76	#include <G4LEAntiXiMinusInelastic.hh>
	77	#include <G4LEDeuteronInelastic.hh>
	78	#include <G4LETritonInelastic.hh>
	79	#include <G4LEAlphaInelastic.hh>
	80	#include <G4LEOmegaMinusInelastic.hh>
	81	#include <G4LEAntiOmegaMinusInelastic.hh>
	82
83	// High-energy Models
84
85	#include <G4HEPionPlusInelastic.hh>
86	#include <G4HEPionMinusInelastic.hh>
87	#include <G4HEKaonPlusInelastic.hh>
88	#include <G4HEKaonZeroInelastic.hh>
89	#include <G4HEKaonZeroInelastic.hh>
90	#include <G4HEKaonMinusInelastic.hh>
91	#include <G4HEProtonInelastic.hh>
92	#include <G4HEAntiProtonInelastic.hh>
93	#include <G4HENeutronInelastic.hh>
94	#include <G4HEAntiNeutronInelastic.hh>
95	#include <G4HELambdaInelastic.hh>
96	#include <G4HEAntiLambdaInelastic.hh>
97	#include <G4HESigmaPlusInelastic.hh>
98	#include <G4HESigmaMinusInelastic.hh>
99	#include <G4HEAntiSigmaPlusInelastic.hh>
100	#include <G4HEAntiSigmaMinusInelastic.hh>
101	#include <G4HEXiZeroInelastic.hh>
102	#include <G4HEXiMinusInelastic.hh>
103	#include <G4HEAntiXiZeroInelastic.hh>
104	#include <G4HEAntiXiMinusInelastic.hh>
105	#include <G4HEOmegaMinusInelastic.hh>
106	#include <G4HEAntiOmegaMinusInelastic.hh>
107
108	// Stopping processes
109	#include <G4AntiProtonAnnihilationAtRest.hh>
110	#include <G4AntiNeutronAnnihilationAtRest.hh>
111
112	#ifdef TRIUMF_STOP_PIMINUS
113	#include <G4PionMinusAbsorptionAtRest.hh>
114	#else
115	#include <G4PiMinusAbsorptionAtRest.hh>
116	#endif
117	#ifdef TRIUMF_STOP_KMINUS
118	#include <G4KaonMinusAbsorption.hh>
119	#else
120	#include <G4KaonMinusAbsorptionAtRest.hh>
121	#endif
122
11f7ee82	123	class TG4PhysicsConstructorHadron: public TG4VPhysicsConstructor
b94a914b	124	{
696e37fa	125	typedef G4std::vector<G4VProcess*> ProcessVector;
696e37fa	126
b94a914b	127	public:
b94a914b	128	TG4PhysicsConstructorHadron(const G4String& name = "Hadron");
11f7ee82	129	TG4PhysicsConstructorHadron(G4int verboseLevel,
	130	G4bool setEM, G4bool setHadron,
	131	const G4String& name = "Hadron");
696e37fa	132	// --> protected
696e37fa	133	// TG4PhysicsConstructorHadron(const TG4PhysicsConstructorHadron& right);
b94a914b	134	virtual ~TG4PhysicsConstructorHadron();
	135
	136	protected:
696e37fa	137	TG4PhysicsConstructorHadron(const TG4PhysicsConstructorHadron& right);
	138
	139	// operators
	140	TG4PhysicsConstructorHadron& operator=(
	141	const TG4PhysicsConstructorHadron& right);
	142
b94a914b	143	// methods
	144	// construct particle and physics
	145	virtual void ConstructParticle();
	146	virtual void ConstructProcess();
696e37fa	147
	148	// data members
	149	// Elastic Process
cc364059	150	G4HadronElasticProcess fElasticProcess; //hadron elastic process
cc364059	151	G4LElastic* fElasticModel; //elastic model
696e37fa	152
696e37fa	153	// Pi +
cc364059	154	G4PionPlusInelasticProcess fPionPlusInelastic; //pi+ inel process
	155	G4LEPionPlusInelastic* fLEPionPlusModel; //pi+ LE inel model
	156	G4HEPionPlusInelastic* fHEPionPlusModel; //pi+ HE inel model
	157	G4MultipleScattering fPionPlusMult; //pi+ msc
	158	G4hIonisation fPionPlusIonisation; //pi+ ionisation
696e37fa	159
696e37fa	160	// Pi -
cc364059	161	G4PionMinusInelasticProcess fPionMinusInelastic; //pi- inel process
	162	G4LEPionMinusInelastic* fLEPionMinusModel; //pi- LE inel model
	163	G4HEPionMinusInelastic* fHEPionMinusModel; //pi- HE inel model
	164	G4MultipleScattering fPionMinusMult; //pi- msc
	165	G4hIonisation fPionMinusIonisation; //pi- ionisation
696e37fa	166	#ifdef TRIUMF_STOP_PIMINUS
cc364059	167	G4PionMinusAbsorptionAtRest fPionMinusAbsorption; //pi- absorption
696e37fa	168	#else
cc364059	169	G4PiMinusAbsorptionAtRest fPionMinusAbsorption; //pi- absorption
696e37fa	170	#endif
	171
	172	// K +
cc364059	173	G4KaonPlusInelasticProcess fKaonPlusInelastic; //kaon+ inel process
	174	G4LEKaonPlusInelastic* fLEKaonPlusModel; //kaon+ LE inel model
	175	G4HEKaonPlusInelastic* fHEKaonPlusModel; //kaon+ HE inel model
	176	G4MultipleScattering fKaonPlusMult; //kaon+ msc
	177	G4hIonisation fKaonPlusIonisation; //kaon+ ionisation
696e37fa	178
696e37fa	179	// K -
cc364059	180	G4KaonMinusInelasticProcess fKaonMinusInelastic; //kaon- inel process
	181	G4LEKaonMinusInelastic* fLEKaonMinusModel; //kaon- LE inel model
	182	G4HEKaonMinusInelastic* fHEKaonMinusModel; //kaon- HE inel model
	183	G4MultipleScattering fKaonMinusMult; //kaon- msc
	184	G4hIonisation fKaonMinusIonisation; //kaon- ionisation
696e37fa	185	#ifdef TRIUMF_STOP_KMINUS
cc364059	186	G4KaonMinusAbsorption fKaonMinusAbsorption; //kaon- absorption
696e37fa	187	#else
cc364059	188	G4PiMinusAbsorptionAtRest fKaonMinusAbsorption; //kaon- absorption
696e37fa	189	#endif
	190
	191	// K0L
cc364059	192	G4KaonZeroLInelasticProcess fKaonZeroLInelastic; //kaon0 inel process
	193	G4LEKaonZeroLInelastic* fLEKaonZeroLModel; //kaon0 LE inel model
	194	G4HEKaonZeroInelastic* fHEKaonZeroLModel; //kaon0 HE inel model
696e37fa	195
696e37fa	196	// K0S
cc364059	197	G4KaonZeroSInelasticProcess fKaonZeroSInelastic; //kaon0S inel process
	198	G4LEKaonZeroSInelastic* fLEKaonZeroSModel; //kaon0S LE inel model
	199	G4HEKaonZeroInelastic* fHEKaonZeroSModel; //kaon0S HE inel mode
696e37fa	200
696e37fa	201	// Proton
cc364059	202	G4ProtonInelasticProcess fProtonInelastic; //p inel process
	203	G4LEProtonInelastic* fLEProtonModel; //p LE inel model
	204	G4HEProtonInelastic* fHEProtonModel; //p HE inel model
	205	G4MultipleScattering fProtonMult; //p msc
	206	G4hIonisation fProtonIonisation; //p ionisation
696e37fa	207
696e37fa	208	// anti-proton
cc364059	209	G4AntiProtonInelasticProcess fAntiProtonInelastic; //p_bar inel process
	210	G4LEAntiProtonInelastic* fLEAntiProtonModel; //p_bar LE inel model
	211	G4HEAntiProtonInelastic* fHEAntiProtonModel; //p_bar HE inel model
	212	G4MultipleScattering fAntiProtonMult; //p_bar msc
	213	G4hIonisation fAntiProtonIonisation;//p_bar ionisation
	214	G4AntiProtonAnnihilationAtRest fAntiProtonAnnihilation;//p_bar annihilation
696e37fa	215
696e37fa	216	// neutron
cc364059	217	G4NeutronInelasticProcess fNeutronInelastic; //n inel process
	218	G4LENeutronInelastic* fLENeutronModel; //n LE inel model
	219	G4HENeutronInelastic* fHENeutronModel; //n HE inel model
	220	G4HadronFissionProcess fNeutronFission; //n fission
	221	G4LFission* fNeutronFissionModel; //n fission model
	222	G4HadronCaptureProcess fNeutronCapture; //n capture
	223	G4LCapture* fNeutronCaptureModel; //n capture model
696e37fa	224
696e37fa	225	// anti-neutron
cc364059	226	G4AntiNeutronInelasticProcess fAntiNeutronInelastic;//n_bar inel process
	227	G4LEAntiNeutronInelastic* fLEAntiNeutronModel; //n_bar LE inel model
	228	G4HEAntiNeutronInelastic* fHEAntiNeutronModel; //n_bar HE inel model
	229	G4AntiNeutronAnnihilationAtRest fAntiNeutronAnnihilation;//n_bar ionisation
696e37fa	230
696e37fa	231	// Lambda
cc364059	232	G4LambdaInelasticProcess fLambdaInelastic; //lambda inel process
	233	G4LELambdaInelastic* fLELambdaModel; //lambda LE inel model
	234	G4HELambdaInelastic* fHELambdaModel; //lambda HE inel model
696e37fa	235
696e37fa	236	// AntiLambda
cc364059	237	G4AntiLambdaInelasticProcess fAntiLambdaInelastic; //lambda_bar inel process
	238	G4LEAntiLambdaInelastic* fLEAntiLambdaModel; //lambda_bar LE inel model
	239	G4HEAntiLambdaInelastic* fHEAntiLambdaModel; //lambda_bar HE inel model
696e37fa	240
696e37fa	241	// SigmaMinus
cc364059	242	G4SigmaMinusInelasticProcess fSigmaMinusInelastic; //sigma- inel process
	243	G4LESigmaMinusInelastic* fLESigmaMinusModel; //sigma- LE inel model
	244	G4HESigmaMinusInelastic* fHESigmaMinusModel; //sigma- HE inel model
	245	G4MultipleScattering fSigmaMinusMult; //sigma- msc
	246	G4hIonisation fSigmaMinusIonisation;//sigma- ionisation
696e37fa	247
696e37fa	248	// AntiSigmaMinus
cc364059	249	G4AntiSigmaMinusInelasticProcess fAntiSigmaMinusInelastic; //sigma-_bar inel process
	250	G4LEAntiSigmaMinusInelastic* fLEAntiSigmaMinusModel; //sigma-_bar LE inel model
	251	G4HEAntiSigmaMinusInelastic* fHEAntiSigmaMinusModel; //sigma-_bar HE inel model
	252	G4MultipleScattering fAntiSigmaMinusMult; //sigma-_bar msc
	253	G4hIonisation fAntiSigmaMinusIonisation;//sigma-_bar ionisation
696e37fa	254
696e37fa	255	// SigmaPlus
cc364059	256	G4SigmaPlusInelasticProcess fSigmaPlusInelastic; //sigma+ inel process
	257	G4LESigmaPlusInelastic* fLESigmaPlusModel; //sigma+ LE inel model
	258	G4HESigmaPlusInelastic* fHESigmaPlusModel; //sigma+ HE inel model
	259	G4MultipleScattering fSigmaPlusMult; //sigma+ msc
	260	G4hIonisation fSigmaPlusIonisation; //sigma+ ionisation
696e37fa	261
696e37fa	262	// AntiSigmaPlus
cc364059	263	G4AntiSigmaPlusInelasticProcess fAntiSigmaPlusInelastic; //sigma+_bar inel process
	264	G4LEAntiSigmaPlusInelastic* fLEAntiSigmaPlusModel; //sigma+_bar LE inel model
	265	G4HEAntiSigmaPlusInelastic* fHEAntiSigmaPlusModel; //sigma+_bar HE inel model
	266	G4MultipleScattering fAntiSigmaPlusMult; //sigma+_bar msc
	267	G4hIonisation fAntiSigmaPlusIonisation; //sigma+_bar ionisation
696e37fa	268
696e37fa	269	// XiZero
cc364059	270	G4XiZeroInelasticProcess fXiZeroInelastic; //xi0 inel process
	271	G4LEXiZeroInelastic* fLEXiZeroModel; //xi0 LE inel model
	272	G4HEXiZeroInelastic* fHEXiZeroModel; //xi0 HE inel model
696e37fa	273
696e37fa	274	// AntiXiZero
cc364059	275	G4AntiXiZeroInelasticProcess fAntiXiZeroInelastic;//xi0_bar inel process
	276	G4LEAntiXiZeroInelastic* fLEAntiXiZeroModel; //xi0_bar LE inel model
	277	G4HEAntiXiZeroInelastic* fHEAntiXiZeroModel; //xi0_bar HE inel model
696e37fa	278
696e37fa	279	// XiMinus
cc364059	280	G4XiMinusInelasticProcess fXiMinusInelastic; //xi- inel process
	281	G4LEXiMinusInelastic* fLEXiMinusModel; //xi- LE inel model
	282	G4HEXiMinusInelastic* fHEXiMinusModel; //xi- HE inel model
	283	G4MultipleScattering fXiMinusMult; //xi- msc
	284	G4hIonisation fXiMinusIonisation; //xi- ionisation
696e37fa	285
696e37fa	286	// AntiXiMinus
cc364059	287	G4AntiXiMinusInelasticProcess fAntiXiMinusInelastic; //xi-_bar inel process
	288	G4LEAntiXiMinusInelastic* fLEAntiXiMinusModel; //xi-_bar LE inel model
	289	G4HEAntiXiMinusInelastic* fHEAntiXiMinusModel; //xi-_bar HE inel model
	290	G4MultipleScattering fAntiXiMinusMult; //xi-_bar msc
	291	G4hIonisation fAntiXiMinusIonisation;//xi-_bar ionisation
696e37fa	292
696e37fa	293	// OmegaMinus
cc364059	294	G4OmegaMinusInelasticProcess fOmegaMinusInelastic; //omega- inel process
	295	G4LEOmegaMinusInelastic* fLEOmegaMinusModel; //omega- LE inel model
	296	G4HEOmegaMinusInelastic* fHEOmegaMinusModel; //omega- HE inel model
	297	G4MultipleScattering fOmegaMinusMult; //omega- msc
	298	G4hIonisation fOmegaMinusIonisation; //omega- ionisation
696e37fa	299
696e37fa	300	// AntiOmegaMinus
cc364059	301	G4AntiOmegaMinusInelasticProcess fAntiOmegaMinusInelastic; //omega-_bar inel process
	302	G4LEAntiOmegaMinusInelastic* fLEAntiOmegaMinusModel; //omega-_bar LE inel model
	303	G4HEAntiOmegaMinusInelastic* fHEAntiOmegaMinusModel; //omega-_bar HE inel model
	304	G4MultipleScattering fAntiOmegaMinusMult; //omega-_bar msc
	305	G4hIonisation fAntiOmegaMinusIonisation;//omega-_bar ionisation
696e37fa	306
696e37fa	307	// Other
cc364059	308	ProcessVector fOtherProcesses; //other process
696e37fa	309
	310
	311	private:
	312	// methods
11f7ee82	313
	314	// EM processes
	315	void ConstructEMProcessForPionPlus();
	316	void ConstructEMProcessForPionMinus();
	317	void ConstructEMProcessForKaonPlus();
	318	void ConstructEMProcessForKaonMinus();
	319	void ConstructEMProcessForProton();
	320	void ConstructEMProcessForAntiProton();
	321	void ConstructEMProcessForSigmaMinus();
	322	void ConstructEMProcessForAntiSigmaMinus();
	323	void ConstructEMProcessForSigmaPlus();
	324	void ConstructEMProcessForAntiSigmaPlus();
	325	void ConstructEMProcessForXiMinus();
	326	void ConstructEMProcessForAntiXiMinus();
	327	void ConstructEMProcessForOmegaMinus();
	328	void ConstructEMProcessForAntiOmegaMinus();
	329	void ConstructEMProcessForOther();
	330
	331	// hadron processes
	332	void ConstructHadProcessForPionPlus();
	333	void ConstructHadProcessForPionMinus();
	334	void ConstructHadProcessForKaonPlus();
	335	void ConstructHadProcessForKaonMinus();
	336	void ConstructHadProcessForKaonZeroLong();
	337	void ConstructHadProcessForKaonZeroShort();
	338	void ConstructHadProcessForProton();
	339	void ConstructHadProcessForAntiProton();
	340	void ConstructHadProcessForNeutron();
	341	void ConstructHadProcessForAntiNeutron();
	342	void ConstructHadProcessForLambda();
	343	void ConstructHadProcessForAntiLambda();
	344	void ConstructHadProcessForSigmaMinus();
	345	void ConstructHadProcessForAntiSigmaMinus();
	346	void ConstructHadProcessForSigmaPlus();
	347	void ConstructHadProcessForAntiSigmaPlus();
	348	void ConstructHadProcessForXiMinus();
	349	void ConstructHadProcessForAntiXiMinus();
	350	void ConstructHadProcessForXiZero();
	351	void ConstructHadProcessForAntiXiZero();
	352	void ConstructHadProcessForOmegaMinus();
	353	void ConstructHadProcessForAntiOmegaMinus();
	354	void ConstructHadProcessForOther();
	355
	356	// data members
	357	G4bool fSetEM; //if true - EM processes are constructed
	358	G4bool fSetHadron;//if true - hadron processes are constructed
b94a914b	359	};
	360
	361	#endif //TG4_PHYSICS_CONSTRUCTOR_HADRON_H
	362