RE05StackingAction Class Reference

#include <RE05StackingAction.hh>

Inheritance diagram for RE05StackingAction:

Public Member Functions
	RE05StackingAction ()
virtual	~RE05StackingAction ()
virtual G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *aTrack)
virtual void	NewStage ()
virtual void	PrepareNewEvent ()
void	SetNRequestMuon (G4int val)
G4int	GetNRequestMuon () const
void	SetNRequestIsoMuon (G4int val)
G4int	GetNRequestIsoMuon () const
void	SetNIsolation (G4int val)
G4int	GetNIsolation () const
void	SetRoIAngle (G4double val)
G4double	GetRoIAngle () const
Public Member Functions inherited from G4UserStackingAction
	G4UserStackingAction ()
virtual	~G4UserStackingAction ()
void	SetStackManager (G4StackManager *value)

Additional Inherited Members
Protected Attributes inherited from G4UserStackingAction
G4StackManager *	stackManager

Detailed Description

Definition at line 45 of file RE05StackingAction.hh.

Constructor & Destructor Documentation

RE05StackingAction::RE05StackingAction

(

)

Definition at line 45 of file RE05StackingAction.cc.

References python.hepunit::deg.

 : trkHits(0), muonHits(0), stage(0)
{ 
  angRoI = 30.0*deg; 
  reqMuon = 2;
  reqIso = 10;
  theMessenger = new RE05StackingActionMessenger(this);
}

RE05StackingAction::~RE05StackingAction ( )

virtual

Definition at line 54 of file RE05StackingAction.cc.

{ delete theMessenger; }

Member Function Documentation

G4ClassificationOfNewTrack RE05StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 58 of file RE05StackingAction.cc.

References fKill, fSuspend, fUrgent, fWaiting, G4Track::GetDefinition(), G4Track::GetParentID(), G4ParticleDefinition::GetPDGCharge(), G4Track::GetTrackStatus(), G4MuonMinus::MuonMinusDefinition(), and G4MuonPlus::MuonPlusDefinition().

{
  G4ClassificationOfNewTrack classification = fWaiting;
  switch(stage)
  {
  case 0: // Stage 0 : Primary muons only
    if(aTrack->GetParentID()==0)
    {
      G4ParticleDefinition * particleType = aTrack->GetDefinition();
      if((particleType==G4MuonPlus::MuonPlusDefinition())
       ||(particleType==G4MuonMinus::MuonMinusDefinition()))
      { classification = fUrgent; }
    }
    break;

  case 1: // Stage 1 : Charged primaries only
          //           Suspended tracks will be sent to the waiting stack
    if(aTrack->GetParentID()!=0) { break; }
    if(aTrack->GetTrackStatus()==fSuspend) { break; }
    if(aTrack->GetDefinition()->GetPDGCharge()==0.) { break; }
    classification = fUrgent;
    break;

  default: // Stage 2 : Accept all primaries
           //           Accept all secondaries in RoI
           //           Kill secondaries outside RoI
    if(aTrack->GetParentID()==0)
    { 
      classification = fUrgent;
      break;
    }
    if((angRoI<0.)||InsideRoI(aTrack,angRoI))
    { 
      classification = fUrgent;
      break;
    }
    classification = fKill;
  }
  return classification;
}

G4int RE05StackingAction::GetNIsolation ( ) const

inline

Definition at line 76 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::GetCurrentValue().

{ return reqIso; }

G4int RE05StackingAction::GetNRequestIsoMuon ( ) const

inline

Definition at line 74 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::GetCurrentValue().

{ return reqIsoMuon; }

G4int RE05StackingAction::GetNRequestMuon ( ) const

inline

Definition at line 72 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::GetCurrentValue().

{ return reqMuon; }

G4double RE05StackingAction::GetRoIAngle ( ) const

inline

Definition at line 78 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::GetCurrentValue().

{ return angRoI; }

void RE05StackingAction::NewStage ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 134 of file RE05StackingAction.cc.

References CLHEP::Hep3Vector::angle(), G4StackManager::clear(), G4THitsCollection< T >::entries(), G4cerr, G4endl, G4StackManager::ReClassify(), and G4UserStackingAction::stackManager.

{
  stage++;
  G4int nhits;
  if(stage==1)
  {
  // Stage 0->1 : check if at least "reqMuon" hits on muon chamber
  //              otherwise abort current event
    if(!muonHits)
    { muonHits = (RE05MuonHitsCollection*)GetCollection("muonCollection"); }
    if(!muonHits)
    { G4cerr << "muonCollection NOT FOUND" << G4endl;
      return; }
    nhits = muonHits->entries();
////    G4cout << "Stage 0->1 : " << nhits << " hits found in the muon chamber."
////         << G4endl;
    if(nhits<reqMuon)
    { 
      stackManager->clear();
////      G4cout << "++++++++ event aborted" << G4endl;
      return;
    }
    stackManager->ReClassify();
    return;
  }

  else if(stage==2)
  {
  // Stage 1->2 : check the isolation of muon tracks
  //              at least "reqIsoMuon" isolated muons
  //              otherwise abort current event.
  //              Isolation requires "reqIso" or less hits
  //              (including own hits) in the RoI region
  //              in the tracker layers.
    nhits = muonHits->entries();
    if(!trkHits)
    { trkHits = (RE05TrackerHitsCollection*)GetCollection("trackerCollection"); }
    if(!trkHits)
    { G4cerr << "trackerCollection NOT FOUND" << G4endl;
      return; }
    G4int nTrkhits = trkHits->entries();
    G4int isoMuon = 0;
    for(G4int j=0;j<nhits;j++)
    {
      G4ThreeVector hitPos = (*muonHits)[j]->GetPos();
      G4int nhitIn = 0;
      for(G4int jj=0;(jj<nTrkhits)&&(nhitIn<=reqIso);jj++)
      {
        G4ThreeVector trkhitPos = (*trkHits)[jj]->GetPos();
        if(trkhitPos.angle(hitPos)<angRoI) nhitIn++;
      }
      if(nhitIn<=reqIso) isoMuon++;
    }
////    G4cout << "Stage 1->2 : " << isoMuon << " isolated muon found." << G4endl;
    if(isoMuon<reqIsoMuon)
    {
      stackManager->clear();
////      G4cout << "++++++++ event aborted" << G4endl;
      return;
    }
    stackManager->ReClassify();
    return;
  }

  else
  {
  // Other stage change : just re-classify
    stackManager->ReClassify();
  }
}

void RE05StackingAction::PrepareNewEvent ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 205 of file RE05StackingAction.cc.

{ 
  stage = 0; 
  trkHits = 0;
  muonHits = 0;
}

void RE05StackingAction::SetNIsolation ( G4int  val )

inline

Definition at line 75 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::SetNewValue().

{ reqIso = val; }

void RE05StackingAction::SetNRequestIsoMuon ( G4int  val )

inline

Definition at line 73 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::SetNewValue().

{ reqIsoMuon = val; }

void RE05StackingAction::SetNRequestMuon ( G4int  val )

inline

Definition at line 71 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::SetNewValue().

{ reqMuon = val; }

void RE05StackingAction::SetRoIAngle ( G4double  val )

inline

Definition at line 77 of file RE05StackingAction.hh.

Referenced by RE05StackingActionMessenger::SetNewValue().

{ angRoI = val; }

---

The documentation for this class was generated from the following files:

• RE05StackingAction.hh
• RE05StackingAction.cc