PWG2/RESONANCES/AliRsnMiniValue.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 ////////////////////////////////////////////////////////////////////////////////
  17 //
  18 //  This class contains all code which is used to compute any of the values
  19 //  which can be of interest within a resonance analysis. Besides the obvious
  20 //  invariant mass, it allows to compute other utility values on all possible
  21 //  targets, in order to allow a wide spectrum of binning and checks.
  22 //  When needed, this object can also define a binning in the variable which
  23 //  it is required to compute, which is used for initializing axes of output
  24 //  histograms (see AliRsnFunction).
  25 //  The value computation requires this object to be passed the object whose
  26 //  informations will be used. This object can be of any allowed input type
  27 //  (track, pair, event), then this class must inherit from AliRsnTarget.
  28 //  Then, when value computation is attempted, a check on target type is done
  29 //  and computation is successful only if expected target matches that of the
  30 //  passed object.
  31 //  In some cases, the value computation can require a support external object,
  32 //  which must then be passed to this class. It can be of any type inheriting
  33 //  from TObject.
  34 //
  35 //  authors: A. Pulvirenti (alberto.pulvirenti@ct.infn.it)
  36 //           M. Vala (martin.vala@cern.ch)
  37 //
  38 ////////////////////////////////////////////////////////////////////////////////
  39
  40 #include "Riostream.h"
  41
  42 #include "AliLog.h"
  43
  44 #include "AliRsnMiniPair.h"
  45 #include "AliRsnMiniEvent.h"
  46
  47 #include "AliRsnMiniValue.h"
  48
  49 ClassImp(AliRsnMiniValue)
  50
  51 //_____________________________________________________________________________
  52 AliRsnMiniValue::AliRsnMiniValue(EType type, Bool_t useMC) :
  53    TNamed(ValueName(type, useMC), ""),
  54    fType(type),
  55    fUseMCInfo(useMC)
  56 {
  57 //
  58 // Constructor
  59 //
  60 }
  61
  62 //_____________________________________________________________________________
  63 AliRsnMiniValue::AliRsnMiniValue(const AliRsnMiniValue& copy) :
  64    TNamed(copy),
  65    fType(copy.fType),
  66    fUseMCInfo(copy.fUseMCInfo)
  67 {
  68 //
  69 // Copy constructor
  70 //
  71 }
  72
  73 //_____________________________________________________________________________
  74 AliRsnMiniValue& AliRsnMiniValue::operator=(const AliRsnMiniValue& copy)
  75 {
  76 //
  77 // Assignment operator.
  78 // Works like copy constructor.
  79 //
  80
  81    TNamed::operator=(copy);
  82    fType = copy.fType;
  83    fUseMCInfo = copy.fUseMCInfo;
  84
  85    return (*this);
  86 }
  87
  88 //_____________________________________________________________________________
  89 const char* AliRsnMiniValue::TypeName(EType type)
  90 {
  91 //
  92 // This method returns a string to give a name to each possible
  93 // computation value.
  94 //
  95
  96    switch (type) {
  97       case kVz:           return "EventVz";
  98       case kMult:         return "EventMult";
  99       case kPlaneAngle:   return "EventPlane";
 100       case kLeadingPt:    return "EventLeadingPt";
 101       case kPt:           return "Pt";
 102       case kPz:           return "Pz";
 103       case kInvMass:      return "InvMass";
 104       case kInvMassRes:   return "InvMassResolution";
 105       case kEta:          return "Eta";
 106       case kMt:           return "Mt";
 107       case kY:            return "Y";
 108       case kPtRatio:      return "PtRatio";
 109       case kDipAngle:     return "DipAngle";
 110       case kCosThetaStar: return "CosThetaStar";
 111       case kAngleLeading: return "AngleToLeading";
 112       default:            return "Undefined";
 113    }
 114 }
 115
 116 //_____________________________________________________________________________
 117 Float_t AliRsnMiniValue::Eval(AliRsnMiniPair *pair, AliRsnMiniEvent *event)
 118 {
 119 //
 120 // Evaluation of the required value.
 121 // In this implementation, fills the member 4-vectors with data
 122 // coming from the object passed as argument, and then returns the value
 123 //
 124
 125    if (!pair && fType > kEventCuts) {
 126       AliError("Null pair passed!");
 127       return 1E20;
 128    }
 129
 130    // compute value depending on types in the enumeration
 131    // if the type does not match any available choice, or if
 132    // the computation is not doable due to any problem
 133    // (not initialized support object, wrong values, risk of floating point errors)
 134    // the method returng kFALSE and sets the computed value to a meaningless number
 135    switch (fType) {
 136       // ---- event values -------------------------------------------------------------------------
 137       case kVz:
 138          return event->Vz();
 139       case kMult:
 140          return event->Mult();
 141       case kPlaneAngle:
 142          return event->Angle();
 143       case kLeadingPt:
 144          return 0.0;
 145       // ---- pair values --------------------------------------------------------------------------
 146       case kPt:
 147          return pair->Pt(fUseMCInfo);
 148       case kInvMass:
 149          return pair->InvMass(fUseMCInfo);
 150       case kEta:
 151          return pair->Eta(fUseMCInfo);
 152       case kInvMassRes:
 153          return pair->InvMassRes();
 154       case kMt:
 155          return pair->Mt(fUseMCInfo);
 156       case kY:
 157          return pair->Y(fUseMCInfo);
 158       case kPtRatio:
 159          return pair->PtRatio(fUseMCInfo);
 160       case kDipAngle:
 161          return pair->DipAngle(fUseMCInfo);
 162       case kCosThetaStar:
 163          return pair->CosThetaStar(fUseMCInfo);
 164       case kAngleLeading:
 165          AliWarning("This method is not yet implemented");
 166          return 1E20;
 167       default:
 168          AliError("Invalid value type");
 169          return 1E20;
 170    }
 171 }