| 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 /// \file electromagnetic/TestEm2/include/Run.hh
27 /// \brief Definition of the Run class
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
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........oooOO0OOooo........oooOO0OOooo......
45
46 class DetectorConstruction;
47 class PrimaryGeneratorAction;
48
49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51 class Run : public G4Run
52 {
53 public:
54 Run(DetectorConstruction*, PrimaryGeneratorAction*);
55 ~Run() override = default;
56
57 void Merge(const G4Run*) override;
58
59 void InitializePerEvent();
60 void FillPerEvent();
61
62 inline void FillPerTrack(G4double, G4double);
63 inline void FillPerStep(G4double, G4int, G4int);
64
65 inline void AddStep(G4double q);
66
67 void EndOfRun(G4double edep, G4double rms, G4double& limit);
68
69 inline void SetVerbose(G4int val) { fVerbose = val; };
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........oooOO0OOooo........oooOO0OOooo......
106
107 inline void Run::FillPerTrack(G4double charge, G4double trkLength)
108 {
109 if (charge != 0.)
110 fChargTrLength += trkLength;
111 else
112 fNeutrTrLength += trkLength;
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
116
117 inline void Run::FillPerStep(G4double dEstep, G4int Lbin, G4int Rbin)
118 {
119 f_dEdL[Lbin] += dEstep;
120 f_dEdR[Rbin] += dEstep;
121 }
122
123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
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........oooOO0OOooo........oooOO0OOooo......
136
137 #endif
138