G4ionEffectiveCharge Class Reference

#include <G4ionEffectiveCharge.hh>

Collaboration diagram for G4ionEffectiveCharge:

Public Member Functions
	G4ionEffectiveCharge ()
virtual	~G4ionEffectiveCharge ()
G4double	EffectiveChargeSquareRatio (const G4ParticleDefinition *p, const G4Material *material, G4double kineticEnergy)
G4double	EffectiveCharge (const G4ParticleDefinition *p, const G4Material *material, G4double kineticEnergy)

Private Member Functions
G4ionEffectiveCharge &	operator= (const G4ionEffectiveCharge &right)
	G4ionEffectiveCharge (const G4ionEffectiveCharge &)

Private Attributes
G4Pow *	g4pow
const G4ParticleDefinition *	lastPart
const G4Material *	lastMat
G4double	lastKinEnergy
G4double	chargeCorrection
G4double	effCharge
G4double	energyHighLimit
G4double	energyLowLimit
G4double	energyBohr
G4double	massFactor
G4double	minCharge

Static Private Attributes
static const G4double	inveplus = 1.0/CLHEP::eplus

Detailed Description

Definition at line 62 of file G4ionEffectiveCharge.hh.

Constructor & Destructor Documentation

G4ionEffectiveCharge() [1/2]

G4ionEffectiveCharge::G4ionEffectiveCharge

(

)

Definition at line 67 of file G4ionEffectiveCharge.cc.

{
  chargeCorrection = 1.0;
  energyHighLimit  = 20.0*MeV;
  energyLowLimit   = 1.0*keV;
  energyBohr       = 25.*keV;
  massFactor       = amu_c2/(proton_mass_c2*keV);
  minCharge        = 1.0;
  lastPart         = 0;
  lastMat          = 0;
  lastKinEnergy    = 0.0;
  effCharge        = eplus;
  g4pow = G4Pow::GetInstance();
}

Here is the call graph for this function:

~G4ionEffectiveCharge()

G4ionEffectiveCharge::~G4ionEffectiveCharge ( )

virtual

Definition at line 84 of file G4ionEffectiveCharge.cc.

{}

G4ionEffectiveCharge() [2/2]

G4ionEffectiveCharge::G4ionEffectiveCharge ( const G4ionEffectiveCharge &  )

private

Member Function Documentation

EffectiveCharge()

G4double G4ionEffectiveCharge::EffectiveCharge	(	const G4ParticleDefinition *	p,
		const G4Material *	material,
		G4double	kineticEnergy
	)

Definition at line 89 of file G4ionEffectiveCharge.cc.

{
  if(p == lastPart && material == lastMat && kineticEnergy == lastKinEnergy)
    return effCharge;

  lastPart      = p;
  lastMat       = material;
  lastKinEnergy = kineticEnergy;

  G4double mass   = p->GetPDGMass();
  G4double charge = p->GetPDGCharge();
  G4double Zi     = charge*inveplus;

  chargeCorrection = 1.0;
  effCharge = charge;

  // The aproximation of ion effective charge from:
  // J.F.Ziegler, J.P. Biersack, U. Littmark
  // The Stopping and Range of Ions in Matter,
  // Vol.1, Pergamon Press, 1985
  // Fast ions or hadrons
  G4double reducedEnergy = kineticEnergy * proton_mass_c2/mass ;

  //G4cout << "e= " << reducedEnergy << " Zi= " << Zi << "  " 
  //<< material->GetName() << G4endl;

  if(Zi < 1.5 || !material || reducedEnergy > Zi*energyHighLimit ) {
    return charge;
  }
  G4double z    = material->GetIonisation()->GetZeffective();
  reducedEnergy = std::max(reducedEnergy,energyLowLimit);

  // Helium ion case
  if( Zi < 2.5 ) {

    static const G4double c[6] = 
      {0.2865,0.1266,-0.001429,0.02402,-0.01135,0.001475};

    G4double Q = std::max(0.0,G4Log(reducedEnergy*massFactor));
    G4double x = c[0];
    G4double y = 1.0;
    for (G4int i=1; i<6; ++i) {
      y *= Q;
      x += y * c[i] ;
    }
    G4double ex;
    if(x < 0.2) { ex = x * (1 - 0.5*x); }
    else        { ex = 1. - G4Exp(-x); }

    G4double tq = 7.6 - Q;
    G4double tq2= tq*tq;
    G4double tt = ( 0.007 + 0.00005 * z );
    if(tq2 < 0.2) { tt *= (1.0 - tq2 + 0.5*tq2*tq2); }
    else          { tt *= G4Exp(-tq2); }

    effCharge = charge*(1.0 + tt) * std::sqrt(ex);

    // Heavy ion case
  } else {
    
    G4double y;
    G4double zi13 = g4pow->A13(Zi);
    G4double zi23 = zi13*zi13;

    // v1 is ion velocity in vF unit
    G4double eF   = material->GetIonisation()->GetFermiEnergy();
    G4double v1sq = reducedEnergy/eF;
    G4double vFsq = eF/energyBohr;
    G4double vF   = std::sqrt(eF/energyBohr);

    // Faster than Fermi velocity
    if ( v1sq > 1.0 ) {
      y = vF * std::sqrt(v1sq) * ( 1.0 + 0.2/v1sq ) / zi23 ;

      // Slower than Fermi velocity
    } else {
      y = 0.692308 * vF * (1.0 + 0.666666*v1sq + v1sq*v1sq/15.0) / zi23 ;
    }

    G4double q;
    G4double y3 = std::pow(y, 0.3) ;
    // G4cout<<"y= "<<y<<" y3= "<<y3<<" v1= "<<v1<<" vF= "<<vF<<G4endl; 
    q = 1.0 - G4Exp( 0.803*y3 - 1.3167*y3*y3 - 0.38157*y - 0.008983*y*y);   
    q = std::max(q,  minCharge/Zi); 
  
    effCharge = q*charge;
    
    G4double tq = 7.6 - G4Log(reducedEnergy/keV);
    G4double tq2= tq*tq;
    G4double sq = 1.0 + ( 0.18 + 0.0015 * z )*G4Exp(-tq2)/ (Zi*Zi);
 
    //    G4cout << "sq= " << sq << G4endl;

    // Screen length according to
    // J.F.Ziegler and J.M.Manoyan, The stopping of ions in compaunds,
    // Nucl. Inst. & Meth. in Phys. Res. B35 (1988) 215-228.

    G4double lambda = 10.0 * vF *g4pow->A23(1.0 - q)/ (zi13 * (6.0 + q));

    G4double lambda2 = lambda*lambda;

    G4double xx = (0.5/q - 0.5)*G4Log(1.0 + lambda2)/vFsq;

    chargeCorrection = sq * (1.0 + xx);
    
  }
  //  G4cout << "G4ionEffectiveCharge: charge= " << charge << " q= " << q 
  //         << " chargeCor= " << chargeCorrection 
  //           << " e(MeV)= " << kineticEnergy/MeV << G4endl;
  return effCharge;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EffectiveChargeSquareRatio()

G4double G4ionEffectiveCharge::EffectiveChargeSquareRatio ( const G4ParticleDefinition *  p, const G4Material *  material, G4double  kineticEnergy  )

inline

Definition at line 106 of file G4ionEffectiveCharge.hh.

{
  G4double charge = effCharge;
  if( kineticEnergy != lastKinEnergy || material != lastMat || p != lastPart) {
    charge = EffectiveCharge(p,material,kineticEnergy);
  }
  charge *= chargeCorrection*inveplus;

  return charge*charge;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

Member Data Documentation

chargeCorrection

G4double G4ionEffectiveCharge::chargeCorrection

private

Definition at line 93 of file G4ionEffectiveCharge.hh.

effCharge

G4double G4ionEffectiveCharge::effCharge

private

Definition at line 94 of file G4ionEffectiveCharge.hh.

energyBohr

G4double G4ionEffectiveCharge::energyBohr

private

Definition at line 98 of file G4ionEffectiveCharge.hh.

energyHighLimit

G4double G4ionEffectiveCharge::energyHighLimit

private

Definition at line 96 of file G4ionEffectiveCharge.hh.

energyLowLimit

G4double G4ionEffectiveCharge::energyLowLimit

private

Definition at line 97 of file G4ionEffectiveCharge.hh.

g4pow

G4Pow* G4ionEffectiveCharge::g4pow

private

Definition at line 87 of file G4ionEffectiveCharge.hh.

inveplus

const G4double G4ionEffectiveCharge::inveplus = 1.0/CLHEP::eplus

staticprivate

Definition at line 86 of file G4ionEffectiveCharge.hh.

lastKinEnergy

G4double G4ionEffectiveCharge::lastKinEnergy

private

Definition at line 91 of file G4ionEffectiveCharge.hh.

lastMat

const G4Material* G4ionEffectiveCharge::lastMat

private

Definition at line 90 of file G4ionEffectiveCharge.hh.

lastPart

const G4ParticleDefinition* G4ionEffectiveCharge::lastPart

private

Definition at line 89 of file G4ionEffectiveCharge.hh.

massFactor

G4double G4ionEffectiveCharge::massFactor

private

Definition at line 99 of file G4ionEffectiveCharge.hh.

minCharge

G4double G4ionEffectiveCharge::minCharge

private

Definition at line 100 of file G4ionEffectiveCharge.hh.

---

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

• G4ionEffectiveCharge.hh
• G4ionEffectiveCharge.cc