83 fParticleChange(nullptr),
84 currentMaterial(nullptr),
133 pname !=
"deuteron" && pname !=
"triton" &&
134 pname !=
"alpha+" && pname !=
"helium" &&
135 pname !=
"hydrogen") {
isIon =
true; }
187 if(cutEnergy < tmax) {
191 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/energy2;
192 cross = (maxEnergy - cutEnergy)/(cutEnergy*maxEnergy)
193 - beta2*
G4Log(maxEnergy/cutEnergy)/tmax;
195 if( 0.0 <
spin ) { cross += 0.5*(maxEnergy - cutEnergy)/energy2; }
215 (p,kineticEnergy,cutEnergy,maxEnergy);
230 (p,kineticEnergy,cutEnergy,maxEnergy);
249 dedx =
DEDX(material, tkin);
252 if (cutEnergy < tmax) {
289 if(e < 0.0) { e = preKinEnergy*0.5; }
310 if(xmin >= xmax) {
return; }
315 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/energy2;
325 deltaKinEnergy = xmin*xmax/(xmin*(1.0 - rndm[0]) + xmax*rndm[0]);
327 f = 1.0 - beta2*deltaKinEnergy/tmax;
330 G4cout <<
"G4BraggIonModel::SampleSecondary Warning! "
331 <<
"Majorant " << grej <<
" < "
332 << f <<
" for e= " << deltaKinEnergy
337 }
while( grej*rndm[1] >= f );
354 if(cost > 1.0) { cost = 1.0; }
355 G4double sint = sqrt((1.0 - cost)*(1.0 + cost));
359 deltaDirection.
set(sint*cos(phi),sint*sin(phi), cost) ;
367 vdp->push_back(delta);
370 kineticEnergy -= deltaKinEnergy;
372 finalP = finalP.
unit();
435 {9.43672f, 0.54398f, 84.341f, 1.3705f, 57.422f},
436 {67.1503f, 0.41409f, 404.512f, 148.97f, 20.99f},
437 {5.11203f, 0.453f, 36.718f, 50.6f, 28.058f},
438 {61.793f, 0.48445f, 361.537f, 57.889f, 50.674f},
439 {7.83464f, 0.49804f, 160.452f, 3.192f, 0.71922f},
440 {19.729f, 0.52153f, 162.341f, 58.35f, 25.668f},
441 {26.4648f, 0.50112f, 188.913f, 30.079f, 16.509f},
442 {7.8655f, 0.5205f, 63.96f, 51.32f, 67.775f},
443 {8.8965f, 0.5148f, 339.36f, 1.7205f, 0.70423f},
444 {2.959f, 0.53255f, 34.247f, 60.655f, 15.153f},
445 {3.80133f, 0.41590f, 12.9966f, 117.83f, 242.28f} };
447 static const G4double atomicWeight[11] = {
448 101.96128f, 44.0098f, 16.0426f, 28.0536f, 42.0804f,
449 104.1512f, 44.665f, 60.0843f, 18.0152f, 18.0152f, 12.0f};
462 G4double shigh =
G4Log( 1.0 + x3*1000.0 + x4*0.001 ) *x2*1000.0;
463 ionloss = slow*shigh / (slow + shigh) ;
464 ionloss *= sqrt(T*1000.0) ;
470 ionloss = slow*shigh / (slow + shigh) ;
512 { 0.35485f, 0.6456f, 6.01525f, 20.8933f, 4.3515f
513 },{ 0.58f, 0.59f, 6.3f, 130.0f, 44.07f
514 },{ 1.42f, 0.49f, 12.25f, 32.0f, 9.161f
515 },{ 2.206f, 0.51f, 15.32f, 0.25f, 8.995f
517 },{ 3.691f, 0.4128f, 18.48f, 50.72f, 9.0f
518 },{ 3.83523f, 0.42993f,12.6125f, 227.41f, 188.97f
520 },{ 1.9259f, 0.5550f, 27.1513f, 26.0665f, 6.2768f
521 },{ 2.81015f, 0.4759f, 50.0253f, 10.556f, 1.0382f
522 },{ 1.533f, 0.531f, 40.44f, 18.41f, 2.718f
523 },{ 2.303f, 0.4861f, 37.01f, 37.96f, 5.092f
525 },{ 9.894f, 0.3081f, 23.65f, 0.384f, 92.93f
526 },{ 4.3f, 0.47f, 34.3f, 3.3f, 12.74f
527 },{ 2.5f, 0.625f, 45.7f, 0.1f, 4.359f
528 },{ 2.1f, 0.65f, 49.34f, 1.788f, 4.133f
529 },{ 1.729f, 0.6562f, 53.41f, 2.405f, 3.845f
530 },{ 1.402f, 0.6791f, 58.98f, 3.528f, 3.211f
531 },{ 1.117f, 0.7044f, 69.69f, 3.705f, 2.156f
532 },{ 2.291f, 0.6284f, 73.88f, 4.478f, 2.066f
533 },{ 8.554f, 0.3817f, 83.61f, 11.84f, 1.875f
534 },{ 6.297f, 0.4622f, 65.39f, 10.14f, 5.036f
536 },{ 5.307f, 0.4918f, 61.74f, 12.4f, 6.665f
537 },{ 4.71f, 0.5087f, 65.28f, 8.806f, 5.948f
538 },{ 6.151f, 0.4524f, 83.0f, 18.31f, 2.71f
539 },{ 6.57f, 0.4322f, 84.76f, 15.53f, 2.779f
540 },{ 5.738f, 0.4492f, 84.6f, 14.18f, 3.101f
541 },{ 5.013f, 0.4707f, 85.8f, 16.55f, 3.211f
542 },{ 4.32f, 0.4947f, 76.14f, 10.85f, 5.441f
543 },{ 4.652f, 0.4571f, 80.73f, 22.0f, 4.952f
544 },{ 3.114f, 0.5236f, 76.67f, 7.62f, 6.385f
545 },{ 3.114f, 0.5236f, 76.67f, 7.62f, 7.502f
547 },{ 3.114f, 0.5236f, 76.67f, 7.62f, 8.514f
548 },{ 5.746f, 0.4662f, 79.24f, 1.185f, 7.993f
549 },{ 2.792f, 0.6346f, 106.1f, 0.2986f, 2.331f
550 },{ 4.667f, 0.5095f, 124.3f, 2.102f, 1.667f
551 },{ 2.44f, 0.6346f, 105.0f, 0.83f, 2.851f
552 },{ 1.413f, 0.7377f, 147.9f, 1.466f, 1.016f
553 },{ 11.72f, 0.3826f, 102.8f, 9.231f, 4.371f
554 },{ 7.126f, 0.4804f, 119.3f, 5.784f, 2.454f
555 },{ 11.61f, 0.3955f, 146.7f, 7.031f, 1.423f
556 },{ 10.99f, 0.41f, 163.9f, 7.1f, 1.052f
558 },{ 9.241f, 0.4275f, 163.1f, 7.954f, 1.102f
559 },{ 9.276f, 0.418f, 157.1f, 8.038f, 1.29f
560 },{ 3.999f, 0.6152f, 97.6f, 1.297f, 5.792f
561 },{ 4.306f, 0.5658f, 97.99f, 5.514f, 5.754f
562 },{ 3.615f, 0.6197f, 86.26f, 0.333f, 8.689f
563 },{ 5.8f, 0.49f, 147.2f, 6.903f, 1.289f
564 },{ 5.6f, 0.49f, 130.0f, 10.0f, 2.844f
565 },{ 3.55f, 0.6068f, 124.7f, 1.112f, 3.119f
566 },{ 3.6f, 0.62f, 105.8f, 0.1692f, 6.026f
567 },{ 5.4f, 0.53f, 103.1f, 3.931f, 7.767f
569 },{ 3.97f, 0.6459f, 131.8f, 0.2233f, 2.723f
570 },{ 3.65f, 0.64f, 126.8f, 0.6834f, 3.411f
571 },{ 3.118f, 0.6519f, 164.9f, 1.208f, 1.51f
572 },{ 3.949f, 0.6209f, 200.5f, 1.878f, 0.9126f
573 },{ 14.4f, 0.3923f, 152.5f, 8.354f, 2.597f
574 },{ 10.99f, 0.4599f, 138.4f, 4.811f, 3.726f
575 },{ 16.6f, 0.3773f, 224.1f, 6.28f, 0.9121f
576 },{ 10.54f, 0.4533f, 159.3f, 4.832f, 2.529f
577 },{ 10.33f, 0.4502f, 162.0f, 5.132f, 2.444f
578 },{ 10.15f, 0.4471f, 165.6f, 5.378f, 2.328f
580 },{ 9.976f, 0.4439f, 168.0f, 5.721f, 2.258f
581 },{ 9.804f, 0.4408f, 176.2f, 5.675f, 1.997f
582 },{ 14.22f, 0.363f, 228.4f, 7.024f, 1.016f
583 },{ 9.952f, 0.4318f, 233.5f, 5.065f, 0.9244f
584 },{ 9.272f, 0.4345f, 210.0f, 4.911f, 1.258f
585 },{ 10.13f, 0.4146f, 225.7f, 5.525f, 1.055f
586 },{ 8.949f, 0.4304f, 213.3f, 5.071f, 1.221f
587 },{ 11.94f, 0.3783f, 247.2f, 6.655f, 0.849f
588 },{ 8.472f, 0.4405f, 195.5f, 4.051f, 1.604f
589 },{ 8.301f, 0.4399f, 203.7f, 3.667f, 1.459f
591 },{ 6.567f, 0.4858f, 193.0f, 2.65f, 1.66f
592 },{ 5.951f, 0.5016f, 196.1f, 2.662f, 1.589f
593 },{ 7.495f, 0.4523f, 251.4f, 3.433f, 0.8619f
594 },{ 6.335f, 0.4825f, 255.1f, 2.834f, 0.8228f
595 },{ 4.314f, 0.5558f, 214.8f, 2.354f, 1.263f
596 },{ 4.02f, 0.5681f, 219.9f, 2.402f, 1.191f
597 },{ 3.836f, 0.5765f, 210.2f, 2.742f, 1.305f
598 },{ 4.68f, 0.5247f, 244.7f, 2.749f, 0.8962f
599 },{ 2.892f, 0.6204f, 208.6f, 2.415f, 1.416f
601 },{ 2.892f, 0.6204f, 208.6f, 2.415f, 1.416f
603 },{ 4.728f, 0.5522f, 217.0f, 3.091f, 1.386f
604 },{ 6.18f, 0.52f, 170.0f, 4.0f, 3.224f
605 },{ 9.0f, 0.47f, 198.0f, 3.8f, 2.032f
606 },{ 2.324f, 0.6997f, 216.0f, 1.599f, 1.399f
607 },{ 1.961f, 0.7286f, 223.0f, 1.621f, 1.296f
608 },{ 1.75f, 0.7427f, 350.1f, 0.9789f, 0.5507f
609 },{ 10.31f, 0.4613f, 261.2f, 4.738f, 0.9899f
610 },{ 7.962f, 0.519f, 235.7f, 4.347f, 1.313f
611 },{ 6.227f, 0.5645f, 231.9f, 3.961f, 1.379f
612 },{ 5.246f, 0.5947f, 228.6f, 4.027f, 1.432f
614 },{ 5.408f, 0.5811f, 235.7f, 3.961f, 1.358f
615 },{ 5.218f, 0.5828f, 245.0f, 3.838f, 1.25f}
627 G4double shigh =
G4Log( 1.0 + x3*1000.0 + x4*0.001 )* x2*1000.0;
628 ionloss = slow*shigh*sqrt(T*1000.0) / (slow + shigh) ;
634 ionloss = slow*shigh / (slow + shigh) ;
666 const G4double* theAtomicNumDensityVector =
682 }
else if(1 == numberOfElements) {
694 for (
G4int i=0; i<numberOfElements; i++)
696 const G4Element* element = (*theElementVector)[i] ;
698 * theAtomicNumDensityVector[i];
714 static const G4double c[6] = {0.2865, 0.1266, -0.001429,
715 0.02402,-0.01135, 0.001475};
720 for (
G4int i=1; i<6; ++i) {
726 w = 1.0 + (0.007 + 0.00005*
z) *
G4Exp( -w*w ) ;
727 w = 4.0 * (1.0 -
G4Exp(-x)) * w * w ;
void set(double x, double y, double z)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
G4int SelectRandomAtomNumber(const G4Material *)
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
std::vector< ExP01TrackerHit * > a
G4double DEDX(const G4Material *material, G4double kineticEnergy)
static constexpr double MeV
G4double GetParticleCharge(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
void SetHighEnergyLimit(G4double)
G4double StoppingPower(const G4Material *material, G4double kineticEnergy)
G4double EffectiveChargeCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
const G4ParticleDefinition * particle
static constexpr double keV
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
static constexpr double twopi_mc2_rcl2
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy) override
Float_t x1[n_points_granero]
const G4ThreeVector & GetMomentumDirection() const
const G4ElementVector * GetElementVector() const
const G4String & GetParticleName() const
const G4String & GetParticleType() const
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length) override
virtual void flatArray(const int size, double *vect)=0
Hep3Vector & rotateUz(const Hep3Vector &)
void SetProposedMomentumDirection(const G4ThreeVector &dir)
G4double EffectiveChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double GetPDGMass() const
G4double G4Log(G4double x)
void SetParticle(const G4ParticleDefinition *p)
G4ParticleDefinition * theElectron
static constexpr double proton_mass_c2
G4double GetTotNbOfElectPerVolume() const
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy) override
G4ParticleChangeForLoss * GetParticleChangeForLoss()
void SetDeexcitationFlag(G4bool val)
static constexpr double amu_c2
G4EmCorrections * EmCorrections()
G4double GetElectronicDEDX(G4int idx, G4double energy) const
G4ParticleDefinition * GetDefinition() const
static constexpr double cm2
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
G4double theZieglerFactor
static G4ASTARStopping * fASTAR
G4VEmFluctuationModel * GetModelOfFluctuations()
static constexpr double electron_mass_c2
static constexpr double twopi
static constexpr double eV
void SetAngularDistribution(G4VEmAngularDistribution *)
G4double GetMeanExcitationEnergy() const
static G4Electron * Electron()
G4int GetIndex(const G4Material *) const
G4ThreeVector GetMomentum() const
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy) final
std::vector< G4Element * > G4ElementVector
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *couple) override
G4bool HasMaterial(const G4Material *material)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
virtual void SetParticleAndCharge(const G4ParticleDefinition *, G4double q2)
G4bool UseAngularGeneratorFlag() const
virtual ~G4BraggIonModel()
G4IonisParamMat * GetIonisation() const
G4double HeEffChargeSquare(G4double z, G4double kinEnergyInMeV) const
const G4Material * currentMaterial
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy) override
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
static G4LossTableManager * Instance()
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
G4VEmAngularDistribution * GetAngularDistribution()
const G4Material * GetMaterial() const
G4double GetTotalMomentum() const
const G4double * GetAtomicNumDensityVector() const
static constexpr double amu
virtual G4double GetParticleCharge(const G4ParticleDefinition *p, const G4Material *mat, G4double kineticEnergy) override
G4double GetTotNbOfAtomsPerVolume() const
Float_t x2[n_points_geant4]
G4double ElectronicStoppingPower(G4double z, G4double kineticEnergy) const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static constexpr double GeV
G4double GetElectronDensity() const
size_t GetNumberOfElements() const
G4ParticleChangeForLoss * fParticleChange
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4BraggIonModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="BraggIon")
G4double GetDensity() const