52 ConstantResolution(0.), NoiseResolution(0.), SamplingResolution(0.),
53 AveLogAlphah(0.), AveLogTmaxh(0.), SigmaLogAlphah(0.), SigmaLogTmaxh(0.),
54 Rhoh(0.), Alphah(0.), Tmaxh(0.), Betah(0.)
134 G4cout <<
"/********************************************/ " <<
G4endl;
135 G4cout <<
" - GFlashHomoShowerParameterisation::Constructor - " <<
G4endl;
136 G4cout <<
"/********************************************/ " <<
G4endl;
146 Ec = 2.66 * std::pow((
X0 *
Z /
A),1.1);
162 G4Exception(
"GFlashHomoShowerParameterisation::GenerateLongitudinalProfile()",
198 (Correlation1h*Random1 + Correlation2h*Random2) );
200 (Correlation1h*Random1 - Correlation2h*Random2) );
215 G4double LongitudinalStepInX0 = LongitudinalStep /
X0;
218 float x3 =
gam(x1,x2);
226 G4double LongitudinalStepInX0 = LongitudinalStep /
X0;
253 Radius =
Rm *
RadiusCore * std::sqrt( Random2/(1. - Random2) );
257 Radius =
Rm *
RadiusTail * std::sqrt( Random2/(1. - Random2) );
283 WeightCore = p1 * std::exp( (p2-Tau)/p3 - std::exp( (p2-Tau) /p3) );
291 std::exp(k4*(Tau-k2)) );
virtual G4double ParSigLogT1()
void GenerateLongitudinalProfile(G4double Energy)
virtual G4double ParWC4()
virtual G4double ParRT5()
virtual G4double ParRC4()
static constexpr double MeV
virtual G4double ParWC5()
virtual G4double ParAveA2()
void ComputeLongitudinalParameters(G4double y)
virtual G4double ParRC1()
Float_t x1[n_points_granero]
G4double IntegrateNspLongitudinal(G4double LongitudinalStep)
G4double GetRadlen() const
G4double ComputeTau(G4double LongitudinalPosition)
virtual G4double ParSigLogA1()
void GenerateEnergyProfile(G4double y)
void GenerateNSpotProfile(G4double y)
virtual G4double ParRT6()
virtual G4double ParAveT1()
virtual G4double ParAveA1()
virtual G4double ParRC2()
virtual G4double ParRho2()
virtual G4double ParSpotA2()
virtual G4double ParSigLogA2()
~GFlashHomoShowerParameterisation()
static constexpr double g
virtual G4double ParWC1()
virtual G4double ParRC3()
virtual G4double ParRho1()
virtual G4double ParWC3()
virtual G4double ParWC6()
void ComputeRadialParameters(G4double y, G4double Tau)
void PrintMaterial(const G4Material *mat)
virtual G4double ParSpotT2()
GVFlashHomoShowerTuning * thePar
GFlashHomoShowerParameterisation(G4Material *aMat, GVFlashHomoShowerTuning *aPar=0)
G4double gam(G4double x, G4double a) const
ThreeVector shoot(const G4int Ap, const G4int Af)
G4double IntegrateEneLongitudinal(G4double LongitudinalStep)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4double GetEffZ(const G4Material *material)
virtual G4double ParRT4()
virtual G4double ParSpotN1()
virtual G4double ParSpotA1()
virtual G4double ParRT2()
virtual G4double ParSpotT1()
virtual G4double ParSpotN2()
G4GLOB_DLL std::ostream G4cout
void SetMaterial(G4Material *mat)
virtual G4double ParAveA3()
virtual G4double ParRT3()
virtual G4double ParRT1()
G4double GenerateExponential(G4double Energy)
virtual G4double ParSigLogT2()
G4double GenerateRadius(G4int ispot, G4double Energy, G4double LongitudinalPosition)
Float_t x2[n_points_geant4]
static constexpr double cm3
G4double GetEffA(const G4Material *material)
virtual G4double ParWC2()
static constexpr double GeV
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetDensity() const