51 #include "G4HCofThisEvent.hh"
52 #include "G4VHitsCollection.hh"
69 :trackerCollID(-1),calorimeterCollID(-1),
70 anticoincidenceCollID(-1), drawFlag(
"all"),
71 theRunAction(runAction)
112 G4Exception(
"GammaRayTelEventAction::BeginOfEventAction()",
113 "GTR0001",
FatalException,
"Null pointer to Run Action: this should not be");
140 G4cout <<
"Number of tracker hits in this event = " << n_hit <<
G4endl;
142 G4int NStrip, NPlane, IsX;
146 for (
int i=0;i<n_hit;i++)
150 ESil = (*THC)[i]->GetEdepSil();
151 NStrip = (*THC)[i]->GetNStrip();
152 NPlane = (*THC)[i]->GetNSilPlane();
153 IsX = (*THC)[i]->GetPlaneType();
158 (*outFile) << std::setw(7) << event_id <<
" " <<
159 ESil/
keV <<
" " << NStrip <<
160 " " << NPlane <<
" " << IsX <<
" " <<
161 (*THC)[i]->GetPos().x()/
mm <<
" "<<
162 (*THC)[i]->GetPos().y()/
mm <<
" "<<
163 (*THC)[i]->GetPos().z()/
mm <<
" "<<
166 G4cout << std::setw(7) << event_id <<
" " <<
167 ESil/
keV <<
" " << NStrip <<
168 " " << NPlane <<
" " << IsX <<
" " <<
169 (*THC)[i]->GetPos().x()/
mm <<
" "<<
170 (*THC)[i]->GetPos().y()/
mm <<
" "<<
171 (*THC)[i]->GetPos().z()/
mm <<
" "<<
197 (*THC)[i]->GetPos().
y()/
mm,
198 (*THC)[i]->GetPos().
z()/
mm);
221 G4int NStrip, NPlane, IsX;
222 for (
G4int i=0;i<n_digi;i++) {
225 NStrip = (*DC)[i]->GetStripNumber();
226 NPlane = (*DC)[i]->GetPlaneNumber();
227 IsX = (*DC)[i]->GetPlaneType();
229 (*outFile) << std::setw(7) << event_id <<
" " << NStrip <<
230 " " << NPlane <<
" " << IsX <<
" " <<
G4endl;