Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 /// \file EventAction.cc 26 /// \file EventAction.cc 27 /// \brief Implementation of the EventAction c 27 /// \brief Implementation of the EventAction class 28 // 28 // 29 // 29 // >> 30 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 33 33 #include "EventAction.hh" 34 #include "EventAction.hh" 34 35 35 #include "HistoManager.hh" << 36 #include "Run.hh" 36 #include "Run.hh" >> 37 #include "HistoManager.hh" 37 38 38 #include "G4Event.hh" 39 #include "G4Event.hh" 39 #include "G4RunManager.hh" 40 #include "G4RunManager.hh" >> 41 #include "G4SystemOfUnits.hh" 40 #include "G4UnitsTable.hh" 42 #include "G4UnitsTable.hh" 41 43 42 //....oooOO0OOooo........oooOO0OOooo........oo 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 43 45 >> 46 EventAction::EventAction() >> 47 :G4UserEventAction(), >> 48 fEdep1(0.), fEdep2(0.), fWeight1(0.), fWeight2(0.), >> 49 fTime0(-1*s) >> 50 { } >> 51 >> 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 53 >> 54 EventAction::~EventAction() >> 55 { } >> 56 >> 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 58 44 void EventAction::BeginOfEventAction(const G4E 59 void EventAction::BeginOfEventAction(const G4Event*) 45 { 60 { 46 fEdep1 = fEdep2 = fWeight1 = fWeight2 = 0.; 61 fEdep1 = fEdep2 = fWeight1 = fWeight2 = 0.; 47 fTime0 = -1 * s; << 62 fTime0 = -1*s; 48 } 63 } 49 64 50 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 51 66 52 void EventAction::AddEdep(G4int iVol, G4double << 67 void EventAction::AddEdep(G4int iVol, G4double edep, >> 68 G4double time, G4double weight) 53 { 69 { 54 // initialize t0 70 // initialize t0 55 if (fTime0 < 0.) fTime0 = time; 71 if (fTime0 < 0.) fTime0 = time; 56 << 72 57 // out of time window ? 73 // out of time window ? 58 const G4double TimeWindow(1 * microsecond); << 74 const G4double TimeWindow (1*microsecond); 59 if (std::fabs(time - fTime0) > TimeWindow) r 75 if (std::fabs(time - fTime0) > TimeWindow) return; 60 << 76 61 if (iVol == 1) { << 77 if (iVol == 1) { fEdep1 += edep; fWeight1 += edep*weight;} 62 fEdep1 += edep; << 78 if (iVol == 2) { fEdep2 += edep; fWeight2 += edep*weight;} 63 fWeight1 += edep * weight; << 64 } << 65 if (iVol == 2) { << 66 fEdep2 += edep; << 67 fWeight2 += edep * weight; << 68 } << 69 } 79 } 70 80 71 //....oooOO0OOooo........oooOO0OOooo........oo 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 72 82 73 void EventAction::EndOfEventAction(const G4Eve 83 void EventAction::EndOfEventAction(const G4Event*) 74 { 84 { 75 G4AnalysisManager* analysisManager = G4Analy << 85 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); 76 << 86 77 G4double Etot = fEdep1 + fEdep2; << 87 G4double Etot = fEdep1 + fEdep2; 78 G4double Wtot = (fWeight1 + fWeight2) / Etot << 88 G4double Wtot = (fWeight1 + fWeight2)/Etot; 79 << 89 80 // pulse height in target << 90 // pulse height in target 81 // << 91 // 82 if (fEdep1 > 0.) { << 92 if (fEdep1 > 0.) { 83 fWeight1 /= fEdep1; << 93 fWeight1 /= fEdep1; 84 analysisManager->FillH1(0, fEdep1, fWeight << 94 analysisManager->FillH1(0, fEdep1, fWeight1); 85 } << 95 } 86 << 96 87 // pulse height in detector << 97 // pulse height in detector 88 // << 98 // 89 if (fEdep2 > 0.) { << 99 if (fEdep2 > 0.) { 90 fWeight2 /= fEdep2; << 100 fWeight2 /= fEdep2; 91 analysisManager->FillH1(1, fEdep2, fWeight << 101 analysisManager->FillH1(1, fEdep2, fWeight2); 92 } << 102 } 93 << 103 94 // total << 104 // total 95 // << 105 // 96 analysisManager->FillH1(2, Etot, Wtot); << 106 analysisManager->FillH1(2, Etot, Wtot); 97 << 107 98 // threshold in target and detector << 108 // threshold in target and detector 99 const G4double Threshold1(10 * keV), Thresho << 109 const G4double Threshold1(10*keV), Threshold2(10*keV); 100 << 110 101 // coincidence, anti-coincidences << 111 //coincidence, anti-coincidences 102 // << 112 // 103 G4bool coincidence = ((fEdep1 >= Threshold1) << 113 G4bool coincidence = ((fEdep1 >= Threshold1) && (fEdep2 >= Threshold2)); 104 G4bool anti_coincidence1 = ((fEdep1 >= Thres << 114 G4bool anti_coincidence1 = ((fEdep1 >= Threshold1) && (fEdep2 < Threshold2)); 105 G4bool anti_coincidence2 = ((fEdep1 < Thresh << 115 G4bool anti_coincidence2 = ((fEdep1 < Threshold1) && (fEdep2 >= Threshold2)); 106 << 116 107 if (coincidence) analysisManager->FillH1(3, << 117 if (coincidence) analysisManager->FillH1(3, fEdep2, fWeight2); 108 if (anti_coincidence1) analysisManager->Fill << 118 if (anti_coincidence1) analysisManager->FillH1(4, fEdep1, fWeight1); 109 if (anti_coincidence2) analysisManager->Fill << 119 if (anti_coincidence2) analysisManager->FillH1(5, fEdep2, fWeight2); 110 << 120 111 // pass energies to Run << 121 // pass energies to Run 112 // << 122 // 113 Run* run = static_cast<Run*>(G4RunManager::G << 123 Run* run = static_cast<Run*>( 114 << 124 G4RunManager::GetRunManager()->GetNonConstCurrentRun()); 115 run->AddEdep(fEdep1, fEdep2); << 125 >> 126 run->AddEdep (fEdep1, fEdep2); 116 } 127 } 117 128 118 //....oooOO0OOooo........oooOO0OOooo........oo 129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 130 >> 131 119 132