message: Update AliAnalysisTaskNuclei.h
  Updated multiplicity selection
diff: #ifndef AliAnalysisTaskNuclei_cxx
  #define AliAnalysisTaskNuclei_cxx
  +#include "AliAnalysisTaskNuclei.h"
  +#include "TSystem.h"
  +#include "AliAnalysisManager.h"
  +#include "AliAODInputHandler.h"
  +#include "AliMCEventHandler.h"
  +#include "AliAnalysisTaskPIDResponse.h"
  +#include "AliPhysicsSelectionTask.h"
  +#include "AliMultSelectionTask.h"
  +#include "TChain.h"
  +#include "TMacro.h"
  +#include "TROOT.h"
  +#include "AliEventCuts.h"
  +#include "AliAODTrack.h"
  +#include "TNtuple.h"
  +#include "AliAnalysisAlien.h"
  // example of analysis task
  class TH1F;
  @@ -18,7 +33,8 @@ class AliPIDResponse; // main class for PID analysis
  using namespace std;
  -class AliAnalysisTaskNuclei : public AliAnalysisTaskSE {
  +class AliAnalysisTaskNuclei : public AliAnalysisTaskSE
  +{
  public:
  AliAnalysisTaskNuclei();
  AliAnalysisTaskNuclei(const char *name);
  @@ -44,12 +60,13 @@ public:
  void SetUseTOFPidCut(Bool_t useTOFPidCut) { fUseTOFPidCut = useTOFPidCut; }
  void SetMomForTOFanaProt(Float_t momTOFprot) { fMomTOFProt = momTOFprot; }
  void SetMomForTOFanaDeut(Float_t momTOFdeut) { fMomTOFDeut = momTOFdeut; }
  - void SetfYlimitFlag(Float_t ylimit) {fYlimit=ylimit;}
  + void SetfYlimitFlag(Float_t ylimit) { fYlimit = ylimit; }//somehow this doesn't actually set the limit. Unsure why. It is now set in the constructor.
  void SetAnalyseAllParticles(Bool_t doAnalyse) { kAnalyseAllParticles = doAnalyse; } // Call this function with argument kTRUE in the AddTask macro in order to run the analysis for protons and deuterons as well.
  + void SetCentralityFramework(int centralityFramework) { fEventCut.fCentralityFramework = centralityFramework; }
  + void SetCentralityEstimator(int est) { fEstimator = est; }
  void CreateHistosTrack(vector<TH1 *> &histos);
  void FillHistosTrack(vector<TH1 *> &histos, AliAODTrack *track);
  - void SetCentralityEstimator(int est) { fEstimator = est; }
  Float_t GetTOFBeta(AliAODTrack *);
  Float_t GetMass2TOF(Float_t, AliAODTrack *);
  AliEventCuts fEventCut;
  @@ -134,10 +151,23 @@ private:
  Float_t fMomTOFDeut;
  Int_t fEstimator;
  - Float_t fYlimit = -999; //!switch for rapidity cut
  + Float_t fYlimit; //!switch for rapidity cut
  Bool_t kAnalyseAllParticles;
  - enum {kSelectedEvents=0, kINT7selected, kDAQincomplete, kV0timing, kClusterTrackletCut, kVertexZ, kVtxNcontrib, kSPDPileUp, kVtxDisplace, kVtxRes, kNbinsEvent};
  + enum
  + {
  + kSelectedEvents = 0,
  + kINT7selected,
  + kDAQincomplete,
  + kV0timing,
  + kClusterTrackletCut,
  + kVertexZ,
  + kVtxNcontrib,
  + kSPDPileUp,
  + kVtxDisplace,
  + kVtxRes,
  + kNbinsEvent
  + };
  AliAnalysisTaskNuclei(const AliAnalysisTaskNuclei &); // not implemented
  AliAnalysisTaskNuclei &operator=(const AliAnalysisTaskNuclei &); // not implemented
  