#include <G4AtimaFluctuations.hh>

Inheritance diagram for G4AtimaFluctuations:

Public Member Functions
G4double	Dispersion (const G4Material , const G4DynamicParticle , const G4double tcut, const G4double tmax, const G4double length) override

	G4AtimaFluctuations (const G4AtimaFluctuations &)=delete

	G4AtimaFluctuations (const G4String &nam="IonFlucAtima")

const G4String &	GetName () const

void	InitialiseMe (const G4ParticleDefinition *) override

G4AtimaFluctuations &	operator= (const G4AtimaFluctuations &right)=delete

G4double	SampleFluctuations (const G4MaterialCutsCouple , const G4DynamicParticle , const G4double tcut, const G4double tmax, const G4double length, const G4double meanLoss) override

void	SetParticleAndCharge (const G4ParticleDefinition *, G4double q2) override

	~G4AtimaFluctuations () override

Private Member Functions
G4double	EnergyTable_interpolate (const G4double table, G4double xval, const G4double y)

Private Attributes
G4double	atomic_mass_unit

G4double	beta2

G4double	charge

G4double	chargeSquare

G4double	dedx_constant

G4double	domega2dx_constant

G4double	effChargeSquare

G4double	electron_mass

G4double	fine_structure

G4LossTableManager *	fManager

G4Pow *	g4calc

G4double	kineticEnergy

G4double	minLoss

G4double	MLN10

const G4String	name

const G4ParticleDefinition *	particle

G4double	particleMass

Static Private Attributes
static const G4double	element_atomic_weights [110]

static const G4double	ls_X_coefficients_a [110][200]

static const G4double	ls_X_coefficients_ahi [110][200]

static G4double	stepE = 0.0

static G4double	tableE [200] = {0.0}

Detailed Description

Definition at line 58 of file G4AtimaFluctuations.hh.

Constructor & Destructor Documentation

◆ G4AtimaFluctuations() [1/2]

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

explicit

Definition at line 74 of file G4AtimaFluctuations.cc.

  : G4VEmFluctuationModel(nam),
    particle(nullptr),
    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));
    }
  }
}

References atomic_mass_unit, beta2, dedx_constant, domega2dx_constant, electron_mass, fine_structure, g4calc, G4Exp(), G4Pow::GetInstance(), kineticEnergy, MLN10, stepE, and tableE.

◆ ~G4AtimaFluctuations()

G4AtimaFluctuations::~G4AtimaFluctuations ( )

override

Definition at line 104 of file G4AtimaFluctuations.cc.

105{}

◆ G4AtimaFluctuations() [2/2]

G4AtimaFluctuations::G4AtimaFluctuations ( const G4AtimaFluctuations & )

delete

Member Function Documentation

◆ Dispersion()

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

overridevirtual

Implements G4VEmFluctuationModel.

Definition at line 146 of file G4AtimaFluctuations.cc.

{
  kineticEnergy = dp->GetKineticEnergy();
  const G4ParticleDefinition* p = dp->GetDefinition();
  if(particle != p) { InitialiseMe(p); }
 
  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);
}

References G4InuclParticleNames::ap, atomic_mass_unit, anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, beta2, cm, cm3, domega2dx_constant, e6, electron_mass, element_atomic_weights, EnergyTable_interpolate(), fine_structure, g, g4calc, G4Exp(), G4Log(), G4lrint(), G4NistManager::GetAtomicMassAmu(), G4DynamicParticle::GetDefinition(), G4Material::GetDensity(), G4Material::GetIonisation(), G4DynamicParticle::GetKineticEnergy(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4IonisParamMat::GetZeffective(), InitialiseMe(), G4NistManager::Instance(), kineticEnergy, ls_X_coefficients_a, ls_X_coefficients_ahi, G4INCL::Math::max(), G4INCL::Math::min(), particle, G4Pow::powA(), tableE, and G4Pow::Z23().

Referenced by SampleFluctuations().

◆ EnergyTable_interpolate()

G4double G4AtimaFluctuations::EnergyTable_interpolate	(	const G4double *	table,
		G4double	xval,
		const G4double *	y
	)

private

Definition at line 216 of file G4AtimaFluctuations.cc.

                                                                                                            {
  G4double r;
  G4int num=200;
  G4double lxval = G4Log(xval)/MLN10;
  if(xval<table[0] || xval>table[num-1])return 0.0;
  if(xval==table[num-1])return y[num-1];
  G4int i = (G4int)(lxval/stepE);
  i = std::min(std::max(i, 0), num-2);
  G4double linstep = table[i+1] - table[i];
  G4double x = 1.0 - ((xval - table[i])/linstep);
  r = (x*y[i]) + ((1-x)*y[i+1]);
  return r;
}

References G4Log(), G4INCL::Math::max(), G4INCL::Math::min(), MLN10, and stepE.

Referenced by Dispersion().

◆ GetName()

const G4String & G4VEmFluctuationModel::GetName ( ) const

inlineinherited

Definition at line 121 of file G4VEmFluctuationModel.hh.

{
  return name;
}

References G4VEmFluctuationModel::name.

Referenced by G4LossTableManager::Register(), and G4EmConfigurator::SetModelForRegion().

◆ 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;
}

References charge, chargeSquare, effChargeSquare, eplus, G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), particle, and particleMass.

Referenced by Dispersion().

◆ operator=()

G4AtimaFluctuations & G4AtimaFluctuations::operator= ( const G4AtimaFluctuations & right )

delete

◆ SampleFluctuations()

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

overridevirtual

Implements G4VEmFluctuationModel.

Definition at line 121 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,tcut,tmax,length);
  
  CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
  //G4cout << "meanLoss= " << meanLoss << " loss= " << siga << G4endl;
  return G4RandGauss::shoot(rndmEngine,meanLoss,std::sqrt(siga));
}

References Dispersion(), G4MaterialCutsCouple::GetMaterial(), eplot::material, minLoss, and G4INCL::DeJongSpin::shoot().

◆ SetParticleAndCharge()

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

overridevirtual

Reimplemented from G4VEmFluctuationModel.

Definition at line 202 of file G4AtimaFluctuations.cc.

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

References charge, chargeSquare, effChargeSquare, eplus, G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), particle, and particleMass.

Field Documentation

◆ atomic_mass_unit

G4double G4AtimaFluctuations::atomic_mass_unit

private

Definition at line 107 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ beta2

G4double G4AtimaFluctuations::beta2

private

Definition at line 123 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ charge

G4double G4AtimaFluctuations::charge

private

Definition at line 102 of file G4AtimaFluctuations.hh.

Referenced by InitialiseMe(), and SetParticleAndCharge().

◆ chargeSquare

G4double G4AtimaFluctuations::chargeSquare

private

Definition at line 103 of file G4AtimaFluctuations.hh.

Referenced by InitialiseMe(), and SetParticleAndCharge().

◆ dedx_constant

G4double G4AtimaFluctuations::dedx_constant

private

Definition at line 108 of file G4AtimaFluctuations.hh.

Referenced by G4AtimaFluctuations().

◆ domega2dx_constant

G4double G4AtimaFluctuations::domega2dx_constant

private

Definition at line 111 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ effChargeSquare

G4double G4AtimaFluctuations::effChargeSquare

private

Definition at line 104 of file G4AtimaFluctuations.hh.

Referenced by InitialiseMe(), and SetParticleAndCharge().

◆ electron_mass

G4double G4AtimaFluctuations::electron_mass

private

Definition at line 109 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ element_atomic_weights

const G4double G4AtimaFluctuations::element_atomic_weights

staticprivate

Definition at line 117 of file G4AtimaFluctuations.hh.

Referenced by Dispersion().

◆ fine_structure

G4double G4AtimaFluctuations::fine_structure

private

Definition at line 110 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ fManager

G4LossTableManager* G4VEmFluctuationModel::fManager

privateinherited

Definition at line 115 of file G4VEmFluctuationModel.hh.

Referenced by G4VEmFluctuationModel::G4VEmFluctuationModel(), and G4VEmFluctuationModel::~G4VEmFluctuationModel().

◆ g4calc

G4Pow* G4AtimaFluctuations::g4calc

private

Definition at line 99 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ kineticEnergy

G4double G4AtimaFluctuations::kineticEnergy

private

Definition at line 122 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

◆ ls_X_coefficients_a

const G4double G4AtimaFluctuations::ls_X_coefficients_a

staticprivate

Definition at line 115 of file G4AtimaFluctuations.hh.

Referenced by Dispersion().

◆ ls_X_coefficients_ahi

const G4double G4AtimaFluctuations::ls_X_coefficients_ahi

staticprivate

Definition at line 116 of file G4AtimaFluctuations.hh.

Referenced by Dispersion().

◆ minLoss

G4double G4AtimaFluctuations::minLoss

private

Definition at line 120 of file G4AtimaFluctuations.hh.

Referenced by SampleFluctuations().

◆ MLN10

G4double G4AtimaFluctuations::MLN10

private

Definition at line 106 of file G4AtimaFluctuations.hh.

Referenced by EnergyTable_interpolate(), and G4AtimaFluctuations().

◆ name

const G4String G4VEmFluctuationModel::name

privateinherited

Definition at line 114 of file G4VEmFluctuationModel.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), G4VEmFluctuationModel::GetName(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

◆ particle

const G4ParticleDefinition* G4AtimaFluctuations::particle

private

Definition at line 97 of file G4AtimaFluctuations.hh.

Referenced by emcalc_gui.MainWindow::cb_select_particle(), emcalc_gui.MainWindow::cb_show_plot(), Dispersion(), InitialiseMe(), and SetParticleAndCharge().

◆ particleMass

G4double G4AtimaFluctuations::particleMass

private

Definition at line 101 of file G4AtimaFluctuations.hh.

Referenced by InitialiseMe(), and SetParticleAndCharge().

◆ stepE

G4double G4AtimaFluctuations::stepE = 0.0

staticprivate

Definition at line 113 of file G4AtimaFluctuations.hh.

Referenced by EnergyTable_interpolate(), and G4AtimaFluctuations().

◆ tableE

G4double G4AtimaFluctuations::tableE = {0.0}

staticprivate

Definition at line 114 of file G4AtimaFluctuations.hh.

Referenced by Dispersion(), and G4AtimaFluctuations().

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

source/processes/electromagnetic/standard/include/G4AtimaFluctuations.hh
source/processes/electromagnetic/standard/src/G4AtimaFluctuations.cc

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4AtimaFluctuations() [1/2]

◆ ~G4AtimaFluctuations()

◆ G4AtimaFluctuations() [2/2]

Member Function Documentation

◆ Dispersion()

◆ EnergyTable_interpolate()

◆ GetName()

◆ InitialiseMe()

◆ operator=()

◆ SampleFluctuations()

◆ SetParticleAndCharge()

Field Documentation

◆ atomic_mass_unit

◆ beta2

◆ charge

◆ chargeSquare

◆ dedx_constant

◆ domega2dx_constant

◆ effChargeSquare

◆ electron_mass

◆ element_atomic_weights

◆ fine_structure

◆ fManager

◆ g4calc

◆ kineticEnergy

◆ ls_X_coefficients_a

◆ ls_X_coefficients_ahi

◆ minLoss

◆ MLN10

◆ name

◆ particle

◆ particleMass

◆ stepE

◆ tableE