|
Geant4 examples
|
00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // $Id$ 00027 // 00030 00031 #include "B1EventAction.hh" 00032 00033 #include "B1RunAction.hh" 00034 #include "B1SteppingAction.hh" 00035 // use of stepping action to get and reset accumulated energy 00036 00037 #include "G4RunManager.hh" 00038 #include "G4Event.hh" 00039 00040 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00041 00042 B1EventAction* B1EventAction::fgInstance = 0; 00043 00044 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00045 00046 B1EventAction* B1EventAction::Instance() 00047 { 00048 // Static acces function via G4RunManager 00049 00050 return fgInstance; 00051 } 00052 00053 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00054 00055 B1EventAction::B1EventAction() 00056 : G4UserEventAction(), 00057 fPrintModulo(100), 00058 fEnergySum(0.), 00059 fEnergy2Sum(0.) 00060 { 00061 fgInstance = this; 00062 } 00063 00064 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00065 00066 B1EventAction::~B1EventAction() 00067 { 00068 fgInstance = 0; 00069 } 00070 00071 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00072 00073 void B1EventAction::BeginOfEventAction(const G4Event* event) 00074 { 00075 G4int eventNb = event->GetEventID(); 00076 if (eventNb%fPrintModulo == 0) { 00077 G4cout << "\n---> Begin of event: " << eventNb << G4endl; 00078 } 00079 00080 // Reset accounted energy in stepping action 00081 B1SteppingAction::Instance()->Reset(); 00082 } 00083 00084 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00085 00086 void B1EventAction::EndOfEventAction(const G4Event* /*event*/) 00087 { 00088 // accumulate statistics 00089 G4double energy = B1SteppingAction::Instance()->GetEnergy(); 00090 fEnergySum += energy; 00091 fEnergy2Sum += energy*energy; 00092 } 00093 00094 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 00095 00096 void B1EventAction::Reset() 00097 { 00098 //reset cumulative quantities 00099 // 00100 fEnergySum = 0.; 00101 fEnergy2Sum = 0.; 00102 } 00103 00104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
1.7.4