33 #include "ExGflashEventAction.hh"
34 #include "ExGflashHit.hh"
52 fNevent(0),fDtime(0.0),fCalorimeterCollectionId(-1)
70 G4cout<<
" ------ Start ExGflashEventAction ----- "<<
G4endl;
81 G4cout <<
"******************************************";
89 G4cout <<
"******************************************"<<
G4endl;
106 G4cout<<
" " << n_hit<<
" hits are stored in ExGflashHitsCollection "<<
G4endl;
116 G4int hitsincrystal[100];
117 for (
int i=0;i<100;i++) energyincrystal[i]=0.;
118 for (
int i=0;i<100;i++) hitsincrystal[i]=0.;
122 for (
int i=0;i<n_hit;i++)
127 totE += (*THC)[i]->GetEdep()/
GeV;
128 G4int num=(*THC)[i]->GetCrystalNum();
130 energyincrystal[
num]+=(*THC)[i]->GetEdep()/
GeV;
131 hitsincrystal[
num]++;
150 for (
int i=0;i<100;i++)
153 if (max <energyincrystal[i])
155 max=energyincrystal[i];
163 energyincrystal[index]+energyincrystal[index+1]+energyincrystal[index-1]+
164 energyincrystal[index-10]+energyincrystal[index-9]+energyincrystal[index-11]+
165 energyincrystal[index+10]+energyincrystal[index+11]+energyincrystal[index+9];
169 energyincrystal[index]+energyincrystal[index+1]+energyincrystal[index-1]+
170 energyincrystal[index+2]+energyincrystal[index-2]+
171 energyincrystal[index-10]+energyincrystal[index-9]+energyincrystal[index-11]+
172 energyincrystal[index-8]+energyincrystal[index-12]+
173 energyincrystal[index+10]+energyincrystal[index+11]+energyincrystal[index+9]+
174 energyincrystal[index+12]+energyincrystal[index+8];
178 hitsincrystal[index]+hitsincrystal[index+1]+hitsincrystal[index-1]+
179 hitsincrystal[index-10]+hitsincrystal[index-9]+hitsincrystal[index-11]+
180 hitsincrystal[index+10]+hitsincrystal[index+11]+hitsincrystal[index+9];
184 hitsincrystal[index]+hitsincrystal[index+1]+hitsincrystal[index-1]+
185 hitsincrystal[index+2]+hitsincrystal[index-2]+
186 hitsincrystal[index-10]+hitsincrystal[index-9]+hitsincrystal[index-11]+
187 hitsincrystal[index-8]+hitsincrystal[index-12]+
188 hitsincrystal[index+10]+hitsincrystal[index+11]+hitsincrystal[index+9]+
189 hitsincrystal[index+12]+hitsincrystal[index+8];
191 G4cout <<
" e1 " << energyincrystal[index]
192 <<
" e3x3 " << e3x3<<
" GeV e5x5 " <<e5x5 <<
G4endl;
194 G4cout <<
" num1 " << hitsincrystal[index]
195 <<
" num3x3 " << num3x3<<
" num5x5 " <<num5x5 <<
G4endl;
198 G4cout <<
" Total energy deposited in the calorimeter: " << totE <<
" (GeV)" <<
G4endl;
200 G4int n_trajectories = 0;
201 if(trajectoryContainer){ n_trajectories = trajectoryContainer->
entries(); }
202 G4cout <<
" " << n_trajectories <<
" trajectories stored in this event." <<
G4endl;