52 ParAveT2(0.), ParSigLogT1(0.), ParSigLogT2(0.),
53 ParSigLogA1(0.), ParSigLogA2(0.), ParRho1(0.), ParRho2(0.), ParsAveA2(0.),
54 AveLogAlphah(0.), AveLogTmaxh(0.), SigmaLogAlphah(0.), SigmaLogTmaxh(0.),
55 Rhoh(0.), Alphah(0.), Tmaxh(0.), Betah(0.), AveLogAlpha(0.), AveLogTmax(0.),
56 SigmaLogAlpha(0.), SigmaLogTmax(0.), Rho(0.),
Alpha(0.), Tmax(0.), Beta(0.)
136 G4cout <<
"/********************************************/ " <<
G4endl;
137 G4cout <<
" - GFlashSamplingShowerParameterisation::Constructor - " <<
G4endl;
138 G4cout <<
"/********************************************/ " <<
G4endl;
176 G4cout <<
"/************ ComputeZAX0EFFetc ************/" <<
G4endl;
177 G4cout <<
" - GFlashSamplingShowerParameterisation::Material - " <<
G4endl;
212 G4cout <<
"/********************************************/ " <<
G4endl;
222 G4Exception(
"GFlashSamplingShowerParameterisation::GenerateLongitudinalProfile()",
267 (Correlation1*Random1 + Correlation2*Random2) ));
269 (Correlation1*Random1 - Correlation2*Random2) ));
273 (Correlation1h*Random1 + Correlation2h*Random2) );
275 (Correlation1h*Random1 - Correlation2h*Random2) );
304 if(Resolution >0.0 && DEne > 0.00)
310 return DEneFluctuated;
321 float x3 =
gam(x1,x2);
391 WeightCore = p1 * std::exp( (p2-Tau)/p3- std::exp( (p2-Tau) /p3) );
399 + std::exp(k4*(Tau-k2)) );
G4double GetEffZ(const G4Material *material)
virtual G4double ParSigLogA2()
ThreeVector shoot(const G4int Ap, const G4int Af)
virtual G4double ParSpotN2()
virtual G4double ParSpotT1()
void GenerateEnergyProfile(G4double y)
GFlashSamplingShowerParameterisation(G4Material *aMat1, G4Material *aMat2, G4double d1, G4double d2, GFlashSamplingShowerTuning *aPar=0)
virtual G4double ParRho1()
virtual G4double ParWC1()
~GFlashSamplingShowerParameterisation()
virtual G4double ParAveA2()
G4double GetDensity() const
virtual G4double ParRho2()
G4double IntegrateEneLongitudinal(G4double LongitudinalStep)
virtual G4double ParRT3()
virtual G4double ParWC3()
G4double GenerateRadius(G4int ispot, G4double Energy, G4double LongitudinalPosition)
G4double SamplingResolution
virtual G4double ParSigLogT1()
void ComputeLongitudinalParameters(G4double y)
virtual G4double ParRT6()
virtual G4double ParRT5()
G4double ComputeTau(G4double LongitudinalPosition)
virtual G4double ParAveA1()
G4GLOB_DLL std::ostream G4cout
virtual G4double ParRT1()
virtual G4double ParSigLogA1()
virtual G4double ParSigLogT2()
G4double SamplingResolution()
G4double GenerateExponential(G4double Energy)
G4double NoiseResolution()
virtual G4double ParRT2()
G4double GetRadlen() const
G4double ConstantResolution
virtual G4double ParSpotA1()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
virtual G4double ParWC6()
G4double ApplySampling(const G4double DEne, const G4double Energy)
virtual G4double ParWC2()
virtual G4double ParSpotA2()
virtual G4double ParAveA3()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double ConstantResolution()
G4double gam(G4double x, G4double a) const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void GenerateNSpotProfile(G4double y)
void SetMaterial(G4Material *mat1, G4Material *mat2)
virtual G4double ParWC4()
G4double GetEffA(const G4Material *material)
virtual G4double ParSpotN1()
virtual G4double ParRC2()
virtual G4double ParAveT1()
virtual G4double ParRC4()
virtual G4double ParRC3()
virtual G4double ParSpotT2()
virtual G4double ParRT4()
void GenerateLongitudinalProfile(G4double Energy)
GFlashSamplingShowerTuning * thePar
virtual G4double ParRC1()
G4double IntegrateNspLongitudinal(G4double LongitudinalStep)
virtual G4double ParWC5()
void ComputeRadialParameters(G4double y, G4double Tau)