GFlashHomoShowerParameterisation Class Reference

#include <GFlashHomoShowerParameterisation.hh>

Inheritance diagram for GFlashHomoShowerParameterisation:

Public Member Functions
	GFlashHomoShowerParameterisation (G4Material *aMat, GVFlashHomoShowerTuning *aPar=0)
	~GFlashHomoShowerParameterisation ()
void	ComputeRadialParameters (G4double y, G4double Tau)
void	GenerateLongitudinalProfile (G4double Energy)
void	ComputeZAX0EFFetc ()
G4double	IntegrateEneLongitudinal (G4double LongitudinalStep)
G4double	IntegrateNspLongitudinal (G4double LongitudinalStep)
G4double	ComputeTau (G4double LongitudinalPosition)
G4double	GeneratePhi ()
G4double	GenerateRadius (G4int ispot, G4double Energy, G4double LongitudinalPosition)
G4double	GenerateExponential (G4double Energy)
void	SetMaterial (G4Material *mat)
G4double	GetAveR99 ()
G4double	GetAveR90 ()
G4double	GetAveTmx ()
G4double	GetAveT99 ()
G4double	GetAveT90 ()
G4double	GetNspot ()
G4double	GetX0 ()
G4double	GetEc ()
G4double	GetRm ()
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 50 of file GFlashHomoShowerParameterisation.hh.

Constructor & Destructor Documentation

GFlashHomoShowerParameterisation::GFlashHomoShowerParameterisation	(	G4Material *	aMat,
		GVFlashHomoShowerTuning *	aPar = 0
	)

Definition at line 49 of file GFlashHomoShowerParameterisation.cc.

References G4cout, G4endl, GVFlashHomoShowerTuning::ParAveA1(), GVFlashHomoShowerTuning::ParAveA2(), GVFlashHomoShowerTuning::ParAveA3(), GVFlashHomoShowerTuning::ParAveT1(), GVFlashHomoShowerTuning::ParRC1(), GVFlashHomoShowerTuning::ParRC2(), GVFlashHomoShowerTuning::ParRC3(), GVFlashHomoShowerTuning::ParRC4(), GVFlashHomoShowerTuning::ParRho1(), GVFlashHomoShowerTuning::ParRho2(), GVFlashHomoShowerTuning::ParRT1(), GVFlashHomoShowerTuning::ParRT2(), GVFlashHomoShowerTuning::ParRT3(), GVFlashHomoShowerTuning::ParRT4(), GVFlashHomoShowerTuning::ParRT5(), GVFlashHomoShowerTuning::ParRT6(), GVFlashHomoShowerTuning::ParSigLogA1(), GVFlashHomoShowerTuning::ParSigLogA2(), GVFlashHomoShowerTuning::ParSigLogT1(), GVFlashHomoShowerTuning::ParSigLogT2(), GVFlashHomoShowerTuning::ParSpotA1(), GVFlashHomoShowerTuning::ParSpotA2(), GVFlashHomoShowerTuning::ParSpotN1(), GVFlashHomoShowerTuning::ParSpotN2(), GVFlashHomoShowerTuning::ParSpotT1(), GVFlashHomoShowerTuning::ParSpotT2(), GVFlashHomoShowerTuning::ParWC1(), GVFlashHomoShowerTuning::ParWC2(), GVFlashHomoShowerTuning::ParWC3(), GVFlashHomoShowerTuning::ParWC4(), GVFlashHomoShowerTuning::ParWC5(), GVFlashHomoShowerTuning::ParWC6(), GVFlashShowerParameterisation::PrintMaterial(), and SetMaterial().

  : GVFlashShowerParameterisation(),
    ConstantResolution(0.), NoiseResolution(0.), SamplingResolution(0.),
    AveLogAlphah(0.), AveLogTmaxh(0.), SigmaLogAlphah(0.), SigmaLogTmaxh(0.),
    Rhoh(0.), Alphah(0.), Tmaxh(0.), Betah(0.)

{  
  if(!aPar) { thePar = new GVFlashHomoShowerTuning; owning = true; }
  else      { thePar = aPar; owning = false; }

  SetMaterial(aMat);
  PrintMaterial(aMat);

  /********************************************/
  /* Homo Calorimeter                         */
  /********************************************/ 
  // 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
  ParSigLogT1 = thePar->ParSigLogT1();     // Sigma T1 (-1.4 + 1.26 ln y)**-1 
  ParSigLogT2 = thePar->ParSigLogT2();

  // variance of 'alpha'
  //
  ParSigLogA1 = thePar->ParSigLogA1();    // Sigma a (-0.58 + 0.86 ln y)**-1 
  ParSigLogA2 = thePar->ParSigLogA2();

  // correlation alpha%T
  //
  ParRho1     = thePar->ParRho1();   // Rho = 0.705 -0.023 ln y 
  ParRho2     = thePar->ParRho2();

  // 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();

  // Coeff for fluctueted radial  profiles for a uniform media
  //
  ParSpotT1   = thePar->ParSpotT1();     // T_spot = T_hom =(0.698 + 0.00212)
  ParSpotT2   = thePar->ParSpotT2();

  ParSpotA1   = thePar->ParSpotA1();     // a_spot= a_hom (0.639 + 0.00334)
  ParSpotA2   = thePar->ParSpotA2();  

  ParSpotN1   = thePar->ParSpotN1();     // N_Spot 93 * ln(Z) E ** 0.876   
  ParSpotN2   = thePar->ParSpotN2(); 

  // Inits

  NSpot         = 0.00;
  AlphaNSpot    = 0.00;
  TNSpot        = 0.00;
  BetaNSpot     = 0.00;

  RadiusCore    = 0.00;
  WeightCore    = 0.00;
  RadiusTail    = 0.00; 

  G4cout << "/********************************************/ " << G4endl;
  G4cout << "  - GFlashHomoShowerParameterisation::Constructor -  " << G4endl;  
  G4cout << "/********************************************/ " << G4endl;
}

GFlashHomoShowerParameterisation::~GFlashHomoShowerParameterisation

(

)

Definition at line 152 of file GFlashHomoShowerParameterisation.cc.

{
  if(owning) { delete thePar; }
}

Member Function Documentation

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 273 of file GFlashHomoShowerParameterisation.cc.

References python.hepunit::GeV, and GVFlashShowerParameterisation::Z.

Referenced by GenerateRadius().

{
  G4double z1 = ParRC1 + ParRC2* std::log(Energy/GeV)  ;       //ok
  G4double z2 = ParRC3+ParRC4*Z ;                              //ok
  RadiusCore   =  z1 + z2 * Tau  ;                             //ok 

  G4double p1 = ParWC1+ParWC2*Z;                               //ok
  G4double p2 = ParWC3+ParWC4*Z;                               //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*Z;                   // 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
}

G4double GFlashHomoShowerParameterisation::ComputeTau ( G4double  LongitudinalPosition )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 264 of file GFlashHomoShowerParameterisation.cc.

References GVFlashShowerParameterisation::X0.

Referenced by GenerateRadius().

{
  G4double tau = LongitudinalPosition / Tmaxh / X0     //<t> = T* a /(a - 1) 
  * (Alphah-1.00) /Alphah * 
  std::exp(AveLogAlphah)/(std::exp(AveLogAlphah)-1.);  //ok 
  return tau;
}

void GFlashHomoShowerParameterisation::ComputeZAX0EFFetc

(

)

G4double GFlashHomoShowerParameterisation::GenerateExponential ( G4double  Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 295 of file GFlashHomoShowerParameterisation.cc.

References G4INCL::DeJongSpin::shoot(), and GVFlashShowerParameterisation::X0.

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

void GFlashHomoShowerParameterisation::GenerateLongitudinalProfile ( G4double  Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 158 of file GFlashHomoShowerParameterisation.cc.

References GVFlashShowerParameterisation::Ec, FatalException, and G4Exception().

{
  if (material==0) 
  {
    G4Exception("GFlashHomoShowerParameterisation::GenerateLongitudinalProfile()",
                "InvalidSetup", FatalException, "No material initialized!");
  }
  
  G4double y = Energy/Ec;
  ComputeLongitudinalParameters(y);  
  GenerateEnergyProfile(y);
  GenerateNSpotProfile(y);
}

G4double GFlashHomoShowerParameterisation::GeneratePhi

(

)

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 236 of file GFlashHomoShowerParameterisation.cc.

References ComputeRadialParameters(), ComputeTau(), DBL_MAX, G4UniformRand, G4INCL::Math::min(), and GVFlashShowerParameterisation::Rm.

{
  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 = Rm * RadiusCore * std::sqrt( Random2/(1. - Random2) );
  }
  else
  {
    Radius = Rm * RadiusTail * std::sqrt( Random2/(1. - Random2) );
  }   
  Radius =  std::min(Radius,DBL_MAX);
  return Radius;
}

G4double GFlashHomoShowerParameterisation::GetAveR90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 73 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::Rm.

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

G4double GFlashHomoShowerParameterisation::GetAveR99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 72 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::Rm.

{return (3.5 * Rm);}

G4double GFlashHomoShowerParameterisation::GetAveT90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 77 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::X0.

{return (2.5* X0*std::exp( AveLogTmaxh) );}

G4double GFlashHomoShowerParameterisation::GetAveT99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 76 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::X0.

{return (X0 * AveLogTmaxh/(AveLogAlphah-1.00));}

G4double GFlashHomoShowerParameterisation::GetAveTmx ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 75 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::X0.

{return (X0 * std::exp(AveLogTmaxh));}

G4double GFlashHomoShowerParameterisation::GetEc ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 81 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::Ec.

{return Ec;} 

G4double GFlashHomoShowerParameterisation::GetNspot ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 79 of file GFlashHomoShowerParameterisation.hh.

{ return NSpot;}

G4double GFlashHomoShowerParameterisation::GetRm ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 82 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::Rm.

{return Rm;} 

G4double GFlashHomoShowerParameterisation::GetX0 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 80 of file GFlashHomoShowerParameterisation.hh.

References GVFlashShowerParameterisation::X0.

{return X0;}  

G4double GFlashHomoShowerParameterisation::IntegrateEneLongitudinal ( G4double  LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 213 of file GFlashHomoShowerParameterisation.cc.

References GVFlashShowerParameterisation::gam(), and GVFlashShowerParameterisation::X0.

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

G4double GFlashHomoShowerParameterisation::IntegrateNspLongitudinal ( G4double  LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 224 of file GFlashHomoShowerParameterisation.cc.

References GVFlashShowerParameterisation::gam(), and GVFlashShowerParameterisation::X0.

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

void GFlashHomoShowerParameterisation::SetMaterial

(

G4Material *

mat

)

Definition at line 139 of file GFlashHomoShowerParameterisation.cc.

References GVFlashShowerParameterisation::A, python.hepunit::cm3, GVFlashShowerParameterisation::density, GVFlashShowerParameterisation::Ec, g(), G4Material::GetDensity(), GVFlashShowerParameterisation::GetEffA(), GVFlashShowerParameterisation::GetEffZ(), G4Material::GetRadlen(), python.hepunit::MeV, GVFlashShowerParameterisation::Rm, GVFlashShowerParameterisation::X0, and GVFlashShowerParameterisation::Z.

Referenced by GFlashHomoShowerParameterisation().

{
  material= mat;
  Z = GetEffZ(material);
  A = GetEffA(material);
  density = material->GetDensity()/(g/cm3);
  X0  = material->GetRadlen(); 
  Ec      = 2.66 * std::pow((X0 * Z / A),1.1); 
  G4double Es = 21*MeV;
  Rm = X0*Es/Ec;
  // PrintMaterial(); 
}

---

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

• GFlashHomoShowerParameterisation.hh
• GFlashHomoShowerParameterisation.cc