1 #ifndef AliEbyEHigherMomentsTask_cxx
2 #define AliEbyEHigherMomentsTask_cxx
4 //=========================================================================//
6 // Analysis Task for Net-Charge Higher Moment Analysis //
7 // Author: Satyajit Jena || Nirbhay K. Behera //
8 // sjena@cern.ch || nbehera@cern.ch //
11 //=========================================================================//
22 #include "TParticle.h"
23 #include "AliAnalysisTaskSE.h"
25 class AliEbyEHigherMomentsTask: public AliAnalysisTaskSE {
27 AliEbyEHigherMomentsTask( const char *name = "HigherMomentAnalysis");
28 virtual ~AliEbyEHigherMomentsTask();
30 virtual void UserCreateOutputObjects();
31 virtual void UserExec(Option_t *option);
32 virtual void doAODEvent();
33 virtual void doMCAODEvent();
35 virtual void Terminate(Option_t *);
37 void SetVertexDiamond(Double_t vx, Double_t vy, Double_t vz) {fVxMax = vx;fVyMax = vy; fVzMax = vz;}
38 void SetCentralityEstimator(const char* centralityEstimator) { fCentralityEstimator = centralityEstimator;}
39 void SetAnalysisType(const char* analysisType) {fAnalysisType = analysisType;}
40 void SetAODtrackCutBit(Int_t bit){ fAODtrackCutBit = bit;}
41 void SetKinematicsCutsAOD(Double_t ptl, Double_t pth, Double_t eta){
44 fEtaLowerLimit = -1.*eta;
45 fEtaHigherLimit = eta;
48 void SetNumberOfPtBins(Int_t nPtBins){ fNptBins = nPtBins;}
54 Bool_t ProperVertex(AliAODEvent *fAOD) const;
55 Bool_t AcceptTrack(AliAODTrack* track) const;
56 Int_t GetPtBin(Double_t pt);
60 TClonesArray *fArrayMC;
61 TString fAnalysisType; // "MC", "ESD", "AOD"
62 TString fCentralityEstimator; // "V0M","TRK","TKL","ZDC","FMD"
65 Double_t fVxMax; //vxmax
66 Double_t fVyMax;//vymax
67 Double_t fVzMax;//vzmax
68 Double_t fPtLowerLimit;
69 Double_t fPtHigherLimit;
72 Double_t fEtaLowerLimit;
73 Double_t fEtaHigherLimit;
74 Int_t fAODtrackCutBit;//track cut bit from track selection (only used for AODs)
81 THnSparse *fTHnCentNplusNminusCh;
82 THnSparse *fTHnCentNplusNminusChTruth;
83 THnSparse *fPtBinNplusNminusCh;
84 THnSparse *fPtBinNplusNminusChTruth;
88 AliEbyEHigherMomentsTask(const AliEbyEHigherMomentsTask&);
89 AliEbyEHigherMomentsTask& operator = (const AliEbyEHigherMomentsTask&);//Not implimented..
90 ClassDef(AliEbyEHigherMomentsTask, 1);