-#ifndef AliEbyEFluctuationAnalysisTask_cxx\r
-#define AliEbyEFluctuationAnalysisTask_cxx\r
-\r
-// Event by event charge fluctuation analysis\r
-// Authors: Satyajit Jena and Panos Cristakoglou\r
-\r
-class TH1F;\r
-class TH2F;\r
-class TString;\r
-class AliESDEvent;\r
-class AliESDtrackCuts;\r
-\r
-#include "AliAnalysisTaskSE.h"\r
-\r
-const Int_t nCentralityBins = 20;\r
-\r
-class AliEbyEFluctuationAnalysisTask : public AliAnalysisTaskSE {\r
- public:\r
- AliEbyEFluctuationAnalysisTask() : AliAnalysisTaskSE(), fESD(0), fOutputList(0), fHistEventStats(0), fHistCentrality(0), fHistNMultMC(0), fHistNPlusNMinusMC(0), fESDtrackCuts(0), fAnalysisType(0), fAnalysisMode(0), fCentralityEstimator("V0M"), fCentralityBins20(kFALSE), fVxMax(3.0),fVyMax(3.0), fVzMax(10.) {\r
- for(Int_t iBin = 0; iBin < nCentralityBins; iBin++) {\r
- fHistNMult[iBin] = NULL;\r
- fHistNPlusNMinus[iBin] = NULL;\r
- }\r
- }\r
- AliEbyEFluctuationAnalysisTask(const char *name);\r
- virtual ~AliEbyEFluctuationAnalysisTask() {}\r
- \r
- virtual void UserCreateOutputObjects();\r
- virtual void UserExec(Option_t *option);\r
- virtual void Terminate(Option_t *);\r
- \r
- void SetAnalysisCutObject(AliESDtrackCuts *const trackCuts) {\r
- fESDtrackCuts = trackCuts;}\r
- void SetVertexDiamond(Double_t vx, Double_t vy, Double_t vz) {\r
- fVxMax = vx;\r
- fVyMax = vy;\r
- fVzMax = vz;\r
- }\r
- \r
- //Centrality\r
- void SetCentralityEstimator(const char* centralityEstimator) {\r
- fCentralityEstimator = centralityEstimator;}\r
- void SetCentralityBins20() {fCentralityBins20 = kTRUE;}\r
-\r
- void SetAnalysisType(const char* analysisType) {\r
- fAnalysisType = analysisType;}\r
- void SetAnalysisMode(const char* analysisMode) {\r
- fAnalysisMode = analysisMode;}\r
-\r
- private:\r
- AliESDEvent *fESD; //! ESD object\r
- TList *fOutputList; //! Output list\r
- TH1F *fHistEventStats; //!event stats\r
- TH1F *fHistCentrality; //!centrality\r
- TH1F *fHistNMult[nCentralityBins]; //! nmult\r
- TH2F *fHistNPlusNMinus[nCentralityBins];//!nplus vs nminus correlation\r
- TH1F *fHistNMultMC; //!nmult MC\r
- TH2F *fHistNPlusNMinusMC;//!nplus vs nminus correlation\r
-\r
- AliESDtrackCuts *fESDtrackCuts; //ESD track cuts\r
-\r
- TString fAnalysisType;//"MC", "ESD", "AOD"\r
- TString fAnalysisMode;//"TPC", "Global"\r
-\r
- TString fCentralityEstimator;//"V0M","TRK","TKL","ZDC","FMD"\r
- Bool_t fCentralityBins20;//centrality bins of 5% width\r
-\r
- Double_t fVxMax;//vxmax\r
- Double_t fVyMax;//vymax\r
- Double_t fVzMax;//vzmax\r
-\r
- AliEbyEFluctuationAnalysisTask(const AliEbyEFluctuationAnalysisTask&); // not implemented\r
- AliEbyEFluctuationAnalysisTask& operator=(const AliEbyEFluctuationAnalysisTask&); // not implemented\r
- \r
- ClassDef(AliEbyEFluctuationAnalysisTask, 1); // example of analysis\r
-};\r
-\r
-#endif\r
+#ifndef AliEbyEFluctuationAnalysisTask_cxx
+#define AliEbyEFluctuationAnalysisTask_cxx
+
+// Event by event charge fluctuation analysis
+// Authors: Satyajit Jena and Panos Cristakoglou
+
+class TH1F;
+class TH2F;
+class TString;
+class AliESDEvent;
+class AliESDtrackCuts;
+
+#include "AliAnalysisTaskSE.h"
+
+const Int_t nCentralityBins = 20;
+
+class AliEbyEFluctuationAnalysisTask : public AliAnalysisTaskSE {
+ public:
+ AliEbyEFluctuationAnalysisTask() : AliAnalysisTaskSE(), fESD(0), fOutputList(0), fHistEventStats(0), fHistCentrality(0), fHistNMultMC(0), fHistNPlusNMinusMC(0), fESDtrackCuts(0), fAnalysisType(0), fAnalysisMode(0), fCentralityEstimator("V0M"), fCentralityBins20(kFALSE), fVxMax(3.0),fVyMax(3.0), fVzMax(10.) {
+ for(Int_t iBin = 0; iBin < nCentralityBins; iBin++) {
+ fHistNMult[iBin] = NULL;
+ fHistNPlusNMinus[iBin] = NULL;
+ }
+ }
+ AliEbyEFluctuationAnalysisTask(const char *name);
+ virtual ~AliEbyEFluctuationAnalysisTask() {}
+
+ virtual void UserCreateOutputObjects();
+ virtual void UserExec(Option_t *option);
+ virtual void Terminate(Option_t *);
+
+ void SetAnalysisCutObject(AliESDtrackCuts *const trackCuts) {
+ fESDtrackCuts = trackCuts;}
+ void SetVertexDiamond(Double_t vx, Double_t vy, Double_t vz) {
+ fVxMax = vx;
+ fVyMax = vy;
+ fVzMax = vz;
+ }
+
+ //Centrality
+ void SetCentralityEstimator(const char* centralityEstimator) {
+ fCentralityEstimator = centralityEstimator;}
+ void SetCentralityBins20() {fCentralityBins20 = kTRUE;}
+
+ void SetAnalysisType(const char* analysisType) {
+ fAnalysisType = analysisType;}
+ void SetAnalysisMode(const char* analysisMode) {
+ fAnalysisMode = analysisMode;}
+
+ private:
+ AliESDEvent *fESD; //! ESD object
+ TList *fOutputList; //! Output list
+ TH1F *fHistEventStats; //!event stats
+ TH1F *fHistCentrality; //!centrality
+ TH1F *fHistNMult[nCentralityBins]; //! nmult
+ TH2F *fHistNPlusNMinus[nCentralityBins];//!nplus vs nminus correlation
+ TH1F *fHistNMultMC; //!nmult MC
+ TH2F *fHistNPlusNMinusMC;//!nplus vs nminus correlation
+
+ AliESDtrackCuts *fESDtrackCuts; //ESD track cuts
+
+ TString fAnalysisType;//"MC", "ESD", "AOD"
+ TString fAnalysisMode;//"TPC", "Global"
+
+ TString fCentralityEstimator;//"V0M","TRK","TKL","ZDC","FMD"
+ Bool_t fCentralityBins20;//centrality bins of 5% width
+
+ Double_t fVxMax;//vxmax
+ Double_t fVyMax;//vymax
+ Double_t fVzMax;//vzmax
+
+ AliEbyEFluctuationAnalysisTask(const AliEbyEFluctuationAnalysisTask&); // not implemented
+ AliEbyEFluctuationAnalysisTask& operator=(const AliEbyEFluctuationAnalysisTask&); // not implemented
+
+ ClassDef(AliEbyEFluctuationAnalysisTask, 1); // example of analysis
+};
+
+#endif