ROOT logo
//
// All cuts for single pions in phi analysis 2010,
// based on track quality and particle identification
// with TPC and TOF.
// Author: Serguey Kiselev.
//
//

#include <Riostream.h>

#include "AliPID.h"
#include "AliPIDResponse.h"
#include "AliRsnCutDaughterKStar2010PP.h"

ClassImp(AliRsnCutDaughterKStar2010PP)

//__________________________________________________________________________________________________
AliRsnCutDaughterKStar2010PP::AliRsnCutDaughterKStar2010PP(const char *name, AliPID::EParticleType pid) :
   AliRsnCut(name, AliRsnTarget::kDaughter),
   fNoPID(kFALSE),
   fPID(pid),
   fCutQuality(Form("%sQuality", name))
{
//
// Constructor
// Initialize track quality cuts to 2010 defaults
//

   fCutQuality.SetPtRange(0.15, 1E+20);
   fCutQuality.SetEtaRange(-0.8, 0.8);
   fCutQuality.SetDCARPtFormula("0.0182+0.0350/pt^1.01");
   fCutQuality.SetDCAZmax(2.0);
   fCutQuality.SetSPDminNClusters(1);
   fCutQuality.SetITSminNClusters(0);
   fCutQuality.SetITSmaxChi2(1E+20);
   fCutQuality.SetTPCminNClusters(70);
   fCutQuality.SetTPCmaxChi2(4.0);
   fCutQuality.SetRejectKinkDaughters();
   fCutQuality.SetAODTestFilterBit(5);
}

//__________________________________________________________________________________________________
Bool_t AliRsnCutDaughterKStar2010PP::IsSelected(TObject *obj)
{
//
// Global check
//

   // coherence check
   if (!TargetOK(obj)) return kFALSE;

   // check track
   AliVTrack *track = dynamic_cast<AliVTrack *>(fDaughter->GetRef());
   if (!track) return kFALSE;

   // check flags
   if ((track->GetStatus() & AliESDtrack::kTPCin   ) == 0) return kFALSE;
   if ((track->GetStatus() & AliESDtrack::kTPCrefit) == 0) return kFALSE;
   if ((track->GetStatus() & AliESDtrack::kITSrefit) == 0) return kFALSE;

   // quality
   if (!fCutQuality.IsSelected(obj)) return kFALSE;

   // if no PID is required, accept
   if (fNoPID) return kTRUE;

   // check initialization of PID object
   AliPIDResponse *pid = fEvent->GetPIDResponse();
   if (!pid) {
      AliFatal("NULL PID response");
      return kFALSE;
   }

   // check if TOF is matched
   // and computes all values used in the PID cut
   Bool_t   isTOF  = MatchTOF(track);
   Double_t pTPC   = track->GetTPCmomentum();
   Double_t p      = track->P();
   Double_t nsTPC  = TMath::Abs(pid->NumberOfSigmasTPC(track, fPID));
   Double_t nsTOF  = isTOF ? TMath::Abs(pid->NumberOfSigmasTOF(track, fPID)) : 1E20;
   Double_t maxTPC = 1E20;
   Double_t maxTOF = 1E20;

   // applies the cut differently depending on the PID and the momentum
   if (isTOF) {
      // TPC: 5sigma cut for all
      if (nsTPC > 5.0) return kFALSE;
      // TOF: 3sigma below 1.5 GeV, 2sigma above
      if (p < 1.5) maxTOF = 3.0; else maxTOF = 2.0;
      return (nsTOF <= maxTOF);
   } else {
      // TPC:
      // all   below   350         MeV: 5sigma
      // all   between 350 and 500 MeV: 3sigma
      // pions above   500         MeV: 2sigma
      // kaons between 500 and 700 MeV: 2sigma
      // kaons above   700         MeV: rejected
      if (pTPC <= 0.35)
         maxTPC = 5.0;
      else if (pTPC > 0.35 && pTPC <= 0.5)
         maxTPC = 3.0;
      else {
         if (fPID == AliPID::kPion)
            maxTPC = 2.0;
         else if (fPID == AliPID::kKaon) {
            if (pTPC <= 0.7)
               maxTPC = 2.0;
            else
               return kFALSE;
         }
      }
      return (nsTPC <= maxTPC);
   }
}
 AliRsnCutDaughterKStar2010PP.cxx:1
 AliRsnCutDaughterKStar2010PP.cxx:2
 AliRsnCutDaughterKStar2010PP.cxx:3
 AliRsnCutDaughterKStar2010PP.cxx:4
 AliRsnCutDaughterKStar2010PP.cxx:5
 AliRsnCutDaughterKStar2010PP.cxx:6
 AliRsnCutDaughterKStar2010PP.cxx:7
 AliRsnCutDaughterKStar2010PP.cxx:8
 AliRsnCutDaughterKStar2010PP.cxx:9
 AliRsnCutDaughterKStar2010PP.cxx:10
 AliRsnCutDaughterKStar2010PP.cxx:11
 AliRsnCutDaughterKStar2010PP.cxx:12
 AliRsnCutDaughterKStar2010PP.cxx:13
 AliRsnCutDaughterKStar2010PP.cxx:14
 AliRsnCutDaughterKStar2010PP.cxx:15
 AliRsnCutDaughterKStar2010PP.cxx:16
 AliRsnCutDaughterKStar2010PP.cxx:17
 AliRsnCutDaughterKStar2010PP.cxx:18
 AliRsnCutDaughterKStar2010PP.cxx:19
 AliRsnCutDaughterKStar2010PP.cxx:20
 AliRsnCutDaughterKStar2010PP.cxx:21
 AliRsnCutDaughterKStar2010PP.cxx:22
 AliRsnCutDaughterKStar2010PP.cxx:23
 AliRsnCutDaughterKStar2010PP.cxx:24
 AliRsnCutDaughterKStar2010PP.cxx:25
 AliRsnCutDaughterKStar2010PP.cxx:26
 AliRsnCutDaughterKStar2010PP.cxx:27
 AliRsnCutDaughterKStar2010PP.cxx:28
 AliRsnCutDaughterKStar2010PP.cxx:29
 AliRsnCutDaughterKStar2010PP.cxx:30
 AliRsnCutDaughterKStar2010PP.cxx:31
 AliRsnCutDaughterKStar2010PP.cxx:32
 AliRsnCutDaughterKStar2010PP.cxx:33
 AliRsnCutDaughterKStar2010PP.cxx:34
 AliRsnCutDaughterKStar2010PP.cxx:35
 AliRsnCutDaughterKStar2010PP.cxx:36
 AliRsnCutDaughterKStar2010PP.cxx:37
 AliRsnCutDaughterKStar2010PP.cxx:38
 AliRsnCutDaughterKStar2010PP.cxx:39
 AliRsnCutDaughterKStar2010PP.cxx:40
 AliRsnCutDaughterKStar2010PP.cxx:41
 AliRsnCutDaughterKStar2010PP.cxx:42
 AliRsnCutDaughterKStar2010PP.cxx:43
 AliRsnCutDaughterKStar2010PP.cxx:44
 AliRsnCutDaughterKStar2010PP.cxx:45
 AliRsnCutDaughterKStar2010PP.cxx:46
 AliRsnCutDaughterKStar2010PP.cxx:47
 AliRsnCutDaughterKStar2010PP.cxx:48
 AliRsnCutDaughterKStar2010PP.cxx:49
 AliRsnCutDaughterKStar2010PP.cxx:50
 AliRsnCutDaughterKStar2010PP.cxx:51
 AliRsnCutDaughterKStar2010PP.cxx:52
 AliRsnCutDaughterKStar2010PP.cxx:53
 AliRsnCutDaughterKStar2010PP.cxx:54
 AliRsnCutDaughterKStar2010PP.cxx:55
 AliRsnCutDaughterKStar2010PP.cxx:56
 AliRsnCutDaughterKStar2010PP.cxx:57
 AliRsnCutDaughterKStar2010PP.cxx:58
 AliRsnCutDaughterKStar2010PP.cxx:59
 AliRsnCutDaughterKStar2010PP.cxx:60
 AliRsnCutDaughterKStar2010PP.cxx:61
 AliRsnCutDaughterKStar2010PP.cxx:62
 AliRsnCutDaughterKStar2010PP.cxx:63
 AliRsnCutDaughterKStar2010PP.cxx:64
 AliRsnCutDaughterKStar2010PP.cxx:65
 AliRsnCutDaughterKStar2010PP.cxx:66
 AliRsnCutDaughterKStar2010PP.cxx:67
 AliRsnCutDaughterKStar2010PP.cxx:68
 AliRsnCutDaughterKStar2010PP.cxx:69
 AliRsnCutDaughterKStar2010PP.cxx:70
 AliRsnCutDaughterKStar2010PP.cxx:71
 AliRsnCutDaughterKStar2010PP.cxx:72
 AliRsnCutDaughterKStar2010PP.cxx:73
 AliRsnCutDaughterKStar2010PP.cxx:74
 AliRsnCutDaughterKStar2010PP.cxx:75
 AliRsnCutDaughterKStar2010PP.cxx:76
 AliRsnCutDaughterKStar2010PP.cxx:77
 AliRsnCutDaughterKStar2010PP.cxx:78
 AliRsnCutDaughterKStar2010PP.cxx:79
 AliRsnCutDaughterKStar2010PP.cxx:80
 AliRsnCutDaughterKStar2010PP.cxx:81
 AliRsnCutDaughterKStar2010PP.cxx:82
 AliRsnCutDaughterKStar2010PP.cxx:83
 AliRsnCutDaughterKStar2010PP.cxx:84
 AliRsnCutDaughterKStar2010PP.cxx:85
 AliRsnCutDaughterKStar2010PP.cxx:86
 AliRsnCutDaughterKStar2010PP.cxx:87
 AliRsnCutDaughterKStar2010PP.cxx:88
 AliRsnCutDaughterKStar2010PP.cxx:89
 AliRsnCutDaughterKStar2010PP.cxx:90
 AliRsnCutDaughterKStar2010PP.cxx:91
 AliRsnCutDaughterKStar2010PP.cxx:92
 AliRsnCutDaughterKStar2010PP.cxx:93
 AliRsnCutDaughterKStar2010PP.cxx:94
 AliRsnCutDaughterKStar2010PP.cxx:95
 AliRsnCutDaughterKStar2010PP.cxx:96
 AliRsnCutDaughterKStar2010PP.cxx:97
 AliRsnCutDaughterKStar2010PP.cxx:98
 AliRsnCutDaughterKStar2010PP.cxx:99
 AliRsnCutDaughterKStar2010PP.cxx:100
 AliRsnCutDaughterKStar2010PP.cxx:101
 AliRsnCutDaughterKStar2010PP.cxx:102
 AliRsnCutDaughterKStar2010PP.cxx:103
 AliRsnCutDaughterKStar2010PP.cxx:104
 AliRsnCutDaughterKStar2010PP.cxx:105
 AliRsnCutDaughterKStar2010PP.cxx:106
 AliRsnCutDaughterKStar2010PP.cxx:107
 AliRsnCutDaughterKStar2010PP.cxx:108
 AliRsnCutDaughterKStar2010PP.cxx:109
 AliRsnCutDaughterKStar2010PP.cxx:110
 AliRsnCutDaughterKStar2010PP.cxx:111
 AliRsnCutDaughterKStar2010PP.cxx:112
 AliRsnCutDaughterKStar2010PP.cxx:113
 AliRsnCutDaughterKStar2010PP.cxx:114