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 EventAction.cc 27 /// \brief Implementation of the EventAction c 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "EventAction.hh" 34 35 #include "HistoManager.hh" 36 #include "Run.hh" 37 38 #include "G4Event.hh" 39 #include "G4RunManager.hh" 40 #include "G4UnitsTable.hh" 41 42 //....oooOO0OOooo........oooOO0OOooo........oo 43 44 void EventAction::BeginOfEventAction(const G4E 45 { 46 fEdep1 = fEdep2 = fWeight1 = fWeight2 = 0.; 47 fTime0 = -1 * s; 48 } 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 51 52 void EventAction::AddEdep(G4int iVol, G4double 53 { 54 // initialize t0 55 if (fTime0 < 0.) fTime0 = time; 56 57 // out of time window ? 58 const G4double TimeWindow(1 * microsecond); 59 if (std::fabs(time - fTime0) > TimeWindow) r 60 61 if (iVol == 1) { 62 fEdep1 += edep; 63 fWeight1 += edep * weight; 64 } 65 if (iVol == 2) { 66 fEdep2 += edep; 67 fWeight2 += edep * weight; 68 } 69 } 70 71 //....oooOO0OOooo........oooOO0OOooo........oo 72 73 void EventAction::EndOfEventAction(const G4Eve 74 { 75 G4AnalysisManager* analysisManager = G4Analy 76 77 G4double Etot = fEdep1 + fEdep2; 78 G4double Wtot = (fWeight1 + fWeight2) / Etot 79 80 // pulse height in target 81 // 82 if (fEdep1 > 0.) { 83 fWeight1 /= fEdep1; 84 analysisManager->FillH1(0, fEdep1, fWeight 85 } 86 87 // pulse height in detector 88 // 89 if (fEdep2 > 0.) { 90 fWeight2 /= fEdep2; 91 analysisManager->FillH1(1, fEdep2, fWeight 92 } 93 94 // total 95 // 96 analysisManager->FillH1(2, Etot, Wtot); 97 98 // threshold in target and detector 99 const G4double Threshold1(10 * keV), Thresho 100 101 // coincidence, anti-coincidences 102 // 103 G4bool coincidence = ((fEdep1 >= Threshold1) 104 G4bool anti_coincidence1 = ((fEdep1 >= Thres 105 G4bool anti_coincidence2 = ((fEdep1 < Thresh 106 107 if (coincidence) analysisManager->FillH1(3, 108 if (anti_coincidence1) analysisManager->Fill 109 if (anti_coincidence2) analysisManager->Fill 110 111 // pass energies to Run 112 // 113 Run* run = static_cast<Run*>(G4RunManager::G 114 115 run->AddEdep(fEdep1, fEdep2); 116 } 117 118 //....oooOO0OOooo........oooOO0OOooo........oo 119