| 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 //
27 /// \file optical/wls/src/WLSEventAction.cc
28 /// \brief Implementation of the WLSEventActio
29 //
30 //
31
32 #include "WLSEventAction.hh"
33
34 #include "WLSEventActionMessenger.hh"
35 #include "WLSPhotonDetHit.hh"
36 #include "WLSRun.hh"
37 #include "WLSRunAction.hh"
38
39 #include "G4AnalysisManager.hh"
40 #include "G4Event.hh"
41 #include "G4EventManager.hh"
42 #include "G4RunManager.hh"
43 #include "G4SDManager.hh"
44 #include "G4TrajectoryContainer.hh"
45 #include "G4VVisManager.hh"
46 #include "Randomize.hh"
47
48 // Purpose: Accumulates statistics regarding h
49 // in the PhotonDet detector
50
51 //....oooOO0OOooo........oooOO0OOooo........oo
52
53 WLSEventAction::WLSEventAction()
54 {
55 fEventMessenger = new WLSEventActionMessenge
56 }
57
58 //....oooOO0OOooo........oooOO0OOooo........oo
59
60 WLSEventAction::~WLSEventAction()
61 {
62 delete fEventMessenger;
63 }
64
65 //....oooOO0OOooo........oooOO0OOooo........oo
66
67 void WLSEventAction::BeginOfEventAction(const
68 {
69 fNTIR = 0;
70 fNExiting = 0;
71 fEscapedEnd = 0;
72 fEscapedMid = 0;
73 fBounce = 0;
74 fWLSBounce = 0;
75 fClad1Bounce = 0;
76 fClad2Bounce = 0;
77 fReflected = 0;
78 fEscaped = 0;
79 fMirror = 0;
80 }
81
82 //....oooOO0OOooo........oooOO0OOooo........oo
83
84 void WLSEventAction::EndOfEventAction(const G4
85 {
86 // Get Hits from the detector if any
87 G4SDManager* SDman = G4SDManager::GetSDMpoin
88 G4String colName = "PhotonDetHitCollection";
89 fMPPCCollID = SDman->GetCollectionID(colName
90
91 G4HCofThisEvent* HCE = evt->GetHCofThisEvent
92 WLSPhotonDetHitsCollection* mppcHC = nullptr
93
94 // Get the hit collections
95 if (HCE) {
96 if (fMPPCCollID >= 0) {
97 mppcHC = (WLSPhotonDetHitsCollection*)(H
98 }
99 }
100
101 // Get hit information about photons that re
102 G4int n_hit = 0;
103 if (mppcHC) {
104 n_hit = mppcHC->entries();
105 }
106
107 auto analysisManager = G4AnalysisManager::In
108 analysisManager->FillH1(2, mppcHC->entries()
109 for (size_t i = 0; i < mppcHC->entries(); ++
110 auto pdHit = (*mppcHC)[i];
111 analysisManager->FillH1(0, pdHit->GetEnerg
112 analysisManager->FillH1(1, pdHit->GetArriv
113 }
114
115 if (fVerboseLevel > 1) {
116 G4cout << "-------------------------------
117 << " In this event, number of:" <<
118 << " Exiting: " << fNExiting
119 << G4endl << " Escaped End: " <<
120 << " Bounced: " << fBounce <
121 << G4endl << " Clad1 Bounce: " <<
122 << " Clad2 Bounce: " << fClad2Bou
123 << G4endl << " Escaped: " <<
124 << G4endl << " Detector hit: " <<
125 }
126
127 auto run = static_cast<WLSRun*>(G4RunManager
128
129 run->AddTIR(fNTIR);
130 run->AddExiting(fNExiting);
131 run->AddEscapedEnd(fEscapedEnd);
132 run->AddEscapedMid(fEscapedMid);
133 run->AddBounce(fBounce);
134 run->AddClad1Bounce(fClad1Bounce);
135 run->AddClad2Bounce(fClad2Bounce);
136 run->AddReflected(fReflected);
137 run->AddEscaped(fEscaped);
138 run->AddMirror(fMirror);
139 run->AddDetectorHits(n_hit);
140 }
141
142 //....oooOO0OOooo........oooOO0OOooo........oo
143
144 G4int WLSEventAction::GetEventNo()
145 {
146 return fpEventManager->GetConstCurrentEvent(
147 }
148
149 //....oooOO0OOooo........oooOO0OOooo........oo
150
151 void WLSEventAction::SetEventVerbose(G4int lev
152 {
153 fVerboseLevel = level;
154 }
155