| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
>> 23 // $Id: G4PiData.cc,v 1.6 2003/06/16 17:03:08 gunter Exp $
>> 24 // GEANT4 tag $Name: geant4-05-02 $
26 // 25 //
27 // ------------------------------------------- 26 // --------------------------------------------------------------------
28 // by J.P Wellisch, Sun Sep 15 2002. 27 // by J.P Wellisch, Sun Sep 15 2002.
29 28
30 #include "G4PiData.hh" 29 #include "G4PiData.hh"
31 #include "G4SystemOfUnits.hh" <<
32 30
33 ////////////////////////////////////////////// << 31 G4PiData::
34 << 32 G4PiData(const G4double * aT, const G4double * aIn, const G4double * anE, G4int nP)
35 G4PiData::G4PiData(const G4double * aT, const <<
36 const G4double * anE, G4int <<
37 { 33 {
38 for(G4int i = 0; i < nP; ++i ) << 34 G4int i=0;
>> 35 for(i=0; i<nP; i++)
39 { 36 {
40 std::pair<G4double, G4double> x; 37 std::pair<G4double, G4double> x;
41 x.first=aT[i]*millibarn; 38 x.first=aT[i]*millibarn;
42 x.second=aIn[i]*millibarn; 39 x.second=aIn[i]*millibarn;
43 std::pair<G4double, std::pair<G4double, G4 40 std::pair<G4double, std::pair<G4double, G4double > > aP;
44 aP.first=anE[i]*GeV; 41 aP.first=anE[i]*GeV;
45 aP.second=x; 42 aP.second=x;
46 push_back(aP); 43 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 } 44 }
57 } 45 }
58 46
59 ////////////////////////////////////////////// << 47 G4bool G4PiData::
60 << 48 AppliesTo(G4double kineticEnergy)
61 G4double G4PiData::ReactionXSection(G4double k <<
62 { 49 {
63 if (kineticEnergy <= fMinE) { return fInel0; << 50 G4bool result = true;
64 if (kineticEnergy >= fMaxE) { return fInel1; << 51 if(kineticEnergy>back().first) result = false;
>> 52 return result;
>> 53 }
65 54
>> 55 G4double G4PiData::
>> 56 ReactionXSection(G4double kineticEnergy)
>> 57 {
>> 58 G4double result = 0;
66 G4PiData::iterator it=begin(); 59 G4PiData::iterator it=begin();
67 while(it!=end()&&kineticEnergy>(*it).first) << 60 while(it!=end()&&kineticEnergy>(*it).first) {it++;}
68 if(it==end()) << 61 if(it==end()) G4Exception("G4PiData used outside validity range");
69 { <<
70 G4ExceptionDescription ed; <<
71 ed << "This cross section is applied for E <<
72 << " outside allowed energy interval" < <<
73 G4Exception("G4PiData::ReactionXSection", <<
74 } <<
75 if(it==begin()) it++; 62 if(it==begin()) it++;
76 G4double x1,x2,e1,e2; 63 G4double x1,x2,e1,e2;
77 e1=(*(it-1)).first; 64 e1=(*(it-1)).first;
78 x1=(*(it-1)).second.second; 65 x1=(*(it-1)).second.second;
79 e2=(*(it)).first; 66 e2=(*(it)).first;
80 x2=(*(it)).second.second; 67 x2=(*(it)).second.second;
81 return std::max(0., x1 + (kineticEnergy-e1)* << 68 result = std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
>> 69 return result;
82 } 70 }
83 71
84 ////////////////////////////////////////////// << 72 G4double G4PiData::
85 << 73 ElasticXSection(G4double kineticEnergy)
86 G4double G4PiData::ElasticXSection(G4double ki <<
87 { 74 {
88 if (kineticEnergy <= fMinE) { return fTot0 - << 75 G4double result = 0;
89 if (kineticEnergy >= fMaxE) { return fTot1 - <<
90 <<
91 G4PiData::iterator it=begin(); 76 G4PiData::iterator it=begin();
92 while(it!=end()&&kineticEnergy>(*it).first) << 77 while(it!=end()&&kineticEnergy>(*it).first) {it++;}
93 if(it==end()) << 78 if(it==end()) G4Exception("G4PiData used outside validity range");
94 { <<
95 G4ExceptionDescription ed; <<
96 ed << "This cross section is applied for E <<
97 << " outside allowed energy interval" < <<
98 G4Exception("G4PiData::ElasticXSection", " <<
99 } <<
100 if(it==begin()) it++; 79 if(it==begin()) it++;
101 G4double x1,x2,e1,e2; 80 G4double x1,x2,e1,e2;
102 e1=(*(it-1)).first; 81 e1=(*(it-1)).first;
103 x1=(*(it-1)).second.first - (*(it-1)).second 82 x1=(*(it-1)).second.first - (*(it-1)).second.second;
104 e2=(*(it)).first; 83 e2=(*(it)).first;
105 x2=(*(it)).second.first - (*(it)).second.sec 84 x2=(*(it)).second.first - (*(it)).second.second;
106 return std::max(0., x1 + (kineticEnergy-e1)* << 85 result = std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
107 } << 86 return result;
108 <<
109 ////////////////////////////////////////////// <<
110 <<
111 G4double G4PiData::TotalXSection(G4double kine <<
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) <<
118 if(it==end()) <<
119 { <<
120 G4ExceptionDescription ed; <<
121 ed << "This cross section is applied for E <<
122 << " outside allowed energy interval" < <<
123 G4Exception("G4PiData::TotalXSection", "ha <<
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)* <<
132 } 87 }
133 88