Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ShellData Class Reference

#include <G4ShellData.hh>

Public Member Functions

 G4ShellData (G4int minZ=1, G4int maxZ=100, G4bool isOccupancy=false)
 
 ~G4ShellData ()
 
size_t NumberOfShells (G4int Z) const
 
G4int ShellId (G4int Z, G4int shellIndex) const
 
G4double ShellOccupancyProbability (G4int Z, G4int shellIndex) const
 
const std::vector< G4double > & ShellIdVector (G4int Z) const
 
G4double BindingEnergy (G4int Z, G4int shellIndex) const
 
void SetOccupancyData ()
 
void LoadData (const G4String &fileName)
 
void PrintData () const
 
G4int SelectRandomShell (G4int Z) const
 

Detailed Description

Definition at line 52 of file G4ShellData.hh.

Constructor & Destructor Documentation

G4ShellData::G4ShellData ( G4int  minZ = 1,
G4int  maxZ = 100,
G4bool  isOccupancy = false 
)

Definition at line 51 of file G4ShellData.cc.

52  : zMin(minZ), zMax(maxZ), occupancyData(isOccupancy)
53 { }
G4ShellData::~G4ShellData ( )

Definition at line 56 of file G4ShellData.cc.

57 {
58  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos;
59  for (pos = idMap.begin(); pos != idMap.end(); ++pos)
60  {
61  std::vector<G4double>* dataSet = (*pos).second;
62  delete dataSet;
63  }
64 
65  std::map<G4int,G4DataVector*,std::less<G4int> >::iterator pos2;
66  for (pos2 = bindingMap.begin(); pos2 != bindingMap.end(); ++pos2)
67  {
68  G4DataVector* dataSet = (*pos2).second;
69  delete dataSet;
70  }
71 
72  if (occupancyData)
73  {
74  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos3;
75  for (pos3 = occupancyPdfMap.begin(); pos3 != occupancyPdfMap.end(); ++pos3)
76  {
77  std::vector<G4double>* dataSet = (*pos3).second;
78  delete dataSet;
79  }
80  }
81 }
static const G4double pos

Member Function Documentation

G4double G4ShellData::BindingEnergy ( G4int  Z,
G4int  shellIndex 
) const

Definition at line 166 of file G4ShellData.cc.

167 {
168  G4double value = 0.;
169 
170  if (Z >= zMin && Z <= zMax)
171  {
172  std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator pos;
173  pos = bindingMap.find(Z);
174  if (pos!= bindingMap.end())
175  {
176  G4DataVector dataSet = *((*pos).second);
177  G4int nData = dataSet.size();
178  if (shellIndex >= 0 && shellIndex < nData)
179  {
180  value = dataSet[shellIndex];
181  }
182  }
183  }
184  return value;
185 }
int G4int
Definition: G4Types.hh:78
const XML_Char int const XML_Char * value
Definition: expat.h:331
double G4double
Definition: G4Types.hh:76
static const G4double pos

Here is the caller graph for this function:

void G4ShellData::LoadData ( const G4String fileName)

Definition at line 234 of file G4ShellData.cc.

235 {
236  // Build the complete string identifying the file with the data set
237 
238  std::ostringstream ost;
239 
240  ost << fileName << ".dat";
241 
242  G4String name(ost.str());
243 
244  char* path = getenv("G4LEDATA");
245  if (!path)
246  {
247  G4String excep("G4ShellData::LoadData()");
248  G4Exception(excep,"em0006",FatalException,"Please set G4LEDATA");
249  return;
250  }
251 
252  G4String pathString(path);
253  G4String dirFile = pathString + name;
254  std::ifstream file(dirFile);
255  std::filebuf* lsdp = file.rdbuf();
256 
257  if (! (lsdp->is_open()) )
258  {
259 
260  G4String excep = "G4ShellData::LoadData()";
261  G4String msg = "data file: " + dirFile + " not found";
262  G4Exception(excep, "em0003",FatalException, msg );
263  return;
264  }
265 
266  G4double a = 0;
267  G4int k = 1;
268  G4int sLocal = 0;
269 
270  G4int Z = 1;
271  G4DataVector* energies = new G4DataVector;
272  std::vector<G4double>* ids = new std::vector<G4double>;
273 
274  do {
275  file >> a;
276  G4int nColumns = 2;
277  if (a == -1)
278  {
279  if (sLocal == 0)
280  {
281  // End of a shell data set
282  idMap[Z] = ids;
283  bindingMap[Z] = energies;
284  G4int n = ids->size();
285  nShells.push_back(n);
286  // Start of new shell data set
287 
288  ids = new std::vector<G4double>;
289  energies = new G4DataVector;
290  Z++;
291  }
292  sLocal++;
293  if (sLocal == nColumns)
294  {
295  sLocal = 0;
296  }
297  }
298 
299  // moved out of the do-while since might go to a leak.
300  // else if (a == -2)
301  // {
302  // End of file; delete the empty vectors created when encountering the last -1 -1 row
303  // delete energies;
304  // delete ids;
305  //nComponents = components.size();
306  // }
307  else
308  {
309  // 1st column is shell id
310  if(k%nColumns != 0)
311  {
312  ids->push_back(a);
313  k++;
314  }
315  else if (k%nColumns == 0)
316  {
317  // 2nd column is binding energy
318  G4double e = a * MeV;
319  energies->push_back(e);
320  k = 1;
321  }
322  }
323  } while (a != -2); // end of file
324  file.close();
325  delete energies;
326  delete ids;
327 
328  // For Doppler broadening: the data set contains shell occupancy and binding energy for each shell
329  // Build additional map with probability for each shell based on its occupancy
330 
331  if (occupancyData)
332  {
333  // Build cumulative from raw shell occupancy
334 
335  for (G4int ZLocal=zMin; ZLocal <= zMax; ZLocal++)
336  {
337  std::vector<G4double> occupancy = ShellIdVector(ZLocal);
338 
339  std::vector<G4double>* prob = new std::vector<G4double>;
340  G4double scale = 1. / G4double(ZLocal);
341 
342  prob->push_back(occupancy[0] * scale);
343  for (size_t i=1; i<occupancy.size(); i++)
344  {
345  prob->push_back(occupancy[i]*scale + (*prob)[i-1]);
346  }
347  occupancyPdfMap[ZLocal] = prob;
348 
349  /*
350  G4double scale = 1. / G4double(Z);
351  // transform((*prob).begin(),(*prob).end(),(*prob).begin(),bind2nd(multiplies<G4double>(),scale));
352 
353  for (size_t i=0; i<occupancy.size(); i++)
354  {
355  (*prob)[i] *= scale;
356  }
357  */
358  }
359  }
360 }
const XML_Char * name
Definition: expat.h:151
std::vector< ExP01TrackerHit * > a
Definition: ExP01Classes.hh:33
int G4int
Definition: G4Types.hh:78
const std::vector< G4double > & ShellIdVector(G4int Z) const
Definition: G4ShellData.cc:97
const G4int n
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
static constexpr double MeV
Definition: G4SIunits.hh:214
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

size_t G4ShellData::NumberOfShells ( G4int  Z) const

Definition at line 84 of file G4ShellData.cc.

85 {
86  G4int z = Z - 1;
87  G4int n = 0;
88 
89  if (Z>= zMin && Z <= zMax)
90  {
91  n = nShells[z];
92  }
93  return n;
94 }
int G4int
Definition: G4Types.hh:78
const G4int n
tuple z
Definition: test.py:28

Here is the caller graph for this function:

void G4ShellData::PrintData ( void  ) const

Definition at line 187 of file G4ShellData.cc.

188 {
189  for (G4int Z = zMin; Z <= zMax; Z++)
190  {
191  G4cout << "---- Shell data for Z = "
192  << Z
193  << " ---- "
194  << G4endl;
195  G4int nSh = nShells[Z-1];
196  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator posId;
197  posId = idMap.find(Z);
198  std::vector<G4double>* ids = (*posId).second;
199  std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator posE;
200  posE = bindingMap.find(Z);
201  G4DataVector* energies = (*posE).second;
202  for (G4int i=0; i<nSh; i++)
203  {
204  G4int id = (G4int) (*ids)[i];
205  G4double e = (*energies)[i] / keV;
206  G4cout << i << ") ";
207 
208  if (occupancyData)
209  {
210  G4cout << " Occupancy: ";
211  }
212  else
213  {
214  G4cout << " Shell id: ";
215  }
216  G4cout << id << " - Binding energy = "
217  << e << " keV ";
218  if (occupancyData)
219  {
220  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator posOcc;
221  posOcc = occupancyPdfMap.find(Z);
222  std::vector<G4double> probs = *((*posOcc).second);
223  G4double prob = probs[i];
224  G4cout << "- Probability = " << prob;
225  }
226  G4cout << G4endl;
227  }
228  G4cout << "-------------------------------------------------"
229  << G4endl;
230  }
231 }
int G4int
Definition: G4Types.hh:78
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
static constexpr double keV
Definition: G4SIunits.hh:216
G4int G4ShellData::SelectRandomShell ( G4int  Z) const

Definition at line 363 of file G4ShellData.cc.

364 {
365  if (Z < zMin || Z > zMax) {
366 
367  G4Exception("G4ShellData::SelectrandomShell","de0001",FatalErrorInArgument, "Z outside boundaries");
368 
369  }
370  G4int shellIndex = 0;
371  std::vector<G4double> prob = ShellVector(Z);
372  G4double random = G4UniformRand();
373 
374  // std::vector<G4double>::const_iterator pos;
375  // pos = lower_bound(prob.begin(),prob.end(),random);
376 
377  // Binary search the shell with probability less or equal random
378 
379  G4int nShellsLocal = NumberOfShells(Z);
380  G4int upperBound = nShellsLocal;
381 
382  while (shellIndex <= upperBound)
383  {
384  G4int midShell = (shellIndex + upperBound) / 2;
385  if ( random < prob[midShell] )
386  upperBound = midShell - 1;
387  else
388  shellIndex = midShell + 1;
389  }
390  if (shellIndex >= nShellsLocal) shellIndex = nShellsLocal - 1;
391 
392  return shellIndex;
393 }
int G4int
Definition: G4Types.hh:78
#define G4UniformRand()
Definition: Randomize.hh:97
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
size_t NumberOfShells(G4int Z) const
Definition: G4ShellData.cc:84
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ShellData::SetOccupancyData ( )
inline

Definition at line 70 of file G4ShellData.hh.

70 { occupancyData = true; }

Here is the caller graph for this function:

G4int G4ShellData::ShellId ( G4int  Z,
G4int  shellIndex 
) const

Definition at line 121 of file G4ShellData.cc.

122 {
123  G4int n = -1;
124 
125  if (Z >= zMin && Z <= zMax)
126  {
127  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
128  pos = idMap.find(Z);
129  if (pos!= idMap.end())
130  {
131  std::vector<G4double> dataSet = *((*pos).second);
132  G4int nData = dataSet.size();
133  if (shellIndex >= 0 && shellIndex < nData)
134  {
135  n = (G4int) dataSet[shellIndex];
136  }
137  }
138  }
139  return n;
140 }
int G4int
Definition: G4Types.hh:78
const G4int n
static const G4double pos

Here is the caller graph for this function:

const std::vector< G4double > & G4ShellData::ShellIdVector ( G4int  Z) const

Definition at line 97 of file G4ShellData.cc.

98 {
99  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
100  if (Z < zMin || Z > zMax) {
101 
102  G4Exception("G4ShellData::ShellIdVector","de0001",FatalErrorInArgument, "Z outside boundaries");
103 
104 
105  } pos = idMap.find(Z);
106  std::vector<G4double>* dataSet = (*pos).second;
107  return *dataSet;
108 }
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
static const G4double pos

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4ShellData::ShellOccupancyProbability ( G4int  Z,
G4int  shellIndex 
) const

Definition at line 143 of file G4ShellData.cc.

144 {
145  G4double prob = -1.;
146 
147  if (Z >= zMin && Z <= zMax)
148  {
149  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
150  pos = idMap.find(Z);
151  if (pos!= idMap.end())
152  {
153  std::vector<G4double> dataSet = *((*pos).second);
154  G4int nData = dataSet.size();
155  if (shellIndex >= 0 && shellIndex < nData)
156  {
157  prob = dataSet[shellIndex];
158  }
159  }
160  }
161  return prob;
162 }
int G4int
Definition: G4Types.hh:78
double G4double
Definition: G4Types.hh:76
static const G4double pos

The documentation for this class was generated from the following files: