Geant4-11
Public Member Functions | Private Member Functions | Private Attributes
G4LindhardSorensenData Class Reference

#include <G4LindhardSorensenData.hh>

Public Member Functions

 G4LindhardSorensenData ()
 
 G4LindhardSorensenData (const G4LindhardSorensenData &)=delete
 
G4double GetDeltaL (G4int Z, G4double gamma) const
 
G4LindhardSorensenData & operator= (const G4LindhardSorensenData &right)=delete
 
 ~G4LindhardSorensenData ()
 

Private Member Functions

G4double ComputeDeltaL (G4int Z, G4double x) const
 
void InitialiseData ()
 

Private Attributes

G4PhysicsVector * data [9]
 
G4Pow * g4calc
 
G4double xmax
 
G4double xmin
 

Detailed Description

Definition at line 52 of file G4LindhardSorensenData.hh.

Constructor & Destructor Documentation

◆ G4LindhardSorensenData() [1/2]

G4LindhardSorensenData::G4LindhardSorensenData ( )
explicit

Definition at line 56 of file G4LindhardSorensenData.cc.

57{
58 g4calc = G4Pow::GetInstance();
59 InitialiseData();
60}
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:41

References g4calc, G4Pow::GetInstance(), and InitialiseData().

◆ ~G4LindhardSorensenData()

G4LindhardSorensenData::~G4LindhardSorensenData ( )

Definition at line 64 of file G4LindhardSorensenData.cc.

65{
66 for(G4int i=0; i<=LVECT; ++i) { delete data[i]; }
67}
LVECT
const G4int LVECT
Definition: G4LindhardSorensenData.cc:51
G4int
int G4int
Definition: G4Types.hh:85

References data, and LVECT.

◆ G4LindhardSorensenData() [2/2]

G4LindhardSorensenData::G4LindhardSorensenData ( const G4LindhardSorensenData &  )
delete

Member Function Documentation

◆ ComputeDeltaL()

G4double G4LindhardSorensenData::ComputeDeltaL ( G4int  Z,
G4double  x 
) const
private

Definition at line 94 of file G4LindhardSorensenData.cc.

95{
96 G4double y(0.0);
97 if(x < xmin) {
98 G4double x1 = (data[idx])->Energy(1);
99 G4double ymin = (*(data[idx]))[0];
100 G4double y1 = (*(data[idx]))[1];
101 y = ymin + (y1 - ymin)*(x - xmin)/(x1 - xmin);
102 } else if(x >= xmax) {
103 G4double x1 = (data[idx])->Energy(NPOINT-2);
104 G4double ymax = (*(data[idx]))[NPOINT-1];
105 G4double y1 = (*(data[idx]))[NPOINT-2];
106 y = y1 + (ymax - y1)*(x - x1)/(xmax - x1);
107 } else {
108 y = (data[idx])->Value(x);
109 }
110 return y;
111}
NPOINT
const G4int NPOINT
Definition: G4LindhardSorensenData.cc:52
G4double
double G4double
Definition: G4Types.hh:83

References data, NPOINT, xmax, and xmin.

Referenced by GetDeltaL().

◆ GetDeltaL()

G4double G4LindhardSorensenData::GetDeltaL ( G4int  Z,
G4double  gamma 
) const

Definition at line 71 of file G4LindhardSorensenData.cc.

72{
73 G4int idx = 0;
74 for(; idx<LVECT; ++idx) {
75 if(Z < zlist[idx+1]) { break; }
76 }
77 idx = std::min(idx, LVECT);
78
79 G4double x = G4Log(gamma - 1.0);
80 G4double y = ComputeDeltaL(idx, x);
81
82 // interpolation over Z if needed
83 if(idx < LVECT && Z > zlist[idx]) {
84 G4double y1 = ComputeDeltaL(idx+1, x);
85 //G4cout << "idx= " << idx << " x= " << x << " y= " << y << " y1= " << y1 << G4endl;
86 y += (y1 - y)*(Z - zlist[idx])/(G4double)(zlist[idx+1] - zlist[idx]);
87 }
88 y *= g4calc->Z23(Z);
89 return y;
90}
zlist
const G4int zlist[9]
Definition: G4LindhardSorensenData.cc:50
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:226
Z
const G4int Z[17]
Definition: G4WaterStopping.cc:51
G4LindhardSorensenData::ComputeDeltaL
G4double ComputeDeltaL(G4int Z, G4double x) const
Definition: G4LindhardSorensenData.cc:94
G4Pow::Z23
G4double Z23(G4int Z) const
Definition: G4Pow.hh:125
G4INCL::Math::min
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
Definition: G4INCLGlobals.hh:117

References ComputeDeltaL(), g4calc, G4Log(), LVECT, G4INCL::Math::min(), Z, G4Pow::Z23(), and zlist.

Referenced by G4LindhardSorensenIonModel::CorrectionsAlongStep().

◆ InitialiseData()

void G4LindhardSorensenData::InitialiseData ( )
private

Definition at line 115 of file G4LindhardSorensenData.cc.

116{
117 xmin = G4Log(0.02);
118 xmax = G4Log(316.22777);
119
120 const G4double lsdata[9][NPOINT] = {
121 {0.0036181621,
122 0.0042618872, 0.0040786701, 0.0039097273, 0.0041240731, 0.0044311195, // 0-5
123 0.0059858073, 0.0079655897, 0.0089669217, 0.010091248, 0.0096319233, // 5-10
124 0.0085344428, 0.0087326058, 0.009780094, 0.010617094, 0.011277997, // 10-15
125 0.01185287, 0.012242278, 0.012507042, 0.012369698, 0.012202436, // 15-20
126 0.012016446, 0.011793182, 0.011612196, 0.011375017, 0.010624408, // 20-25
127 0.0096685612, 0.0093765113, 0.0091152632, 0.0087458522, 0.0082645153, // 25-30
128 0.0068625676, 0.0046558921, 0.0021660968, -0.00072760644, -0.0062653709,// 30-35
129 -0.015834368, -0.03224412, -0.062296044, -0.10327705, -0.16545368},
130 {-0.024837796,
131 -0.018943357, -0.013251703, -0.007646936, -0.003751624, 7.4932758e-05, // 0-5
132 0.0027996278, 0.0057534909, 0.0078441157, 0.010341958, 0.012419648, // 5-10
133 0.01491866, 0.016403769, 0.018106955, 0.019249728, 0.020704966, // 10-15
134 0.021636244, 0.022526456, 0.023429501, 0.024251631, 0.024727856, // 15-20
135 0.024997337, 0.025225657, 0.025335138, 0.025474553, 0.025306356, // 20-25
136 0.025048151, 0.02467553, 0.02416496, 0.023519917, 0.022699646, // 25-30
137 0.020902526, 0.017645158, 0.013219058, 0.0059989118, -0.0042559395, // 30-35
138 -0.020223315, -0.044480728, -0.074942757, -0.1108863, -0.15439805},
139 {-0.04948514,
140 -0.039259418, -0.029641846, -0.022112859, -0.0152972, -0.0077213168, // 0-5
141 -0.0031032447, 0.0023506153, 0.0058573025, 0.0097554723, 0.012692892, // 5-10
142 0.016404575, 0.018642388, 0.021492082, 0.023328047, 0.025358859, // 10-15
143 0.026839018, 0.028169751, 0.029240668, 0.030035325, 0.030689491, // 15-20
144 0.031113343, 0.031518392, 0.031657748, 0.031794607, 0.031672646, // 20-25
145 0.031507453, 0.031024952, 0.030410522, 0.029266343, 0.027729316, // 25-30
146 0.025132583, 0.02141173, 0.016127626, 0.008039769, -0.0046793613, // 30-35
147 -0.02196456, -0.043877567, -0.070852217, -0.09979099, -0.13187648},
148 {-0.07502957,
149 -0.063927817, -0.053895181, -0.043621797, -0.033329203, -0.024765223, // 0-5
150 -0.015996122, -0.0084995741, -0.00017499271, 0.0064195209, 0.01117532,// 5-10
151 0.016956954, 0.02155737, 0.025213602, 0.029869549, 0.031980316, // 10-15
152 0.034664781, 0.036995048, 0.038244843, 0.039836367, 0.040930356, // 15-20
153 0.041677336, 0.042211864, 0.042472572, 0.042617755, 0.042393662, // 20-25
154 0.042049893, 0.041341375, 0.040408112, 0.03867031, 0.036457345, // 25-30
155 0.032742289, 0.027897265, 0.020305954, 0.0099602019, -0.0036658833, // 30-35
156 -0.020344557, -0.039700638, -0.059781744, -0.080649004, -0.10263347},
157 {-0.081885964,
158 -0.072049323, -0.06299877, -0.053400445, -0.043830227, -0.03415406, // 0-5
159 -0.024723298, -0.015706543, -0.0070305988, 0.0016621818, 0.0093560784, // 5-10
160 0.017338568, 0.023286722, 0.028930541, 0.034048359, 0.037750635, // 10-15
161 0.041512165, 0.04362858, 0.046323682, 0.048044285, 0.049379817, // 15-20
162 0.050501289, 0.050963981, 0.051463496, 0.051292532, 0.051074821, // 20-25
163 0.050350498, 0.049370817, 0.047660517, 0.045247985, 0.04188952, // 25-30
164 0.037009647, 0.030171627, 0.021186013, 0.0092064517, -0.0051358689, // 30-35
165 -0.01977885, -0.036354277, -0.053187271, -0.069851489, -0.087776477},
166 {-0.083068958,
167 -0.074302116, -0.065040071, -0.056500699, -0.047189921, -0.03781489, // 0-5
168 -0.028666422, -0.0192974, -0.0097949279, -0.0011649946, 0.0071301285, // 5-10
169 0.01548354, 0.023483416, 0.030515545, 0.036162102, 0.040602642, // 10-15
170 0.045011738, 0.048259036, 0.051136495, 0.053524791, 0.054971064, // 15-20
171 0.056247898, 0.056729496, 0.057135059, 0.056923775, 0.05657171, // 20-25
172 0.055504504, 0.054134086, 0.051812977, 0.048857213, 0.044525663, // 25-30
173 0.038593442, 0.03027527, 0.020452704, 0.0082978722, -0.0045556908, // 30-35
174 -0.018750622, -0.033642805, -0.048788529, -0.064509221, -0.080215679},
175 {-0.083061344,
176 -0.075180377, -0.066630961, -0.058114277, -0.049303999, -0.040618218, // 0-5
177 -0.031608369, -0.022553882, -0.012699373, -0.0033434259, 0.0051129584, // 5-10
178 0.013817276, 0.023132109, 0.030886148, 0.037688311, 0.043478106, // 10-15
179 0.048602814, 0.052502168, 0.055969817, 0.058142903, 0.060860174, // 15-20
180 0.061732795, 0.062844014, 0.062911637, 0.06278651, 0.061898673, // 20-25
181 0.060500452, 0.058323088, 0.055297604, 0.051242701, 0.04596475, // 25-30
182 0.038863441, 0.030302781, 0.019055839, 0.0072837125, -0.0047516889, // 30-35
183 -0.017733076, -0.031547911, -0.04523347, -0.059727542, -0.074202114},
184 {-0.082449782,
185 -0.07492693, -0.067041495, -0.059614338, -0.051061455, -0.042587329, // 0-5
186 -0.034242658, -0.025296495, -0.016157818, -0.0064198954, 0.0033148315, // 5-10
187 0.012927383, 0.021914897, 0.03053429, 0.03854699, 0.045293871, // 10-15
188 0.05141546, 0.056473799, 0.060050334, 0.063292547, 0.065215264, // 15-20
189 0.066961692, 0.067369593, 0.067623539, 0.066825274, 0.065740573, // 20-25
190 0.063673411, 0.061033708, 0.057399807, 0.052544441, 0.04583072, // 25-30
191 0.037931379, 0.028569029, 0.017544895, 0.0068751554, -0.0038105519, // 30-35
192 -0.015233268, -0.026956839, -0.038975048, -0.051387898, -0.064344384},
193 {-0.081232852,
194 -0.074976912, -0.068421001, -0.060951614, -0.053624488, -0.045755024, // 0-5
195 -0.037708151, -0.029246587, -0.02018635, -0.01073862, -0.00085387172, // 5-10
196 0.0090831897, 0.018970381, 0.02761289, 0.036215831, 0.044355404, // 10-15
197 0.051067631, 0.057143777, 0.061348214, 0.065462308, 0.06715594, // 15-20
198 0.069155687, 0.069845259, 0.069828188, 0.069170195, 0.067585183, // 20-25
199 0.065486015, 0.062437523, 0.058111347, 0.052410482, 0.045597673, // 25-30
200 0.037571853, 0.028019555, 0.016564385, 0.0060254421, -0.0043254115, // 30-35
201 -0.015675051, -0.027353574, -0.03931666, -0.050325729, -0.06174447}};
202
203 for(G4int i=0; i<=LVECT; ++i) {
204 data[i] = new G4PhysicsLinearVector(xmin, xmax, NPOINT-1, true);
205 for(std::size_t j=0; j<NPOINT; ++j) {
206 data[i]->PutValue(j, lsdata[i][j]);
207 }
208 data[i]->FillSecondDerivatives();
209 }
210}
G4PhysicsVector::PutValue
void PutValue(const std::size_t index, const G4double value)
G4PhysicsVector::FillSecondDerivatives
void FillSecondDerivatives(const G4SplineType=G4SplineType::Base, const G4double dir1=0.0, const G4double dir2=0.0)
Definition: G4PhysicsVector.cc:205

References data, G4PhysicsVector::FillSecondDerivatives(), G4Log(), LVECT, NPOINT, G4PhysicsVector::PutValue(), xmax, and xmin.

Referenced by G4LindhardSorensenData().

◆ operator=()

G4LindhardSorensenData & G4LindhardSorensenData::operator= ( const G4LindhardSorensenData &  right)
delete

Field Documentation

◆ data

G4PhysicsVector* G4LindhardSorensenData::data[9]
private

Definition at line 76 of file G4LindhardSorensenData.hh.

Referenced by ComputeDeltaL(), InitialiseData(), and ~G4LindhardSorensenData().

◆ g4calc

G4Pow* G4LindhardSorensenData::g4calc
private

Definition at line 77 of file G4LindhardSorensenData.hh.

Referenced by G4LindhardSorensenData(), and GetDeltaL().

◆ xmax

G4double G4LindhardSorensenData::xmax
private

Definition at line 74 of file G4LindhardSorensenData.hh.

Referenced by ComputeDeltaL(), and InitialiseData().

◆ xmin

G4double G4LindhardSorensenData::xmin
private

Definition at line 73 of file G4LindhardSorensenData.hh.

Referenced by ComputeDeltaL(), and InitialiseData().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3