PWGLF/totEt/AliAnalysisEtSelectorEmcal.cxx

   1 //_________________________________________________________________________
   2 //  Utility Class for transverse energy studies
   3 //  Selection class for EMCAL
   4 //
   5 //*-- Authors: Oystein Djuvsland (Bergen)
   6 //_________________________________________________________________________
   7
   8
   9 #include "AliAnalysisEtSelectorEmcal.h"
  10 #include "AliAnalysisEtCuts.h"
  11 #include "TParticle.h"
  12 void AliAnalysisEtSelectorEmcal::SetEvent(const AliESDEvent* event)
  13 {
  14     fEvent = event;
  15     if(!fInitialized) Init(event);
  16 }
  17
  18 AliAnalysisEtSelectorEmcal::AliAnalysisEtSelectorEmcal(AliAnalysisEtCuts* cuts):AliAnalysisEtSelector(cuts)
  19 ,fInitialized(kFALSE)
  20 {
  21
  22 }
  23
  24 AliAnalysisEtSelectorEmcal::~AliAnalysisEtSelectorEmcal()
  25 {
  26
  27 }
  28
  29 void AliAnalysisEtSelectorEmcal::Init()
  30 {
  31     AliAnalysisEtSelector::Init();
  32 }
  33
  34 Int_t AliAnalysisEtSelectorEmcal::Init(const AliESDEvent* event)
  35 {
  36
  37     AliAnalysisEtSelector::Init(event);
  38     Printf("Initializing selector for run: %d", event->GetRunNumber());
  39     fInitialized = kTRUE;
  40     return 0;
  41 }
  42
  43 TRefArray* AliAnalysisEtSelectorEmcal::GetClusters()
  44 {
  45
  46   if(!fClusterArray) fClusterArray = new TRefArray;
  47
  48   if(fClusterArray)
  49   {
  50     fEvent->GetEMCALClusters(fClusterArray);
  51   }
  52   else
  53   {
  54     Printf("Could not initialize cluster array");
  55   }
  56   return fClusterArray;
  57 }
  58
  59 Bool_t AliAnalysisEtSelectorEmcal::CutMinEnergy(const AliESDCaloCluster& cl) const
  60 {
  61   return cl.E() > fCuts->GetReconstructedEmcalClusterEnergyCut();
  62 }
  63
  64 Bool_t AliAnalysisEtSelectorEmcal::CutMinEnergy(const TParticle& p) const
  65 {
  66     return p.Energy() > fCuts->GetReconstructedEmcalClusterEnergyCut();
  67 }
  68
  69 Bool_t AliAnalysisEtSelectorEmcal::CutDistanceToBadChannel(const AliESDCaloCluster& ) const
  70 {
  71     return kTRUE;
  72 }
  73
  74 Bool_t AliAnalysisEtSelectorEmcal::CutTrackMatching(const AliESDCaloCluster& ) const
  75 {
  76     return kTRUE;
  77 }
  78
  79 Bool_t AliAnalysisEtSelectorEmcal::CutGeometricalAcceptance(const TParticle& part) const
  80 {
  81   return TMath::Abs(part.Eta()) < fCuts->GetGeometryEmcalEtaAccCut()
  82           && part.Phi() < fCuts->GetGeometryEmcalPhiAccMaxCut()*TMath::Pi()/180.
  83           && part.Phi() > fCuts->GetGeometryEmcalPhiAccMinCut()*TMath::Pi()/180.;
  84 }
  85
  86 Bool_t AliAnalysisEtSelectorEmcal::CutGeometricalAcceptance(const AliVTrack& part) const
  87 {
  88   return TMath::Abs(part.Eta()) < fCuts->GetGeometryEmcalEtaAccCut()
  89           && part.Phi() < fCuts->GetGeometryEmcalPhiAccMaxCut()*TMath::Pi()/180.
  90           && part.Phi() > fCuts->GetGeometryEmcalPhiAccMinCut()*TMath::Pi()/180.;
  91 }
  92
  93
  94
  95
  96
  97