GFlashSamplingShowerParameterisation Class Reference

#include <GFlashSamplingShowerParameterisation.hh>

Inheritance diagram for GFlashSamplingShowerParameterisation:

Public Member Functions
	GFlashSamplingShowerParameterisation (G4Material *aMat1, G4Material *aMat2, G4double d1, G4double d2, GFlashSamplingShowerTuning *aPar=0)
	~GFlashSamplingShowerParameterisation ()
void	ComputeRadialParameters (G4double y, G4double Tau)
void	GenerateLongitudinalProfile (G4double Energy)
void	ComputeZAX0EFFetc ()
G4double	IntegrateEneLongitudinal (G4double LongitudinalStep)
G4double	IntegrateNspLongitudinal (G4double LongitudinalStep)
G4double	ComputeTau (G4double LongitudinalPosition)
void	SetMaterial (G4Material *mat1, G4Material *mat2)
G4double	GeneratePhi ()
G4double	GenerateRadius (G4int ispot, G4double Energy, G4double LongitudinalPosition)
G4double	GenerateExponential (G4double Energy)
G4double	GetAveR99 ()
G4double	GetAveR90 ()
G4double	GetAveTmx ()
G4double	GetAveT99 ()
G4double	GetAveT90 ()
G4double	GetNspot ()
G4double	GetX0 ()
G4double	GetEc ()
G4double	GetRm ()
G4double	GetEffDensity ()
G4double	ApplySampling (const G4double DEne, const G4double Energy)
Public Member Functions inherited from GVFlashShowerParameterisation
	GVFlashShowerParameterisation ()
virtual	~GVFlashShowerParameterisation ()
G4double	GeneratePhi ()
G4double	GetEffZ (const G4Material *material)
G4double	GetEffA (const G4Material *material)
G4double	gam (G4double x, G4double a) const
void	PrintMaterial (const G4Material *mat)

Additional Inherited Members
Protected Attributes inherited from GVFlashShowerParameterisation
GVFlashHomoShowerTuning *	thePar
G4double	density
G4double	A
G4double	Z
G4double	X0
G4double	Ec
G4double	Rm
G4double	NSpot

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 180 of file GFlashSamplingShowerParameterisation.cc.

{
   delete thePar;
}

Member Function Documentation

ApplySampling()

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

Definition at line 352 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;
    G4double x1 = DEne / Resolution;
    G4double x2 = 1.0 / Resolution;
    DEneFluctuated = G4RandGamma::shoot(x1, x2);
  }
  return DEneFluctuated;
}

ComputeRadialParameters()

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 437 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
}

Referenced by GenerateRadius().

ComputeTau()

G4double GFlashSamplingShowerParameterisation::ComputeTau ( G4double LongitudinalPosition )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 427 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;
}

Referenced by GenerateRadius().

ComputeZAX0EFFetc()

void GFlashSamplingShowerParameterisation::ComputeZAX0EFFetc

(

)

Definition at line 212 of file GFlashSamplingShowerParameterisation.cc.

{
  G4cout << "/************ ComputeZAX0EFFetc ************/" << G4endl;
  G4cout << "  - GFlashSamplingShowerParameterisation::Material -  " << G4endl;
 
  const G4double Es = 21.2 * MeV;
 
  // 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;
}

Referenced by GFlashSamplingShowerParameterisation().

GenerateExponential()

G4double GFlashSamplingShowerParameterisation::GenerateExponential ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 464 of file GFlashSamplingShowerParameterisation.cc.

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

GenerateLongitudinalProfile()

void GFlashSamplingShowerParameterisation::GenerateLongitudinalProfile ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 257 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);
}

GeneratePhi()

G4double GFlashSamplingShowerParameterisation::GeneratePhi

(

)

GenerateRadius()

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 400 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;
}

GetAveR90()

G4double GFlashSamplingShowerParameterisation::GetAveR90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 73 of file GFlashSamplingShowerParameterisation.hh.

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

GetAveR99()

G4double GFlashSamplingShowerParameterisation::GetAveR99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 72 of file GFlashSamplingShowerParameterisation.hh.

{ return (3.5 * Rmeff); }

GetAveT90()

G4double GFlashSamplingShowerParameterisation::GetAveT90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 77 of file GFlashSamplingShowerParameterisation.hh.

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

GetAveT99()

G4double GFlashSamplingShowerParameterisation::GetAveT99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 76 of file GFlashSamplingShowerParameterisation.hh.

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

GetAveTmx()

G4double GFlashSamplingShowerParameterisation::GetAveTmx ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 75 of file GFlashSamplingShowerParameterisation.hh.

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

GetEc()

G4double GFlashSamplingShowerParameterisation::GetEc ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 81 of file GFlashSamplingShowerParameterisation.hh.

{ return Eceff; }

GetEffDensity()

G4double GFlashSamplingShowerParameterisation::GetEffDensity ( )

inline

Definition at line 83 of file GFlashSamplingShowerParameterisation.hh.

{ return Rhoeff; }

GetNspot()

G4double GFlashSamplingShowerParameterisation::GetNspot ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 79 of file GFlashSamplingShowerParameterisation.hh.

{ return NSpot; }

GetRm()

G4double GFlashSamplingShowerParameterisation::GetRm ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 82 of file GFlashSamplingShowerParameterisation.hh.

{ return Rmeff; }

GetX0()

G4double GFlashSamplingShowerParameterisation::GetX0 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 80 of file GFlashSamplingShowerParameterisation.hh.

{ return X0eff; }

IntegrateEneLongitudinal()

G4double GFlashSamplingShowerParameterisation::IntegrateEneLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 375 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;
}

IntegrateNspLongitudinal()

G4double GFlashSamplingShowerParameterisation::IntegrateNspLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 388 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;
}

SetMaterial()

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

Definition at line 187 of file GFlashSamplingShowerParameterisation.cc.

{
  const G4double Es = 21.2 * 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);
  // Ec1 = 610.0 * MeV / (Z1 + 1.24);
  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);
  // Ec2 = 610.0 * MeV / (Z2 + 1.24);
  Rm2 = X02 * Es / Ec2;
  // PrintMaterial();
}

Referenced by GFlashSamplingShowerParameterisation().

---

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

• GFlashSamplingShowerParameterisation.hh
• GFlashSamplingShowerParameterisation.cc