| 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 runAndEvent/RE01/src/RE01EventAction
27 /// \brief Implementation of the RE01EventActi
28 //
29 //
30
31 #include "RE01EventAction.hh"
32
33 #include "RE01CalorimeterHit.hh"
34 #include "RE01TrackerHit.hh"
35 #include "RE01Trajectory.hh"
36
37 #include "G4Event.hh"
38 #include "G4EventManager.hh"
39 #include "G4HCofThisEvent.hh"
40 #include "G4PrimaryParticle.hh"
41 #include "G4PrimaryVertex.hh"
42 #include "G4SDManager.hh"
43 #include "G4SystemOfUnits.hh"
44 #include "G4TrajectoryContainer.hh"
45 #include "G4UImanager.hh"
46 #include "G4VHitsCollection.hh"
47 #include "G4VVisManager.hh"
48 #include "G4ios.hh"
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51 RE01EventAction::RE01EventAction() : G4UserEve
52 {
53 ;
54 }
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57 RE01EventAction::~RE01EventAction()
58 {
59 ;
60 }
61
62 //....oooOO0OOooo........oooOO0OOooo........oo
63 void RE01EventAction::BeginOfEventAction(const
64 {
65 G4SDManager* SDman = G4SDManager::GetSDMpoin
66 if (fTrackerCollID < 0 || fCalorimeterCollID
67 G4String colNam;
68 fTrackerCollID = SDman->GetCollectionID(co
69 fCalorimeterCollID = SDman->GetCollectionI
70 }
71 }
72
73 //....oooOO0OOooo........oooOO0OOooo........oo
74 void RE01EventAction::EndOfEventAction(const G
75 {
76 G4cout << ">>> Summary of Event " << evt->Ge
77 if (evt->GetNumberOfPrimaryVertex() == 0) {
78 G4cout << "Event is empty." << G4endl;
79 return;
80 }
81
82 if (fTrackerCollID < 0 || fCalorimeterCollID
83
84 G4HCofThisEvent* HCE = evt->GetHCofThisEvent
85 RE01TrackerHitsCollection* THC = 0;
86 RE01CalorimeterHitsCollection* CHC = 0;
87 if (HCE) {
88 THC = (RE01TrackerHitsCollection*)(HCE->Ge
89 CHC = (RE01CalorimeterHitsCollection*)(HCE
90 }
91
92 if (THC) {
93 int n_hit = THC->entries();
94 G4cout << G4endl;
95 G4cout << "Tracker hits "
96 << "-------------------------------
97 G4cout << n_hit << " hits are stored in RE
98 if (fpEventManager->GetVerboseLevel() > 0)
99 G4cout << "List of hits in tracker" << G
100 for (int i = 0; i < n_hit; i++) {
101 (*THC)[i]->Print();
102 }
103 }
104 }
105 if (CHC) {
106 int n_hit = CHC->entries();
107 G4cout << G4endl;
108 G4cout << "Calorimeter hits "
109 << "-------------------------------
110 G4cout << n_hit << " hits are stored in RE
111 G4double totE = 0;
112 for (int i = 0; i < n_hit; i++) {
113 totE += (*CHC)[i]->GetEdep();
114 }
115 G4cout << " Total energy deposition in
116 }
117
118 // get number of stored trajectories
119 G4TrajectoryContainer* trajectoryContainer =
120 G4int n_trajectories = 0;
121 if (trajectoryContainer) n_trajectories = tr
122 // extract the trajectories and print them o
123 G4cout << G4endl;
124 G4cout << "Trajectories in tracker "
125 << "---------------------------------
126 if (fpEventManager->GetVerboseLevel() > 0) {
127 for (G4int i = 0; i < n_trajectories; i++)
128 RE01Trajectory* trj = (RE01Trajectory*)(
129 trj->ShowTrajectory();
130 }
131 }
132
133 G4cout << G4endl;
134 G4cout << "Primary particles "
135 << "---------------------------------
136 G4int n_vertex = evt->GetNumberOfPrimaryVert
137 for (G4int iv = 0; iv < n_vertex; iv++) {
138 G4PrimaryVertex* pv = evt->GetPrimaryVerte
139 G4cout << G4endl;
140 G4cout << "Primary vertex " << G4ThreeVect
141 << " at t = " << (pv->GetT0()) /
142 if (fpEventManager->GetVerboseLevel() > 0)
143 G4PrimaryParticle* pp = pv->GetPrimary()
144 while (pp) {
145 PrintPrimary(pp, 0);
146 pp = pp->GetNext();
147 }
148 }
149 }
150 }
151
152 //....oooOO0OOooo........oooOO0OOooo........oo
153 void RE01EventAction::PrintPrimary(G4PrimaryPa
154 {
155 for (G4int ii = 0; ii <= ind; ii++) {
156 G4cout << " ";
157 }
158 G4cout << "==PDGcode " << pp->GetPDGcode() <
159 if (pp->GetG4code() != 0) {
160 G4cout << "(" << pp->GetG4code()->GetParti
161 }
162 else {
163 G4cout << "is not defined in G4";
164 }
165 G4cout << " " << pp->GetMomentum() / GeV <<
166 if (pp->GetTrackID() < 0) {
167 G4cout << G4endl;
168 }
169 else {
170 G4cout << ">>> G4Track ID " << pp->GetTrac
171 }
172
173 G4PrimaryParticle* daughter = pp->GetDaughte
174 while (daughter) {
175 PrintPrimary(daughter, ind + 1);
176 daughter = daughter->GetNext();
177 }
178 }
179