#include <GFlashSamplingShowerParameterisation.hh>

Inheritance diagram for GFlashSamplingShowerParameterisation:

Public Member Functions
G4double	ApplySampling (const G4double DEne, const G4double Energy)

void	ComputeRadialParameters (G4double y, G4double Tau)

G4double	ComputeTau (G4double LongitudinalPosition)

void	ComputeZAX0EFFetc ()

G4double	gam (G4double x, G4double a) const

G4double	GenerateExponential (G4double Energy)

void	GenerateLongitudinalProfile (G4double Energy)

G4double	GeneratePhi ()

G4double	GenerateRadius (G4int ispot, G4double Energy, G4double LongitudinalPosition)

G4double	GetAveR90 ()

G4double	GetAveR99 ()

G4double	GetAveT90 ()

G4double	GetAveT99 ()

G4double	GetAveTmx ()

G4double	GetEc ()

G4double	GetEffA (const G4Material *material)

G4double	GetEffZ (const G4Material *material)

G4double	GetNspot ()

G4double	GetRm ()

G4double	GetX0 ()

	GFlashSamplingShowerParameterisation (G4Material aMat1, G4Material aMat2, G4double d1, G4double d2, GFlashSamplingShowerTuning *aPar=0)

G4double	IntegrateEneLongitudinal (G4double LongitudinalStep)

G4double	IntegrateNspLongitudinal (G4double LongitudinalStep)

void	PrintMaterial (const G4Material *mat)

void	SetMaterial (G4Material mat1, G4Material mat2)

	~GFlashSamplingShowerParameterisation ()

Protected Attributes
G4double	A

G4double	density

G4double	Ec

G4double	Rm

G4double	X0

G4double	Z

Private Member Functions
void	ComputeLongitudinalParameters (G4double y)

void	GenerateEnergyProfile (G4double y)

void	GenerateNSpotProfile (G4double y)

Private Attributes
G4double	A1

G4double	A2

G4double	Aeff

G4double	Alpha

G4double	Alphah

G4double	AlphaNSpot

G4double	AveLogAlpha

G4double	AveLogAlphah

G4double	AveLogTmax

G4double	AveLogTmaxh

G4double	Beta

G4double	Betah

G4double	BetaNSpot

G4double	ConstantResolution

G4double	d1

G4double	d2

G4double	density1

G4double	density2

G4double	Ec1

G4double	Ec2

G4double	Eceff

G4double	ehat

MyGamma *	fGamma

G4double	Fs

G4Material *	material1

G4Material *	material2

G4double	NoiseResolution

G4double	NSpot

G4double	ParAveA1

G4double	ParAveA2

G4double	ParAveA3

G4double	ParAveT1

G4double	ParAveT2

G4double	ParRC1

G4double	ParRC2

G4double	ParRC3

G4double	ParRC4

G4double	ParRho1

G4double	ParRho2

G4double	ParRT1

G4double	ParRT2

G4double	ParRT3

G4double	ParRT4

G4double	ParRT5

G4double	ParRT6

G4double	ParsAveA1

G4double	ParsAveA2

G4double	ParsAveT1

G4double	ParsAveT2

G4double	ParSigLogA1

G4double	ParSigLogA2

G4double	ParSigLogT1

G4double	ParSigLogT2

G4double	ParsRC1

G4double	ParsRC2

G4double	ParsRho1

G4double	ParsRho2

G4double	ParsRT1

G4double	ParsRT2

G4double	ParsSigLogA1

G4double	ParsSigLogA2

G4double	ParsSigLogT1

G4double	ParsSigLogT2

G4double	ParsSpotA1

G4double	ParsSpotA2

G4double	ParsSpotN1

G4double	ParsSpotN2

G4double	ParsSpotT1

G4double	ParsSpotT2

G4double	ParsWC1

G4double	ParsWC2

G4double	ParWC1

G4double	ParWC2

G4double	ParWC3

G4double	ParWC4

G4double	ParWC5

G4double	ParWC6

G4double	RadiusCore

G4double	RadiusTail

G4double	Rho

G4double	Rhoeff

G4double	Rhoh

G4double	Rm1

G4double	Rm2

G4double	Rmeff

G4double	SamplingResolution

G4double	SigmaLogAlpha

G4double	SigmaLogAlphah

G4double	SigmaLogTmax

G4double	SigmaLogTmaxh

GFlashSamplingShowerTuning *	thePar

G4double	Tmax

G4double	Tmaxh

G4double	TNSpot

G4double	WeightCore

G4double	X01

G4double	X02

G4double	X0eff

G4double	Z1

G4double	Z2

G4double	Zeff

Detailed Description

Definition at line 49 of file GFlashSamplingShowerParameterisation.hh.

Constructor & Destructor Documentation

◆ GFlashSamplingShowerParameterisation()

GFlashSamplingShowerParameterisation::GFlashSamplingShowerParameterisation	(	G4Material *	aMat1,
		G4Material *	aMat2,
		G4double	d1,
		G4double	d2,
		GFlashSamplingShowerTuning *	aPar = `0`
	)

Parameters

aPar
aMat1	passive material,
dd1	- passive layer thickness
aMat2	active material,
dd2	- acive layer thickness

Definition at line 46 of file GFlashSamplingShowerParameterisation.cc.

  : GVFlashShowerParameterisation(),
    ParAveT2(0.), ParSigLogT1(0.), ParSigLogT2(0.),
    ParSigLogA1(0.), ParSigLogA2(0.), ParRho1(0.), ParRho2(0.), ParsAveA2(0.),
    AveLogAlphah(0.), AveLogTmaxh(0.), SigmaLogAlphah(0.), SigmaLogTmaxh(0.),
    Rhoh(0.), Alphah(0.), Tmaxh(0.), Betah(0.), AveLogAlpha(0.), AveLogTmax(0.),
    SigmaLogAlpha(0.), SigmaLogTmax(0.), Rho(0.), Alpha(0.), Tmax(0.), Beta(0.)
{  
  if(!aPar) {
    thePar = new GFlashSamplingShowerTuning;
  } else {
    thePar = aPar;
  }
 
  SetMaterial(aMat1,aMat2 );
  d1=dd1;
  d2=dd2;
 
  // Longitudinal Coefficients for a homogenious calo
 
  // shower max
  ParAveT1    = thePar->ParAveT1();   // ln (ln y -0.812)
  ParAveA1    = thePar->ParAveA1();   // ln a (0.81 + (0.458 + 2.26/Z)ln y)
  ParAveA2    = thePar->ParAveA2();
  ParAveA3    = thePar->ParAveA3();
  // Variance of shower max sampling
  ParSigLogT1 = thePar->ParsSigLogT1();  // Sigma T1 (-1.4 + 1.26 ln y)**-1 --> bug : these two lines were missing,
  ParSigLogT2 = thePar->ParsSigLogT2();  // leaving ParSigLogT1, ParSigLogT2 as 0.0
  // variance of 'alpha'
  ParSigLogA1 = thePar->ParSigLogA1();   // Sigma a (-0.58 + 0.86 ln y)**-1 --> bug : these two lines were missing
  ParSigLogA2 = thePar->ParSigLogA2();   // leaving ParSigLogA1 ParSigLogAé as 0.0
  // correlation alpha%T
  ParRho1     = thePar->ParRho1();       // Rho = 0.705 -0.023 ln y  --> bug : these two lines were missing,
  ParRho2     = thePar->ParRho2();       // leaving ParRho1 and ParRho2 being 0.0
  // Sampling
  ParsAveT1   = thePar->ParsAveT1();  // T_sam = log(exp( log T_hom) + t1*Fs-1 + t2*(1-ehat));
  ParsAveT2   = thePar->ParsAveT2();
  ParsAveA1   = thePar->ParsAveA1();  
  // Variance of shower max sampling
  ParsSigLogT1 = thePar->ParsSigLogT1();    // Sigma T1 (-2.5 + 1.25 ln y)**-1 --> bug ParSigLogT1() was called instead of ParsSigLogT1(); Same for T2.
  ParsSigLogT2 = thePar->ParsSigLogT2();
  // variance of 'alpha'
  ParsSigLogA1 = thePar->ParsSigLogA1();    // Sigma a (-0.82 + 0.79 ln y)**-1 --> bug ParSigLogA1() was called instead of ParsSigLogA1(); Same for A2
  ParsSigLogA2 = thePar->ParsSigLogA2();
  // correlation alpha%T
  ParsRho1     = thePar->ParsRho1();   // Rho = 0.784 -0.023 ln y --> bug was using ParRho1() and ParRho2()
  ParsRho2     = thePar->ParsRho2();
 
  // Radial Coefficients
  // r_C (tau)= z_1 +z_2 tau
  // r_t (tau)= k1 (std::exp (k3(tau -k2 ))+std::exp (k_4 (tau- k_2))))
  ParRC1 =   thePar->ParRC1();  // z_1 = 0.0251 + 0.00319 ln E
  ParRC2 =   thePar->ParRC2();
  ParRC3 =   thePar->ParRC3();  // z_2 = 0.1162 + - 0.000381 Z
  ParRC4 =   thePar->ParRC4();
 
  ParWC1 = thePar->ParWC1();
  ParWC2 = thePar->ParWC2();
  ParWC3 = thePar->ParWC3();
  ParWC4 = thePar->ParWC4();
  ParWC5 = thePar->ParWC5(); 
  ParWC6 = thePar->ParWC6();
  ParRT1 = thePar->ParRT1();
  ParRT2 = thePar->ParRT2();
  ParRT3 = thePar->ParRT3();
  ParRT4 = thePar->ParRT4(); 
  ParRT5 = thePar->ParRT5();
  ParRT6 = thePar->ParRT6();
 
  //additional sampling parameter
  ParsRC1= thePar->ParsRC1();
  ParsRC2= thePar->ParsRC2();
  ParsWC1= thePar->ParsWC1();
  ParsWC2=  thePar->ParsWC2(); 
  ParsRT1= thePar->ParsRT1();
  ParsRT2= thePar->ParsRT2();
 
  // Coeff for fluctuedted radial  profiles for a sampling media
  ParsSpotT1   = thePar->ParSpotT1();     // T_spot = T_hom =(0.698 + 0.00212)
  ParsSpotT2   = thePar->ParSpotT2();
  ParsSpotA1   = thePar->ParSpotA1();     // a_spot= a_hom (0.639 + 0.00334)
  ParsSpotA2   = thePar->ParSpotA2();    
  ParsSpotN1   = thePar->ParSpotN1();     // N_Spot 93 * ln(Z) E ** 0.876   
  ParsSpotN2   = thePar->ParSpotN2(); 
  SamplingResolution  = thePar->SamplingResolution();
  ConstantResolution  = thePar->ConstantResolution(); 
  NoiseResolution     = thePar->NoiseResolution();
 
  // Inits
  NSpot         = 0.00;
  AlphaNSpot    = 0.00;
  TNSpot        = 0.00;
  BetaNSpot     = 0.00;
  RadiusCore    = 0.00;
  WeightCore    = 0.00;
  RadiusTail    = 0.00;   
  ComputeZAX0EFFetc();
 
  G4cout << "/********************************************/ " << G4endl;
  G4cout << "  - GFlashSamplingShowerParameterisation::Constructor -  " << G4endl;
  G4cout << "/********************************************/ " << G4endl;  
}

◆ ~GFlashSamplingShowerParameterisation()

GFlashSamplingShowerParameterisation::~GFlashSamplingShowerParameterisation ( )

Definition at line 154 of file GFlashSamplingShowerParameterisation.cc.

{
   delete thePar;
}

References thePar.

Member Function Documentation

◆ ApplySampling()

G4double GFlashSamplingShowerParameterisation::ApplySampling	(	const G4double	DEne,
		const G4double	Energy
	)

Definition at line 322 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double DEneFluctuated = DEne;
  G4double Resolution     = std::pow(SamplingResolution,2);
 
  //       +pow(NoiseResolution,2)/  //@@@@@@@@ FIXME 
  //                         Energy*(1.*MeV)+
  //                         pow(ConstantResolution,2)*
  //                          Energy/(1.*MeV);
 
  if(Resolution >0.0 && DEne > 0.00)
  {
    G4float x1=DEne/Resolution;
    G4float x2 = G4RandGamma::shoot(x1, 1.0)*Resolution;     
    DEneFluctuated=x2;
  }
  return DEneFluctuated;
}

References SamplingResolution, and G4INCL::DeJongSpin::shoot().

◆ ComputeLongitudinalParameters()

void GFlashSamplingShowerParameterisation::ComputeLongitudinalParameters ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 246 of file GFlashSamplingShowerParameterisation.cc.

{
  AveLogTmaxh    = std::log(  std::max(  ParAveT1 + std::log(y),                           0.1  )  );        // ok
  AveLogAlphah   = std::log(  std::max(  ParAveA1 + (ParAveA2+ParAveA3/Zeff)*std::log(y),  0.1  )  );        // ok
  // hom
  SigmaLogTmaxh  = std::min(  0.5,  1.00/( ParSigLogT1 + ParSigLogT2*std::log(y) )  );                       // ok
  SigmaLogAlphah = std::min(  0.5,  1.00/( ParSigLogA1 + ParSigLogA2*std::log(y) )  );                       // ok
  Rhoh           = ParRho1 + ParRho2*std::log(y);                                                             //ok
  // if sampling 
  AveLogTmax    = std::max(  0.1,  std::log(std::exp(AveLogTmaxh)  + ParsAveT1/Fs + ParsAveT2*(1-ehat))  );  // ok
  AveLogAlpha   = std::max(  0.1,  std::log(std::exp(AveLogAlphah) + ParsAveA1/Fs)                       );  // ok
  //
  SigmaLogTmax  = std::min(  0.5,  1.00 / (ParsSigLogT1 + ParsSigLogT2*std::log(y))                      );  // ok
  SigmaLogAlpha = std::min(  0.5,  1.00 / (ParsSigLogA1 + ParsSigLogA2*std::log(y))                      );  // ok
  Rho           = ParsRho1 + ParsRho2*std::log(y);                                                           // ok
 
 
  if (0) {
    G4cout << " y            = " << y << G4endl;
    G4cout << " std::log(std::exp(AveLogTmaxh)  + ParsAveT1/Fs + ParsAveT2*(1-ehat)) = "
       << " std::log(" << std::exp(AveLogTmaxh) << " + " << ParsAveT1/Fs << " + " << ParsAveT2*(1-ehat) << ") = "
       << " std::log(" << std::exp(AveLogTmaxh) << " + " << ParsAveT1 << "/" << Fs << " + " << ParsAveT2 << "*" << (1-ehat) << ") = "
       << " std::log(" << std::exp(AveLogTmaxh)  + ParsAveT1/Fs + ParsAveT2*(1-ehat) << ")" << G4endl;
    G4cout << " AveLogTmaxh    " <<  AveLogTmaxh     << G4endl;
    G4cout << " AveLogAlphah   " <<  AveLogAlphah    << G4endl;
    G4cout << " SigmaLogTmaxh  " <<  SigmaLogTmaxh   << G4endl;
    G4cout << " 1.00/( ParSigLogT1 + ParSigLogT2*std::log(y)  ) = " <<  1.00 << "/" << ( ParSigLogT1 + ParSigLogT2*std::log(y)  ) << " = "
       << 1.00 << "/" << "(" << ParSigLogT1 << " + " <<  ParSigLogT2*std::log(y) << " )  = "
       << 1.00 << "/" << "(" << ParSigLogT1 << " + " <<  ParSigLogT2 << "*" << std::log(y) << " ) "
       << G4endl;
    G4cout << " SigmaLogAlphah " <<  SigmaLogAlphah  << G4endl;
    G4cout << " Rhoh           " <<  Rhoh            << G4endl;
    G4cout << " AveLogTmax     " <<  AveLogTmax      << G4endl;
    G4cout << " AveLogAlpha    " <<  AveLogAlpha     << G4endl;
    G4cout << " SigmaLogTmax   " <<  SigmaLogTmax    << G4endl;
    G4cout << " SigmaLogAlpha  " <<  SigmaLogAlpha   << G4endl;
    G4cout << " Rho            " <<  Rho             << G4endl;
  }
}

References AveLogAlpha, AveLogAlphah, AveLogTmax, AveLogTmaxh, ehat, Fs, G4cout, G4endl, G4INCL::Math::max(), G4INCL::Math::min(), ParAveA1, ParAveA2, ParAveA3, ParAveT1, ParRho1, ParRho2, ParsAveA1, ParsAveT1, ParsAveT2, ParSigLogA1, ParSigLogA2, ParSigLogT1, ParSigLogT2, ParsRho1, ParsRho2, ParsSigLogA1, ParsSigLogA2, ParsSigLogT1, ParsSigLogT2, Rho, Rhoh, SigmaLogAlpha, SigmaLogAlphah, SigmaLogTmax, SigmaLogTmaxh, and Zeff.

Referenced by GenerateLongitudinalProfile().

◆ ComputeRadialParameters()

void GFlashSamplingShowerParameterisation::ComputeRadialParameters	(	G4double	y,
		G4double	Tau
	)

virtual

Implements GVFlashShowerParameterisation.

Definition at line 411 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double z1 = ParRC1 + ParRC2* std::log(Energy/GeV);         //ok
  G4double z2 = ParRC3+ParRC4*Zeff;                            //ok
  RadiusCore  =  z1 + z2 * Tau;                                //ok 
  G4double p1 = ParWC1+ParWC2*Zeff;                            //ok
  G4double p2 = ParWC3+ParWC4*Zeff;                            //ok
  G4double p3 = ParWC5+ParWC6*std::log(Energy/GeV);            //ok
  WeightCore   =  p1 * std::exp( (p2-Tau)/p3-  std::exp( (p2-Tau) /p3) ); //ok
  
  G4double k1 = ParRT1+ParRT2*Zeff;                // ok
  G4double k2 = ParRT3;                            // ok
  G4double k3 = ParRT4;                            // ok
  G4double k4 = ParRT5+ParRT6* std::log(Energy/GeV);    // ok
  
  RadiusTail   = k1*(std::exp(k3*(Tau-k2))
               + std::exp(k4*(Tau-k2)) );            //ok
 
  // sampling calorimeter  
 
  RadiusCore   = RadiusCore + ParsRC1*(1-ehat) + ParsRC2/Fs*std::exp(-Tau); //ok
  WeightCore   = WeightCore + (1-ehat)
                            * (ParsWC1+ParsWC2/Fs * std::exp(-std::pow((Tau-1.),2))); //ok
  RadiusTail   = RadiusTail + (1-ehat)* ParsRT1+ ParsRT2/Fs *std::exp(-Tau);     //ok  
}

References ehat, Fs, GeV, ParRC1, ParRC2, ParRC3, ParRC4, ParRT1, ParRT2, ParRT3, ParRT4, ParRT5, ParRT6, ParsRC1, ParsRC2, ParsRT1, ParsRT2, ParsWC1, ParsWC2, ParWC1, ParWC2, ParWC3, ParWC4, ParWC5, ParWC6, RadiusCore, RadiusTail, WeightCore, and Zeff.

Referenced by GenerateRadius().

◆ ComputeTau()

G4double GFlashSamplingShowerParameterisation::ComputeTau ( G4double LongitudinalPosition )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 400 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double tau = LongitudinalPosition / Tmax/ X0eff     //<t> = T* a /(a - 1) 
                 * (Alpha-1.00) /Alpha
                 * std::exp(AveLogAlpha)/(std::exp(AveLogAlpha)-1.);  //ok 
  return tau;
}

References Alpha, AveLogAlpha, Tmax, and X0eff.

Referenced by GenerateRadius().

◆ ComputeZAX0EFFetc()

void GFlashSamplingShowerParameterisation::ComputeZAX0EFFetc ( )

Definition at line 185 of file GFlashSamplingShowerParameterisation.cc.

{
  G4cout << "/************ ComputeZAX0EFFetc ************/" << G4endl;
  G4cout << "  - GFlashSamplingShowerParameterisation::Material -  " << G4endl;
 
  G4double Es = 21*MeV; //constant
 
  // material and geometry parameters for a sampling calorimeter
  G4double denominator = (d1*density1 + d2*density2);
  G4double W1  = (d1*density1) / denominator;
  G4double W2  = (d2*density2) / denominator;
  Zeff   = ( W1*Z1 ) + ( W2*Z2 );    //X0*Es/Ec;
  Aeff   = ( W1*A1 ) + ( W2*A2 );
  Rhoeff = ( ( d1*density1 ) + ( d2*density2 ) )  / ( d1 + d2 ); // --> was G4double ( d2  + d1 );
  X0eff  = (W1 * Rhoeff) / (X01 * density1) + (W2 * Rhoeff) / (X02 * density2 );
  X0eff  = 1./ X0eff;
  Rmeff =  1./  (  (  ((W1*Ec1)/X01) +  ((W2*Ec2)/X02)  )  /  Es  ) ;
  Eceff =  X0eff * (   (W1*Ec1)/X01  +   (W2*Ec2)/X02   );
  Fs =  X0eff/(d1+d2);// --> was G4double ((d1/mm )+(d2/mm) ); Can't understand if dividing by mm makes sense... looks weird.
  ehat = (  1. / ( 1 + 0.007*(Z1- Z2) )  );
 
  G4cout << "W1= " << W1 << G4endl;
  G4cout << "W2= " << W2 << G4endl;
  G4cout << "effective quantities Zeff = "<<Zeff<< G4endl;
  G4cout << "effective quantities Aeff = "<<Aeff<< G4endl;
  G4cout << "effective quantities Rhoeff = "<<Rhoeff/g *cm3<<" g/cm3"  << G4endl;
  G4cout << "effective quantities X0eff = "<<X0eff/cm <<" cm" << G4endl;
 
  X0eff = X0eff * Rhoeff;
 
  G4cout << "effective quantities X0eff = "<<X0eff/g*cm2 <<" g/cm2" << G4endl;  
  X0eff = X0eff /Rhoeff;
  G4cout << "effective quantities RMeff = "<<Rmeff/cm<<"  cm" << G4endl; 
  Rmeff = Rmeff* Rhoeff;
  G4cout << "effective quantities RMeff = "<<Rmeff/g *cm2<<" g/cm2" << G4endl;  
  Rmeff = Rmeff/ Rhoeff;
  G4cout << "effective quantities Eceff = "<<Eceff/MeV<< " MeV"<< G4endl;  
  G4cout << "effective quantities Fs = "<<Fs<<G4endl;
  G4cout << "effective quantities ehat = "<<ehat<<G4endl;
  G4cout << "/********************************************/ " <<G4endl; 
}

References A1, A2, Aeff, cm, cm2, cm3, d1, d2, density1, density2, Ec1, Ec2, Eceff, ehat, Fs, g, G4cout, G4endl, MeV, Rhoeff, Rmeff, X01, X02, X0eff, Z1, Z2, and Zeff.

Referenced by GFlashSamplingShowerParameterisation().

◆ gam()

G4double GVFlashShowerParameterisation::gam	(	G4double	x,
		G4double	a
	)		const

inherited

Definition at line 121 of file GVFlashShowerParameterisation.cc.

{
  return fGamma->Gamma(a, x); 
}

References GVFlashShowerParameterisation::fGamma, and MyGamma::Gamma().

Referenced by GFlashHomoShowerParameterisation::IntegrateEneLongitudinal(), IntegrateEneLongitudinal(), GFlashHomoShowerParameterisation::IntegrateNspLongitudinal(), and IntegrateNspLongitudinal().

◆ GenerateEnergyProfile()

void GFlashSamplingShowerParameterisation::GenerateEnergyProfile ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 288 of file GFlashSamplingShowerParameterisation.cc.

{ 
  G4double Correlation1  = std::sqrt( (1+Rho )/2 );
  G4double Correlation2  = std::sqrt( (1-Rho )/2 );
  G4double Correlation1h = std::sqrt( (1+Rhoh)/2 );
  G4double Correlation2h = std::sqrt( (1-Rhoh)/2 );
  G4double Random1 = G4RandGauss::shoot();
  G4double Random2 = G4RandGauss::shoot();
 
  Tmax  = std::max(  1.,   std::exp( AveLogTmax  + SigmaLogTmax  * (Correlation1*Random1 + Correlation2*Random2) ) );
  Alpha = std::max(  1.1,  std::exp( AveLogAlpha + SigmaLogAlpha * (Correlation1*Random1 - Correlation2*Random2) ) );
  Beta  = (Alpha-1.00)/Tmax;
  //Parameters for Enenrgy Profile including correaltion and sigmas  
  Tmaxh  = std::exp( AveLogTmaxh  + SigmaLogTmaxh  *
  (Correlation1h*Random1 + Correlation2h*Random2) );  
  Alphah = std::exp( AveLogAlphah + SigmaLogAlphah *
  (Correlation1h*Random1 - Correlation2h*Random2) );
  Betah  = (Alphah-1.00)/Tmaxh;
}

References Alpha, Alphah, AveLogAlpha, AveLogAlphah, AveLogTmax, AveLogTmaxh, Beta, Betah, G4INCL::Math::max(), Rho, Rhoh, G4INCL::DeJongSpin::shoot(), SigmaLogAlpha, SigmaLogAlphah, SigmaLogTmax, SigmaLogTmaxh, Tmax, and Tmaxh.

Referenced by GenerateLongitudinalProfile().

◆ GenerateExponential()

G4double GFlashSamplingShowerParameterisation::GenerateExponential ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 440 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double ParExp1 =  9./7.*X0eff;
  G4double random  = -ParExp1*G4RandExponential::shoot() ;
  return random;
}

References G4INCL::DeJongSpin::shoot(), and X0eff.

◆ GenerateLongitudinalProfile()

void GFlashSamplingShowerParameterisation::GenerateLongitudinalProfile ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 229 of file GFlashSamplingShowerParameterisation.cc.

{
  if ((material1==0) || (material2 ==0))
  {
    G4Exception("GFlashSamplingShowerParameterisation::GenerateLongitudinalProfile()",
                "InvalidSetup", FatalException, "No material initialized!");
  }  
  G4double y = Energy/Eceff;
  ComputeLongitudinalParameters(y);  
  GenerateEnergyProfile(y);
  GenerateNSpotProfile(y);
}

References ComputeLongitudinalParameters(), Eceff, FatalException, G4Exception(), GenerateEnergyProfile(), GenerateNSpotProfile(), material1, and material2.

◆ GenerateNSpotProfile()

void GFlashSamplingShowerParameterisation::GenerateNSpotProfile ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 310 of file GFlashSamplingShowerParameterisation.cc.

{
  TNSpot     = Tmaxh *  (ParsSpotT1+ParsSpotT2*Zeff); //ok.
  TNSpot     = std::max(0.5,Tmaxh * (ParsSpotT1+ParsSpotT2*Zeff));
  AlphaNSpot = Alphah * (ParsSpotA1+ParsSpotA2*Zeff);     
  BetaNSpot  = (AlphaNSpot-1.00)/TNSpot;           // ok 
  NSpot      = ParsSpotN1 /SamplingResolution * std::pow(y*Eceff/GeV,ParsSpotN2 );
}

References Alphah, AlphaNSpot, BetaNSpot, Eceff, GeV, G4INCL::Math::max(), NSpot, ParsSpotA1, ParsSpotA2, ParsSpotN1, ParsSpotN2, ParsSpotT1, ParsSpotT2, SamplingResolution, Tmaxh, TNSpot, and Zeff.

Referenced by GenerateLongitudinalProfile().

◆ GeneratePhi()

G4double GFlashSamplingShowerParameterisation::GeneratePhi ( )

◆ GenerateRadius()

G4double GFlashSamplingShowerParameterisation::GenerateRadius	(	G4int	ispot,
		G4double	Energy,
		G4double	LongitudinalPosition
	)

virtual

Implements GVFlashShowerParameterisation.

Definition at line 370 of file GFlashSamplingShowerParameterisation.cc.

{
  if(ispot < 1) 
  {
    // Determine lateral parameters in the middle of the step.
    // They depend on energy & position along step
    //
    G4double Tau = ComputeTau(LongitudinalPosition);
    ComputeRadialParameters(Energy,Tau);  
  }
  
  G4double Radius;
  G4double Random1 = G4UniformRand();
  G4double Random2 = G4UniformRand(); 
  if(Random1  <WeightCore) //WeightCore = p < w_i  
  {
    Radius = Rmeff * RadiusCore * std::sqrt( Random2/(1. - Random2) );
  }
  else
  {
    Radius = Rmeff * RadiusTail * std::sqrt( Random2/(1. - Random2) );
  }   
  Radius =  std::min(Radius,DBL_MAX);
  return Radius;
}

References ComputeRadialParameters(), ComputeTau(), DBL_MAX, G4UniformRand, G4INCL::Math::min(), RadiusCore, RadiusTail, Rmeff, and WeightCore.

◆ GetAveR90()

G4double GFlashSamplingShowerParameterisation::GetAveR90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 76 of file GFlashSamplingShowerParameterisation.hh.

76{return (1.5 * Rmeff);} //ok

References Rmeff.

◆ GetAveR99()

G4double GFlashSamplingShowerParameterisation::GetAveR99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 75 of file GFlashSamplingShowerParameterisation.hh.

75{return (3.5 * Rmeff);}

References Rmeff.

◆ GetAveT90()

G4double GFlashSamplingShowerParameterisation::GetAveT90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 80 of file GFlashSamplingShowerParameterisation.hh.

80{return (2.5* X0eff* std::exp( AveLogTmax));}

References AveLogTmax, and X0eff.

◆ GetAveT99()

G4double GFlashSamplingShowerParameterisation::GetAveT99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 79 of file GFlashSamplingShowerParameterisation.hh.

79{return (X0eff*AveLogTmax/(AveLogAlpha-1.00));}

References AveLogAlpha, AveLogTmax, and X0eff.

◆ GetAveTmx()

G4double GFlashSamplingShowerParameterisation::GetAveTmx ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 78 of file GFlashSamplingShowerParameterisation.hh.

78{return (X0eff*std::exp(AveLogTmax));}

References AveLogTmax, and X0eff.

◆ GetEc()

G4double GFlashSamplingShowerParameterisation::GetEc ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 84 of file GFlashSamplingShowerParameterisation.hh.

84{return Eceff;}

References Eceff.

◆ GetEffA()

G4double GVFlashShowerParameterisation::GetEffA ( const G4Material * material )

inherited

Definition at line 82 of file GVFlashShowerParameterisation.cc.

{
  // Returns A or effective A=sum(pi*Ai) (if compound/mixture)
  // of given material
  //
  G4double a = 0.;
  G4int nofElements = mat->GetNumberOfElements();
  if (nofElements > 1) {
    for (G4int i=0; i<nofElements; i++) {
      G4double aOfElement = mat->GetElement(i)->GetA()/(g/mole);
      G4double massFraction = mat->GetFractionVector()[i];     
      a += aOfElement*massFraction;
    }
  }
  else { 
    a = mat->GetA()/(g/mole);
  }
  return a;
}

References g, G4Element::GetA(), G4Material::GetA(), G4Material::GetElement(), G4Material::GetFractionVector(), G4Material::GetNumberOfElements(), and mole.

Referenced by GFlashHomoShowerParameterisation::SetMaterial(), and SetMaterial().

◆ GetEffZ()

G4double GVFlashShowerParameterisation::GetEffZ ( const G4Material * material )

inherited

Definition at line 59 of file GVFlashShowerParameterisation.cc.

{
  // Returns Z or effective Z=sum(pi*Zi) (if compound/mixture)
  // of given material
  //
  G4double z = 0.;
  G4int nofElements = mat->GetNumberOfElements();
  if (nofElements > 1) 
  {
    for (G4int i=0; i<nofElements; i++) {
      G4double zOfElement = mat->GetElement(i)->GetZ();
      G4double massFraction = mat->GetFractionVector()[i];
      // cout << mat->GetElement(i)->GetName()
      //      <<" Z= "<<zOfElement << " , Fraction= "<<massFraction <<endl;
      z += zOfElement*massFraction;
    }
  }
  else { 
    z = mat->GetZ(); 
  }  
  return z;
}

References G4Material::GetElement(), G4Material::GetFractionVector(), G4Material::GetNumberOfElements(), G4Element::GetZ(), and G4Material::GetZ().

Referenced by GFlashHomoShowerParameterisation::SetMaterial(), and SetMaterial().

◆ GetNspot()

G4double GFlashSamplingShowerParameterisation::GetNspot ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 82 of file GFlashSamplingShowerParameterisation.hh.

82{return NSpot;}

References NSpot.

◆ GetRm()

G4double GFlashSamplingShowerParameterisation::GetRm ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 85 of file GFlashSamplingShowerParameterisation.hh.

85{return Rmeff;}

References Rmeff.

◆ GetX0()

G4double GFlashSamplingShowerParameterisation::GetX0 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 83 of file GFlashSamplingShowerParameterisation.hh.

83{return X0eff;}

References X0eff.

◆ IntegrateEneLongitudinal()

G4double GFlashSamplingShowerParameterisation::IntegrateEneLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 344 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double LongitudinalStepInX0 = LongitudinalStep / X0eff;
  G4float x1= Betah*LongitudinalStepInX0;
  G4float x2= Alphah;
  float x3 =  gam(x1,x2);
  G4double DEne=x3;
  return DEne;
}

References Alphah, Betah, GVFlashShowerParameterisation::gam(), and X0eff.

◆ IntegrateNspLongitudinal()

G4double GFlashSamplingShowerParameterisation::IntegrateNspLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 357 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double LongitudinalStepInX0 = LongitudinalStep / X0eff; 
  G4float x1 = BetaNSpot*LongitudinalStepInX0;
  G4float x2 = AlphaNSpot;
  G4float x3 =  gam(x1,x2);
  G4double DNsp = x3;
  return DNsp;
}

References AlphaNSpot, BetaNSpot, GVFlashShowerParameterisation::gam(), and X0eff.

◆ PrintMaterial()

void GVFlashShowerParameterisation::PrintMaterial ( const G4Material * mat )

inherited

Definition at line 102 of file GVFlashShowerParameterisation.cc.

{
  G4cout<<"/********************************************/ " << G4endl;
  G4cout<<"  - GVFlashShowerParameterisation::Material -  " << G4endl;
  G4cout<<"        Material : " << mat->GetName()  << G4endl;
  G4cout<<"   Z  = " << Z      << G4endl;
  G4cout<<"   A  = " << A      << G4endl;
  G4cout<<"   X0 = " << X0/cm  << " cm" << G4endl;
  G4cout<<"   Rm = " << Rm/cm  << " cm" << G4endl;
  G4cout<<"   Ec = " << Ec/MeV << " MeV"<< G4endl;
  G4cout<<"/********************************************/ " << G4endl; 
}

References GVFlashShowerParameterisation::A, cm, GVFlashShowerParameterisation::Ec, G4cout, G4endl, G4Material::GetName(), MeV, GVFlashShowerParameterisation::Rm, GVFlashShowerParameterisation::X0, and GVFlashShowerParameterisation::Z.

Referenced by GFlashHomoShowerParameterisation::GFlashHomoShowerParameterisation().

◆ SetMaterial()

void GFlashSamplingShowerParameterisation::SetMaterial	(	G4Material *	mat1,
		G4Material *	mat2
	)

Definition at line 161 of file GFlashSamplingShowerParameterisation.cc.

{
  G4double Es = 21*MeV;
  material1= mat1;
  Z1 = GetEffZ(material1);
  A1 = GetEffA(material1);
  density1 = material1->GetDensity();
  X01  = material1->GetRadlen();   
  Ec1      = 2.66 * std::pow((X01 * Z1 / A1),1.1); 
  Rm1 = X01*Es/Ec1;
 
  material2= mat2;
  Z2 = GetEffZ(material2);
  A2 = GetEffA(material2);
  density2 = material2->GetDensity();
  X02  = material2->GetRadlen();   
  Ec2      = 2.66 * std::pow((X02 * Z2 / A2),1.1); 
  Rm2 = X02*Es/Ec2; 
  // PrintMaterial(); 
}