G4BohrFluctuations Class Reference

#include <G4BohrFluctuations.hh>

Inheritance diagram for G4BohrFluctuations:

Collaboration diagram for G4BohrFluctuations:

Public Member Functions
	G4BohrFluctuations (const G4String &nam="BohrFluc")
virtual	~G4BohrFluctuations ()
G4double	SampleFluctuations (const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double, G4double, G4double)
G4double	Dispersion (const G4Material *, const G4DynamicParticle *, G4double, G4double)
void	InitialiseMe (const G4ParticleDefinition *)
Public Member Functions inherited from G4VEmFluctuationModel
	G4VEmFluctuationModel (const G4String &nam)
virtual	~G4VEmFluctuationModel ()
virtual void	SetParticleAndCharge (const G4ParticleDefinition *, G4double q2)
const G4String &	GetName () const

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

Private Attributes
const G4ParticleDefinition *	particle
G4double	particleMass
G4double	chargeSquare
G4double	minNumberInteractionsBohr
G4double	minFraction
G4double	xmin
G4double	minLoss
G4double	kineticEnergy
G4double	beta2

Detailed Description

Definition at line 58 of file G4BohrFluctuations.hh.

Constructor & Destructor Documentation

G4BohrFluctuations() [1/2]

G4BohrFluctuations::G4BohrFluctuations

(

const G4String &

nam = "BohrFluc"

)

Definition at line 64 of file G4BohrFluctuations.cc.

 :G4VEmFluctuationModel(nam),
  particle(0),
  minNumberInteractionsBohr(2.0),
  minFraction(0.2),
  xmin(0.2),
  minLoss(0.001*eV)
{
  particleMass   = proton_mass_c2;
  chargeSquare   = 1.0;
  kineticEnergy  = 0.0;
  beta2          = 0.0;
}

~G4BohrFluctuations()

G4BohrFluctuations::~G4BohrFluctuations ( )

virtual

Definition at line 80 of file G4BohrFluctuations.cc.

{}

G4BohrFluctuations() [2/2]

G4BohrFluctuations::G4BohrFluctuations ( const G4BohrFluctuations &  )

private

Member Function Documentation

Dispersion()

G4double G4BohrFluctuations::Dispersion ( const G4Material *  material, const G4DynamicParticle *  dp, G4double  tmax, G4double  length  )

virtual

Implements G4VEmFluctuationModel.

Definition at line 150 of file G4BohrFluctuations.cc.

{
  if(!particle) { InitialiseMe(dp->GetDefinition()); }

  G4double electronDensity = material->GetElectronDensity();
  kineticEnergy = dp->GetKineticEnergy();
  G4double etot = kineticEnergy + particleMass;
  beta2 = kineticEnergy*(kineticEnergy + 2.0*particleMass)/(etot*etot);
  G4double siga  = (1.0/beta2 - 0.5) * twopi_mc2_rcl2 * tmax * length
                 * electronDensity * chargeSquare;

  return siga;
}

Here is the call graph for this function:

Here is the caller graph for this function:

InitialiseMe()

void G4BohrFluctuations::InitialiseMe ( const G4ParticleDefinition *  part )

virtual

Reimplemented from G4VEmFluctuationModel.

Definition at line 85 of file G4BohrFluctuations.cc.

{
  particle       = part;
  particleMass   = part->GetPDGMass();
  G4double q     = part->GetPDGCharge()/eplus;
  chargeSquare   = q*q;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

SampleFluctuations()

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

virtual

Implements G4VEmFluctuationModel.

Definition at line 96 of file G4BohrFluctuations.cc.

{
  if(meanLoss <= minLoss) { return meanLoss; }
  const G4Material* material = couple->GetMaterial();
  G4double siga = Dispersion(material,dp,tmax,length);
  G4double loss = meanLoss;

  G4double navr = meanLoss*meanLoss/siga;
  //G4cout << "### meanLoss= " << meanLoss << "  navr= " << navr << G4endl;
  if (navr >= minNumberInteractionsBohr) {
 
    // Increase fluctuations for big fractional energy loss
    if ( meanLoss > minFraction*kineticEnergy ) {
      G4double gam = (kineticEnergy - meanLoss)/particleMass + 1.0;
      G4double b2  = 1.0 - 1.0/(gam*gam);
      if(b2 < xmin*beta2) b2 = xmin*beta2;
      G4double x   = b2/beta2;
      G4double x3  = 1.0/(x*x*x);
      siga *= 0.25*(1.0 + x)*(x3 + (1.0/b2 - 0.5)/(1.0/beta2 - 0.5) );
    }
    siga = sqrt(siga);
    G4double twomeanLoss = meanLoss + meanLoss;
    //G4cout << "siga= " << siga << " 2edp= " << twomeanLoss <<G4endl;

    if(twomeanLoss < siga) {
      G4double x;
      do {
        loss = twomeanLoss*G4UniformRand();
        x = (loss - meanLoss)/siga;
        // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
      } while (1.0 - 0.5*x*x < G4UniformRand());
    } else {
      do {
        loss = G4RandGauss::shoot(meanLoss,siga);
        // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
      } while (0.0 > loss || loss > twomeanLoss);
    }

  // Poisson fluctuations
  } else {
    G4double n    = (G4double)(G4Poisson(navr));
    loss = meanLoss*n/navr;
  }
  //G4cout << "loss= " << loss << G4endl;

  return loss;
}

Here is the call graph for this function:

Member Data Documentation

beta2

G4double G4BohrFluctuations::beta2

private

Definition at line 97 of file G4BohrFluctuations.hh.

chargeSquare

G4double G4BohrFluctuations::chargeSquare

private

Definition at line 89 of file G4BohrFluctuations.hh.

kineticEnergy

G4double G4BohrFluctuations::kineticEnergy

private

Definition at line 96 of file G4BohrFluctuations.hh.

minFraction

G4double G4BohrFluctuations::minFraction

private

Definition at line 92 of file G4BohrFluctuations.hh.

minLoss

G4double G4BohrFluctuations::minLoss

private

Definition at line 94 of file G4BohrFluctuations.hh.

minNumberInteractionsBohr

G4double G4BohrFluctuations::minNumberInteractionsBohr

private

Definition at line 91 of file G4BohrFluctuations.hh.

particle

const G4ParticleDefinition* G4BohrFluctuations::particle

private

Definition at line 86 of file G4BohrFluctuations.hh.

particleMass

G4double G4BohrFluctuations::particleMass

private

Definition at line 88 of file G4BohrFluctuations.hh.

xmin

G4double G4BohrFluctuations::xmin

private

Definition at line 93 of file G4BohrFluctuations.hh.

---

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

• G4BohrFluctuations.hh
• G4BohrFluctuations.cc