Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4PenelopeSamplingData Class Reference

#include <G4PenelopeSamplingData.hh>

Public Member Functions

 G4PenelopeSamplingData (G4int npoints=150)
 
 ~G4PenelopeSamplingData ()
 
void AddPoint (G4double x0, G4double pac0, G4double a0, G4double b0, size_t ITTL0, size_t ITTU0)
 
size_t GetNumberOfStoredPoints ()
 
void Clear ()
 
void DumpTable ()
 
G4double GetX (size_t index)
 
G4double GetPAC (size_t index)
 
G4double GetA (size_t index)
 
G4double GetB (size_t index)
 
G4double SampleValue (G4double rndm)
 

Detailed Description

Definition at line 49 of file G4PenelopeSamplingData.hh.

Constructor & Destructor Documentation

G4PenelopeSamplingData::G4PenelopeSamplingData ( G4int  npoints = 150)

Definition at line 37 of file G4PenelopeSamplingData.cc.

37  :
38  np(nPoints)
39 {
40  //create vectors
41  x = new G4DataVector();
42  pac = new G4DataVector();
43  a = new G4DataVector();
44  b = new G4DataVector();
45  ITTL = new std::vector<size_t>;
46  ITTU = new std::vector<size_t>;
47 }
G4PenelopeSamplingData::~G4PenelopeSamplingData ( )

Definition at line 50 of file G4PenelopeSamplingData.cc.

51 {
52  if (x) delete x;
53  if (pac) delete pac;
54  if (a) delete a;
55  if (b) delete b;
56  if (ITTL) delete ITTL;
57  if (ITTU) delete ITTU;
58 }

Member Function Documentation

void G4PenelopeSamplingData::AddPoint ( G4double  x0,
G4double  pac0,
G4double  a0,
G4double  b0,
size_t  ITTL0,
size_t  ITTU0 
)

Definition at line 97 of file G4PenelopeSamplingData.cc.

99 {
100  x->push_back(x0);
101  pac->push_back(pac0);
102  a->push_back(a0);
103  b->push_back(b0);
104  ITTL->push_back(ITTL0);
105  ITTU->push_back(ITTU0);
106 
107  //check how many points we do have now
108  size_t nOfPoints = GetNumberOfStoredPoints();
109 
110  if (nOfPoints > ((size_t)np))
111  {
112  G4cout << "G4PenelopeSamplingData::AddPoint() " << G4endl;
113  G4cout << "WARNING: Up to now there are " << nOfPoints << " points in the table" << G4endl;
114  G4cout << "while the anticipated (declared) number is " << np << G4endl;
115  }
116  return;
117 }
const G4double a0
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

void G4PenelopeSamplingData::Clear ( )

Definition at line 79 of file G4PenelopeSamplingData.cc.

80 {
81  if (x) delete x;
82  if (pac) delete pac;
83  if (a) delete a;
84  if (b) delete b;
85  if (ITTL) delete ITTL;
86  if (ITTU) delete ITTU;
87  //create vectors
88  x = new G4DataVector();
89  pac = new G4DataVector();
90  a = new G4DataVector();
91  b = new G4DataVector();
92  ITTL = new std::vector<size_t>;
93  ITTU = new std::vector<size_t>;
94 }
void G4PenelopeSamplingData::DumpTable ( )

Definition at line 120 of file G4PenelopeSamplingData.cc.

121 {
122 
123  G4cout << "*************************************************************************" << G4endl;
124  G4cout << GetNumberOfStoredPoints() << " points" << G4endl;
125  G4cout << "*************************************************************************" << G4endl;
126  for (size_t i=0;i<GetNumberOfStoredPoints();i++)
127  {
128  G4cout << i << " " << (*x)[i] << " " << (*pac)[i] << " " << (*a)[i] << " " <<
129  (*b)[i] << " " << (*ITTL)[i] << " " << (*ITTU)[i] << G4endl;
130  }
131  G4cout << "*************************************************************************" << G4endl;
132 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

G4double G4PenelopeSamplingData::GetA ( size_t  index)

Definition at line 153 of file G4PenelopeSamplingData.cc.

154 {
155  if (index < a->size())
156  return (*a)[index];
157  else
158  return 0;
159 }
G4double G4PenelopeSamplingData::GetB ( size_t  index)

Definition at line 162 of file G4PenelopeSamplingData.cc.

163 {
164  if (index < b->size())
165  return (*b)[index];
166  else
167  return 0;
168 }
size_t G4PenelopeSamplingData::GetNumberOfStoredPoints ( )

Definition at line 61 of file G4PenelopeSamplingData.cc.

62 {
63  size_t points = x->size();
64 
65  //check everything is all right
66  if (pac->size() != points || a->size() != points ||
67  b->size() != points || ITTL->size() != points ||
68  ITTU->size() != points)
69  {
71  ed << "Data vectors look to have different dimensions !" << G4endl;
72  G4Exception("G4PenelopeSamplingData::GetNumberOfStoredPoints()","em2040",
73  FatalException,ed);
74  }
75  return points;
76 }
std::ostringstream G4ExceptionDescription
Definition: globals.hh:76
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4PenelopeSamplingData::GetPAC ( size_t  index)

Definition at line 144 of file G4PenelopeSamplingData.cc.

145 {
146  if (index < pac->size())
147  return (*pac)[index];
148  else
149  return 0;
150 }
G4double G4PenelopeSamplingData::GetX ( size_t  index)

Definition at line 135 of file G4PenelopeSamplingData.cc.

136 {
137  if (index < x->size())
138  return (*x)[index];
139  else
140  return 0;
141 }

Here is the caller graph for this function:

G4double G4PenelopeSamplingData::SampleValue ( G4double  rndm)

Definition at line 171 of file G4PenelopeSamplingData.cc.

172 {
173  //One passes here a random number in (0,1).
174  //Notice: it possible that is between (0,b) with b<1
175  size_t points = GetNumberOfStoredPoints();
176 
177  size_t itn = (size_t) (maxRand*(points-1));
178  size_t i = (*ITTL)[itn];
179  size_t j = (*ITTU)[itn];
180 
181  while ((j-i) > 1)
182  {
183  size_t k = (i+j)/2;
184  if (maxRand > (*pac)[k])
185  i = k;
186  else
187  j = k;
188  }
189 
190  //Sampling from the rational inverse cumulative distribution
191  G4double result = 0;
192 
193  G4double rr = maxRand - (*pac)[i];
194  if (rr > 1e-16)
195  {
196  G4double d = (*pac)[i+1]-(*pac)[i];
197  result = (*x)[i]+
198  ((1.0+(*a)[i]+(*b)[i])*d*rr/
199  (d*d+((*a)[i]*d+(*b)[i]*rr)*rr))*((*x)[i+1]-(*x)[i]);
200  }
201  else
202  result = (*x)[i];
203 
204  return result;
205 }
G4double G4ParticleHPJENDLHEData::G4double result
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:


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