G4UrbanFluctuation Class Reference

#include <G4UrbanFluctuation.hh>

Inheritance diagram for G4UrbanFluctuation:

Public Member Functions
	G4UrbanFluctuation (const G4String &nam="UrbanFluc")
	~G4UrbanFluctuation () override
G4UrbanFluctuation &	operator= (const G4UrbanFluctuation &right)=delete
	G4UrbanFluctuation (const G4UrbanFluctuation &)=delete
Public Member Functions inherited from G4UniversalFluctuation
	G4UniversalFluctuation (const G4String &nam="UniFluc")
	~G4UniversalFluctuation () override
G4double	SampleFluctuations (const G4MaterialCutsCouple *, const G4DynamicParticle *, const G4double, const G4double, const G4double, const G4double) override
G4double	Dispersion (const G4Material *, const G4DynamicParticle *, const G4double, const G4double, const G4double) override
void	InitialiseMe (const G4ParticleDefinition *) override
void	SetParticleAndCharge (const G4ParticleDefinition *, G4double q2) override
G4UniversalFluctuation &	operator= (const G4UniversalFluctuation &right)=delete
	G4UniversalFluctuation (const G4UniversalFluctuation &)=delete
Public Member Functions inherited from G4VEmFluctuationModel
	G4VEmFluctuationModel (const G4String &nam)
virtual	~G4VEmFluctuationModel ()
virtual G4double	SampleFluctuations (const G4MaterialCutsCouple *, const G4DynamicParticle *, const G4double tcut, const G4double tmax, const G4double length, const G4double meanLoss)=0
virtual G4double	Dispersion (const G4Material *, const G4DynamicParticle *, const G4double tcut, const G4double tmax, const 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

Protected Member Functions
G4double	SampleGlandz (CLHEP::HepRandomEngine *rndm, const G4Material *, const G4double tcut) override
Protected Member Functions inherited from G4UniversalFluctuation
virtual G4double	SampleGlandz (CLHEP::HepRandomEngine *rndm, const G4Material *, const G4double tcut)
void	AddExcitation (CLHEP::HepRandomEngine *rndm, const G4double ax, const G4double ex, G4double &eav, G4double &eloss, G4double &esig2)
void	SampleGauss (CLHEP::HepRandomEngine *rndm, const G4double eav, const G4double esig2, G4double &eloss)

Additional Inherited Members
Protected Attributes inherited from G4UniversalFluctuation
G4double	particleMass = 0.0
G4double	m_Inv_particleMass = DBL_MAX
G4double	m_massrate = DBL_MAX
G4double	chargeSquare = 1.0
G4double	ipotFluct = 0.0
G4double	ipotLogFluct = 0.0
G4double	e0 = 0.0
G4double	minNumberInteractionsBohr = 10.0
G4double	minLoss
G4double	nmaxCont = 8.0
G4double	rate = 0.56
G4double	fw = 4.0
G4double	a0 = 42.0
G4double	w2 = 0.0
G4double	meanLoss = 0.0
const G4ParticleDefinition *	particle = nullptr
G4double *	rndmarray = nullptr
G4int	sizearray = 30

Detailed Description

Definition at line 54 of file G4UrbanFluctuation.hh.

Constructor & Destructor Documentation

G4UrbanFluctuation() [1/2]

G4UrbanFluctuation::G4UrbanFluctuation ( const G4String &  nam = "UrbanFluc" )

explicit

Definition at line 55 of file G4UrbanFluctuation.cc.

 : G4UniversalFluctuation(nam)
{}

~G4UrbanFluctuation()

G4UrbanFluctuation::~G4UrbanFluctuation ( )

overridedefault

G4UrbanFluctuation() [2/2]

G4UrbanFluctuation::G4UrbanFluctuation ( const G4UrbanFluctuation &  )

delete

Member Function Documentation

operator=()

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

delete

SampleGlandz()

G4double G4UrbanFluctuation::SampleGlandz ( CLHEP::HepRandomEngine *  rndm, const G4Material *  material, const G4double  tcut  )

overrideprotectedvirtual

Reimplemented from G4UniversalFluctuation.

Definition at line 65 of file G4UrbanFluctuation.cc.

{
  if (material != lastMaterial) {
    auto ioni = material->GetIonisation();
    f1Fluct = ioni->GetF1fluct();
    f2Fluct = ioni->GetF2fluct();
    e1Fluct = ioni->GetEnergy1fluct();
    e2Fluct = ioni->GetEnergy2fluct();
    e1LogFluct = ioni->GetLogEnergy1fluct();
    e2LogFluct = ioni->GetLogEnergy2fluct();
    esmall = 0.5*std::sqrt(e0*ipotFluct);  
    lastMaterial = material;   
  }
 
  G4double a1(0.0), a2(0.0), a3(0.0);
  G4double loss = 0.0;
  G4double e1 = e1Fluct;
  G4double e2 = e2Fluct;
 
  if(tcut > ipotFluct) {
    if(w2 > ipotLogFluct)  {
      if(w2 > e2LogFluct) {
    const G4double C = meanLoss*(1.-rate)/(w2-ipotLogFluct);
    a1 = C*f1Fluct*(w2-e1LogFluct)/e1Fluct;
    a2 = C*f2Fluct*(w2-e2LogFluct)/e2Fluct;
      } else {
    a1 = meanLoss*(1.-rate)/e1;
      }
      if(a1 < a0) {
        const G4double fwnow = 0.5+(fw-0.5)*std::sqrt(a1/a0);
        a1 /= fwnow;
        e1 *= fwnow;
      } else {
        a1 /= fw;
        e1 *= fw;
      }
    }   
  }
 
  const G4double w1 = tcut/e0;
  a3 = rate*meanLoss*(tcut-e0)/(e0*tcut*G4Log(w1));
  if(a1+a2 <= 0.) { a3 /= rate; }
 
  //'nearly' Gaussian fluctuation if a1>nmaxCont&&a2>nmaxCont&&a3>nmaxCont  
  G4double emean = 0.;
  G4double sig2e = 0.;
 
  // excitation of type 1
  if(a1 > 0.0) { AddExcitation(rndmEngineF, a1, e1, emean, loss, sig2e); }
 
  // excitation of type 2
  if(a2 > 0.0) { AddExcitation(rndmEngineF, a2, e2, emean, loss, sig2e); }
 
  if(sig2e > 0.0) { SampleGauss(rndmEngineF, emean, sig2e, loss); }
 
  // ionisation 
  if(a3 > 0.) {
    emean = 0.;
    sig2e = 0.;
    G4double p3 = a3;
    G4double alfa = 1.;
    if(a3 > nmaxCont) {
      alfa = w1*(nmaxCont+a3)/(w1*nmaxCont+a3);
      const G4double alfa1  = alfa*G4Log(alfa)/(alfa-1.);
      const G4double namean = a3*w1*(alfa-1.)/((w1-1.)*alfa);
      emean += namean*e0*alfa1;
      sig2e += e0*e0*namean*(alfa-alfa1*alfa1);
      p3 -= namean;
    }
 
    const G4double w3 = alfa*e0;
    if(tcut > w3) {
      const G4double w = (tcut-w3)/tcut;
      const G4int nnb = (G4int)G4Poisson(p3);
      if(nnb > 0) {
        if(nnb > sizearray) {
          sizearray = nnb;
          delete [] rndmarray;
          rndmarray = new G4double[nnb];
        }
        rndmEngineF->flatArray(nnb, rndmarray);
        for (G4int k=0; k<nnb; ++k) { loss += w3/(1.-w*rndmarray[k]); }
      }
    }
    if(sig2e > 0.0) { SampleGauss(rndmEngineF, emean, sig2e, loss); }
  }
  //G4cout << "### loss=" << loss << G4endl;
  return loss;
}

---

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

• G4UrbanFluctuation.hh
• G4UrbanFluctuation.cc