#include <GFlashHomoShowerParameterisation.hh>

Inheritance diagram for GFlashHomoShowerParameterisation:

Public Member Functions
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 ()

	GFlashHomoShowerParameterisation (G4Material aMat, GVFlashHomoShowerTuning aPar=0)

G4double	IntegrateEneLongitudinal (G4double LongitudinalStep)

G4double	IntegrateNspLongitudinal (G4double LongitudinalStep)

void	PrintMaterial (const G4Material *mat)

void	SetMaterial (G4Material *mat)

	~GFlashHomoShowerParameterisation ()

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	Alphah

G4double	AlphaNSpot

G4double	AveLogAlphah

G4double	AveLogTmaxh

G4double	Betah

G4double	BetaNSpot

G4double	ConstantResolution

MyGamma *	fGamma

G4Material *	material

G4double	NoiseResolution

G4double	NSpot

G4double	ParAveA1

G4double	ParAveA2

G4double	ParAveA3

G4double	ParAveT1

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	ParSigLogA1

G4double	ParSigLogA2

G4double	ParSigLogT1

G4double	ParSigLogT2

G4double	ParSpotA1

G4double	ParSpotA2

G4double	ParSpotN1

G4double	ParSpotN2

G4double	ParSpotT1

G4double	ParSpotT2

G4double	ParWC1

G4double	ParWC2

G4double	ParWC3

G4double	ParWC4

G4double	ParWC5

G4double	ParWC6

G4double	RadiusCore

G4double	RadiusTail

G4double	Rhoh

G4double	SamplingResolution

G4double	SigmaLogAlphah

G4double	SigmaLogTmaxh

GVFlashHomoShowerTuning *	thePar

G4double	Tmaxh

G4double	TNSpot

G4double	WeightCore

Detailed Description

Definition at line 49 of file GFlashHomoShowerParameterisation.hh.

Constructor & Destructor Documentation

◆ GFlashHomoShowerParameterisation()

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

Definition at line 47 of file GFlashHomoShowerParameterisation.cc.

  : 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;
  } else {
    thePar = aPar;
  }
 
  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::~GFlashHomoShowerParameterisation ( )

Definition at line 154 of file GFlashHomoShowerParameterisation.cc.

{
  delete thePar;
}

References thePar.

Member Function Documentation

◆ ComputeLongitudinalParameters()

void GFlashHomoShowerParameterisation::ComputeLongitudinalParameters ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 175 of file GFlashHomoShowerParameterisation.cc.

{
  AveLogTmaxh  = std::log(ParAveT1 + std::log(y));
    //ok  <ln T hom>
  AveLogAlphah = std::log(ParAveA1 + (ParAveA2+ParAveA3/Z)*std::log(y));
    //ok  <ln alpha hom> 
 
  SigmaLogTmaxh  = 1.00/( ParSigLogT1 + ParSigLogT2*std::log(y)) ;
    //ok sigma (ln T hom)
  SigmaLogAlphah = 1.00/( ParSigLogA1 + ParSigLogA2*std::log(y));
    //ok sigma (ln alpha hom)
  Rhoh           = ParRho1+ParRho2*std::log(y);        //ok             
}

References AveLogAlphah, AveLogTmaxh, ParAveA1, ParAveA2, ParAveA3, ParAveT1, ParRho1, ParRho2, ParSigLogA1, ParSigLogA2, ParSigLogT1, ParSigLogT2, Rhoh, SigmaLogAlphah, SigmaLogTmaxh, and GVFlashShowerParameterisation::Z.

Referenced by GenerateLongitudinalProfile().

◆ ComputeRadialParameters()

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 274 of file GFlashHomoShowerParameterisation.cc.

{
  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
}

References GeV, ParRC1, ParRC2, ParRC3, ParRC4, ParRT1, ParRT2, ParRT3, ParRT4, ParRT5, ParRT6, ParWC1, ParWC2, ParWC3, ParWC4, ParWC5, ParWC6, RadiusCore, RadiusTail, WeightCore, and GVFlashShowerParameterisation::Z.

Referenced by GenerateRadius().

◆ ComputeTau()

G4double GFlashHomoShowerParameterisation::ComputeTau ( G4double LongitudinalPosition )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 265 of file GFlashHomoShowerParameterisation.cc.

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

References Alphah, AveLogAlphah, Tmaxh, and GVFlashShowerParameterisation::X0.

Referenced by GenerateRadius().

◆ ComputeZAX0EFFetc()

void GFlashHomoShowerParameterisation::ComputeZAX0EFFetc ( )

◆ 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 IntegrateEneLongitudinal(), GFlashSamplingShowerParameterisation::IntegrateEneLongitudinal(), IntegrateNspLongitudinal(), and GFlashSamplingShowerParameterisation::IntegrateNspLongitudinal().

◆ GenerateEnergyProfile()

void GFlashHomoShowerParameterisation::GenerateEnergyProfile ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 189 of file GFlashHomoShowerParameterisation.cc.

{ 
  G4double Correlation1h = std::sqrt((1+Rhoh)/2);
  G4double Correlation2h = std::sqrt((1-Rhoh)/2);
 
  G4double Random1 = G4RandGauss::shoot();
  G4double Random2 = G4RandGauss::shoot();
 
  // 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 Alphah, AveLogAlphah, AveLogTmaxh, Betah, Rhoh, G4INCL::DeJongSpin::shoot(), SigmaLogAlphah, SigmaLogTmaxh, and Tmaxh.

Referenced by GenerateLongitudinalProfile().

◆ GenerateExponential()

G4double GFlashHomoShowerParameterisation::GenerateExponential ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 296 of file GFlashHomoShowerParameterisation.cc.

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

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

◆ GenerateLongitudinalProfile()

void GFlashHomoShowerParameterisation::GenerateLongitudinalProfile ( G4double Energy )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 159 of file GFlashHomoShowerParameterisation.cc.

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

References ComputeLongitudinalParameters(), GVFlashShowerParameterisation::Ec, FatalException, G4Exception(), GenerateEnergyProfile(), GenerateNSpotProfile(), and material.

◆ GenerateNSpotProfile()

void GFlashHomoShowerParameterisation::GenerateNSpotProfile ( G4double y )

privatevirtual

Implements GVFlashShowerParameterisation.

Definition at line 206 of file GFlashHomoShowerParameterisation.cc.

{
  TNSpot     = Tmaxh *  (ParSpotT1+ParSpotT2*Z);   // ok
  AlphaNSpot = Alphah * (ParSpotA1+ParSpotA2*Z);   
  BetaNSpot  = (AlphaNSpot-1.00)/TNSpot;           // ok 
  NSpot      = ParSpotN1 * std::log(Z)*std::pow((y*Ec)/GeV,ParSpotN2 ); // ok
}

References Alphah, AlphaNSpot, BetaNSpot, GVFlashShowerParameterisation::Ec, GeV, NSpot, ParSpotA1, ParSpotA2, ParSpotN1, ParSpotN2, ParSpotT1, ParSpotT2, Tmaxh, TNSpot, and GVFlashShowerParameterisation::Z.

Referenced by GenerateLongitudinalProfile().

◆ GeneratePhi()

G4double GFlashHomoShowerParameterisation::GeneratePhi ( )

◆ GenerateRadius()

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

virtual

Implements GVFlashShowerParameterisation.

Definition at line 237 of file GFlashHomoShowerParameterisation.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 = Rm * RadiusCore * std::sqrt( Random2/(1. - Random2) );
  }
  else
  {
    Radius = Rm * 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, GVFlashShowerParameterisation::Rm, and WeightCore.

◆ GetAveR90()

G4double GFlashHomoShowerParameterisation::GetAveR90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 72 of file GFlashHomoShowerParameterisation.hh.

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

References GVFlashShowerParameterisation::Rm.

◆ GetAveR99()

G4double GFlashHomoShowerParameterisation::GetAveR99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 71 of file GFlashHomoShowerParameterisation.hh.

71{return (3.5 * Rm);}

References GVFlashShowerParameterisation::Rm.

◆ GetAveT90()

G4double GFlashHomoShowerParameterisation::GetAveT90 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 76 of file GFlashHomoShowerParameterisation.hh.

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

References AveLogTmaxh, and GVFlashShowerParameterisation::X0.

◆ GetAveT99()

G4double GFlashHomoShowerParameterisation::GetAveT99 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 75 of file GFlashHomoShowerParameterisation.hh.

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

References AveLogAlphah, AveLogTmaxh, and GVFlashShowerParameterisation::X0.

◆ GetAveTmx()

G4double GFlashHomoShowerParameterisation::GetAveTmx ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 74 of file GFlashHomoShowerParameterisation.hh.

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

References AveLogTmaxh, and GVFlashShowerParameterisation::X0.

◆ GetEc()

G4double GFlashHomoShowerParameterisation::GetEc ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 80 of file GFlashHomoShowerParameterisation.hh.

80{return Ec;}

References GVFlashShowerParameterisation::Ec.

◆ 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 SetMaterial(), and GFlashSamplingShowerParameterisation::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 SetMaterial(), and GFlashSamplingShowerParameterisation::SetMaterial().

◆ GetNspot()

G4double GFlashHomoShowerParameterisation::GetNspot ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 78 of file GFlashHomoShowerParameterisation.hh.

78{ return NSpot;}

References NSpot.

◆ GetRm()

G4double GFlashHomoShowerParameterisation::GetRm ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 81 of file GFlashHomoShowerParameterisation.hh.

81{return Rm;}

References GVFlashShowerParameterisation::Rm.

◆ GetX0()

G4double GFlashHomoShowerParameterisation::GetX0 ( )

inlinevirtual

Implements GVFlashShowerParameterisation.

Definition at line 79 of file GFlashHomoShowerParameterisation.hh.

79{return X0;}

References GVFlashShowerParameterisation::X0.

◆ IntegrateEneLongitudinal()

G4double GFlashHomoShowerParameterisation::IntegrateEneLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 214 of file GFlashHomoShowerParameterisation.cc.

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

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

◆ IntegrateNspLongitudinal()

G4double GFlashHomoShowerParameterisation::IntegrateNspLongitudinal ( G4double LongitudinalStep )

virtual

Implements GVFlashShowerParameterisation.

Definition at line 225 of file GFlashHomoShowerParameterisation.cc.

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

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

◆ 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().

◆ SetMaterial()

void GFlashHomoShowerParameterisation::SetMaterial ( G4Material * mat )

Definition at line 141 of file GFlashHomoShowerParameterisation.cc.

{
  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(); 
}