numuCCBwd::FindMuonCandidateAction Class Reference

Find Muon candidate. More...

#include <numuCCBwdSelection.hxx>

Inheritance diagram for numuCCBwd::FindMuonCandidateAction:

Public Member Functions
bool	Apply (AnaEventC &event, ToyBoxB &box) const
StepBase *	MakeClone ()
	MANDATORY FUNCTIONS !!!
Public Member Functions inherited from StepBase
void	SetTitle (const std::string &title)
void	SetType (TypeEnum type)
	Set the type of this step.
void	SetBreak (bool b)
	Set whether the cut sequence should be stopped when this cut fails.
void	SetIndex (Int_t index)
	Set the index of this step.
const std::string &	Title () const
	Return the title of this step.
TypeEnum	Type () const
	Return the type of this step.
bool	Break () const
	Return true if the cut sequence should be stopped when a cut fails.
Int_t	Index () const
	Return the index of this step.
std::string	ConvertType () const
	Convert into a string the type of this step.
void	DisableInBranch (Int_t branch)
	Disable step in a given branch.
void	EnableInBranch (Int_t branch)
	Enable step in a given branch.
bool	IsDisabledInBranch (Int_t branch) const
	Returns true if the step is disabled in the specified branch.
const std::vector< StepBase * > &	GetNextSteps () const
	Returns the vector of next steps.
void	RemoveNextStep (const std::string &title, Int_t ID=0)
	Remove step with a given title (and branch ID) from next steps.
void	AddBranchUniqueID (UInt_t ibranch)
void	RemoveBranchUniqueID (UInt_t branch)
	Remove a branch unique ID from this step.
const std::vector< UInt_t > &	GetBranchUniqueIDs () const
	Returns the vector of branche unique IDs associated to this step.
void	AddNextStep (StepBase *step, Int_t branch=-1)
void	GetNextSteps (std::vector< StepBase *> &steps, bool withSplits=false)
	Get the vector of next steps.
BranchStatus	GetNextStepsInBranch (std::vector< StepBase *> &steps, Int_t branch)
	Get the vector of next steps in a given branch, returns a status code.
void	ValidateStepRecursive () const
	check that the step and all its suns are correct (it was branch unique IDs defined)
	ClassDef (StepBase, 1)

Additional Inherited Members
Public Types inherited from StepBase
enum	TypeEnum { kCut =0, kAction, kUnknown }
	Enumerator describing the values that _type member can take.
enum	BranchStatus { BranchOK =0, BranchEmpty, NoBranches, BranchOutOfRange }
Protected Attributes inherited from StepBase
TypeEnum	_type
	The type of this step: cut or action.
bool	_break
	true if the cut sequence should be stopped when a cut fails
std::string	_title
	The title of this step, which is used by the DrawingTools.
std::vector< bool >	_disabledInBranch
	Is the step disabled in a given branch.
std::vector< StepBase * >	_nextSteps
	The vector of next steps.
std::vector< UInt_t >	_branchUniqueIDs
	Branch unique IDs associated to this step.
Int_t	_index
	Index of the step in the selection.

Detailed Description

Find Muon candidate.

Definition at line 150 of file numuCCBwdSelection.hxx.

Member Function Documentation

Apply()

bool numuCCBwd::FindMuonCandidateAction::Apply ( AnaEventC &  event, ToyBoxB &  box  ) const

virtual

The return value should specify whether the step was passed. For actions the return value is ignored

Reimplemented from StepBase.

Definition at line 292 of file numuCCBwdSelection.cxx.

                                                                                 {
//**************************************************

  /*
  Can't use the usual find leading track action since it requires the upstream and not the downstream end to be inside the FGD FV
  Muon candidate is taken as the highest mom "positive" track whose downstream end is inside the FGD FV, >18 TPC hits needed
  <=> HMPT FGD FV end <=> HMNT FGD FV start     
  */
  
  // Cast the ToyBox to the appropriate type
  ToyBoxTracker& box = *static_cast<ToyBoxTracker*>(&boxB); 

  box.MainTrack = NULL;

  EventBoxB* EventBox = event.EventBoxes[EventBoxId::kEventBoxTracker];
  AnaRecObjectC** selTracksNS = EventBox->RecObjectsInGroup[EventBoxTracker::kTracksWithTPC];
  int nTPC = EventBox->nRecObjectsInGroup[EventBoxTracker::kTracksWithTPC];
  
  //Sort according to TPC mom TODO : put it in a correction
  std::vector<AnaTrackB*> selTracksV;
  for (Int_t i=0;i<nTPC; ++i){
    selTracksV.push_back( static_cast<AnaTrackB*>(selTracksNS[i]) );
  }
  std::sort(selTracksV.begin(), selTracksV.end(), HMFirst);
  
  
  //Copy vector in an array
  AnaTrackB** selTracks = new AnaTrackB*[nTPC];
  std::copy(selTracksV.begin(), selTracksV.end(), selTracks);
  
  
  
  //loop over tpc tracks
  for (Int_t i=0;i<nTPC; ++i){
    AnaTrackB* track = selTracks[i];    

    if( numuCCBwd::EndsInFGDFV(*track) != box.DetectorFV ) continue;

    //apply the quality cut (>18 hits in TPC)
    if ( ( box.DetectorFV == SubDetId::kFGD1 && !numuCCBwd::TrackQualityCut(*track, SubDetId::kTPC1 ) ) 
      || ( box.DetectorFV == SubDetId::kFGD2 && !numuCCBwd::TrackQualityCut(*track, SubDetId::kTPC2 ) ) ){
      continue;
    }


    //Require a "positive" charge (implies a negative one in the backward hypothesis)
    if( track->Charge == 1 ){
      //Take the first (HM) track which satifsies the previous requirements
      box.MainTrack = track;
      
      delete[] selTracks;
      return true;
    }
  }
  
  delete[] selTracks;
  return false;
  
}

---

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

• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheSelections/v3r12/src/numuCCBwdSelection.hxx
• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheSelections/v3r12/src/numuCCBwdSelection.cxx