numuCC4piMultiPi::FindEcalPhotonsAction Class Reference

Inheritance diagram for numuCC4piMultiPi::FindEcalPhotonsAction:

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)

Public Attributes
int	MostUpstreamLayerHitCut
bool	prod5Cut

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

Definition at line 189 of file numuCC4piMultiPiSelection.hxx.

Member Function Documentation

Apply()

bool FindEcalPhotonsAction::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 222 of file numuCC4piMultiPiSelection.cxx.

                                                                                      {
//*********************************************************************
    // Find whether there is a pi0 track in the ECal or not 
    // In order to apply a pi0 veto cut later 
    ToyBoxCC4piMultiPi * cc4piMultiPibox = static_cast<ToyBoxCC4piMultiPi*>(&box);
    SubDetId::SubDetEnum det = static_cast<SubDetId::SubDetEnum>(box.DetectorFV);

    // Get all Ecal objects
    EventBoxTracker::RecObjectGroupEnum groupID;
    groupID = EventBoxTracker::kTracksWithECal;
    EventBoxB* EventBox = event.EventBoxes[EventBoxId::kEventBoxTracker];
    
    AnaTrackB *highestObject = NULL;
    double higherEnergyObject = 0;

    //loop over all ECal objects
    for(int i = 0; i < EventBox->nRecObjectsInGroup[groupID]; i++ ) {
        AnaTrackB *allECALObjects = static_cast<AnaTrackB*>(EventBox->RecObjectsInGroup[groupID][i]);
        if(!(allECALObjects->nECALSegments)) continue;
        // Check for isolated ecal object
        if(!(anaUtils::TrackUsesOnlyDet(*allECALObjects,SubDetId::kDSECAL) ||
             anaUtils::TrackUsesOnlyDet(*allECALObjects,SubDetId::kTopTECAL) ||
             anaUtils::TrackUsesOnlyDet(*allECALObjects,SubDetId::kBottomTECAL) ||
             anaUtils::TrackUsesOnlyDet(*allECALObjects,SubDetId::kLeftTECAL) ||
             anaUtils::TrackUsesOnlyDet(*allECALObjects,SubDetId::kRightTECAL))
            )  continue;

        AnaECALParticleB* ecalBase = allECALObjects->ECALSegments[0];    
        AnaTrackB* object =  allECALObjects; 
        cc4piMultiPibox->ISOEcal.push_back(object);   
        
        // Find the most energetic object in ECal 
        if(ecalBase->EMEnergy > higherEnergyObject){
            higherEnergyObject = ecalBase->EMEnergy;
            highestObject = object;
        } 
    }// end of loop on allECALTracks

    // Save the highest energetic ecal object in the box and the pi0 track
    if (highestObject){
        cc4piMultiPibox->HObject  = highestObject;
        // Check if the highest energetic ecal object is shower-like
        if( numuCC4piMultiPiUtils::ECALPi0Selection(event, box, highestObject, MostUpstreamLayerHitCut, det) ) {
            cc4piMultiPibox->Pi0Ecaltrack  = highestObject;
            cc4piMultiPibox->nPi0Ecaltracks = cc4piMultiPibox->nPi0Ecaltracks+1; 
        }
    }  
    return true;
}

---

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

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