| 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 #ifndef PAR04EVENTACTION_HH
27 #define PAR04EVENTACTION_HH
28
29 #include "G4Timer.hh" // for G4Timer
30 #include "G4UserEventAction.hh" // for G4UserEventAction
31
32 #include <G4Types.hh> // for G4int, G4double
33 #include <vector> // for vector
34 class G4Event;
35 class Par04DetectorConstruction;
36 class Par04ParallelFullWorld;
37
38 /**
39 * @brief Event action class for hits' analysis.
40 *
41 * Analysis of single-particle events and developed showers in the detector.
42 * At the end of the event basic variables are calculated and saved in the
43 * histograms.
44 * Additionally ntuple with cell energies and IDs (in cylindrical coordinates) is stored.
45 *
46 */
47
48 class Par04EventAction : public G4UserEventAction
49 {
50 public:
51 Par04EventAction(Par04DetectorConstruction* aDetector, Par04ParallelFullWorld* aParallel);
52 virtual ~Par04EventAction();
53
54 /// Timer is started
55 virtual void BeginOfEventAction(const G4Event* aEvent) final;
56 /// Hits collection is retrieved, analysed, and histograms are filled.
57 virtual void EndOfEventAction(const G4Event* aEvent) final;
58 inline std::vector<G4double>& GetCalEdep() { return fCalEdep; }
59 inline std::vector<G4int>& GetCalRho() { return fCalRho; }
60 inline std::vector<G4int>& GetCalPhi() { return fCalPhi; }
61 inline std::vector<G4int>& GetCalZ() { return fCalZ; }
62 inline std::vector<G4double>& GetPhysicalCalEdep() { return fCalPhysicalEdep; }
63 inline std::vector<G4int>& GetPhysicalCalLayer() { return fCalPhysicalLayer; }
64 inline std::vector<G4int>& GetPhysicalCalSlice() { return fCalPhysicalSlice; }
65 inline std::vector<G4int>& GetPhysicalCalRow() { return fCalPhysicalRow; }
66 void StartTimer();
67 void StopTimer();
68
69 private:
70 /// ID of a hit collection to analyse
71 G4int fHitCollectionID = -1;
72 G4int fPhysicalFullHitCollectionID = -1;
73 G4int fPhysicalFastHitCollectionID = -1;
74 /// Timer measurement from Geant4
75 G4Timer fTimer;
76 /// Pointer to detector construction to retrieve (once) the detector
77 /// dimensions and size of readout
78 Par04DetectorConstruction* fDetector = nullptr;
79 Par04ParallelFullWorld* fParallel = nullptr;
80 /// Size of cell along Z axis
81 G4double fCellSizeZ = 0;
82 /// Size of cell along radius of cylinder
83 G4double fCellSizeRho = 0;
84 /// Size of cell in azimuthal angle
85 G4double fCellSizePhi = 0;
86 /// Number of readout cells along radius
87 G4int fCellNbRho = 0;
88 /// Number of readout cells in azimuthal angle
89 G4int fCellNbPhi = 0;
90 /// Number of readout cells along z axis
91 G4int fCellNbZ = 0;
92 /// Number of physical readout layers
93 G4int fPhysicalNbLayers = 0;
94 /// Number of physical readout slices
95 G4int fPhysicalNbSlices = 0;
96 /// Number of physical readout rows
97 G4int fPhysicalNbRows = 0;
98 /// Cell energy deposits to be stored in ntuple
99 std::vector<G4double> fCalEdep;
100 /// Cell ID of radius to be stored in ntuple
101 std::vector<G4int> fCalRho;
102 /// Cell ID of azimuthal angle to be stored in ntuple
103 std::vector<G4int> fCalPhi;
104 /// Cell ID of z axis to be stored in ntuple
105 std::vector<G4int> fCalZ;
106 /// Physical cell energy deposits to be stored in ntuple
107 std::vector<G4double> fCalPhysicalEdep;
108 /// Physical layer ID to be stored in ntuple
109 std::vector<G4int> fCalPhysicalLayer;
110 /// Physical slice ID to be stored in ntuple
111 std::vector<G4int> fCalPhysicalSlice;
112 /// Physical row ID to be stored in ntuple
113 std::vector<G4int> fCalPhysicalRow;
114 };
115
116 #endif /* PAR04EVENTACTION_HH */
117