84 minLoss(0.001*CLHEP::
eV)
118 if(meanLoss <=
minLoss)
return meanLoss;
145 siga *= 0.25*(1.0 +
x)*(x3 + (1.0/b2 - 0.5)/(1.0/beta2 - 0.5) );
149 G4double twomeanLoss = meanLoss + meanLoss;
159 }
while (0.0 > loss || twomeanLoss < loss);
162 }
else if(sn > 0.1) {
169 loss = twomeanLoss*rndmEngine->
flat();
194 G4double siga = (1. -
beta2*0.5)*tmax*length*electronDensity*
214 /(tmax*(1.0 -
beta2));
215 if(fac_cut > 0.01 && fac > 0.01) {
242 {-0.3291, -0.8312, 0.2460, -1.0220},
243 {-0.5615, -0.5898, 0.5205, -0.7258},
244 {-0.5280, -0.4981, 0.5519, -0.5865},
245 {-0.5125, -0.4625, 0.5660, -0.5190},
246 {-0.5127, -0.8595, 0.5626, -0.8721},
247 {-0.5174, -1.1930, 0.5565, -1.1980},
248 {-0.5179, -1.1850, 0.5560, -1.2070},
249 {-0.5209, -0.9355, 0.5590, -1.0250},
250 {-0.5255, -0.7766, 0.5720, -0.9412},
252 {-0.5776, -0.6665, 0.6598, -0.8484},
253 {-0.6013, -0.6045, 0.7321, -0.7671},
254 {-0.5781, -0.5518, 0.7605, -0.6919},
255 {-0.5587, -0.4981, 0.7835, -0.6195},
256 {-0.5466, -0.4656, 0.7978, -0.5771},
257 {-0.5406, -0.4690, 0.8031, -0.5718},
258 {-0.5391, -0.5061, 0.8024, -0.5974},
259 {-0.5380, -0.6483, 0.7962, -0.6970},
260 {-0.5355, -0.7722, 0.7962, -0.7839},
261 {-0.5329, -0.7720, 0.7988, -0.7846},
263 {-0.5335, -0.7671, 0.7984, -0.7933},
264 {-0.5324, -0.7612, 0.7998, -0.8031},
265 {-0.5305, -0.7300, 0.8031, -0.7990},
266 {-0.5307, -0.7178, 0.8049, -0.8216},
267 {-0.5248, -0.6621, 0.8165, -0.7919},
268 {-0.5180, -0.6502, 0.8266, -0.7986},
269 {-0.5084, -0.6408, 0.8396, -0.8048},
270 {-0.4967, -0.6331, 0.8549, -0.8093},
271 {-0.4861, -0.6508, 0.8712, -0.8432},
272 {-0.4700, -0.6186, 0.8961, -0.8132},
274 {-0.4545, -0.5720, 0.9227, -0.7710},
275 {-0.4404, -0.5226, 0.9481, -0.7254},
276 {-0.4288, -0.4778, 0.9701, -0.6850},
277 {-0.4199, -0.4425, 0.9874, -0.6539},
278 {-0.4131, -0.4188, 0.9998, -0.6332},
279 {-0.4089, -0.4057, 1.0070, -0.6218},
280 {-0.4039, -0.3913, 1.0150, -0.6107},
281 {-0.3987, -0.3698, 1.0240, -0.5938},
282 {-0.3977, -0.3608, 1.0260, -0.5852},
283 {-0.3972, -0.3600, 1.0260, -0.5842},
285 {-0.3985, -0.3803, 1.0200, -0.6013},
286 {-0.3985, -0.3979, 1.0150, -0.6168},
287 {-0.3968, -0.3990, 1.0160, -0.6195},
288 {-0.3971, -0.4432, 1.0050, -0.6591},
289 {-0.3944, -0.4665, 1.0010, -0.6825},
290 {-0.3924, -0.5109, 0.9921, -0.7235},
291 {-0.3882, -0.5158, 0.9947, -0.7343},
292 {-0.3838, -0.5125, 0.9999, -0.7370},
293 {-0.3786, -0.4976, 1.0090, -0.7310},
294 {-0.3741, -0.4738, 1.0200, -0.7155},
296 {-0.3969, -0.4496, 1.0320, -0.6982},
297 {-0.3663, -0.4297, 1.0430, -0.6828},
298 {-0.3630, -0.4120, 1.0530, -0.6689},
299 {-0.3597, -0.3964, 1.0620, -0.6564},
300 {-0.3555, -0.3809, 1.0720, -0.6454},
301 {-0.3525, -0.3607, 1.0820, -0.6289},
302 {-0.3505, -0.3465, 1.0900, -0.6171},
303 {-0.3397, -0.3570, 1.1020, -0.6384},
304 {-0.3314, -0.3552, 1.1130, -0.6441},
305 {-0.3235, -0.3531, 1.1230, -0.6498},
307 {-0.3150, -0.3483, 1.1360, -0.6539},
308 {-0.3060, -0.3441, 1.1490, -0.6593},
309 {-0.2968, -0.3396, 1.1630, -0.6649},
310 {-0.2935, -0.3225, 1.1760, -0.6527},
311 {-0.2797, -0.3262, 1.1940, -0.6722},
312 {-0.2704, -0.3202, 1.2100, -0.6770},
313 {-0.2815, -0.3227, 1.2480, -0.6775},
314 {-0.2880, -0.3245, 1.2810, -0.6801},
315 {-0.3034, -0.3263, 1.3270, -0.6778},
316 {-0.2936, -0.3215, 1.3430, -0.6835},
318 {-0.3282, -0.3200, 1.3980, -0.6650},
319 {-0.3260, -0.3070, 1.4090, -0.6552},
320 {-0.3511, -0.3074, 1.4470, -0.6442},
321 {-0.3501, -0.3064, 1.4500, -0.6442},
322 {-0.3490, -0.3027, 1.4550, -0.6418},
323 {-0.3487, -0.3048, 1.4570, -0.6447},
324 {-0.3478, -0.3074, 1.4600, -0.6483},
325 {-0.3501, -0.3283, 1.4540, -0.6669},
326 {-0.3494, -0.3373, 1.4550, -0.6765},
327 {-0.3485, -0.3373, 1.4570, -0.6774},
329 {-0.3462, -0.3300, 1.4630, -0.6728},
330 {-0.3462, -0.3225, 1.4690, -0.6662},
331 {-0.3453, -0.3094, 1.4790, -0.6553},
332 {-0.3844, -0.3134, 1.5240, -0.6412},
333 {-0.3848, -0.3018, 1.5310, -0.6303},
334 {-0.3862, -0.2955, 1.5360, -0.6237},
335 {-0.4262, -0.2991, 1.5860, -0.6115},
336 {-0.4278, -0.2910, 1.5900, -0.6029},
337 {-0.4303, -0.2817, 1.5940, -0.5927},
338 {-0.4315, -0.2719, 1.6010, -0.5829},
340 {-0.4359, -0.2914, 1.6050, -0.6010},
341 {-0.4365, -0.2982, 1.6080, -0.6080},
342 {-0.4253, -0.3037, 1.6120, -0.6150},
343 {-0.4335, -0.3245, 1.6160, -0.6377},
344 {-0.4307, -0.3292, 1.6210, -0.6447},
345 {-0.4284, -0.3204, 1.6290, -0.6380},
346 {-0.4227, -0.3217, 1.6360, -0.6438}
350 if( 0 > iz ) { iz = 0; }
351 else if(95 < iz ) { iz = 95; }
360 if(ss < slim) { s1 = 1.0/slim; }
362 else if(s1*ss < 1.0) { s1 = 1.0/
ss; }
373 {0.1014, 0.3700, 0.9642, 3.987},
374 {0.1955, 0.6941, 2.522, 1.040},
375 {0.05058, 0.08975, 0.1419, 10.80},
376 {0.05009, 0.08660, 0.2751, 3.787},
377 {0.01273, 0.03458, 0.3951, 3.812}
406 if(y <= 0.2) x *= (y*(1.0 - 0.5*
y));
410 y = energy - b[i][1];
412 G4double s2 = factor * x * b[i][0] / (y*y + x*
x);
432 else f *=
G4Log(4.0*eF/I);
virtual void InitialiseMe(const G4ParticleDefinition *) final
std::vector< ExP01TrackerHit * > a
G4double A13(G4double A) const
static constexpr double MeV
virtual void InitialiseMe(const G4ParticleDefinition *) override
static constexpr double keV
static constexpr double twopi_mc2_rcl2
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double, G4double, G4double) override
virtual G4double Dispersion(const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length) override
G4double GetPDGCharge() const
virtual void SetParticleAndCharge(const G4ParticleDefinition *, G4double q2) override
G4double Factor(const G4Material *, G4double Zeff)
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss) override
G4double GetPDGMass() const
G4double G4Log(G4double x)
G4UniversalFluctuation uniFluct
static constexpr double proton_mass_c2
virtual ~G4IonFluctuations()
G4double expA(G4double A) const
static constexpr double amu_c2
static G4Pow * GetInstance()
G4double GetFermiEnergy() const
G4double powA(G4double A, G4double y) const
static constexpr double electron_mass_c2
virtual void SetParticleAndCharge(const G4ParticleDefinition *, G4double q2) final
static constexpr double eV
G4double GetMeanExcitationEnergy() const
ThreeVector shoot(const G4int Ap, const G4int Af)
const G4ParticleDefinition * particle
static constexpr double eplus
static const G4double fac
G4IonisParamMat * GetIonisation() const
G4double GetZeffective() const
G4double RelativisticFactor(const G4Material *, G4double Zeff)
G4double GetKineticEnergy() const
const G4Material * GetMaterial() const
G4IonFluctuations(const G4String &nam="IonFluc")
G4double GetElectronDensity() const
size_t GetNumberOfElements() const