27 #include "ExN04StackingAction.hh"
31 #include "G4HCofThisEvent.hh"
36 #include "ExN04StackingActionMessenger.hh"
41 : trkHits(0), muonHits(0), stage(0)
50 {
delete theMessenger; }
84 if((angRoI<0.)||InsideRoI(aTrack,angRoI))
89 classification =
fKill;
91 return classification;
105 for(
G4int i=0;i<nhits;i++)
109 if(angl<ang) {
return true; }
124 return HCE->
GetHC(colID);
143 G4cout <<
"Stage 0->1 : " << nhits <<
" hits found in the muon chamber."
171 for(
G4int j=0;j<nhits;j++)
175 for(
G4int jj=0;(jj<nTrkhits)&&(nhitIn<=reqIso);jj++)
178 if(trkhitPos.
angle(hitPos)<angRoI) nhitIn++;
180 if(nhitIn<=reqIso) isoMuon++;
182 G4cout <<
"Stage 1->2 : " << isoMuon <<
" isolated muon found." <<
G4endl;
183 if(isoMuon<reqIsoMuon)