30 #include "ExN04StackingAction.hh"
34 #include "G4HCofThisEvent.hh"
39 #include "ExN04StackingActionMessenger.hh"
43 : trkHits(0), muonHits(0), stage(0)
52 {
delete theMessenger; }
86 if((angRoI<0.)||InsideRoI(aTrack,angRoI))
91 classification =
fKill;
93 return classification;
107 for(
G4int i=0;i<nhits;i++)
111 if(angl<ang) {
return true; }
126 return HCE->
GetHC(colID);
145 G4cout <<
"Stage 0->1 : " << nhits <<
" hits found in the muon chamber."
173 for(
G4int j=0;j<nhits;j++)
177 for(
G4int jj=0;(jj<nTrkhits)&&(nhitIn<=reqIso);jj++)
180 if(trkhitPos.
angle(hitPos)<angRoI) nhitIn++;
182 if(nhitIn<=reqIso) isoMuon++;
184 G4cout <<
"Stage 1->2 : " << isoMuon <<
" isolated muon found." <<
G4endl;
185 if(isoMuon<reqIsoMuon)