Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4QHadronElasticDataSet Class Reference

#include <G4QHadronElasticDataSet.hh>

+ Inheritance diagram for G4QHadronElasticDataSet:

Public Member Functions

 G4QHadronElasticDataSet (const G4String &dataSetName="CHIPSElasticXS")
 
 ~G4QHadronElasticDataSet ()
 
G4bool IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
 
G4double GetIsoCrossSection (const G4DynamicParticle *P, G4int Z, G4int A, const G4Isotope *, const G4Element *, const G4Material *)
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
void DumpPhysicsTable (const G4ParticleDefinition &)
 
void Description () const
 
- Public Member Functions inherited from G4VCrossSectionDataSet
 G4VCrossSectionDataSet (const G4String &nam="")
 
virtual ~G4VCrossSectionDataSet ()
 
virtual G4bool IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
 
virtual G4bool IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
 
G4double GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
 
G4double ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
 
virtual G4double GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
 
virtual G4double GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
 
virtual G4Isotope * SelectIsotope (const G4Element *, G4double kinEnergy)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void DumpPhysicsTable (const G4ParticleDefinition &)
 
virtual void CrossSectionDescription (std::ostream &) const
 
void SetVerboseLevel (G4int value)
 
G4double GetMinKinEnergy () const
 
void SetMinKinEnergy (G4double value)
 
G4double GetMaxKinEnergy () const
 
void SetMaxKinEnergy (G4double value)
 
const G4String & GetName () const
 

Additional Inherited Members

- Protected Member Functions inherited from G4VCrossSectionDataSet
void SetName (const G4String &)
 
- Protected Attributes inherited from G4VCrossSectionDataSet
G4int verboseLevel
 

Detailed Description

Definition at line 78 of file G4QHadronElasticDataSet.hh.

Constructor & Destructor Documentation

◆ G4QHadronElasticDataSet()

G4QHadronElasticDataSet::G4QHadronElasticDataSet ( const G4String &  dataSetName = "CHIPSElasticXS")

Definition at line 49 of file G4QHadronElasticDataSet.cc.

50 : G4VCrossSectionDataSet(dataSetName)
51{
52 //Isotopes = G4QIsotope::Get(); // Pointer to the G4QIsotopes singleton
53 Description();
54}

◆ ~G4QHadronElasticDataSet()

G4QHadronElasticDataSet::~G4QHadronElasticDataSet ( )
inline

Definition at line 83 of file G4QHadronElasticDataSet.hh.

83{};

Member Function Documentation

◆ BuildPhysicsTable()

void G4QHadronElasticDataSet::BuildPhysicsTable ( const G4ParticleDefinition &  )
inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 88 of file G4QHadronElasticDataSet.hh.

88{};

◆ Description()

void G4QHadronElasticDataSet::Description ( ) const

Definition at line 56 of file G4QHadronElasticDataSet.cc.

57{
58 char* dirName = getenv("G4PhysListDocDir");
59 if (dirName) {
60 std::ofstream outFile;
61 G4String outFileName = GetName() + ".html";
62 G4String pathName = G4String(dirName) + "/" + outFileName;
63
64 outFile.open(pathName);
65 outFile << "<html>\n";
66 outFile << "<head>\n";
67
68 outFile << "<title>Description of CHIPSElasticXS</title>\n";
69 outFile << "</head>\n";
70 outFile << "<body>\n";
71
72 outFile << "CHIPSElasticXS provides hadron-nuclear elastic cross\n"
73 << "sections for all hadrons at all energies. These cross\n"
74 << "sections represent parameterizations developed by M. Kossov.\n";
75
76 outFile << "</body>\n";
77 outFile << "</html>\n";
78 outFile.close();
79 }
80}
G4VCrossSectionDataSet::GetName
const G4String & GetName() const
Definition: G4VCrossSectionDataSet.hh:211

Referenced by G4QHadronElasticDataSet().

◆ DumpPhysicsTable()

void G4QHadronElasticDataSet::DumpPhysicsTable ( const G4ParticleDefinition &  )
inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4QHadronElasticDataSet.hh.

89{};

◆ GetIsoCrossSection()

G4double G4QHadronElasticDataSet::GetIsoCrossSection ( const G4DynamicParticle *  P,
G4int  Z,
G4int  A,
const G4Isotope *  ,
const G4Element *  ,
const G4Material *   
)
virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 182 of file G4QHadronElasticDataSet.cc.

185{
186 G4ParticleDefinition* particle = Pt->GetDefinition();
187 G4double Momentum=Pt->GetTotalMomentum();
188 G4VQCrossSection* CSmanager=0;
189 G4VQCrossSection* CSmanager2=0;
190
191 G4int pPDG=0;
192 if(particle == G4Neutron::Neutron())
193 {
194 CSmanager=G4QNeutronElasticCrossSection::GetPointer();
195 pPDG=2112;
196 }
197 else if(particle == G4Proton::Proton())
198 {
199 CSmanager=G4QProtonElasticCrossSection::GetPointer();
200 pPDG=2212;
201 }
202 else if(particle == G4PionMinus::PionMinus())
203 {
204 CSmanager=G4QPionMinusElasticCrossSection::GetPointer();
205 pPDG=-211;
206 }
207 else if(particle == G4PionPlus::PionPlus())
208 {
209 CSmanager=G4QPionPlusElasticCrossSection::GetPointer();
210 pPDG=211;
211 }
212 else if(particle == G4KaonMinus::KaonMinus())
213 {
214 CSmanager=G4QKaonMinusElasticCrossSection::GetPointer();
215 pPDG=-321;
216 }
217 else if(particle == G4KaonPlus::KaonPlus())
218 {
219 CSmanager=G4QKaonPlusElasticCrossSection::GetPointer();
220 pPDG=321;
221 }
222 else if(particle == G4KaonZeroLong::KaonZeroLong() ||
223 particle == G4KaonZeroShort::KaonZeroShort() ||
224 particle == G4KaonZero::KaonZero() ||
225 particle == G4AntiKaonZero::AntiKaonZero() )
226 {
227 CSmanager=G4QKaonMinusElasticCrossSection::GetPointer();
228 CSmanager2=G4QKaonPlusElasticCrossSection::GetPointer();
229 if(G4UniformRand() > 0.5) pPDG= 321;
230 else pPDG=-321;
231 }
232 else if(particle == G4Lambda::Lambda())
233 {
234 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
235 pPDG=3122;
236 }
237 else if(particle == G4SigmaPlus::SigmaPlus())
238 {
239 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
240 pPDG=3222;
241 }
242 else if(particle == G4SigmaMinus::SigmaMinus())
243 {
244 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
245 pPDG=3112;
246 }
247 else if(particle == G4SigmaZero::SigmaZero())
248 {
249 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
250 pPDG=3212;
251 }
252 else if(particle == G4XiMinus::XiMinus())
253 {
254 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
255 pPDG=3312;
256 }
257 else if(particle == G4XiZero::XiZero())
258 {
259 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
260 pPDG=3322;
261 }
262 else if(particle == G4OmegaMinus::OmegaMinus())
263 {
264 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
265 pPDG=3334;
266 }
267 else if(particle == G4AntiNeutron::AntiNeutron())
268 {
269 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
270 pPDG=-2112;
271 }
272 else if(particle == G4AntiProton::AntiProton())
273 {
274 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
275 pPDG=-2212;
276 }
277 else if(particle == G4AntiLambda::AntiLambda())
278 {
279 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
280 pPDG=-3122;
281 }
282 else if(particle == G4AntiSigmaPlus::AntiSigmaPlus())
283 {
284 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
285 pPDG=-3222;
286 }
287 else if(particle == G4AntiSigmaMinus::AntiSigmaMinus())
288 {
289 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
290 pPDG=-3112;
291 }
292 else if(particle == G4AntiSigmaZero::AntiSigmaZero())
293 {
294 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
295 pPDG=-3212;
296 }
297 else if(particle == G4AntiXiMinus::AntiXiMinus())
298 {
299 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
300 pPDG=-3312;
301 }
302 else if(particle == G4AntiXiZero::AntiXiZero())
303 {
304 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
305 pPDG=-3322;
306 }
307 else if(particle == G4AntiOmegaMinus::AntiOmegaMinus())
308 {
309 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
310 pPDG=-3334;
311 }
312 else
313 {
314 G4cerr << "-ERROR-G4QHadronElasticDataSet::GetIsoZACrossSection: PDG="
315 << particle->GetPDGEncoding() << " isn't supported by CHIPS" << G4endl;
316 //throw G4HadronicException(__FILE__, __LINE__,
317 // "G4QHadronElasticDataSet::GetIsoZACrossSection: Particle isn't supported by CHIPS");
318 return 0;
319 }
320 G4int N=A-Z;
321 G4double CSI=CSmanager->GetCrossSection(true, Momentum, Z, N, pPDG); // CS(j,i) basic
322 if(CSmanager2) CSI = (CSI + CSmanager2->GetCrossSection(true, Momentum, Z, N, pPDG))/2;
323 return CSI;
324}
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cerr
G4DLLIMPORT std::ostream G4cerr
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
G4AntiKaonZero::AntiKaonZero
static G4AntiKaonZero * AntiKaonZero()
Definition: G4AntiKaonZero.cc:103
G4AntiLambda::AntiLambda
static G4AntiLambda * AntiLambda()
Definition: G4AntiLambda.cc:108
G4AntiNeutron::AntiNeutron
static G4AntiNeutron * AntiNeutron()
Definition: G4AntiNeutron.cc:103
G4AntiOmegaMinus::AntiOmegaMinus
static G4AntiOmegaMinus * AntiOmegaMinus()
Definition: G4AntiOmegaMinus.cc:110
G4AntiProton::AntiProton
static G4AntiProton * AntiProton()
Definition: G4AntiProton.cc:93
G4AntiSigmaMinus::AntiSigmaMinus
static G4AntiSigmaMinus * AntiSigmaMinus()
Definition: G4AntiSigmaMinus.cc:106
G4AntiSigmaPlus::AntiSigmaPlus
static G4AntiSigmaPlus * AntiSigmaPlus()
Definition: G4AntiSigmaPlus.cc:108
G4AntiSigmaZero::AntiSigmaZero
static G4AntiSigmaZero * AntiSigmaZero()
Definition: G4AntiSigmaZero.cc:99
G4AntiXiMinus::AntiXiMinus
static G4AntiXiMinus * AntiXiMinus()
Definition: G4AntiXiMinus.cc:105
G4AntiXiZero::AntiXiZero
static G4AntiXiZero * AntiXiZero()
Definition: G4AntiXiZero.cc:106
G4KaonMinus::KaonMinus
static G4KaonMinus * KaonMinus()
Definition: G4KaonMinus.cc:113
G4KaonPlus::KaonPlus
static G4KaonPlus * KaonPlus()
Definition: G4KaonPlus.cc:113
G4KaonZeroLong::KaonZeroLong
static G4KaonZeroLong * KaonZeroLong()
Definition: G4KaonZeroLong.cc:115
G4KaonZeroShort::KaonZeroShort
static G4KaonZeroShort * KaonZeroShort()
Definition: G4KaonZeroShort.cc:104
G4KaonZero::KaonZero
static G4KaonZero * KaonZero()
Definition: G4KaonZero.cc:104
G4Lambda::Lambda
static G4Lambda * Lambda()
Definition: G4Lambda.cc:108
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4OmegaMinus::OmegaMinus
static G4OmegaMinus * OmegaMinus()
Definition: G4OmegaMinus.cc:110
G4PionMinus::PionMinus
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:98
G4PionPlus::PionPlus
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:98
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
G4QProtonElasticCrossSection::GetPointer
static G4VQCrossSection * GetPointer()
Definition: G4QProtonElasticCrossSection.cc:143
G4SigmaMinus::SigmaMinus
static G4SigmaMinus * SigmaMinus()
Definition: G4SigmaMinus.cc:106
G4SigmaPlus::SigmaPlus
static G4SigmaPlus * SigmaPlus()
Definition: G4SigmaPlus.cc:108
G4SigmaZero::SigmaZero
static G4SigmaZero * SigmaZero()
Definition: G4SigmaZero.cc:99
G4VQCrossSection::GetCrossSection
virtual G4double GetCrossSection(G4bool, G4double, G4int, G4int, G4int pPDG=0)
Definition: G4VQCrossSection.hh:99
G4XiMinus::XiMinus
static G4XiMinus * XiMinus()
Definition: G4XiMinus.cc:106
G4XiZero::XiZero
static G4XiZero * XiZero()
Definition: G4XiZero.cc:106

◆ IsIsoApplicable()

G4bool G4QHadronElasticDataSet::IsIsoApplicable ( const G4DynamicParticle *  Pt,
G4int  Z,
G4int  A,
const G4Element *  elm = 0,
const G4Material *  mat = 0 
)
virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 82 of file G4QHadronElasticDataSet.cc.

84{
85 G4ParticleDefinition* particle = Pt->GetDefinition();
86 if (particle == G4Neutron::Neutron() ) return true; // @@ isotopes?
87 else if (particle == G4Proton::Proton() ) return true;
88 else if (particle == G4PionMinus::PionMinus() ) return true;
89 else if (particle == G4PionPlus::PionPlus() ) return true;
90 else if (particle == G4KaonPlus::KaonPlus() ) return true;
91 else if (particle == G4KaonMinus::KaonMinus() ) return true;
92 else if (particle == G4KaonZeroLong::KaonZeroLong() ) return true;
93 else if (particle == G4KaonZeroShort::KaonZeroShort() ) return true;
94 else if (particle == G4Lambda::Lambda() ) return true;
95 else if (particle == G4SigmaPlus::SigmaPlus() ) return true;
96 else if (particle == G4SigmaMinus::SigmaMinus() ) return true;
97 else if (particle == G4SigmaZero::SigmaZero() ) return true;
98 else if (particle == G4XiMinus::XiMinus() ) return true;
99 else if (particle == G4XiZero::XiZero() ) return true;
100 else if (particle == G4OmegaMinus::OmegaMinus() ) return true;
101 else if (particle == G4AntiNeutron::AntiNeutron() ) return true;
102 else if (particle == G4AntiProton::AntiProton() ) return true;
103 else if (particle == G4AntiLambda::AntiLambda() ) return true;
104 else if (particle == G4AntiSigmaPlus::AntiSigmaPlus() ) return true;
105 else if (particle == G4AntiSigmaMinus::AntiSigmaMinus() ) return true;
106 else if (particle == G4AntiSigmaZero::AntiSigmaZero() ) return true;
107 else if (particle == G4AntiXiMinus::AntiXiMinus() ) return true;
108 else if (particle == G4AntiXiZero::AntiXiZero() ) return true;
109 else if (particle == G4AntiOmegaMinus::AntiOmegaMinus() ) return true;
110 return false;
111}
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
The documentation for this class was generated from the following files: