30 #include "ExN04StackingAction.hh"
34 #include "G4HCofThisEvent.hh"
39 #include "ExN04StackingActionMessenger.hh"
44 : trkHits(0), muonHits(0), stage(0)
53 {
delete theMessenger; }
87 if((angRoI<0.)||InsideRoI(aTrack,angRoI))
92 classification =
fKill;
94 return classification;
108 for(
G4int i=0;i<nhits;i++)
112 if(angl<ang) {
return true; }
127 return HCE->
GetHC(colID);
146 G4cout <<
"Stage 0->1 : " << nhits <<
" hits found in the muon chamber."
174 for(
G4int j=0;j<nhits;j++)
178 for(
G4int jj=0;(jj<nTrkhits)&&(nhitIn<=reqIso);jj++)
181 if(trkhitPos.
angle(hitPos)<angRoI) nhitIn++;
183 if(nhitIn<=reqIso) isoMuon++;
185 G4cout <<
"Stage 1->2 : " << isoMuon <<
" isolated muon found." <<
G4endl;
186 if(isoMuon<reqIsoMuon)