| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file electromagnetic/TestEm2/include/Run.
27 /// \brief Definition of the Run class
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #ifndef Run_h
34 #define Run_h 1
35
36 #include "DetectorConstruction.hh"
37
38 #include "G4AnalysisManager.hh"
39 #include "G4Run.hh"
40
41 #include <vector>
42 typedef std::vector<G4double> MyVector;
43
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 class DetectorConstruction;
47 class PrimaryGeneratorAction;
48
49 //....oooOO0OOooo........oooOO0OOooo........oo
50
51 class Run : public G4Run
52 {
53 public:
54 Run(DetectorConstruction*, PrimaryGenerato
55 ~Run() override = default;
56
57 void Merge(const G4Run*) override;
58
59 void InitializePerEvent();
60 void FillPerEvent();
61
62 inline void FillPerTrack(G4double, G4doubl
63 inline void FillPerStep(G4double, G4int, G
64
65 inline void AddStep(G4double q);
66
67 void EndOfRun(G4double edep, G4double rms,
68
69 inline void SetVerbose(G4int val) { fVerbo
70
71 private:
72 void Reset();
73
74 DetectorConstruction* fDet = nullptr;
75 PrimaryGeneratorAction* fKin = nullptr;
76
77 G4int f_nLbin = kMaxBin;
78 MyVector f_dEdL;
79 MyVector fSumELongit;
80 MyVector fSumE2Longit;
81 MyVector fSumELongitCumul;
82 MyVector fSumE2LongitCumul;
83
84 G4int f_nRbin = kMaxBin;
85 MyVector f_dEdR;
86 MyVector fSumERadial;
87 MyVector fSumE2Radial;
88 MyVector fSumERadialCumul;
89 MyVector fSumE2RadialCumul;
90
91 G4double fChargTrLength = 0.;
92 G4double fSumChargTrLength = 0.;
93 G4double fSum2ChargTrLength = 0.;
94
95 G4double fNeutrTrLength = 0.;
96 G4double fSumNeutrTrLength = 0.;
97 G4double fSum2NeutrTrLength = 0.;
98
99 G4double fChargedStep = 0.;
100 G4double fNeutralStep = 0.;
101
102 G4int fVerbose = 0;
103 };
104
105 //....oooOO0OOooo........oooOO0OOooo........oo
106
107 inline void Run::FillPerTrack(G4double charge,
108 {
109 if (charge != 0.)
110 fChargTrLength += trkLength;
111 else
112 fNeutrTrLength += trkLength;
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oo
116
117 inline void Run::FillPerStep(G4double dEstep,
118 {
119 f_dEdL[Lbin] += dEstep;
120 f_dEdR[Rbin] += dEstep;
121 }
122
123 //....oooOO0OOooo........oooOO0OOooo........oo
124
125 inline void Run::AddStep(G4double q)
126 {
127 if (q == 0.0) {
128 fNeutralStep += 1.0;
129 }
130 else {
131 fChargedStep += 1.0;
132 }
133 }
134
135 //....oooOO0OOooo........oooOO0OOooo........oo
136
137 #endif
138