PWG2/RESONANCES/AliRsnCut.cxx

   1 //
   2 // *** Class AliRsnCut ***
   3 //
   4 // Cut base class: all other cuts inherit from it.
   5 // The 'core' of the class is the method "IsSelected()" which
   6 // must be overloaded by any specific cut implementation.
   7 //
   8 // This class provides some default instruments to check values
   9 // agains a reference or an allowed range, in order to permit
  10 // a unique way to execute such kind of checks.
  11 //
  12 // authors: Alberto Pulvirenti (alberto.pulvirenti@ct.infn.it)
  13 //          Martin Vala (martin.vala@cern.ch)
  14 //
  15
  16 #include "AliRsnDaughter.h"
  17 #include "AliRsnMother.h"
  18 #include "AliRsnEvent.h"
  19
  20 #include "AliRsnCut.h"
  21
  22 ClassImp(AliRsnCut)
  23
  24 //______________________________________________________________________________
  25 AliRsnCut::AliRsnCut(const char *name, RSNTARGET target) :
  26   AliRsnTarget(name, target),
  27   fMinI(0),
  28   fMaxI(0),
  29   fMinD(0.),
  30   fMaxD(0.),
  31   fCutValueI(0),
  32   fCutValueD(0.0),
  33   fCutResult(kTRUE)
  34 {
  35 //
  36 // Default constructor.
  37 //
  38 }
  39
  40 //______________________________________________________________________________
  41 AliRsnCut::AliRsnCut
  42 (const char *name, RSNTARGET target, Int_t imin, Int_t imax, Double_t dmin, Double_t dmax) :
  43   AliRsnTarget(name, target),
  44   fMinI(imin),
  45   fMaxI(imax),
  46   fMinD(dmin),
  47   fMaxD(dmax),
  48   fCutValueI(0),
  49   fCutValueD(0.0),
  50   fCutResult(kTRUE)
  51 {
  52 //
  53 // Constructor with arguments.
  54 // This is provided to allow a quick setting of all data members.
  55 //
  56 }
  57
  58 //______________________________________________________________________________
  59 AliRsnCut::AliRsnCut
  60 (const char *name, RSNTARGET target, Double_t dmin, Double_t dmax, Int_t imin, Int_t imax) :
  61   AliRsnTarget(name, target),
  62   fMinI(imin),
  63   fMaxI(imax),
  64   fMinD(dmin),
  65   fMaxD(dmax),
  66   fCutValueI(0),
  67   fCutValueD(0.0),
  68   fCutResult(kTRUE)
  69 {
  70 //
  71 // Constructor with arguments.
  72 // This is provided to allow a quick setting of all data members.
  73 //
  74 }
  75
  76 //______________________________________________________________________________
  77 AliRsnCut::AliRsnCut(const AliRsnCut& copy) :
  78   AliRsnTarget(copy),
  79   fMinI(copy.fMinI),
  80   fMaxI(copy.fMaxI),
  81   fMinD(copy.fMinD),
  82   fMaxD(copy.fMaxD),
  83   fCutValueI(copy.fCutValueI),
  84   fCutValueD(copy.fCutValueD),
  85   fCutResult(copy.fCutResult)
  86 {
  87 //
  88 // Copy constructor.
  89 // Don't duplicate memory occupancy for pointer
  90 //
  91 }
  92
  93 //______________________________________________________________________________
  94 AliRsnCut& AliRsnCut::operator=(const AliRsnCut& copy)
  95 {
  96 //
  97 // Assignment operator.
  98 // Don't duplicate memory occupancy for pointer
  99 //
 100
 101   AliRsnTarget::operator=(copy);
 102
 103   fMinI      = copy.fMinI;
 104   fMaxI      = copy.fMaxI;
 105   fMinD      = copy.fMinD;
 106   fMaxD      = copy.fMaxD;
 107   fCutValueI = copy.fCutValueI;
 108   fCutValueD = copy.fCutValueD;
 109
 110   return (*this);
 111 }
 112
 113 //______________________________________________________________________________
 114 Bool_t AliRsnCut::IsSelected(TObject* /*object*/)
 115 {
 116 //
 117 // Virtual cut-checking method.
 118 // In this implementation, it does nothing, and all classes
 119 // inheriting from this, should provide a proper implementation
 120 // which must return kTRUE if the cut is passed, and kFALSE otherwise.
 121 //
 122
 123   AliWarning("This virtual function must be implemented properly");
 124   return kTRUE;
 125 }
 126
 127 //______________________________________________________________________________
 128 Bool_t AliRsnCut::OkValueI()
 129 {
 130 //
 131 // This method is used when the cut consists in comparing the cut value
 132 // with a reference integer value to which it must be equal.
 133 //
 134
 135   // eval result
 136   fCutResult = (fCutValueI == fMinI);
 137
 138   // print debug message
 139   AliDebug(AliLog::kDebug + 2, "=== CUT DEBUG ====================================");
 140   AliDebug(AliLog::kDebug + 2, Form("Cut name     : %s", GetName()));
 141   AliDebug(AliLog::kDebug + 2, Form("Checked value: %d", fCutValueI));
 142   AliDebug(AliLog::kDebug + 2, Form("Cut value    : %d", fMinI));
 143   AliDebug(AliLog::kDebug + 2, Form("Cut result   : %s", (fCutResult ? "PASSED" : "NOT PASSED")));
 144   AliDebug(AliLog::kDebug + 2, "=== END CUT DEBUG ================================");
 145
 146   return fCutResult;
 147 }
 148
 149 //______________________________________________________________________________
 150 Bool_t AliRsnCut::OkValueD()
 151 {
 152 //
 153 // This method is used when the cut consists in comparing the cut value
 154 // with a reference double value to which it must be equal (or at least, almost).
 155 //
 156
 157   // eval result
 158   fCutResult = (TMath::Abs(fCutValueD - fMinD) < 1E-6);
 159
 160   // print debug message
 161   AliDebug(AliLog::kDebug + 2, "=== CUT DEBUG ====================================");
 162   AliDebug(AliLog::kDebug + 2, Form("Cut name     : %s", GetName()));
 163   AliDebug(AliLog::kDebug + 2, Form("Checked value: %f", fCutValueD));
 164   AliDebug(AliLog::kDebug + 2, Form("Cut value    : %f", fMinD));
 165   AliDebug(AliLog::kDebug + 2, Form("Cut result   : %s", (fCutResult ? "PASSED" : "NOT PASSED")));
 166   AliDebug(AliLog::kDebug + 2, "=== END CUT DEBUG ================================");
 167
 168   return fCutResult;
 169 }
 170
 171 //______________________________________________________________________________
 172 Bool_t AliRsnCut::OkRangeI()
 173 {
 174 //
 175 // This method is used when the cut consists in an allowed range
 176 // where the cut value must be included to pass the cut.
 177 // Then, the cut result is kTRUE if the cut value is inside this range.
 178 //
 179
 180   // eval result
 181   fCutResult = ((fCutValueI >= fMinI) && (fCutValueI <= fMaxI));
 182
 183   // print debug message
 184   AliDebug(AliLog::kDebug + 2, "=== CUT DEBUG ====================================");
 185   AliDebug(AliLog::kDebug + 2, Form("Cut name     : %s", GetName()));
 186   AliDebug(AliLog::kDebug + 2, Form("Checked value: %d", fCutValueI));
 187   AliDebug(AliLog::kDebug + 2, Form("Cut range    : %d , %d", fMinI, fMaxI));
 188   AliDebug(AliLog::kDebug + 2, Form("Cut result   : %s", (fCutResult ? "PASSED" : "NOT PASSED")));
 189   AliDebug(AliLog::kDebug + 2, "=== END CUT DEBUG ================================");
 190
 191   return fCutResult;
 192 }
 193
 194 //______________________________________________________________________________
 195 Bool_t AliRsnCut::OkRangeD()
 196 {
 197 //
 198 // This method is used when the cut consists in an allowed range
 199 // where the cut value must be included to pass the cut.
 200 // Then, the cut result is kTRUE if the cut value is inside this range.
 201 //
 202
 203   // eval result
 204   fCutResult = ((fCutValueD >= fMinD) && (fCutValueD <= fMaxD));
 205
 206   // print debug message
 207   AliDebug(AliLog::kDebug + 2, "=== CUT DEBUG ====================================");
 208   AliDebug(AliLog::kDebug + 2, Form("Cut name     : %s", GetName()));
 209   AliDebug(AliLog::kDebug + 2, Form("Checked value: %f", fCutValueD));
 210   AliDebug(AliLog::kDebug + 2, Form("Cut range    : %f , %f", fMinD, fMaxD));
 211   AliDebug(AliLog::kDebug + 2, Form("Cut result   : %s", (fCutResult ? "PASSED" : "NOT PASSED")));
 212   AliDebug(AliLog::kDebug + 2, "=== END CUT DEBUG ================================");
 213
 214   return fCutResult;
 215 }
 216
 217 //______________________________________________________________________________
 218 void AliRsnCut::Print(Option_t*) const
 219 {
 220 //
 221 // Override TObject::Print() method,
 222 // and print some useful info about the cut general parameters.
 223 //
 224
 225   AliInfo("=== CUT DETAILS ====================================");
 226   AliInfo(Form("Cut name     : [%s]", GetName()));
 227   AliInfo(Form("Cut target   : [%s]", GetTargetTypeName()));
 228   AliInfo(Form("Cut edges [D]: [%f - %f]", fMinD, fMaxD));
 229   AliInfo(Form("Cut edges [I]: [%d - %d]", fMinI, fMaxI));
 230   AliInfo("====================================================");
 231 }