Added AliPPVsMultUtils to AliAnalysisUtils

ANALYSIS/AliAnalysisUtils.cxx

@@ -9,6 +9,7 @@
 #include <TMatrixDSym.h>
 #include <TMath.h>
 #include "AliVMultiplicity.h"
+#include "AliPPVsMultUtils.h"
 
 #include "AliAnalysisUtils.h"
 
@@ -33,7 +34,8 @@ AliAnalysisUtils::AliAnalysisUtils():TObject(),
   fnSigmaPlpDiamZSPD(5.),
   fUseSPDCutInMultBins(kFALSE),
   fASPDCvsTCut(65.),
-  fBSPDCvsTCut(4.)
+  fBSPDCvsTCut(4.),
+  fPPVsMultUtils(0x0)
 {
   // Default contructor
 }
@@ -265,3 +267,13 @@ Double_t AliAnalysisUtils::GetWDist(const AliVVertex* v0, const AliVVertex* v1)
 
 }
 
+//______________________________________________________________________
+Float_t AliAnalysisUtils::GetMultiplicityPercentile(AliVEvent *event, TString lMethod ){
+  if(!fPPVsMultUtils)
+    fPPVsMultUtils=new AliPPVsMultUtils();
+  if( (event->InheritsFrom("AliAODEvent")) || (event->InheritsFrom("AliESDEvent")) ) return fPPVsMultUtils->GetMultiplicityPercentile(event,lMethod);
+  else {
+    AliFatal("Event is neither of AOD nor ESD type"); 
+    return -999.;
+  }
+}

ANALYSIS/AliAnalysisUtils.h

@@ -15,6 +15,7 @@
 
 class AliVEvent;
 class AliVVertex;
+class AliPPVsMultUtils;
 
 class AliAnalysisUtils : public TObject {
 
@@ -59,6 +60,9 @@ class AliAnalysisUtils : public TObject {
   void SetASPDCvsTCut(Float_t a) { fASPDCvsTCut = a; }
   void SetBSPDCvsTCut(Float_t b) { fBSPDCvsTCut = b; }
 
+  //multiplicity selection in pp
+  Float_t GetMultiplicityPercentile(AliVEvent *event, TString lMethod = "V0M");
+
  private:
 
   Bool_t fisAOD; // flag for AOD:1 or ESD:0
@@ -85,8 +89,10 @@ class AliAnalysisUtils : public TObject {
 
   Float_t fASPDCvsTCut; // constant for the linear cut in SPD clusters vs tracklets
   Float_t fBSPDCvsTCut; // slope for the linear cut in SPD  clusters vs tracklets
-
-  ClassDef(AliAnalysisUtils,2) // base helper class
+
+  AliPPVsMultUtils *fPPVsMultUtils; //! multiplicity selection in pp
+
+  ClassDef(AliAnalysisUtils,3) // base helper class
 };
 #endif