| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 //
29 // ---------------------------------------------------------------
30 // GEANT 4 class header file
31 //
32 // GFlashHomoShowerParameterisation
33 //
34 // Class description:
35 //
36 // GFlash homogeneous shower parameterisation.
37
38 // Author: Joanna Weng - 02.2004
39 // ---------------------------------------------------------------
40 #ifndef GFlashHomoShowerParameterisation_h
41 #define GFlashHomoShowerParameterisation_h 1
42
43 #include "globals.hh"
44 #include "GVFlashHomoShowerTuning.hh"
45 #include "GVFlashShowerParameterisation.hh"
46
47 class G4Material;
48
49 class GFlashHomoShowerParameterisation : public GVFlashShowerParameterisation
50 {
51 public: // with description
52 GFlashHomoShowerParameterisation(G4Material* aMat, GVFlashHomoShowerTuning* aPar = 0);
53 ~GFlashHomoShowerParameterisation();
54
55 void ComputeRadialParameters(G4double y, G4double Tau);
56 void GenerateLongitudinalProfile(G4double Energy);
57 void ComputeZAX0EFFetc();
58
59 G4double IntegrateEneLongitudinal(G4double LongitudinalStep);
60 G4double IntegrateNspLongitudinal(G4double LongitudinalStep);
61 G4double ComputeTau(G4double LongitudinalPosition);
62
63 G4double GeneratePhi();
64 G4double GenerateRadius(G4int ispot, G4double Energy, G4double LongitudinalPosition);
65 G4double GenerateExponential(G4double Energy);
66 void SetMaterial(G4Material* mat);
67
68 inline G4double GetAveR99() { return (3.5 * Rm); }
69 inline G4double GetAveR90() { return (1.5 * Rm); } // ok
70
71 inline G4double GetAveTmx() {return (X0 * std::exp(AveLogTmaxh));}
72 inline G4double GetAveT99() {return (X0 * AveLogTmaxh/(AveLogAlphah-1.00));}
73 inline G4double GetAveT90() {return (2.5* X0*std::exp( AveLogTmaxh) );}
74
75 inline G4double GetNspot(){ return NSpot;}
76 inline G4double GetX0() { return X0; }
77 inline G4double GetEc() { return Ec; }
78 inline G4double GetRm() { return Rm; }
79
80 private:
81 G4Material* material;
82
83 // Resolution
84 G4double ConstantResolution;
85 G4double NoiseResolution;
86 G4double SamplingResolution;
87
88 // parametrization parameters
89 GVFlashHomoShowerTuning* thePar;
90
91 // Cashed parameters:
92 // Longitudinal Coefficients for a homogeneous calo
93 G4double ParAveT1;
94 G4double ParAveA1, ParAveA2, ParAveA3;
95 G4double ParSigLogT1, ParSigLogT2;
96 G4double ParSigLogA1, ParSigLogA2;
97 G4double ParRho1, ParRho2;
98
99 void ComputeLongitudinalParameters(G4double y);
100 void GenerateEnergyProfile(G4double y);
101 void GenerateNSpotProfile(G4double y);
102
103 // Radial Coefficients
104 G4double ParRC1, ParRC2, ParRC3, ParRC4;
105 G4double ParWC1, ParWC2, ParWC3;
106 G4double ParWC4, ParWC5, ParWC6;
107 G4double ParRT1, ParRT2, ParRT3, ParRT4;
108 G4double ParRT5, ParRT6;
109
110 // Spot multiplicity Coefficients
111 G4double ParSpotT1, ParSpotT2, ParSpotA1, ParSpotA2;
112 G4double ParSpotN1, ParSpotN2;
113
114 // PARAMETRISATION variables (Energy & position dependent)
115 // Longitudinal
116 // homogeneous
117 G4double AveLogAlphah, AveLogTmaxh;
118 G4double SigmaLogAlphah, SigmaLogTmaxh;
119 G4double Rhoh;
120 G4double Alphah, Tmaxh, Betah;
121
122 // Multiplicity
123 G4double NSpot, AlphaNSpot, TNSpot, BetaNSpot;
124
125 // Radial
126 G4double RadiusCore, WeightCore, RadiusTail;
127 };
128
129 #endif
130
131