Geant4_10
G4NucleonNuclearCrossSection.hh
Go to the documentation of this file.
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 // author: Vladimir.Grichine@cern.ch
27 //
28 // Implements data from: Barashenkov V.S., Nucleon-Nucleus Cross Section,
29 // Preprint JINR P2-89-770, p. 12, Dubna 1989 (scanned version from KEK)
30 // Based on G. Folger version of G4PiNuclearCrossSection class
31 //
32 // Modified:
33 // 05.03.07 V.Ivanchenko - add IfZAApplicable, remove "debug"
34 // 06.03.07 V.Ivanchenko - add GetElasticCrossSection for combined dataset
35 //
36 
37 #ifndef G4NucleonNuclearCrossSection_h
38 #define G4NucleonNuclearCrossSection_h
39 
41 #include "G4ParticleDefinition.hh"
42 
43 #include "globals.hh"
44 #include "G4PiData.hh"
45 #include "G4HadTmpUtil.hh"
46 
48 {
49 public:
50 
53 
54  static const char* Default_Name() {return "G4NucleonNuclearCrossSection";}
55 
56  virtual G4bool IsElementApplicable(const G4DynamicParticle* aParticle,
57  G4int Z,
58  const G4Material* mat = 0);
59 
60  virtual G4double GetElementCrossSection(const G4DynamicParticle* aParticle,
61  G4int Z,
62  const G4Material* mat = 0);
63 
64  virtual void CrossSectionDescription(std::ostream&) const;
65 
66  inline G4double GetElasticCrossSection(const G4DynamicParticle* aParticle,
67  G4int Z);
68 
69  inline G4double GetTotalXsc() { return fTotalXsc; };
70  inline G4double GetElasticXsc(){ return fElasticXsc; };
71 
72 private:
73 
74  G4double Interpolate(G4int Z1, G4int Z2, G4int Z, G4double x1, G4double x2);
75 
76 // add Hydrogen from PDG group.
77 
78 static const G4double e1[44];
79 
80 static const G4double he_m_t[44];
81 static const G4double he_m_in[44];
82 static const G4double he_p_in[44];
83 
84 static const G4double be_m_t[44];
85 static const G4double be_m_in[44];
86 static const G4double be_p_in[44];
87 
88 static const G4double c_m_t[44];
89 static const G4double c_m_in[44];
90 static const G4double c_p_in[44];
91 
92 
93 static const G4double e2[44];
94 
95 static const G4double n_m_t[44];
96 static const G4double n_m_in[44];
97 static const G4double n_p_in[44];
98 
99 static const G4double o_m_t[44];
100 static const G4double o_m_in[44];
101 static const G4double o_p_in[44];
102 
103 static const G4double na_m_t[44];
104 static const G4double na_m_in[44];
105 static const G4double na_p_in[44];
106 
107 
108 static const G4double e3[45];
109 
110 // static const G4double e3_1[31];
111 
112 static const G4double al_m_t[45];
113 static const G4double al_m_in[45];
114 static const G4double al_p_in[45];
115 
116 static const G4double si_m_t[45];
117 static const G4double si_m_in[45];
118 static const G4double si_p_in[45];
119 
120 static const G4double ca_m_t[45];
121 static const G4double ca_m_in[45];
122 static const G4double ca_p_in[45];
123 
124 
125 static const G4double e4[47];
126 
127 static const G4double fe_m_t[47];
128 static const G4double fe_m_in[47];
129 static const G4double fe_p_in[47];
130 
131 static const G4double cu_m_t[47];
132 static const G4double cu_m_in[47];
133 static const G4double cu_p_in[47];
134 
135 static const G4double mo_m_t[47];
136 static const G4double mo_m_in[47];
137 static const G4double mo_p_in[47];
138 
139 
140 static const G4double e5[48];
141 
142 static const G4double cd_m_t[48];
143 static const G4double cd_m_in[48];
144 static const G4double cd_p_in[48];
145 
146 static const G4double sn_m_t[48];
147 static const G4double sn_m_in[48];
148 static const G4double sn_p_in[48];
149 
150 static const G4double w_m_t[48];
151 static const G4double w_m_in[48];
152 static const G4double w_p_in[48];
153 
154 static const G4double e6[46];
155 
156 // static const G4double e7[46];
157 
158 static const G4double pb_m_t[46];
159 static const G4double pb_m_in[46];
160 static const G4double pb_p_in[46];
161 
162 static const G4double u_m_t[46];
163 static const G4double u_m_in[46];
164 static const G4double u_p_in[46];
165 
166 // vectors for treatment
167 
168 std::vector< G4int > theZ;
169 std::vector< G4PiData* > thePipData;
170 std::vector< G4PiData* > thePimData;
171 
172  // cross sections
173 
174  G4double fTotalXsc;
175  G4double fElasticXsc;
176 
177  // particles
178  const G4ParticleDefinition* theProton;
179  const G4ParticleDefinition* theNeutron;
180 
181 };
182 
183 inline
185  const G4DynamicParticle* dp, G4int Z)
186 {
187  GetElementCrossSection(dp, Z);
188  return fElasticXsc;
189 }
190 
191 #endif
Double_t Z2
Definition: plot.C:266
Double_t x2[nxs]
Definition: Style.C:19
virtual G4bool IsElementApplicable(const G4DynamicParticle *aParticle, G4int Z, const G4Material *mat=0)
virtual G4double GetElementCrossSection(const G4DynamicParticle *aParticle, G4int Z, const G4Material *mat=0)
virtual void CrossSectionDescription(std::ostream &) const
int G4int
Definition: G4Types.hh:78
Float_t mat
Definition: plot.C:40
G4double GetElasticCrossSection(const G4DynamicParticle *aParticle, G4int Z)
Float_t Z
Definition: plot.C:39
bool G4bool
Definition: G4Types.hh:79
Double_t x1[nxs]
Definition: Style.C:18
double G4double
Definition: G4Types.hh:76
Double_t Z1
Definition: plot.C:266