35 #include "G4HCofThisEvent.hh"
40 #include "ExN04StackingActionMessenger.hh"
42 #include "ExN04StackingAction.hh"
46 : trkHits(0), muonHits(0), stage(0)
96 classification =
fKill;
98 return classification;
113 for(
G4int i=0;i<nhits;i++)
117 if(angl<ang) {
return true; }
133 return HCE->
GetHC(colID);
153 G4cout <<
"Stage 0->1 : " << nhits <<
" hits found in the muon chamber."
182 for(
G4int j=0;j<nhits;j++)
186 for(
G4int jj=0;(jj<nTrkhits)&&(nhitIn<=
reqIso);jj++)
191 if(nhitIn<=
reqIso) isoMuon++;
193 G4cout <<
"Stage 1->2 : " << isoMuon <<
" isolated muon found." <<
G4endl;
G4int GetCollectionID(G4String colName)
static G4RunManager * GetRunManager()
G4VHitsCollection * GetCollection(G4String colName)
G4ClassificationOfNewTrack
virtual ~ExN04StackingAction()
G4double GetPDGCharge() const
G4ParticleDefinition * GetDefinition() const
static G4MuonMinus * MuonMinusDefinition()
G4HCofThisEvent * GetHCofThisEvent() const
double angle(const Hep3Vector &) const
ExN04MuonHitsCollection * muonHits
static constexpr double deg
G4bool InsideRoI(const G4Track *aTrack, G4double ang)
const G4ThreeVector & GetPosition() const
virtual G4ClassificationOfNewTrack ClassifyNewTrack(const G4Track *aTrack)
G4GLOB_DLL std::ostream G4cerr
G4VHitsCollection * GetHC(G4int i)
const G4Event * GetCurrentEvent() const
static G4SDManager * GetSDMpointer()
G4StackManager * stackManager
ExN04StackingActionMessenger * theMessenger
G4GLOB_DLL std::ostream G4cout
static G4MuonPlus * MuonPlusDefinition()
G4TrackStatus GetTrackStatus() const
G4int GetParentID() const
virtual void PrepareNewEvent()
ExN04TrackerHitsCollection * trkHits