| Geant4 Cross Reference |
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 //
27 // -------------------------------------------------------------------
28 //
29 // GEANT4 Class header file
30 //
31 //
32 // File name: G4ComponentSAIDTotalXS
33 //
34 // Authors: G.Folger, V.Ivanchenko, D.Wright
35 //
36 // Modifications:
37 //
38
39 //
40 // Class Description
41 // Total, elastic and inelastic hadron/nucleon cross sections
42 // from SAID database, G4SAIDXSDATA environment variable
43 // should be defined
44 // Class Description - End
45
46 #ifndef G4ComponentSAIDTotalXS_h
47 #define G4ComponentSAIDTotalXS_h 1
48
49 #include "G4VComponentCrossSection.hh"
50 #include "G4ParticleDefinition.hh"
51 #include "globals.hh"
52 #include "G4Threading.hh"
53
54 enum G4SAIDCrossSectionType
55 {
56 saidUnknown = 0,
57 saidPP = 1,
58 saidNP = 2,
59 saidPIPP = 3,
60 saidPINP = 4,
61 saidPINP_PI0N = 5,
62 saidPINP_ETAN = 6,
63 saidGP_PI0P = 7,
64 saidGP_PIPN = 8,
65 saidGN_PINP = 9,
66 saidGN_PI0N = 10,
67 saidGP_ETAP = 11,
68 saidGP_ETAPP = 12,
69 numberOfSaidXS = 13
70 };
71
72 class G4PhysicsVector;
73
74 class G4ComponentSAIDTotalXS : public G4VComponentCrossSection
75 {
76 public: //with description
77
78 G4ComponentSAIDTotalXS();
79
80 virtual ~G4ComponentSAIDTotalXS();
81
82 virtual
83 G4double GetTotalElementCrossSection(const G4ParticleDefinition*,
84 G4double kinEnergy,
85 G4int /*Z*/, G4double /*N*/);
86
87 virtual
88 G4double GetTotalIsotopeCrossSection(const G4ParticleDefinition*,
89 G4double kinEnergy,
90 G4int /*Z*/, G4int /*N*/);
91
92 virtual
93 G4double GetInelasticElementCrossSection(const G4ParticleDefinition*,
94 G4double kinEnergy,
95 G4int /*Z*/, G4double /*N*/);
96
97 virtual
98 G4double GetInelasticIsotopeCrossSection(const G4ParticleDefinition*,
99 G4double kinEnergy,
100 G4int /*Z*/, G4int /*N*/);
101
102 virtual
103 G4double GetElasticElementCrossSection(const G4ParticleDefinition*,
104 G4double kinEnergy,
105 G4int /*Z*/, G4double /*N*/);
106
107 virtual
108 G4double GetElasticIsotopeCrossSection(const G4ParticleDefinition*,
109 G4double kinEnergy,
110 G4int /*Z*/, G4int /*N*/);
111
112 G4double GetChargeExchangeCrossSection(const G4ParticleDefinition* prim,
113 const G4ParticleDefinition* sec,
114 G4double kinEnergy,
115 G4int /*Z*/, G4int /*N*/);
116
117 virtual void Description(std::ostream&) const final;
118
119 private:
120
121 G4SAIDCrossSectionType GetType(const G4ParticleDefinition* prim,
122 const G4ParticleDefinition* sec,
123 G4int Z, G4int N);
124
125 void Initialise(G4SAIDCrossSectionType tp);
126
127 void ReadData(G4int index, G4PhysicsVector*,
128 const G4String&, const G4String&);
129
130 void PrintWarning(const G4ParticleDefinition* prim,
131 const G4ParticleDefinition* sec,
132 G4int /*Z*/, G4int /*N*/,
133 const G4String&, const G4String&);
134
135 G4ComponentSAIDTotalXS & operator=(const G4ComponentSAIDTotalXS &right);
136 G4ComponentSAIDTotalXS(const G4ComponentSAIDTotalXS&);
137
138 static const G4String fnames[numberOfSaidXS];
139 G4PhysicsVector* elastdata[numberOfSaidXS];
140 G4PhysicsVector* inelastdata[numberOfSaidXS];
141
142 #ifdef G4MULTITHREADED
143 static G4Mutex saidXSMutex;
144 #endif
145
146 };
147
148 #endif
149