Updated working version of the ebye charge fluctuation code

[u/mrichter/AliRoot.git] / PWG2 / EBYE / Fluctuations / AliEbyEFluctuationAnalysisTask.h

#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
  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