| 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 //
27 /// \file RE05/src/RE05EventAction.cc
28 /// \brief Implementation of the RE05EventAction class
29 //
30
31 #include "RE05EventAction.hh"
32
33 #include "RE05CalorimeterHit.hh"
34 #include "RE05MuonHit.hh"
35 #include "RE05TrackerHit.hh"
36
37 #include "G4Event.hh"
38 #include "G4EventManager.hh"
39 #include "G4HCofThisEvent.hh"
40 #include "G4SDManager.hh"
41 #include "G4SystemOfUnits.hh"
42 #include "G4Trajectory.hh"
43 #include "G4TrajectoryContainer.hh"
44 #include "G4UImanager.hh"
45 #include "G4VHitsCollection.hh"
46 #include "G4VVisManager.hh"
47 #include "G4ios.hh"
48
49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51 RE05EventAction::RE05EventAction()
52 : G4UserEventAction(), fTrackerCollID(-1), fCalorimeterCollID(-1), fMuonCollID(-1)
53 {}
54
55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
56
57 RE05EventAction::~RE05EventAction() {}
58
59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60
61 void RE05EventAction::BeginOfEventAction(const G4Event*)
62 {
63 G4SDManager* SDman = G4SDManager::GetSDMpointer();
64 if (fTrackerCollID < 0 || fCalorimeterCollID < 0 || fMuonCollID < 0) {
65 G4String colNam;
66 fTrackerCollID = SDman->GetCollectionID(colNam = "trackerCollection");
67 fCalorimeterCollID = SDman->GetCollectionID(colNam = "calCollection");
68 fMuonCollID = SDman->GetCollectionID(colNam = "muonCollection");
69 }
70 }
71
72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73
74 void RE05EventAction::EndOfEventAction(const G4Event* evt)
75 {
76 G4cout << ">>> Event " << evt->GetEventID() << G4endl;
77
78 if (fTrackerCollID < 0 || fCalorimeterCollID < 0 || fMuonCollID < 0) return;
79
80 G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
81 RE05TrackerHitsCollection* THC = 0;
82 RE05CalorimeterHitsCollection* CHC = 0;
83 RE05MuonHitsCollection* MHC = 0;
84 if (HCE) {
85 THC = (RE05TrackerHitsCollection*)(HCE->GetHC(fTrackerCollID));
86 CHC = (RE05CalorimeterHitsCollection*)(HCE->GetHC(fCalorimeterCollID));
87 MHC = (RE05MuonHitsCollection*)(HCE->GetHC(fMuonCollID));
88 }
89
90 if (THC) {
91 int n_hit = THC->entries();
92 G4cout << " " << n_hit << " hits are stored in RE05TrackerHitsCollection." << G4endl;
93 }
94 if (CHC) {
95 int n_hit = CHC->entries();
96 G4cout << " " << n_hit << " hits are stored in RE05CalorimeterHitsCollection." << G4endl;
97 G4double totE = 0;
98 for (int i = 0; i < n_hit; i++) {
99 totE += (*CHC)[i]->GetEdep();
100 }
101 G4cout << " Total energy deposition in calorimeter : " << totE / GeV << " (GeV)" << G4endl;
102 }
103 if (MHC) {
104 int n_hit = MHC->entries();
105 G4cout << " " << n_hit << " hits are stored in RE05MuonHitsCollection." << G4endl;
106 }
107 }
108
109 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
110