ExN04StackingAction Class Reference

#include <ExN04StackingAction.hh>

Inheritance diagram for ExN04StackingAction:

Collaboration 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
	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
Public Member Functions inherited from G4UserStackingAction
	G4UserStackingAction ()
virtual	~G4UserStackingAction ()
void	SetStackManager (G4StackManager *value)

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

Private Attributes
ExN04TrackerHitsCollection *	trkHits
ExN04MuonHitsCollection *	muonHits
ExN04StackingActionMessenger *	theMessenger
G4int	stage
G4int	reqMuon
G4int	reqIsoMuon
G4int	reqIso
G4double	angRoI

Additional Inherited Members
Protected Attributes inherited from G4UserStackingAction
G4StackManager *	stackManager

Detailed Description

Definition at line 44 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

Constructor & Destructor Documentation

ExN04StackingAction() [1/2]

ExN04StackingAction::ExN04StackingAction

(

)

Definition at line 45 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

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

~ExN04StackingAction() [1/2]

ExN04StackingAction::~ExN04StackingAction ( )

virtual

Definition at line 55 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

{ delete theMessenger; }

ExN04StackingAction() [2/2]

ExN04StackingAction::ExN04StackingAction

(

)

~ExN04StackingAction() [2/2]

virtual ExN04StackingAction::~ExN04StackingAction ( )

virtual

Member Function Documentation

ClassifyNewTrack() [1/2]

G4ClassificationOfNewTrack ExN04StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 60 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

{
  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;
}

Here is the call graph for this function:

ClassifyNewTrack() [2/2]

virtual G4ClassificationOfNewTrack ExN04StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Reimplemented from G4UserStackingAction.

GetCollection() [1/2]

G4VHitsCollection* ExN04StackingAction::GetCollection ( G4String  colName )

private

GetCollection() [2/2]

G4VHitsCollection * ExN04StackingAction::GetCollection ( G4String  colName )

private

Definition at line 124 of file eventgenerator/HepMC/HepMCEx01/src/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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNIsolation() [1/2]

G4int ExN04StackingAction::GetNIsolation ( ) const

inline

Definition at line 58 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ return reqIso; }

Here is the caller graph for this function:

GetNIsolation() [2/2]

G4int ExN04StackingAction::GetNIsolation ( ) const

inline

Definition at line 76 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ return reqIso; }

GetNRequestIsoMuon() [1/2]

G4int ExN04StackingAction::GetNRequestIsoMuon ( ) const

inline

Definition at line 56 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ return reqIsoMuon; }

Here is the caller graph for this function:

GetNRequestIsoMuon() [2/2]

G4int ExN04StackingAction::GetNRequestIsoMuon ( ) const

inline

Definition at line 74 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ return reqIsoMuon; }

GetNRequestMuon() [1/2]

G4int ExN04StackingAction::GetNRequestMuon ( ) const

inline

Definition at line 54 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ return reqMuon; }

Here is the caller graph for this function:

GetNRequestMuon() [2/2]

G4int ExN04StackingAction::GetNRequestMuon ( ) const

inline

Definition at line 72 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ return reqMuon; }

GetRoIAngle() [1/2]

G4double ExN04StackingAction::GetRoIAngle ( ) const

inline

Definition at line 60 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ return angRoI; }

Here is the call graph for this function:

Here is the caller graph for this function:

GetRoIAngle() [2/2]

G4double ExN04StackingAction::GetRoIAngle ( ) const

inline

Definition at line 78 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ return angRoI; }

InsideRoI() [1/2]

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

private

InsideRoI() [2/2]

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

private

Definition at line 102 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

{
  if(!muonHits)
  { muonHits = (ExN04MuonHitsCollection*)GetCollection("muonCollection"); }
  if(!muonHits)
  { G4cerr << "muonCollection NOT FOUND" << G4endl;
    return true; }

  G4int nhits = muonHits->entries();

  const G4ThreeVector trPos = aTrack->GetPosition();
  for(G4int i=0;i<nhits;i++)
  {
    G4ThreeVector muHitPos = (*muonHits)[i]->GetPos();
    G4double angl = muHitPos.angle(trPos);
    if(angl<ang) { return true; }
  }

  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

NewStage() [1/2]

void ExN04StackingAction::NewStage ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 139 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

{
  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 = (ExN04MuonHitsCollection*)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 =
      (ExN04TrackerHitsCollection*)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();
  }
}

Here is the call graph for this function:

NewStage() [2/2]

virtual void ExN04StackingAction::NewStage ( )

virtual

Reimplemented from G4UserStackingAction.

PrepareNewEvent() [1/2]

void ExN04StackingAction::PrepareNewEvent ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 212 of file eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc.

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

PrepareNewEvent() [2/2]

virtual void ExN04StackingAction::PrepareNewEvent ( )

virtual

Reimplemented from G4UserStackingAction.

SetNIsolation() [1/2]

void ExN04StackingAction::SetNIsolation ( G4int  val )

inline

Definition at line 57 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ reqIso = val; }

Here is the caller graph for this function:

SetNIsolation() [2/2]

void ExN04StackingAction::SetNIsolation ( G4int  val )

inline

Definition at line 75 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ reqIso = val; }

SetNRequestIsoMuon() [1/2]

void ExN04StackingAction::SetNRequestIsoMuon ( G4int  val )

inline

Definition at line 55 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ reqIsoMuon = val; }

Here is the caller graph for this function:

SetNRequestIsoMuon() [2/2]

void ExN04StackingAction::SetNRequestIsoMuon ( G4int  val )

inline

Definition at line 73 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ reqIsoMuon = val; }

SetNRequestMuon() [1/2]

void ExN04StackingAction::SetNRequestMuon ( G4int  val )

inline

Definition at line 53 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ reqMuon = val; }

Here is the caller graph for this function:

SetNRequestMuon() [2/2]

void ExN04StackingAction::SetNRequestMuon ( G4int  val )

inline

Definition at line 71 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ reqMuon = val; }

SetRoIAngle() [1/2]

void ExN04StackingAction::SetRoIAngle ( G4double  val )

inline

Definition at line 59 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

{ angRoI = val; }

Here is the caller graph for this function:

SetRoIAngle() [2/2]

void ExN04StackingAction::SetRoIAngle ( G4double  val )

inline

Definition at line 77 of file parallel/TopC/ParN04/include/ExN04StackingAction.hh.

{ angRoI = val; }

Member Data Documentation

angRoI

G4double ExN04StackingAction::angRoI

private

Definition at line 74 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

muonHits

ExN04MuonHitsCollection * ExN04StackingAction::muonHits

private

Definition at line 67 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

reqIso

G4int ExN04StackingAction::reqIso

private

Definition at line 73 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

reqIsoMuon

G4int ExN04StackingAction::reqIsoMuon

private

Definition at line 72 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

reqMuon

G4int ExN04StackingAction::reqMuon

private

Definition at line 71 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

stage

G4int ExN04StackingAction::stage

private

Definition at line 70 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

theMessenger

ExN04StackingActionMessenger * ExN04StackingAction::theMessenger

private

Definition at line 68 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

trkHits

ExN04TrackerHitsCollection * ExN04StackingAction::trkHits

private

Definition at line 66 of file eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh.

---

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

• eventgenerator/HepMC/HepMCEx01/include/ExN04StackingAction.hh
• eventgenerator/HepMC/HepMCEx01/src/ExN04StackingAction.cc