3 #ifndef DNDETANALYSIS_H
4 #define DNDETANALYSIS_H
6 // ------------------------------------------------------
8 // Class for dn/deta analysis
10 // ------------------------------------------------------
13 // - more documentation
14 // - add debug statements
15 // - add more histograms
16 // - add functionality to set the bin sizes
17 // - figure out correct way to treat the errors
18 // - add functionality to make dn/deta for different mult classes?
21 #include "AlidNdEtaCorrection.h"
22 #include "AliPWG0Helper.h"
30 class dNdEtaAnalysis : public TNamed
33 enum { kVertexBinning = 1+2 }; // the first is for the whole vertex range, the others divide the vertex range
36 dNdEtaAnalysis(const Char_t* name, const Char_t* title, AliPWG0Helper::AnalysisMode analysisMode = (AliPWG0Helper::AnalysisMode) (AliPWG0Helper::kTPC | AliPWG0Helper::kFieldOn));
37 virtual ~dNdEtaAnalysis();
39 dNdEtaAnalysis(const dNdEtaAnalysis &c);
40 dNdEtaAnalysis &operator=(const dNdEtaAnalysis &c);
41 virtual void Copy(TObject &c) const;
43 void FillTrack(Float_t vtx, Float_t eta, Float_t pt);
44 void FillEvent(Float_t vtx, Float_t n);
45 void FillTriggeredEvent(Float_t n);
47 void Finish(AlidNdEtaCorrection* correction, Float_t ptCut, AlidNdEtaCorrection::CorrectionType correctionType, const char* tag = "", Int_t backgroundSubtraction = 0, TH1* combinatoricsCorrection = 0);
49 void DrawHistograms(Bool_t simple = kFALSE);
50 void LoadHistograms(const Char_t* dir = 0);
51 void SaveHistograms();
53 virtual Long64_t Merge(TCollection* list);
55 AliCorrection* GetData() const { return fData; }
57 TH1F* GetPtHistogram() const { return fPtDist; }
59 TH1F* GetdNdEtaHistogram(Int_t i = 0) const { return fdNdEta[i]; }
60 TH1F* GetdNdEtaPtCutOffCorrectedHistogram(Int_t i = 0) const { return fdNdEtaPtCutOffCorrected[i]; }
62 void SetAnalysisMode(AliPWG0Helper::AnalysisMode mode) { fAnalysisMode = mode; }
63 AliPWG0Helper::AnalysisMode GetAnalysisMode() { return fAnalysisMode; }
65 void SetVtxMin(Float_t vtxMin) {fvtxMin = vtxMin;}
66 void SetVtxMax(Float_t vtxMax) {fvtxMax = vtxMax;}
69 Float_t GetVtxMin(Float_t eta);
71 AliCorrection* fData; // we store the data in an AliCorrection
72 TH1F* fMult; // (raw) multiplicity distribution of triggered events
74 TH1F* fPtDist; // pt distribution
76 TH1F* fdNdEta[kVertexBinning]; // dndeta results for different vertex bins (0 = full range)
77 TH1F* fdNdEtaPtCutOffCorrected[kVertexBinning]; // dndeta results for different vertex bins (0 = full range), pt cut off corrected
79 AliPWG0Helper::AnalysisMode fAnalysisMode; // detector (combination) used for the analysis
80 TString fTag; // tag saved that describes the applied correction
82 Float_t fvtxMin; // mininum vtx on which to integrate
83 Float_t fvtxMax; // maximum vtx on which to intgrate
85 ClassDef(dNdEtaAnalysis, 3)