--- /dev/null
+#ifndef ALIANALYSISCENTRALCUTESD_H
+#define ALIANALYSISCENTRALCUTESD_H
+
+
+/*
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
+ * See cxx source for full Copyright notice
+ * $Id$
+ */
+
+// ---------------------------------------------------
+// ESD particle level cuts for azimuthal isotropic
+// expansion in highly central collisions analysis
+// author: Cristian Andrei
+// acristian@niham.nipne.ro
+// ---------------------------------------------------
+
+
+#include <TString.h>
+#include <TPDGCode.h>
+
+#include "AliAnalysisCuts.h"
+
+
+class TObject;
+class TList;
+class TF1;
+
+class AliESDtrack;
+
+
+class AliAnalysisCentralCutESD: public AliAnalysisCuts {
+public:
+ AliAnalysisCentralCutESD(const char *name="AliAnalysisCentralCutESD", const char *title="ESD_cuts");
+ virtual ~AliAnalysisCentralCutESD();
+
+ Bool_t IsSelected(TObject* obj);
+ Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+
+ void SetPartType(PDG_t type) {fReqPID = kTRUE; fPartType = type;}
+ void SetPIDtype(TString type) {fPIDtype = type;}
+ void SetPriorFunctions(Bool_t pfunc) {fPriorsFunc = pfunc;}
+
+ void SetReqIsCharged() {fReqCharge = kTRUE;}
+
+private:
+ AliAnalysisCentralCutESD(const AliAnalysisCentralCutESD &ref);
+ AliAnalysisCentralCutESD& operator=(const AliAnalysisCentralCutESD &ref);
+
+ Bool_t fReqPID; //kTRUE -> run the PID
+
+ Bool_t fReqCharge; //kTRUE -> only charged particles selected
+
+ PDG_t fPartType; //can be kPi*, kK* or kProton
+
+ TString fPIDtype; //PID method -> can be Custom or Bayesian
+ Bool_t fPriorsFunc; // if kTRUE -> use function priors
+
+ Double_t fPartPriors[10]; //prior probabilities
+ TF1 *fElectronFunction; //momentum dependence of the prior probs
+ TF1 *fMuonFunction; //momentum dependence of the prior probs
+ TF1 *fPionFunction; //momentum dependence of the prior probs
+ TF1 *fKaonFunction; //momentum dependence of the prior probs
+ TF1 *fProtonFunction; //momentum dependence of the prior probs
+
+ Bool_t IsA(AliESDtrack *track, PDG_t type);
+
+ Bool_t IsCharged (AliESDtrack* const track) const;
+
+ Double_t GetPriors(Int_t i, Double_t P);
+
+ ClassDef(AliAnalysisCentralCutESD, 1);
+};
+
+#endif