Loading...
Searching...
No Matches
G4PiData Class Reference
#include <G4PiData.hh>
Inheritance diagram for G4PiData:
Public Member Functions | |
| G4PiData (const G4double *aTotal, const G4double *aInelastic, const G4double *anEnergy, G4int nPoints) | |
| ~G4PiData ()=default | |
| G4double | ReactionXSection (G4double kineticEnergy) |
| G4double | ElasticXSection (G4double kineticEnergy) |
| G4double | TotalXSection (G4double kineticEnergy) |
| G4PiData & | operator= (const G4PiData &right)=delete |
| G4PiData (const G4PiData &)=delete | |
Detailed Description
Definition at line 44 of file G4PiData.hh.
Constructor & Destructor Documentation
◆ G4PiData() [1/2]
|
explicit |
Definition at line 35 of file G4PiData.cc.
37{
39 {
40 std::pair<G4double, G4double> x;
41 x.first=aT[i]*millibarn;
42 x.second=aIn[i]*millibarn;
43 std::pair<G4double, std::pair<G4double, G4double > > aP;
44 aP.first=anE[i]*GeV;
45 aP.second=x;
46 push_back(aP);
47 if (i == 0) {
48 fMinE = aP.first;
49 fTot0 = x.first;
50 fInel0 = x.second;
51 } else if (i+1 == nP) {
52 fMaxE = aP.first;
53 fTot1 = x.first;
54 fInel1 = x.second;
55 }
56 }
57}
Referenced by G4PiData(), and operator=().
◆ ~G4PiData()
|
default |
◆ G4PiData() [2/2]
|
delete |
Member Function Documentation
◆ ElasticXSection()
Definition at line 86 of file G4PiData.cc.
87{
88 if (kineticEnergy <= fMinE) { return fTot0 - fInel0; }
89 if (kineticEnergy >= fMaxE) { return fTot1 - fInel1; }
90
91 G4PiData::iterator it=begin();
92 while(it!=end()&&kineticEnergy>(*it).first) {it++;} /* Loop checking, 08.01.2016, W. Pokorski */
93 if(it==end())
94 {
95 G4ExceptionDescription ed;
96 ed << "This cross section is applied for E(MeV)= " << kineticEnergy
99 }
100 if(it==begin()) it++;
101 G4double x1,x2,e1,e2;
102 e1=(*(it-1)).first;
103 x1=(*(it-1)).second.first - (*(it-1)).second.second;
104 e2=(*(it)).first;
105 x2=(*(it)).second.first - (*(it)).second.second;
106 return std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
107}
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition G4Exception.cc:59
◆ operator=()
◆ ReactionXSection()
Definition at line 61 of file G4PiData.cc.
62{
63 if (kineticEnergy <= fMinE) { return fInel0; }
64 if (kineticEnergy >= fMaxE) { return fInel1; }
65
66 G4PiData::iterator it=begin();
67 while(it!=end()&&kineticEnergy>(*it).first) {++it;} /* Loop checking, 08.01.2016, W. Pokorski */
68 if(it==end())
69 {
70 G4ExceptionDescription ed;
71 ed << "This cross section is applied for E(MeV)= " << kineticEnergy
74 }
75 if(it==begin()) it++;
76 G4double x1,x2,e1,e2;
77 e1=(*(it-1)).first;
78 x1=(*(it-1)).second.second;
79 e2=(*(it)).first;
80 x2=(*(it)).second.second;
81 return std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
82}
◆ TotalXSection()
Definition at line 111 of file G4PiData.cc.
112{
113 if (kineticEnergy <= fMinE) { return fTot0; }
114 if (kineticEnergy >= fMaxE) { return fTot1; }
115
116 G4PiData::iterator it=begin();
117 while(it!=end()&&kineticEnergy>(*it).first) {it++;} /* Loop checking, 08.01.2016, W. Pokorski */
118 if(it==end())
119 {
120 G4ExceptionDescription ed;
121 ed << "This cross section is applied for E(MeV)= " << kineticEnergy
124 }
125 if(it==begin()) it++;
126 G4double x1,x2,e1,e2;
127 e1=(*(it-1)).first;
128 x1=(*(it-1)).second.first;
129 e2=(*(it)).first;
130 x2=(*(it)).second.first;
131 return std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
132}
The documentation for this class was generated from the following files:
