AliEMCalPtTaskTrackSelectionESD.cxx

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998-2014, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
#include <TClonesArray.h>
#include <TObjArray.h>
#include "AliEMCalPtTaskTrackSelectionESD.h"
#include <memory>

#include "AliESDEvent.h"
#include "AliESDtrack.h"
#include "AliESDtrackCuts.h"
#include "AliLog.h"
#include "AliPicoTrack.h"
#include "AliVCuts.h"

ClassImp(EMCalTriggerPtAnalysis::AliEMCalPtTaskTrackSelectionESD)

namespace EMCalTriggerPtAnalysis {

AliEMCalPtTaskTrackSelectionESD::AliEMCalPtTaskTrackSelectionESD():
    AliEMCalPtTaskVTrackSelection()
{
}

AliEMCalPtTaskTrackSelectionESD::AliEMCalPtTaskTrackSelectionESD(AliESDtrackCuts* cuts):
    AliEMCalPtTaskVTrackSelection()
{
  this->AddTrackCuts(cuts);
}


TObjArray* AliEMCalPtTaskTrackSelectionESD::GetAcceptedTracks(const TClonesArray* const tracks) {
  if(!fListOfTracks) fListOfTracks = new TObjArray;
  else fListOfTracks->Clear();
  for(TIter trackIter = TIter(tracks).Begin(); trackIter != TIter::End(); ++trackIter){
    if(IsTrackAccepted(static_cast<AliVTrack *>(*trackIter))) fListOfTracks->Add(*trackIter);
  }
  return fListOfTracks;
}

TObjArray* AliEMCalPtTaskTrackSelectionESD::GetAcceptedTracks(const AliVEvent* const event) {
  if(!fListOfTracks) fListOfTracks = new TObjArray;
  else fListOfTracks->Clear();
  const AliESDEvent *esd = dynamic_cast<const AliESDEvent *>(event);
  if(!esd){
    AliError("Event not of type AliESDEvent");
    return fListOfTracks;
  }
  for(int itrk = 0; itrk < esd->GetNumberOfTracks(); itrk++){
    AliESDtrack *trk = static_cast<AliESDtrack *>(esd->GetTrack(itrk));
    if(IsTrackAccepted(trk)) fListOfTracks->AddLast(trk);
  }
  return fListOfTracks;
}

bool AliEMCalPtTaskTrackSelectionESD::IsTrackAccepted(AliVTrack* const trk) {
  if(!fListOfCuts) return kTRUE;
  AliESDtrack *esdt = dynamic_cast<AliESDtrack *>(trk);
  if(!esdt){
    AliPicoTrack *picoTrack = dynamic_cast<AliPicoTrack *>(trk);
    if(picoTrack)
      esdt = dynamic_cast<AliESDtrack*>(picoTrack->GetTrack());
    else{
      AliError("Neither Pico nor ESD track");
      return kFALSE;
    }
  }

  bool result = true;
  for(TIter cutIter = TIter(fListOfCuts).Begin(); cutIter != TIter::End(); ++cutIter){
    if(!(static_cast<AliVCuts *>(*cutIter))->IsSelected(esdt)) result = false;
  }
  return result;
}

} /* namespace EMCalTriggerPtAnalysis */