| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //
28 // Author: Maria Grazia Pia (Maria.Grazia.Pia@
29 //
30 // History:
31 // -----------
32 // 31 Jul 2001 MGP Created
33 //
34 // -------------------------------------------
35
36 #include "G4RDDopplerProfile.hh"
37 #include "G4DataVector.hh"
38 #include "G4SystemOfUnits.hh"
39 #include "G4RDVEMDataSet.hh"
40 #include "G4RDEMDataSet.hh"
41 #include "G4RDCompositeEMDataSet.hh"
42 #include "G4RDVDataSetAlgorithm.hh"
43 #include "G4RDLogLogInterpolation.hh"
44
45 #include <fstream>
46 #include <sstream>
47 #include "Randomize.hh"
48
49 // The following deprecated header is included
50 //#include "function.h"
51
52 // Constructor
53
54 G4RDDopplerProfile::G4RDDopplerProfile(G4int m
55 : zMin(minZ), zMax(maxZ)
56 {
57 nBiggs = 31;
58
59 LoadBiggsP("/doppler/p-biggs");
60
61 for (G4int Z=zMin; Z<zMax+1; Z++)
62 {
63 LoadProfile("/doppler/profile",Z);
64 }
65 }
66
67 // Destructor
68 G4RDDopplerProfile::~G4RDDopplerProfile()
69 {
70 std::map<G4int,G4RDVEMDataSet*,std::less<G4i
71 for (pos = profileMap.begin(); pos != profil
72 {
73 G4RDVEMDataSet* dataSet = (*pos).second;
74 delete dataSet;
75 dataSet = 0;
76 }
77 }
78
79
80 size_t G4RDDopplerProfile::NumberOfProfiles(G4
81 {
82 G4int n = 0;
83 if (Z>= zMin && Z <= zMax) n = nShells[Z-1];
84 return n;
85 }
86
87
88 const G4RDVEMDataSet* G4RDDopplerProfile::Prof
89 {
90 std::map<G4int,G4RDVEMDataSet*,std::less<G4i
91 if (Z < zMin || Z > zMax)
92 G4Exception("G4RDDopplerProfile::Profiles(
93 FatalException, "Z outside bou
94 pos = profileMap.find(Z);
95 G4RDVEMDataSet* dataSet = (*pos).second;
96 return dataSet;
97 }
98
99
100 const G4RDVEMDataSet* G4RDDopplerProfile::Prof
101 {
102 const G4RDVEMDataSet* profis = Profiles(Z);
103 const G4RDVEMDataSet* profi = profis->GetCom
104 return profi;
105 }
106
107
108 void G4RDDopplerProfile::PrintData() const
109 {
110 for (G4int Z=zMin; Z<zMax; Z++)
111 {
112 const G4RDVEMDataSet* profis = Profiles(
113 profis->PrintData();
114 }
115 }
116
117
118 void G4RDDopplerProfile::LoadBiggsP(const G4St
119 {
120 std::ostringstream ost;
121 ost << fileName << ".dat";
122 G4String name(ost.str());
123
124 const char* path = G4FindDataDir("G4LEDATA")
125 if (!path)
126 {
127 G4String excep("G4LEDATA environment var
128 G4Exception("G4RDDopplerProfile::LoadBig
129 "InvalidSetup", FatalExcepti
130 }
131
132 G4String pathString(path);
133 G4String dirFile = pathString + name;
134 std::ifstream file(dirFile);
135 std::filebuf* lsdp = file.rdbuf();
136
137 if (! (lsdp->is_open()) )
138 {
139 G4String s1("Data file: ");
140 G4String s2(" not found");
141 G4String excep = s1 + dirFile + s2;
142 G4Exception("G4RDDopplerProfile::LoadBig
143 "DataNotFound", FatalExcepti
144 }
145
146 G4double p;
147 while(!file.eof())
148 {
149 file >> p;
150 biggsP.push_back(p);
151 }
152
153 // Make sure that the number of data loaded
154 if (biggsP.size() != nBiggs)
155 G4Exception("G4RDDopplerProfile::LoadBiggs
156 FatalException, "Number of mom
157 }
158
159
160 void G4RDDopplerProfile::LoadProfile(const G4S
161 {
162 std::ostringstream ost;
163 ost << fileName << "-" << Z << ".dat";
164 G4String name(ost.str());
165
166 const char* path = G4FindDataDir("G4LEDATA")
167 if (!path)
168 {
169 G4String excep("G4LEDATA environment var
170 G4Exception("G4RDDopplerProfile::LoadPro
171 "InvalidSetup", FatalExcepti
172 }
173
174 G4String pathString(path);
175 G4String dirFile = pathString + name;
176 std::ifstream file(dirFile);
177 std::filebuf* lsdp = file.rdbuf();
178
179 if (! (lsdp->is_open()) )
180 {
181 G4String s1("Data file: ");
182 G4String s2(" not found");
183 G4String excep = s1 + dirFile + s2;
184 G4Exception("G4RDDopplerProfile::LoadPro
185 "DataNotFound", FatalExcepti
186 }
187
188 G4double p;
189 G4int nShell = 0;
190
191 // Create CompositeDataSet for the current Z
192 G4RDVDataSetAlgorithm* interpolation = new G
193 G4RDVEMDataSet* dataSetForZ = new G4RDCompos
194
195 while (!file.eof())
196 {
197 nShell++;
198 G4DataVector* profi = new G4DataVector;
199 G4DataVector* biggs = new G4DataVector;
200
201 // Read in profile data for the current
202 for (size_t i=0; i<nBiggs; i++)
203 {
204 file >> p;
205 profi->push_back(p);
206 biggs->push_back(biggsP[i]);
207 // if (i == 16) G4cout << "profile = "
208 }
209
210 // Create G4RDEMDataSet for the current
211 G4RDVDataSetAlgorithm* algo = interpolat
212 G4RDVEMDataSet* dataSet = new G4RDEMData
213
214 // Add current shell profile component t
215 dataSetForZ->AddComponent(dataSet);
216 }
217
218 // Fill in number of shells for the current
219 nShells.push_back(nShell);
220
221 profileMap[Z] = dataSetForZ;
222 }
223
224
225 G4double G4RDDopplerProfile::RandomSelectMome
226 {
227 G4double value = 0.;
228 const G4RDVEMDataSet* profis = Profiles(Z);
229 value = profis->RandomSelect(shellIndex);
230 return value;
231 }
232