G4AtimaFluctuations Class Reference

#include <G4AtimaFluctuations.hh>

Inheritance diagram for G4AtimaFluctuations:

Public Member Functions
	G4AtimaFluctuations (const G4String &nam="IonFlucAtima")
virtual	~G4AtimaFluctuations ()
virtual G4double	SampleFluctuations (const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss) override
virtual G4double	Dispersion (const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length) override
virtual void	InitialiseMe (const G4ParticleDefinition *) override
virtual void	SetParticleAndCharge (const G4ParticleDefinition *, G4double q2) override
Public Member Functions inherited from G4VEmFluctuationModel
	G4VEmFluctuationModel (const G4String &nam)
virtual	~G4VEmFluctuationModel ()
virtual G4double	SampleFluctuations (const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss)=0
virtual G4double	Dispersion (const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length)=0
virtual void	InitialiseMe (const G4ParticleDefinition *)
virtual void	SetParticleAndCharge (const G4ParticleDefinition *, G4double q2)
const G4String &	GetName () const
G4VEmFluctuationModel &	operator= (const G4VEmFluctuationModel &right)=delete
	G4VEmFluctuationModel (const G4VEmFluctuationModel &)=delete

Detailed Description

Definition at line 58 of file G4AtimaFluctuations.hh.

Constructor & Destructor Documentation

G4AtimaFluctuations()

G4AtimaFluctuations::G4AtimaFluctuations ( const G4String &  nam = "IonFlucAtima" )

explicit

Definition at line 74 of file G4AtimaFluctuations.cc.

  : G4VEmFluctuationModel(nam),
    particle(0),
    particleMass(CLHEP::proton_mass_c2),
    charge(1.0),
    chargeSquare(1.0),
    effChargeSquare(1.0),
    minLoss(0.001*CLHEP::eV)
{
  g4calc = G4Pow::GetInstance();
  kineticEnergy = 0.0;
  beta2 = 0.0;
  MLN10 = 2.30258509299;
  atomic_mass_unit = 931.4940954; // MeV/c^2
  dedx_constant = 0.3070749187;  //4*pi*Na*me*c^2*r_e^2  //MeV cm^2
  electron_mass = 0.510998928;   // MeV/c^2
  fine_structure = 1.0/137.035999139;
  domega2dx_constant = dedx_constant*electron_mass;  //4*pi*Na*me*c^2*r_e^2  //MeV^2 cm^2
  if(tableE[0] == 0.0) {
    G4double logmin = 0.;
    G4double logmax = 5.;
    stepE = (logmax-logmin)/(G4double)(199);
    for(G4int i=0; i<200; ++i){
      tableE[i] = G4Exp(MLN10*(logmin + ((G4double)i)*stepE));
    }
  }
}

~G4AtimaFluctuations()

G4AtimaFluctuations::~G4AtimaFluctuations ( )

virtual

Definition at line 104 of file G4AtimaFluctuations.cc.

{}

Member Function Documentation

Dispersion()

G4double G4AtimaFluctuations::Dispersion ( const G4Material *  mat, const G4DynamicParticle *  dp, G4double  tmax, G4double  length  )

overridevirtual

Implements G4VEmFluctuationModel.

Definition at line 146 of file G4AtimaFluctuations.cc.

{
  kineticEnergy = dp->GetKineticEnergy();
  const G4ParticleDefinition* p = dp->GetDefinition();
 
  G4double ap = p->GetPDGMass()/atomic_mass_unit;
  G4double zp = p->GetPDGCharge();
  G4double ep = kineticEnergy/ap;// ep in MeV/u
  G4double zt = mat->GetIonisation()->GetZeffective();
  G4double at = G4NistManager::Instance()->GetAtomicMassAmu(G4lrint(zt));
 
  G4int z = G4lrint(zp);
  if(z > 109) { z = 109; }
 
  G4double gamma=1.0 + ep/atomic_mass_unit;
  beta2=1.0-1.0/(gamma*gamma);
  G4double beta = std::sqrt(beta2);
  //z_eff_Pierce_Blann(Z, beta);
  G4double zp_eff = zp*(1.0-G4Exp(-0.95*137.035999139*beta/g4calc->Z23(z)));
  //
  G4double f = domega2dx_constant*zp_eff*zp_eff*zt/at;
  //
  G4double cor =
    24.89 * g4calc->powA(zt,1.2324)/(electron_mass*1e6 * beta2)*
    G4Log( 2.0*electron_mass*1e6*beta2/(33.05*g4calc->powA(zt,1.6364)));
  cor = std::max(cor, 0.0 );
 
  //Lindhard corrections
  if(ep<tableE[0])ep = tableE[0];
  G4double da = (ap - element_atomic_weights[z])/element_atomic_weights[z];
 
  G4double v3 = EnergyTable_interpolate(tableE,ep,ls_X_coefficients_a[z-1]);
  G4double v4 = EnergyTable_interpolate(tableE,ep,ls_X_coefficients_ahi[z-1]);
  G4double dif = v4 - v3;
  G4double X = v3+(dif*da/0.05);
  X *= gamma*gamma;
  //
  G4double sse = 0.;
  if(ep<30.0){
    //Energy straggling Firsov
    G4double factor = 4.8184e-3*g4calc->powA(zp+zt,8.0/3.0)/at;
    sse = std::min(f*(X+cor), factor*beta2/fine_structure/fine_structure);
  }else{
   // sse = f*(X+cor); 
    sse = f*X;//09/10/2018
  }
  //
  return sse*length/(cm)*mat->GetDensity()/(g/cm3);
}

Referenced by SampleFluctuations().

InitialiseMe()

void G4AtimaFluctuations::InitialiseMe ( const G4ParticleDefinition *  part )

overridevirtual

Reimplemented from G4VEmFluctuationModel.

Definition at line 109 of file G4AtimaFluctuations.cc.

{
  particle       = part;
  particleMass   = part->GetPDGMass();
  charge         = part->GetPDGCharge()/eplus;
  chargeSquare   = charge*charge;
  effChargeSquare= chargeSquare;
  uniFluct.InitialiseMe(part);
}

SampleFluctuations()

G4double G4AtimaFluctuations::SampleFluctuations ( const G4MaterialCutsCouple *  couple, const G4DynamicParticle *  dp, G4double  tmax, G4double  length, G4double  meanLoss  )

overridevirtual

Implements G4VEmFluctuationModel.

Definition at line 122 of file G4AtimaFluctuations.cc.

{
  //  G4cout << "### meanLoss= " << meanLoss << G4endl;
  if(meanLoss <= minLoss) return meanLoss;
 
 // G4cout << "G4AtimaFluctuations::SampleFluctuations E(MeV)= " 
 //    << dp->GetKineticEnergy()
 //    << "  Elim(MeV)= " << parameter*charge*particleMass 
 //    << " " << parameter << " " << charge << " " << particleMass << G4endl; 
 
  const G4Material* material = couple->GetMaterial();
  G4double siga = Dispersion(material,dp,tmax,length);
  
  CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
  //G4cout << "meanLoss= " << meanLoss << " loss= " << siga << G4endl;
  return G4RandGauss::shoot(rndmEngine,meanLoss,std::sqrt(siga));
}

SetParticleAndCharge()

void G4AtimaFluctuations::SetParticleAndCharge ( const G4ParticleDefinition *  part, G4double  q2  )

overridevirtual

Reimplemented from G4VEmFluctuationModel.

Definition at line 201 of file G4AtimaFluctuations.cc.

{
  if(part != particle) {
    particle       = part;
    particleMass   = part->GetPDGMass();
    charge         = part->GetPDGCharge()/eplus;
    chargeSquare   = charge*charge;
  }
  effChargeSquare  = q2;
  uniFluct.SetParticleAndCharge(part, q2);
}

---

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

• G4AtimaFluctuations.hh
• G4AtimaFluctuations.cc