source: trunk/source/processes/hadronic/cross_sections/include/G4NucleonNuclearCrossSection.hh @ 1340

Last change on this file since 1340 was 1340, checked in by garnier, 14 years ago

update ti head
File size: 5.8 KB
Line 
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
40#include "G4VCrossSectionDataSet.hh"
41#include "G4ParticleDefinition.hh"
42
43#include "globals.hh"
44#include "G4PiData.hh"
45#include "G4HadTmpUtil.hh"
46
47class G4NucleonNuclearCrossSection : public G4VCrossSectionDataSet
48{
49public:
50 
51  G4NucleonNuclearCrossSection();
52  virtual ~G4NucleonNuclearCrossSection();
53
54  virtual G4bool IsApplicable(const G4DynamicParticle* aParticle,
55                              const G4Element* anElement);
56
57  virtual G4bool IsIsoApplicable(const G4DynamicParticle* aParticle,
58                                 G4int Z, G4int A);
59
60  G4double GetCrossSection(const G4DynamicParticle* aParticle, 
61                           const G4Element* anElement,
62                           G4double T=0.);
63
64  G4double GetZandACrossSection(const G4DynamicParticle* aParticle, 
65                                G4int ZZ, G4int AA, G4double T=0. );
66
67  G4double GetElasticCrossSection(const G4DynamicParticle* aParticle, 
68                                  G4int ZZ, G4int AA);
69
70  G4double GetTotalXsc()  { return fTotalXsc;   };
71  G4double GetElasticXsc(){ return fElasticXsc; };
72
73
74  void BuildPhysicsTable(const G4ParticleDefinition&) {}
75  void DumpPhysicsTable(const G4ParticleDefinition&) {}
76 
77private:
78
79  G4double Interpolate(G4int Z1, G4int Z2, G4int Z, G4double x1, G4double x2);
80
81// add Hydrogen from PDG group.
82
83static const G4double e1[44];
84
85static const G4double he_m_t[44];
86static const G4double he_m_in[44];
87static const G4double he_p_in[44];
88
89static const G4double be_m_t[44];
90static const G4double be_m_in[44];
91static const G4double be_p_in[44];
92
93static const G4double c_m_t[44];
94static const G4double c_m_in[44];
95static const G4double c_p_in[44];
96
97
98static const G4double e2[44];
99
100static const G4double n_m_t[44];
101static const G4double n_m_in[44];
102static const G4double n_p_in[44];
103
104static const G4double o_m_t[44];
105static const G4double o_m_in[44];
106static const G4double o_p_in[44];
107
108static const G4double na_m_t[44];
109static const G4double na_m_in[44];
110static const G4double na_p_in[44];
111
112
113static const G4double e3[45];
114
115// static const G4double e3_1[31];
116
117static const G4double al_m_t[45];
118static const G4double al_m_in[45];
119static const G4double al_p_in[45];
120
121static const G4double si_m_t[45];
122static const G4double si_m_in[45];
123static const G4double si_p_in[45];
124
125static const G4double ca_m_t[45];
126static const G4double ca_m_in[45];
127static const G4double ca_p_in[45];
128
129
130static const G4double e4[47];
131
132static const G4double fe_m_t[47];
133static const G4double fe_m_in[47];
134static const G4double fe_p_in[47];
135
136static const G4double cu_m_t[47];
137static const G4double cu_m_in[47];
138static const G4double cu_p_in[47];
139
140static const G4double mo_m_t[47];
141static const G4double mo_m_in[47];
142static const G4double mo_p_in[47];
143
144
145static const G4double e5[48];
146
147static const G4double cd_m_t[48];
148static const G4double cd_m_in[48];
149static const G4double cd_p_in[48];
150
151static const G4double sn_m_t[48];
152static const G4double sn_m_in[48];
153static const G4double sn_p_in[48];
154
155static const G4double w_m_t[48];
156static const G4double w_m_in[48];
157static const G4double w_p_in[48];
158
159
160
161
162static const G4double e6[46];
163
164
165// static const G4double e7[46];
166
167static const G4double pb_m_t[46];
168static const G4double pb_m_in[46];
169static const G4double pb_p_in[46];
170
171static const G4double u_m_t[46];
172static const G4double u_m_in[46];
173static const G4double u_p_in[46];
174
175// vectors for treatment
176
177std::vector< G4int >     theZ;
178std::vector< G4PiData* > thePipData;
179std::vector< G4PiData* > thePimData;
180
181  // cross sections
182
183  G4double fTotalXsc;
184  G4double fElasticXsc;
185
186  // particles
187  const G4ParticleDefinition* theProton;
188  const G4ParticleDefinition* theNeutron;
189
190};
191
192inline
193G4double G4NucleonNuclearCrossSection::GetElasticCrossSection(
194         const G4DynamicParticle* dp, G4int ZZ, G4int AA)
195{
196  GetZandACrossSection(dp, ZZ, AA);
197  return fElasticXsc;
198}
199
200#endif
Note: See TracBrowser for help on using the repository browser.