ExN04StackingAction Class Reference

#include <Doxymodules_eventgenerator.h>

Inheritance diagram for ExN04StackingAction:

Public Member Functions
	ExN04StackingAction ()
virtual	~ExN04StackingAction ()
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

Private Member Functions
G4bool	InsideRoI (const G4Track *aTrack, G4double ang)
G4VHitsCollection *	GetCollection (G4String colName)

Private Attributes
ExN04TrackerHitsCollection *	fTrkHits
ExN04MuonHitsCollection *	fMuonHits
ExN04StackingActionMessenger *	fMessenger
G4int	fStage
G4int	fReqMuon
G4int	fReqIsoMuon
G4int	fReqIso
G4double	fAngRoI

Detailed Description

Definition at line 48 of file Doxymodules_eventgenerator.h.

Constructor & Destructor Documentation

ExN04StackingAction()

ExN04StackingAction::ExN04StackingAction

(

)

Definition at line 44 of file ExN04StackingAction.cc.

 : G4UserStackingAction(),
   fTrkHits(0), fMuonHits(0), fStage(0)
{
  fAngRoI = 30.0*deg;
  fReqMuon = 2;
  fReqIso = 10;
  fMessenger = new ExN04StackingActionMessenger(this);
}

~ExN04StackingAction()

ExN04StackingAction::~ExN04StackingAction ( )

virtual

Definition at line 55 of file ExN04StackingAction.cc.

{ delete fMessenger; }

Member Function Documentation

ClassifyNewTrack()

G4ClassificationOfNewTrack ExN04StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Definition at line 60 of file ExN04StackingAction.cc.

{
  G4ClassificationOfNewTrack classification = fWaiting;
  switch(fStage)
  {
  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((fAngRoI<0.)||InsideRoI(aTrack,fAngRoI))
    {
      classification = fUrgent;
      break;
    }
    classification = fKill;
  }
  return classification;
}

NewStage()

void ExN04StackingAction::NewStage ( )

virtual

Definition at line 139 of file ExN04StackingAction.cc.

{
  fStage++;
  G4int nhits;
  if(fStage==1)
  {
  // Stage 0->1 : check if at least "fReqMuon" hits on muon chamber
  //              otherwise abort current event
    if(!fMuonHits)
    { fMuonHits = (ExN04MuonHitsCollection*)GetCollection("muonCollection"); }
    if(!fMuonHits)
    { G4cerr << "muonCollection NOT FOUND" << G4endl;
      return; }
    nhits = fMuonHits->entries();
    G4cout << "Stage 0->1 : " << nhits << " hits found in the muon chamber."
         << G4endl;
    if(nhits<fReqMuon)
    {
      stackManager->clear();
      G4cout << "++++++++ event aborted" << G4endl;
      return;
    }
    stackManager->ReClassify();
    return;
  }
 
  else if(fStage==2)
  {
  // Stage 1->2 : check the isolation of muon tracks
  //              at least "fReqIsoMuon" isolated muons
  //              otherwise abort current event.
  //              Isolation requires "fReqIso" or less hits
  //              (including own hits) in the RoI region
  //              in the tracker layers.
    nhits = fMuonHits->entries();
    if(!fTrkHits)
    { fTrkHits =
      (ExN04TrackerHitsCollection*)GetCollection("trackerCollection"); }
    if(!fTrkHits)
    { G4cerr << "trackerCollection NOT FOUND" << G4endl;
      return; }
    G4int nTrkhits = fTrkHits->entries();
    G4int isoMuon = 0;
    for(G4int j=0;j<nhits;j++)
    {
      G4ThreeVector hitPos = (*fMuonHits)[j]->GetPos();
      G4int nhitIn = 0;
      for(G4int jj=0;(jj<nTrkhits)&&(nhitIn<=fReqIso);jj++)
      {
        G4ThreeVector trkhitPos = (*fTrkHits)[jj]->GetPos();
        if(trkhitPos.angle(hitPos)<fAngRoI) nhitIn++;
      }
      if(nhitIn<=fReqIso) isoMuon++;
    }
    G4cout << "Stage 1->2 : " << isoMuon << " isolated muon found." << G4endl;
    if(isoMuon<fReqIsoMuon)
    {
      stackManager->clear();
      G4cout << "++++++++ event aborted" << G4endl;
      return;
    }
    stackManager->ReClassify();
    return;
  }
 
  else
  {
  // Other fStage change : just re-classify
    stackManager->ReClassify();
  }
}

PrepareNewEvent()

void ExN04StackingAction::PrepareNewEvent ( )

virtual

Definition at line 212 of file ExN04StackingAction.cc.

{
  fStage = 0;
  fTrkHits = 0;
  fMuonHits = 0;
}

SetNRequestMuon()

void ExN04StackingAction::SetNRequestMuon ( G4int  val )

inline

Definition at line 52 of file ExN04StackingAction.hh.

{ fReqMuon = val; }

GetNRequestMuon()

G4int ExN04StackingAction::GetNRequestMuon ( ) const

inline

Definition at line 53 of file ExN04StackingAction.hh.

{ return fReqMuon; }

SetNRequestIsoMuon()

void ExN04StackingAction::SetNRequestIsoMuon ( G4int  val )

inline

Definition at line 54 of file ExN04StackingAction.hh.

{ fReqIsoMuon = val; }

GetNRequestIsoMuon()

G4int ExN04StackingAction::GetNRequestIsoMuon ( ) const

inline

Definition at line 55 of file ExN04StackingAction.hh.

{ return fReqIsoMuon; }

SetNIsolation()

void ExN04StackingAction::SetNIsolation ( G4int  val )

inline

Definition at line 56 of file ExN04StackingAction.hh.

{ fReqIso = val; }

GetNIsolation()

G4int ExN04StackingAction::GetNIsolation ( ) const

inline

Definition at line 57 of file ExN04StackingAction.hh.

{ return fReqIso; }

SetRoIAngle()

void ExN04StackingAction::SetRoIAngle ( G4double  val )

inline

Definition at line 58 of file ExN04StackingAction.hh.

{ fAngRoI = val; }

GetRoIAngle()

G4double ExN04StackingAction::GetRoIAngle ( ) const

inline

Definition at line 59 of file ExN04StackingAction.hh.

{ return fAngRoI; }

InsideRoI()

G4bool ExN04StackingAction::InsideRoI ( const G4Track *  aTrack, G4double  ang  )

private

Definition at line 102 of file ExN04StackingAction.cc.

{
  if(!fMuonHits)
  { fMuonHits = (ExN04MuonHitsCollection*)GetCollection("muonCollection"); }
  if(!fMuonHits)
  { G4cerr << "muonCollection NOT FOUND" << G4endl;
    return true; }
 
  G4int nhits = fMuonHits->entries();
 
  const G4ThreeVector trPos = aTrack->GetPosition();
  for(G4int i=0;i<nhits;i++)
  {
    G4ThreeVector muHitPos = (*fMuonHits)[i]->GetPos();
    G4double angl = muHitPos.angle(trPos);
    if(angl<ang) { return true; }
  }
 
  return false;
}

GetCollection()

G4VHitsCollection * ExN04StackingAction::GetCollection ( G4String  colName )

private

Definition at line 124 of file ExN04StackingAction.cc.

{
  G4SDManager* SDMan = G4SDManager::GetSDMpointer();
  G4RunManager* runMan = G4RunManager::GetRunManager();
  int colID = SDMan->GetCollectionID(colName);
  if(colID>=0)
  {
    const G4Event* currentEvent = runMan->GetCurrentEvent();
    G4HCofThisEvent* HCE = currentEvent->GetHCofThisEvent();
    return HCE->GetHC(colID);
  }
  return 0;
}

Member Data Documentation

fTrkHits

ExN04TrackerHitsCollection* ExN04StackingAction::fTrkHits

private

Definition at line 65 of file ExN04StackingAction.hh.

fMuonHits

ExN04MuonHitsCollection* ExN04StackingAction::fMuonHits

private

Definition at line 66 of file ExN04StackingAction.hh.

fMessenger

ExN04StackingActionMessenger* ExN04StackingAction::fMessenger

private

Definition at line 67 of file ExN04StackingAction.hh.

fStage

G4int ExN04StackingAction::fStage

private

Definition at line 69 of file ExN04StackingAction.hh.

fReqMuon

G4int ExN04StackingAction::fReqMuon

private

Definition at line 70 of file ExN04StackingAction.hh.

fReqIsoMuon

G4int ExN04StackingAction::fReqIsoMuon

private

Definition at line 71 of file ExN04StackingAction.hh.

fReqIso

G4int ExN04StackingAction::fReqIso

private

Definition at line 72 of file ExN04StackingAction.hh.

fAngRoI

G4double ExN04StackingAction::fAngRoI

private

Definition at line 73 of file ExN04StackingAction.hh.

---

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

• .doxygen/Doxymodules_eventgenerator.h
• extended/eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh
• extended/eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc