PWG2/RESONANCES/AliRsnCutPrimaryVertex.cxx

   1 //
   2 // Class AliRsnCutPrimaryVertex
   3 //
   4 // General implementation of a single cut strategy, which can be:
   5 // - a value contained in a given interval  [--> IsBetween()   ]
   6 // - a value equal to a given reference     [--> MatchesValue()]
   7 //
   8 // In all cases, the reference value(s) is (are) given as data members
   9 // and each kind of cut requires a given value type (Int, UInt, Double),
  10 // but the cut check procedure is then automatized and chosen thanks to
  11 // an enumeration of the implemented cut types.
  12 // At the end, the user (or any other point which uses this object) has
  13 // to use the method IsSelected() to check if this cut has been passed.
  14 //
  15 // authors: Martin Vala (martin.vala@cern.ch)
  16 //          Alberto Pulvirenti (alberto.pulvirenti@ct.infn.it)
  17 //
  18
  19 // #include "AliLog.h"
  20 // #include "AliESDEvent.h"
  21 // #include "AliESDVertex.h"
  22 //
  23 // #include "AliRsnEvent.h"
  24 // #include "AliRsnCutPrimaryVertex.h"
  25
  26 // #include "Riostream.h"
  27 // #include "TMath.h"
  28 //
  29 // #include "AliLog.h"
  30 // #include "AliESDEvent.h"
  31 // #include "AliESDVertex.h"
  32 //
  33 // #include "AliRsnEvent.h"
  34 #include "AliRsnCutPrimaryVertex.h"
  35
  36 ClassImp(AliRsnCutPrimaryVertex)
  37
  38 //_________________________________________________________________________________________________
  39 AliRsnCutPrimaryVertex::AliRsnCutPrimaryVertex() :
  40     AliRsnCut()
  41 {
  42 //
  43 // Default constructor.
  44 //
  45 }
  46
  47 //_________________________________________________________________________________________________
  48 AliRsnCutPrimaryVertex::AliRsnCutPrimaryVertex
  49 (const char *name, Int_t nContributors) :
  50     AliRsnCut(name, 0, nContributors)
  51 {
  52 //
  53 // Main constructor.
  54 // the cut range is outside the interval 0 - min number of contributors.
  55 //
  56 }
  57
  58 //_________________________________________________________________________________________________
  59 Bool_t AliRsnCutPrimaryVertex::IsSelected(AliRsnCut::ETarget /*tgt*/, AliRsnDaughter*/*const track*/) const
  60 {
  61 //
  62 // Cut checker.
  63 //
  64   AliWarning("Cannot apply this cut to particles");
  65   return kTRUE;
  66 }
  67
  68 //_________________________________________________________________________________________________
  69 Bool_t AliRsnCutPrimaryVertex::IsSelected(AliRsnCut::ETarget, AliRsnPairParticle*/*const pair*/) const
  70 {
  71 //
  72 // Cut checker
  73 //
  74
  75   AliWarning("Cannot apply this cut to pairs");
  76   return kTRUE;
  77 }
  78
  79 //_________________________________________________________________________________________________
  80 Bool_t AliRsnCutPrimaryVertex::IsSelected(AliRsnCut::ETarget, AliRsnEvent*event)
  81 {
  82 //
  83 // Cut checker
  84 //
  85
  86   // retrieve ESD event
  87   AliESDEvent *esd = dynamic_cast<AliESDEvent*>(event->GetRef());
  88   if (!esd) {
  89     AliDebug(AliLog::kDebug+2, "NO ESD");
  90     return kTRUE;
  91   }
  92
  93   // check primary vertex and eventually fill step 1
  94   // if a vertex with tracks was successfully reconstructed,
  95   // it is used for computing DCA;
  96   // otherwise, the one computed with SPD is used.
  97   // This is known from the "Status" parameter of the vertex itself.
  98   const AliESDVertex *v = esd->GetPrimaryVertex();
  99   if (!v->GetStatus()) v = esd->GetPrimaryVertexSPD();
 100   if (!v->GetStatus()) {
 101     AliDebug(AliLog::kDebug+2, "Bad vertex status");
 102     return kFALSE;
 103   }
 104
 105   fCutValueI = v->GetNContributors();
 106
 107   return !OkRange();
 108 }
 109
 110 //_________________________________________________________________________________________________
 111 Bool_t AliRsnCutPrimaryVertex::IsSelected(AliRsnCut::ETarget, AliRsnEvent*/*ev1*/, AliRsnEvent*/*ev2*/) const
 112 {
 113 //
 114 // Cut checker
 115 //
 116
 117   AliWarning("Cannot apply this cut to event mixing");
 118   return kTRUE;
 119 }