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)) );
virtual G4double ParSigLogT1()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double NoiseResolution()
virtual G4double ParWC4()
virtual G4double ParRT5()
virtual G4double ParRC4()
static constexpr double MeV
virtual G4double ParWC5()
virtual G4double ParAveA2()
virtual G4double ParRC1()
static constexpr double mm
Float_t x1[n_points_granero]
G4double GetRadlen() const
virtual G4double ParSigLogA1()
G4double GenerateRadius(G4int ispot, G4double Energy, G4double LongitudinalPosition)
virtual G4double ParRT6()
G4double IntegrateEneLongitudinal(G4double LongitudinalStep)
virtual G4double ParAveT1()
void GenerateEnergyProfile(G4double y)
virtual G4double ParAveA1()
virtual G4double ParRC2()
virtual G4double ParRho2()
virtual G4double ParSpotA2()
G4double ApplySampling(const G4double DEne, const G4double Energy)
~GFlashSamplingShowerParameterisation()
virtual G4double ParSigLogA2()
G4double SamplingResolution
static constexpr double g
virtual G4double ParWC1()
virtual G4double ParRC3()
static constexpr double cm2
virtual G4double ParRho1()
G4double ConstantResolution
GFlashSamplingShowerParameterisation(G4Material *aMat1, G4Material *aMat2, G4double d1, G4double d2, GFlashSamplingShowerTuning *aPar=0)
virtual G4double ParWC3()
virtual G4double ParWC6()
void GenerateLongitudinalProfile(G4double Energy)
GFlashSamplingShowerTuning * thePar
virtual G4double ParSpotT2()
G4double ComputeTau(G4double LongitudinalPosition)
G4double ConstantResolution()
G4double gam(G4double x, G4double a) const
ThreeVector shoot(const G4int Ap, const G4int Af)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4double GetEffZ(const G4Material *material)
G4double SamplingResolution()
void GenerateNSpotProfile(G4double y)
virtual G4double ParRT4()
void ComputeRadialParameters(G4double y, G4double Tau)
virtual G4double ParSpotN1()
virtual G4double ParSpotA1()
virtual G4double ParRT2()
virtual G4double ParSpotT1()
virtual G4double ParSpotN2()
static constexpr double cm
G4GLOB_DLL std::ostream G4cout
virtual G4double ParAveA3()
virtual G4double ParRT3()
G4double GenerateExponential(G4double Energy)
virtual G4double ParRT1()
void SetMaterial(G4Material *mat1, G4Material *mat2)
virtual G4double ParSigLogT2()
Float_t x2[n_points_geant4]
static constexpr double cm3
G4double GetEffA(const G4Material *material)
virtual G4double ParWC2()
static constexpr double GeV
G4double IntegrateNspLongitudinal(G4double LongitudinalStep)
void ComputeLongitudinalParameters(G4double y)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetDensity() const