ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RE06Run.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RE06Run.hh
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 //
28 //
29 //
30 
31 #ifndef RE06Run_h
32 #define RE06Run_h 1
33 
34 #include "globals.hh"
35 #include "G4Run.hh"
36 
37 #include "G4THitsMap.hh"
38 
39 class G4Event;
40 
41 class RE06Run : public G4Run
42 {
43 public:
44  RE06Run();
45  virtual ~RE06Run();
46 
47  virtual void RecordEvent(const G4Event*);
48  virtual void Merge(const G4Run*);
49 
50  G4double GetTotalE(G4int i) const { return GetTotal(fMapSum[i][0]); }
51  G4double GetNGamma(G4int i) const { return GetTotal(fMapSum[i][1]); }
52  G4double GetNElectron(G4int i) const { return GetTotal(fMapSum[i][2]); }
53  G4double GetNPositron(G4int i) const { return GetTotal(fMapSum[i][3]); }
54  G4double GetTotalL(G4int i) const { return GetTotal(fMapSum[i][4]); }
55  G4double GetNStep(G4int i) const { return GetTotal(fMapSum[i][5]); }
56 
57  G4double GetEMinGamma(G4int i) const { return FindMinimum(fMapMin[i][0]);}
58  G4double GetEMinElectron(G4int i) const { return FindMinimum(fMapMin[i][1]);}
59  G4double GetEMinPositron(G4int i) const { return FindMinimum(fMapMin[i][2]);}
60 
62  {
63  G4double* p = fMapPara[i][j][k];
64  if(p) return *p;
65  return 0.;
66  }
67 
68 private:
69  G4double GetTotal(const G4THitsMap<G4double> &map) const;
70  G4double FindMinimum(const G4THitsMap<G4double> &map) const;
71 
72  // Maps for accumulation
73  // fMapSum[i][j]
74  // i = 0 : Calor-A_abs j = 0 : total eDep
75  // i = 1 : Calor-A_gap j = 1 : number of gamma
76  // i = 2 : Calor-B_abs j = 2 : number of electron
77  // i = 3 : Calor-B_gap j = 3 : number of positron
78  // i = 4 : Calor-C_abs j = 4 : total step length for e+/e-
79  // i = 5 : Calor-C_gap j = 5 : total number of steps for e+/e-
82 
83  // Maps for minimum value
84  // i = 0 : Calor-A_abs j = 0 : minimum kinE at generation for gamma
85  // i = 1 : Calor-A_gap j = 1 : minimum kinE at generation for electron
86  // i = 2 : Calor-B_abs j = 2 : minimum kinE at generation for positron
87  // i = 3 : Calor-B_gap
88  // i = 4 : Calor-C_abs
89  // i = 5 : Calor-C_gap
92 
93  // Maps for accumulation in parallel world
94  // fMapPara[i][j]
95  // i = 0 : Calor-AP_para j = 0 : total eDep
96  // i = 1 : Calor-BP_para j = 1 : number of gamma
97  // i = 2 : Calor-CP_para j = 2 : number of electron
98  // j = 3 : number of positron
99  // j = 4 : total step length for e+/e-
100  // j = 5 : total number of steps for e+/e-
103 
104 };
105 
106 #endif
107