ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
EventAction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file EventAction.cc
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 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collaboration publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // and papers
31 // M. Batmunkh et al. J Radiat Res Appl Sci 8 (2015) 498-507
32 // O. Belov et al. Physica Medica 32 (2016) 1510-1520
33 // The Geant4-DNA web site is available at http://geant4-dna.org
34 //
35 // -------------------------------------------------------------------
36 // November 2016
37 // -------------------------------------------------------------------
38 //
39 //
42 
43 #include "G4Event.hh"
44 #include "Randomize.hh"
45 #include "EventAction.hh"
46 #include "RunAction.hh"
47 #include "G4SystemOfUnits.hh"
48 #include "G4ios.hh"
49 
50 #include "Analysis.hh"
51 #include "Run.hh"
52 #include "G4RunManager.hh"
53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
54 
55 EventAction::EventAction(RunAction* run)
56 :fRunAction(run)
57 {
58  remove ("OutputPerEvent.out");
59 }
60 
61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
62 
63 EventAction::~EventAction()
64 {
65 
66 }
67 
68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
69 
70 void EventAction::BeginOfEventAction(const G4Event* evt)
71 {
72  G4int evtNb = evt->GetEventID();
73  fRunAction->SetNumEvent(evtNb);
74  //
75  fRunAction->SetEdepALL(0);
76  fRunAction->SetEdepMedium(0);
77  fRunAction->SetEdepSlice(0);
78  fRunAction->SetEdepNeuron(0);
79  fRunAction->SetEdepSoma(0);
80  fRunAction->SetEdepDend(0);
81  fRunAction->SetEdepAxon(0);
82 }
83 
84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
85 
86 void EventAction::EndOfEventAction(const G4Event* evt)
87 {
88  //Get Postion and Momentum of primary
89  // beam index of per track
90  G4int evtNb = evt->GetEventID();
91  // const G4ThreeVector& pos = evt->GetPrimaryVertex()->GetPosition();
92  // const G4ThreeVector& mom = evt->GetPrimaryVertex()
93  // ->GetPrimary()->GetMomentum();
94 
95  Run* run = static_cast<Run*>(
96  G4RunManager::GetRunManager()->GetNonConstCurrentRun());
97  run->AddEdepALL(fRunAction->GetEdepALL());
98  run->AddEdepMedium(fRunAction->GetEdepMedium());
99  run->AddEdepSlice(fRunAction->GetEdepSlice());
100  run->AddEdepNeuron(fRunAction->GetEdepNeuron());
101  run->AddEdepSoma(fRunAction->GetEdepSoma());
102  run->AddEdepDend(fRunAction->GetEdepDend());
103  run->AddEdepAxon(fRunAction->GetEdepAxon());
104  // to calculate LET
105  /*
106  if (fRunAction->GetEdepALL() > 0.)
107  {
108  G4double LET = fRunAction->GetEdepALL()
109  / ftrackLength() ;
110  run->AddPrimaryLET(LET);
111  }
112  */
113 /* G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
114  analysisManager->FillNtupleDColumn(4,0,fRunAction->GetEdepALL()/keV );
115  analysisManager->FillNtupleDColumn(4,1,fRunAction->GetEdepMedium()/keV );
116  analysisManager->FillNtupleDColumn(4,2,
117  (fRunAction->GetEdepSlice()+fRunAction->GetEdepNeuron())/keV);
118  analysisManager->FillNtupleDColumn(4,3,fRunAction->GetEdepNeuron()/keV );
119  analysisManager->AddNtupleRow(4);
120 */
121  std::ofstream OutputEdep("OutputPerEvent.out", std::ios::app);
122  OutputEdep<< evtNb+1 << '\t' << " " // event number
123  // edep in all volume
124  << fRunAction->GetEdepALL()/keV << '\t' << " "
125  // outside bounding box
126  << fRunAction->GetEdepMedium()/keV << '\t' << " "
127  << (fRunAction->GetEdepSlice()+fRunAction->GetEdepNeuron())/keV
128  << '\t' << " "// inside Bounding Slice
129  << fRunAction->GetEdepNeuron()/keV << '\t' << " "
130  << fRunAction->GetEdepSoma()/keV << '\t' << " "
131  << fRunAction->GetEdepDend()/keV << '\t' << " "
132  << fRunAction->GetEdepAxon()/keV << '\t' << " "
133  << G4endl;
134 
135  /*
136  if (fRunAction->GetEdepNeuron() > 0.)
137  {
138  std::ofstream WriteDataInNeurN("0-Neur-EventN.out", std::ios::app);
139  WriteDataInNeurN << evtNb+1 << '\t' << " " //
140  << pos.x()/um << '\t' << " "
141  << pos.y()/um << '\t' << " "
142  << pos.z()/um << '\t' << " "
143  << G4endl;
144  }
145  */
146 }