00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049 #ifndef G4BGGPionInelasticXS_h
00050 #define G4BGGPionInelasticXS_h
00051
00052 #include "globals.hh"
00053 #include "G4VCrossSectionDataSet.hh"
00054 #include "G4ParticleDefinition.hh"
00055 #include "G4Element.hh"
00056 #include "G4HadTmpUtil.hh"
00057
00058
00059 class G4GlauberGribovCrossSection;
00060 class G4UPiNuclearCrossSection;
00061 class G4HadronNucleonXsc;
00062
00063 class G4ComponentSAIDTotalXS;
00064 class G4Material;
00065 class G4Element;
00066 class G4Isotope;
00067
00068 class G4BGGPionInelasticXS : public G4VCrossSectionDataSet
00069 {
00070 public:
00071
00072 G4BGGPionInelasticXS (const G4ParticleDefinition*);
00073
00074 virtual ~G4BGGPionInelasticXS();
00075
00076 virtual
00077 G4bool IsElementApplicable(const G4DynamicParticle*, G4int Z,
00078 const G4Material* mat = 0);
00079
00080 virtual
00081 G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A,
00082 const G4Element* elm = 0,
00083 const G4Material* mat = 0);
00084
00085 virtual
00086 G4double GetElementCrossSection(const G4DynamicParticle*, G4int Z,
00087 const G4Material* mat = 0);
00088
00089 virtual
00090 G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A,
00091 const G4Isotope* iso = 0,
00092 const G4Element* elm = 0,
00093 const G4Material* mat = 0);
00094
00095 virtual
00096 void BuildPhysicsTable(const G4ParticleDefinition&);
00097
00098 virtual void CrossSectionDescription(std::ostream&) const;
00099
00100 private:
00101
00102 G4double CoulombFactor(G4double kinEnergy, G4int Z);
00103
00104 G4BGGPionInelasticXS & operator=(const G4BGGPionInelasticXS &right);
00105 G4BGGPionInelasticXS(const G4BGGPionInelasticXS&);
00106
00107 G4double fGlauberEnergy;
00108 G4double fLowEnergy;
00109 G4double fSAIDHighEnergyLimit;
00110 G4double theGlauberFac[93];
00111 G4double theCoulombFac[93];
00112 G4int theA[93];
00113
00114 const G4ParticleDefinition* particle;
00115 const G4ParticleDefinition* theProton;
00116 G4GlauberGribovCrossSection* fGlauber;
00117 G4UPiNuclearCrossSection* fPion;
00118 G4HadronNucleonXsc* fHadron;
00119
00120 G4ComponentSAIDTotalXS* fSAID;
00121 G4bool isPiplus;
00122 G4bool isInitialized;
00123 };
00124
00125 #endif