Added a value for computin absolute resolution instead of relative one
authorpulvir <pulvir@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 2 Aug 2011 17:20:45 +0000 (17:20 +0000)
committerpulvir <pulvir@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 2 Aug 2011 17:20:45 +0000 (17:20 +0000)
PWG2/RESONANCES/AliRsnMiniPair.h
PWG2/RESONANCES/AliRsnMiniValue.cxx
PWG2/RESONANCES/AliRsnMiniValue.h

index 58413e4..7b3c9a2 100644 (file)
@@ -35,6 +35,7 @@ public:
    Double_t        Eta(Bool_t mc)            const  {return fSum[ID(mc)].Eta();}
    Double_t        InvMass(Bool_t mc)        const  {return fSum[ID(mc)].M();}
    Double_t        InvMassRes()              const;            
+   Double_t        InvMassDiff()             const;            
    Double_t        Mt(Bool_t mc)             const  {return fRef[ID(mc)].Mt();}
    Double_t        Y(Bool_t mc)              const  {return fRef[ID(mc)].Rapidity();}
    Double_t        PtRatio(Bool_t mc)        const;
@@ -77,6 +78,17 @@ inline Double_t AliRsnMiniPair::InvMassRes() const
    return (fSum[0].M() - fSum[1].M()) / fSum[1].M();
 }
 
+inline Double_t AliRsnMiniPair::InvMassDiff() const
+{
+//
+// Return invariant mass resolution
+//
+
+   if (fSum[1].M() <= 0.0) return 1E20;
+   
+   return (fSum[0].M() - fSum[1].M());
+}
+
 inline Double_t AliRsnMiniPair::PtRatio(Bool_t mc) const
 {
 //
index 894bf56..5bc310a 100644 (file)
@@ -102,6 +102,7 @@ const char* AliRsnMiniValue::TypeName(EType type)
       case kPz:           return "Pz";
       case kInvMass:      return "InvMass";
       case kInvMassRes:   return "InvMassResolution";
+      case kInvMassDiff:  return "InvMassDifference";
       case kEta:          return "Eta";
       case kMt:           return "Mt";
       case kY:            return "Y";
@@ -151,6 +152,8 @@ Float_t AliRsnMiniValue::Eval(AliRsnMiniPair *pair, AliRsnMiniEvent *event)
          return pair->Eta(fUseMCInfo);
       case kInvMassRes:
          return pair->InvMassRes();
+      case kInvMassDiff:
+         return pair->InvMassDiff();
       case kMt:
          return pair->Mt(fUseMCInfo);
       case kY:
index 0184938..684247a 100644 (file)
@@ -26,6 +26,7 @@ public:
       kPz,            // pair longitudinal momentum
       kInvMass,       // pair invariant mass (with reconstructed momenta)
       kInvMassRes,    // pair invariant mass resolution
+      kInvMassDiff,   // pair invariant mass difference (MC - reconstructed)
       kEta,           // pair pseudo-rapidity
       kMt,            // pair transverse mass (need a reference mass)
       kY,             // pair rapidity (need a reference mass)