LArStackingAction Class Reference

#include "LArStackingAction.h"

Inheritance diagram for LArStackingAction:

Public Member Functions
	LArStackingAction (int)
virtual	~LArStackingAction ()
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 Attributes
G4int	fStack

Private Member Functions
std::string	InsideTPC (const G4Track *aTrack)

Private Attributes
G4int	fstage
G4int	freqMuon
G4int	freqIsoMuon
G4int	freqIso
G4double	fangRoI

Detailed Description

Definition at line 40 of file LArStackingAction.h.

Constructor & Destructor Documentation

LArStackingAction::LArStackingAction

(

int

)

Definition at line 52 of file LArStackingAction.cxx.

References fStack.

 : fstage(0)
 , freqMuon(2)
 , freqIsoMuon(0)
 , freqIso(10)
 , fangRoI(30.*CLHEP::deg)
{ 
  //theMessenger = new LArStackingActionMessenger(this);
  fStack = dum;
    // Positive values effect action in this routine. Negative values
    // effect action in G4BadIdeaAction.

}

LArStackingAction::~LArStackingAction ( )

virtual

Definition at line 66 of file LArStackingAction.cxx.

{ //delete theMessenger; 
}

Member Function Documentation

G4ClassificationOfNewTrack LArStackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Definition at line 71 of file LArStackingAction.cxx.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), fStack, fstage, geo::GeometryCore::GetLArTPCVolumeName(), and InsideTPC().

{
  G4ClassificationOfNewTrack classification = fWaiting;
  art::ServiceHandle<geo::Geometry> geom;
  TString volName(InsideTPC(aTrack));
  Double_t buffer = 500; // Keep muNucl neutrals within 5m (for now) of larTPC.

  // These 3 for now for investigation in gdb.
  //int pdg = aTrack->GetDefinition()->GetPDGEncoding();
  int ppdg = aTrack->GetParentID();
  TString process("NA");
  if (ppdg) process = (TString)aTrack->GetCreatorProcess()->GetProcessName();


  switch(fstage){
  case 0: // Fstage 0 : Primary muons only
    if (aTrack->GetParentID()==0)
    {
      G4ParticleDefinition *particleType = aTrack->GetDefinition();
      if( ((particleType==G4MuonPlus::MuonPlusDefinition())
          || (particleType==G4MuonMinus::MuonMinusDefinition())
           )
          && !volName.Contains("unknown")
          ){ 
        classification = fUrgent; 
      }
    }
    if (volName.Contains("unknown"))   classification = fKill;
    break;
    
  case 1: // Stage 1 : K0,Lambda,n's made urgent here.
          //           Suspended tracks will be sent to the waiting stack
    if(aTrack->GetTrackStatus()==fSuspend) { break; }

    if ((aTrack->GetDefinition()->GetPDGEncoding()==2112 || aTrack->GetDefinition()->GetPDGEncoding()==130 || aTrack->GetDefinition()->GetPDGEncoding()==310 || aTrack->GetDefinition()->GetPDGEncoding()==311 || aTrack->GetDefinition()->GetPDGEncoding()==3122 ) && (aTrack->GetParentID()==1) && !volName.Contains("unknown"))
      {
        
        const G4ThreeVector tr4Pos = aTrack->GetPosition();
        // G4 returns positions in mm, have to convert to cm for LArSoft coordinate systems
        const TVector3 trPos(tr4Pos.x()/CLHEP::cm,tr4Pos.y()/CLHEP::cm,tr4Pos.z()/CLHEP::cm);
        //double locNeut = trPos.Mag();
        classification = fUrgent; 
        // std::cout << "LArStackingAction: DetHalfWidth, Height, FullLength: " << geom->DetHalfWidth() << ", " << geom->DetHalfHeight() << ", " << geom->DetLength() << std::endl;
        
        if (
            trPos.X() < (geom->DetHalfWidth()*2.0 + buffer) && trPos.X() > (-buffer) &&
            trPos.Y() < (geom->DetHalfHeight()*2.0 + buffer) && trPos.Y() > (-geom->DetHalfHeight()*2.0 - buffer) &&
            trPos.Z() < (geom->DetLength() + buffer) && trPos.Z() > (-buffer) 
            )

          { classification = fUrgent; break; }
        // These tracks need to be "scored" cuz every now and then they
        // might get to the LAr.
        else
          { classification = fKill; break; }
        
        
      }
    
        //    if(aTrack->GetDefinition()->GetPDGCharge()==0.) { break; }
    if (volName.Contains("unknown"))   classification = fKill;
    break;

  default: 
    // Track all other Primaries. Accept all secondaries in TPC.
    // Kill muon ionization electrons outside TPC
    // ignore primaries since they have no creator process
    
    if(aTrack->GetParentID() == 0 && !volName.Contains("unknown")){ 
      classification = fUrgent;
      break;
    }
       
    if(volName.Contains(geom->GetLArTPCVolumeName()) && aTrack->GetParentID()!=0)
      { 
        classification = fUrgent;
        if (fStack & 0x4 && 
            aTrack->GetCreatorProcess()->GetProcessName().contains("muIoni")
            ) 
          {
            classification = fKill;
          }
        break;
      }
    else if (volName.Contains("unknown") ){
      classification = fKill;
      break;
    }
    // Leave this here, even though I claim we've Killed these in stage 2.
    if(aTrack->GetDefinition()->GetPDGEncoding()==11 
       && aTrack->GetCreatorProcess()->GetProcessName().contains("muIoni") )
      {
        classification = fKill;
        break;
      }
    // For now, kill every other thing, no matter where it is.
    classification = fKill;

  } // end switch

  return classification;
}

G4int LArStackingAction::GetNIsolation ( ) const

inline

Definition at line 76 of file LArStackingAction.h.

References freqIso.

{ return freqIso;     }

G4int LArStackingAction::GetNRequestIsoMuon ( ) const

inline

Definition at line 74 of file LArStackingAction.h.

References freqIsoMuon.

{ return freqIsoMuon; }

G4int LArStackingAction::GetNRequestMuon ( ) const

inline

Definition at line 72 of file LArStackingAction.h.

References freqMuon.

{ return freqMuon;    }

G4double LArStackingAction::GetRoIAngle ( ) const

inline

Definition at line 78 of file LArStackingAction.h.

References fangRoI.

{ return fangRoI;     }

std::string LArStackingAction::InsideTPC ( const G4Track *  aTrack )

private

Definition at line 175 of file LArStackingAction.cxx.

References geo::GeometryCore::VolumeName().

Referenced by ClassifyNewTrack().

{

  art::ServiceHandle<geo::Geometry> geom;
  const G4ThreeVector tr4Pos = aTrack->GetPosition();

  // G4 returns positions in mm, have to convert to cm for LArSoft coordinate systems
  const TVector3 trPos(tr4Pos.x()/CLHEP::cm,tr4Pos.y()/CLHEP::cm,tr4Pos.z()/CLHEP::cm);

  const std::string volName(geom->VolumeName(trPos));

  return volName;
}

void LArStackingAction::NewStage ( )

virtual

Definition at line 190 of file LArStackingAction.cxx.

References fstage.

{

  // Here when Urgent stack is empty. Waiting stack about to be made Urgent,
  // upon saying ReClassify().
  fstage++;

    // I yanked the ExN04's use here of stackManager->clear(), which clears stack
    // and prepares to end the event. Think I may wanna do something like this if
    // muon has been tracked, its doca is large, there are no hit voxels in the TPC,
    // and I'm in further need of optimization.

  if(fstage==1){
  // Stage 0->1 : check if at least "reqMuon" hits on muon chamber
  //              otherwise abort current event

    stackManager->ReClassify();
    return;
  }

  else if(fstage==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.
    stackManager->ReClassify();
    return;
  }

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

void LArStackingAction::PrepareNewEvent ( )

virtual

Definition at line 227 of file LArStackingAction.cxx.

References fstage.

{ 
  fstage = 0; 
  //trkHits = 0;
  //muonHits = 0;
}

void LArStackingAction::SetNIsolation ( G4int  val )

inline

Definition at line 75 of file LArStackingAction.h.

{ freqIso = val;      }

void LArStackingAction::SetNRequestIsoMuon ( G4int  val )

inline

Definition at line 73 of file LArStackingAction.h.

{ freqIsoMuon = val;  }

void LArStackingAction::SetNRequestMuon ( G4int  val )

inline

Definition at line 71 of file LArStackingAction.h.

{ freqMuon = val;     }

void LArStackingAction::SetRoIAngle ( G4double  val )

inline

Definition at line 77 of file LArStackingAction.h.

{ fangRoI = val;      }

Member Data Documentation

G4double LArStackingAction::fangRoI

private

Definition at line 66 of file LArStackingAction.h.

Referenced by GetRoIAngle().

G4int LArStackingAction::freqIso

private

Definition at line 65 of file LArStackingAction.h.

Referenced by GetNIsolation().

G4int LArStackingAction::freqIsoMuon

private

Definition at line 64 of file LArStackingAction.h.

Referenced by GetNRequestIsoMuon().

G4int LArStackingAction::freqMuon

private

Definition at line 63 of file LArStackingAction.h.

Referenced by GetNRequestMuon().

G4int LArStackingAction::fStack

Definition at line 51 of file LArStackingAction.h.

Referenced by ClassifyNewTrack(), and LArStackingAction().

G4int LArStackingAction::fstage

private

Definition at line 62 of file LArStackingAction.h.

Referenced by ClassifyNewTrack(), NewStage(), and PrepareNewEvent().

---

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

• larsim/v06_53_00/source/larsim/LArG4/LArStackingAction.h
• larsim/v06_53_00/source/larsim/LArG4/LArStackingAction.cxx