File indexing completed on 2025-10-26 08:45:24
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040 #ifndef GVFlashShowerParameterisation_h
0041 #define GVFlashShowerParameterisation_h 1
0042
0043 #include "globals.hh"
0044 #include "GVFlashHomoShowerTuning.hh"
0045
0046 class MyGamma;
0047 class G4Material;
0048
0049 class GVFlashShowerParameterisation
0050 {
0051 public:
0052 GVFlashShowerParameterisation();
0053 virtual ~GVFlashShowerParameterisation();
0054
0055 virtual void ComputeRadialParameters(G4double y, G4double Tau) = 0;
0056 virtual void GenerateLongitudinalProfile(G4double Energy) = 0;
0057 virtual G4double IntegrateEneLongitudinal(G4double LongitudinalStep) = 0;
0058 virtual G4double IntegrateNspLongitudinal(G4double LongitudinalStep) = 0;
0059 virtual G4double ComputeTau(G4double LongitudinalPosition) = 0;
0060 virtual G4double GenerateRadius(G4int ispot, G4double Energy,
0061 G4double LongitudinalPosition) = 0;
0062 virtual void ComputeLongitudinalParameters(G4double y) = 0;
0063 virtual void GenerateEnergyProfile(G4double y) = 0;
0064 virtual void GenerateNSpotProfile(G4double y) = 0;
0065 virtual G4double GenerateExponential(G4double Energy) = 0;
0066
0067 virtual G4double GetAveR99() = 0;
0068 virtual G4double GetAveR90() = 0;
0069
0070 virtual G4double GetAveTmx() = 0;
0071 virtual G4double GetAveT99() = 0;
0072 virtual G4double GetAveT90() = 0;
0073
0074 virtual G4double GetNspot() = 0;
0075 virtual G4double GetX0() = 0;
0076 virtual G4double GetEc() = 0;
0077 virtual G4double GetRm() = 0;
0078
0079 G4double GeneratePhi();
0080 G4double GetEffZ(const G4Material* material);
0081 G4double GetEffA(const G4Material* material);
0082 G4double gam(G4double x, G4double a) const;
0083 void PrintMaterial(const G4Material* mat);
0084
0085 protected:
0086 GVFlashHomoShowerTuning* thePar;
0087
0088 G4double density, A, Z, X0, Ec, Rm;
0089
0090 G4double NSpot;
0091
0092 private:
0093 MyGamma* fGamma;
0094 };
0095
0096 #endif