source: trunk/source/processes/hadronic/cross_sections/include/G4GlauberGribovCrossSection.hh@ 1036

Last change on this file since 1036 was 962, checked in by garnier, 17 years ago

update processes
File size: 8.8 KB
RevLine 
[819]1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26// Calculation of the total, elastic and inelastic cross-sections
27// based on parametrisations of (proton, pion, kaon, photon) nucleon
28// cross-sections and the hadron-nucleous cross-section model in
29// the framework of Glauber-Gribov approach
30//
31//
32//
33//
34//
35// 17.07.06 V. Grichine - first implementation
36// 22.01.07 V.Ivanchenko - add interface with Z and A
37// 05.03.07 V.Ivanchenko - add IfZAApplicable
38// 06.03.07 V.Ivanchenko - add GetElasticGlauberGribov and GetElasticGlauberGribov
39// for combined dataset
40//
41//
42
43#ifndef G4GlauberGribovCrossSection_h
44#define G4GlauberGribovCrossSection_h
45
46#include "globals.hh"
47#include "G4Proton.hh"
48#include "G4Nucleus.hh"
49
50#include "G4VCrossSectionDataSet.hh"
51
52class G4ParticleDefinition;
53
54class G4GlauberGribovCrossSection : public G4VCrossSectionDataSet
55{
56public:
57
58 G4GlauberGribovCrossSection ();
59 virtual ~G4GlauberGribovCrossSection ();
60
61 virtual
62 G4bool IsApplicable(const G4DynamicParticle* aDP, const G4Element*);
63
64 virtual
65 G4bool IsZAApplicable(const G4DynamicParticle* aDP, G4double Z, G4double A);
66
67 virtual
68 G4double GetCrossSection(const G4DynamicParticle*,
69 const G4Element*,
70 G4double aTemperature = 0.0);
71
72 virtual
73 G4double GetIsoZACrossSection(const G4DynamicParticle*,
74 G4double Z, G4double A,
75 G4double aTemperature = 0.0);
76
77 virtual
78 void BuildPhysicsTable(const G4ParticleDefinition&)
79 {}
80
81 virtual
82 void DumpPhysicsTable(const G4ParticleDefinition&)
83 {G4cout << "G4GlauberGribovCrossSection: uses Glauber-Gribov formula"<<G4endl;}
84
85 G4double GetRatioSD(const G4DynamicParticle*, G4double At, G4double Zt);
86 G4double GetRatioQE(const G4DynamicParticle*, G4double At, G4double Zt);
87
88 G4double GetHadronNucleonXsc(const G4DynamicParticle*, const G4Element*);
89 G4double GetHadronNucleonXsc(const G4DynamicParticle*, G4double At, G4double Zt);
90
91 G4double GetHadronNucleonXscPDG(const G4DynamicParticle*, const G4Element*);
92 G4double GetHadronNucleonXscPDG(const G4DynamicParticle*, G4double At, G4double Zt);
93
94 G4double GetHadronNucleonXscNS(const G4DynamicParticle*, const G4Element*);
95 G4double GetHadronNucleonXscNS(const G4DynamicParticle*,G4double At, G4double Zt);
96
97 G4double GetHNinelasticXsc(const G4DynamicParticle*, const G4Element*);
98 G4double GetHNinelasticXsc(const G4DynamicParticle*, G4double At, G4double Zt);
99 G4double GetHNinelasticXscVU(const G4DynamicParticle*, G4double At, G4double Zt);
100 G4double GetHadronNucleonXscMK(const G4DynamicParticle* aParticle,
101 const G4ParticleDefinition* nucleon );
102
103 G4double CalculateEcmValue ( const G4double , const G4double , const G4double );
104
105 G4double CalcMandelstamS( const G4double , const G4double , const G4double );
106
107 G4double GetElasticGlauberGribov(const G4DynamicParticle*,G4double Z, G4double A);
108 G4double GetInelasticGlauberGribov(const G4DynamicParticle*,G4double Z, G4double A);
109
110 G4double GetTotalGlauberGribovXsc() { return fTotalXsc; };
111 G4double GetElasticGlauberGribovXsc() { return fElasticXsc; };
112 G4double GetInelasticGlauberGribovXsc(){ return fInelasticXsc; };
113 G4double GetProductionGlauberGribovXsc(){ return fProductionXsc; };
114 G4double GetDiffractionGlauberGribovXsc(){ return fDiffractionXsc; };
115 G4double GetRadiusConst() { return fRadiusConst; };
116
117 G4double GetNucleusRadius(const G4DynamicParticle*, const G4Element*);
118 G4double GetNucleusRadius(G4double At);
119
[962]120 inline G4double GetParticleBarCorTot( const G4ParticleDefinition* theParticle, G4double Z );
121 inline G4double GetParticleBarCorIn( const G4ParticleDefinition* theParticle, G4double Z );
122
123 inline void SetEnergyLowerLimit(G4double E ){fLowerLimit=E;};
124
[819]125private:
126
127 const G4double fUpperLimit;
[962]128 G4double fLowerLimit;
[819]129 const G4double fRadiusConst;
130
[962]131 static const G4double fNeutronBarCorrectionTot[93];
132 static const G4double fNeutronBarCorrectionIn[93];
133
134 static const G4double fProtonBarCorrectionTot[93];
135 static const G4double fProtonBarCorrectionIn[93];
136
137 static const G4double fPionPlusBarCorrectionTot[93];
138 static const G4double fPionPlusBarCorrectionIn[93];
139
140 static const G4double fPionMinusBarCorrectionTot[93];
141 static const G4double fPionMinusBarCorrectionIn[93];
142
[819]143 G4double fTotalXsc, fElasticXsc, fInelasticXsc, fProductionXsc, fDiffractionXsc;
144 G4double fHadronNucleonXsc;
145
146 G4ParticleDefinition* theGamma;
147 G4ParticleDefinition* theProton;
148 G4ParticleDefinition* theNeutron;
149 G4ParticleDefinition* theAProton;
150 G4ParticleDefinition* theANeutron;
151 G4ParticleDefinition* thePiPlus;
152 G4ParticleDefinition* thePiMinus;
153 G4ParticleDefinition* thePiZero;
154 G4ParticleDefinition* theKPlus;
155 G4ParticleDefinition* theKMinus;
156 G4ParticleDefinition* theK0S;
157 G4ParticleDefinition* theK0L;
158 G4ParticleDefinition* theL;
159 G4ParticleDefinition* theAntiL;
160 G4ParticleDefinition* theSPlus;
161 G4ParticleDefinition* theASPlus;
162 G4ParticleDefinition* theSMinus;
163 G4ParticleDefinition* theASMinus;
164 G4ParticleDefinition* theS0;
165 G4ParticleDefinition* theAS0;
166 G4ParticleDefinition* theXiMinus;
167 G4ParticleDefinition* theXi0;
168 G4ParticleDefinition* theAXiMinus;
169 G4ParticleDefinition* theAXi0;
170 G4ParticleDefinition* theOmega;
171 G4ParticleDefinition* theAOmega;
172 G4ParticleDefinition* theD;
173 G4ParticleDefinition* theT;
174 G4ParticleDefinition* theA;
175 G4ParticleDefinition* theHe3;
176
177};
178
[962]179////////////////////////////////////////////////////////////////
180//
181// Inlines
182
[819]183inline
184G4double G4GlauberGribovCrossSection::GetElasticGlauberGribov(
185 const G4DynamicParticle* dp, G4double Z, G4double A)
186{
187 GetIsoZACrossSection(dp, Z, A);
188 return fElasticXsc;
189}
190
[962]191/////////////////////////////////////////////////////////////////
192
[819]193inline
194G4double G4GlauberGribovCrossSection::GetInelasticGlauberGribov(
195 const G4DynamicParticle* dp, G4double Z, G4double A)
196{
197 GetIsoZACrossSection(dp, Z, A);
198 return fInelasticXsc;
199}
200
[962]201/////////////////////////////////////////////////////////////////////
202//
203// return correction at Tkin = 90*GeV GG -> Barashenkov tot xsc, when it
204// is available, else return 1.0
205
206
207inline G4double G4GlauberGribovCrossSection::GetParticleBarCorTot(
208 const G4ParticleDefinition* theParticle, G4double Z )
209{
210 G4int iZ = G4int(Z);
211
212 if( iZ >= 2 && iZ <= 92)
213 {
214 if( theParticle == theProton ) return fProtonBarCorrectionTot[iZ];
215 else if( theParticle == theNeutron) return fNeutronBarCorrectionTot[iZ];
216 else if( theParticle == thePiPlus ) return fPionPlusBarCorrectionTot[iZ];
217 else if( theParticle == thePiMinus) return fPionMinusBarCorrectionTot[iZ];
218 else return 1.0;
219 }
220 else return 1.0;
221}
222
223/////////////////////////////////////////////////////////////////////
224//
225// return correction at Tkin = 90*GeV GG -> Barashenkov in xsc, when it
226// is available, else return 1.0
227
228
229inline G4double G4GlauberGribovCrossSection::GetParticleBarCorIn(
230 const G4ParticleDefinition* theParticle, G4double Z )
231{
232 G4int iZ = G4int(Z);
233
234 if( iZ >= 2 && iZ <= 92)
235 {
236 if( theParticle == theProton ) return fProtonBarCorrectionIn[iZ];
237 else if( theParticle == theNeutron) return fNeutronBarCorrectionIn[iZ];
238 else if( theParticle == thePiPlus ) return fPionPlusBarCorrectionIn[iZ];
239 else if( theParticle == thePiMinus) return fPionMinusBarCorrectionIn[iZ];
240 else return 1.0;
241 }
242 else return 1.0;
243}
244
[819]245#endif
Note: See TracBrowser for help on using the repository browser.