| 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 "G4DopplerProfile.hh"
37 #include "G4DataVector.hh"
38 #include "G4SystemOfUnits.hh"
39 #include "G4VEMDataSet.hh"
40 #include "G4EMDataSet.hh"
41 #include "G4CompositeEMDataSet.hh"
42 #include "G4VDataSetAlgorithm.hh"
43 #include "G4LogLogInterpolation.hh"
44
45 #include <fstream>
46 #include <sstream>
47 #include "Randomize.hh"
48
49 //....oooOO0OOooo........oooOO0OOooo........oo
50
51 G4DopplerProfile::G4DopplerProfile(G4int minZ,
52 : zMin(minZ), zMax(maxZ)
53 {
54 LoadBiggsP("/doppler/p-biggs");
55
56 for (G4int Z=zMin; Z<zMax+1; Z++)
57 {
58 LoadProfile("/doppler/profile",Z);
59 }
60 }
61
62 //....oooOO0OOooo........oooOO0OOooo........oo
63
64 // Destructor
65 G4DopplerProfile::~G4DopplerProfile()
66 {
67 for (auto& pos : profileMap)
68 {
69 G4VEMDataSet* dataSet = pos.second;
70 delete dataSet;
71 dataSet = nullptr;
72 }
73 }
74
75 //....oooOO0OOooo........oooOO0OOooo........oo
76
77 size_t G4DopplerProfile::NumberOfProfiles(G4in
78 {
79 G4int n = 0;
80 if (Z>= zMin && Z <= zMax) n = nShells[Z-1];
81 return n;
82 }
83
84 //....oooOO0OOooo........oooOO0OOooo........oo
85
86 const G4VEMDataSet* G4DopplerProfile::Profiles
87 {
88 if (Z < zMin || Z > zMax)
89 G4Exception("G4DopplerProfile::Profiles",
90 "em1005",FatalException,"Z outside bou
91 auto pos = profileMap.find(Z);
92 G4VEMDataSet* dataSet = (*pos).second;
93 return dataSet;
94 }
95
96 //....oooOO0OOooo........oooOO0OOooo........oo
97
98 const G4VEMDataSet* G4DopplerProfile::Profile(
99 {
100 const G4VEMDataSet* profis = Profiles(Z);
101 const G4VEMDataSet* profi = profis->GetCompo
102 return profi;
103 }
104
105 //....oooOO0OOooo........oooOO0OOooo........oo
106
107 void G4DopplerProfile::PrintData() const
108 {
109 for (G4int Z=zMin; Z<zMax; Z++)
110 {
111 const G4VEMDataSet* profis = Profiles(Z)
112 profis->PrintData();
113 }
114 }
115
116 //....oooOO0OOooo........oooOO0OOooo........oo
117
118 void G4DopplerProfile::LoadBiggsP(const G4Stri
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 G4Exception("G4DopplerProfile::LoadBiggs
128 "em0006",FatalException,"G4LEDATA envi
129 return;
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("G4DopplerProfile::LoadBiggs
143 "em0003",FatalException,excep);
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("G4DopplerProfile::LoadBiggsP"
156 "em1006",FatalException,"Number of mom
157 }
158
159 //....oooOO0OOooo........oooOO0OOooo........oo
160
161 void G4DopplerProfile::LoadProfile(const G4Str
162 {
163 std::ostringstream ost;
164 ost << fileName << "-" << Z << ".dat";
165 G4String name(ost.str());
166
167 const char* path = G4FindDataDir("G4LEDATA")
168 if (!path)
169 {
170 G4String excep("G4LEDATA environment var
171 G4Exception("G4DopplerProfile::LoadProfi
172 "em0006",FatalException,excep);
173 return;
174 }
175
176 G4String pathString(path);
177 G4String dirFile = pathString + name;
178 std::ifstream file(dirFile);
179 std::filebuf* lsdp = file.rdbuf();
180
181 if (! (lsdp->is_open()) )
182 {
183 G4String s1("data file: ");
184 G4String s2(" not found");
185 G4String excep = s1 + dirFile + s2;
186 G4Exception("G4DopplerProfile::LoadProfi
187 "em0003",FatalException,excep);
188 }
189
190 G4double p;
191 G4int nShell = 0;
192
193 // Create CompositeDataSet for the current Z
194 G4VDataSetAlgorithm* interpolation = new G4L
195 G4VEMDataSet* dataSetForZ = new G4CompositeE
196
197 while (!file.eof())
198 {
199 nShell++;
200 G4DataVector* profi = new G4DataVector;
201 G4DataVector* biggs = new G4DataVector;
202
203 // Read in profile data for the current
204 for (size_t i=0; i<nBiggs; i++)
205 {
206 file >> p;
207 profi->push_back(p);
208 biggs->push_back(biggsP[i]);
209 }
210
211 // Create G4EMDataSet for the current sh
212 G4VDataSetAlgorithm* algo = interpolatio
213 G4VEMDataSet* dataSet = new G4EMDataSet(
214
215 // Add current shell profile component t
216 dataSetForZ->AddComponent(dataSet);
217 }
218
219 // Fill in number of shells for the current
220 nShells.push_back(nShell);
221
222 profileMap[Z] = dataSetForZ;
223 }
224
225 //....oooOO0OOooo........oooOO0OOooo........oo
226
227 G4double G4DopplerProfile::RandomSelectMoment
228 {
229 G4double value = 0.;
230 const G4VEMDataSet* profis = Profiles(Z);
231 value = profis->RandomSelect(shellIndex);
232 return value;
233 }
234