]>
Commit | Line | Data |
---|---|---|
d66d46f7 | 1 | /* |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. | |
3 | * See cxx source for full Copyright notice | |
4 | * $Id$ | |
5 | */ | |
6 | ||
7 | /********************************************************** | |
8 | * In this class azimuthal correlators in mixed harmonics * | |
9 | * are implemented in terms of Q-vectors. This approach * | |
10 | * doesn't require evaluation of nested loops. This class * | |
11 | * can be used to: * | |
12 | * * | |
13 | * a) Extract subdominant harmonics (like v1 and v4); * | |
5c078806 | 14 | * b) Study flow of two-particle resonances; * |
15 | * c) Study strong parity violation. * | |
d66d46f7 | 16 | * * |
17 | * Author: Ante Bilandzic (abilandzic@gmail.com) * | |
18 | *********************************************************/ | |
19 | ||
20 | #ifndef ALIFLOWANALYSISWITHMIXEDHARMONICS_H | |
21 | #define ALIFLOWANALYSISWITHMIXEDHARMONICS_H | |
22 | ||
23 | #include "AliFlowCommonConstants.h" // needed as include | |
24 | #include "TMatrixD.h" | |
25 | ||
929098e4 | 26 | class TDirectoryFile; |
d66d46f7 | 27 | class TList; |
28 | class TFile; | |
929098e4 | 29 | class TH1F; |
30 | class TH1D; | |
d66d46f7 | 31 | class TH2; |
67a6b178 | 32 | class TH2D; |
d66d46f7 | 33 | class TProfile; |
34 | class TProfile2D; | |
35 | ||
36 | class AliFlowEventSimple; | |
37 | class AliFlowCommonHist; | |
38 | class AliFlowCommonHistResults; | |
39 | ||
40 | //================================================================================================================ | |
41 | ||
42 | class AliFlowAnalysisWithMixedHarmonics | |
43 | { | |
44 | public: | |
45 | AliFlowAnalysisWithMixedHarmonics(); | |
46 | virtual ~AliFlowAnalysisWithMixedHarmonics(); | |
47 | // 0.) Methods called in the constructor: | |
5c078806 | 48 | virtual void InitializeArrays(); |
d66d46f7 | 49 | // 1.) Method Init() and methods called within Init(): |
50 | virtual void Init(); | |
51 | virtual void CrossCheckSettings(); | |
52 | virtual void AccessConstants(); | |
53 | virtual void BookAndNestAllLists(); | |
54 | virtual void BookProfileHoldingSettings(); | |
55 | virtual void BookCommonHistograms(); | |
56 | virtual void BookAllEventByEventQuantities(); | |
57 | virtual void BookAllAllEventQuantities(); | |
58 | virtual void BookAndFillWeightsHistograms(); | |
5c078806 | 59 | // 2.) Method Make() and methods called within Make(): |
d66d46f7 | 60 | virtual void Make(AliFlowEventSimple *anEvent); |
61 | virtual void CheckPointersUsedInMake(); | |
62 | virtual void Calculate3pCorrelator(); | |
63 | virtual void CalculateNonIsotropicTerms(); | |
5c078806 | 64 | virtual void CalculateDifferential3pCorrelator(); |
d66d46f7 | 65 | virtual void ResetEventByEventQuantities(); |
66 | // 3.) Method Finish() and methods called within Finish(): | |
67 | virtual void Finish(); | |
68 | virtual void CheckPointersUsedInFinish(); | |
69 | virtual void AccessSettings(); | |
d66d46f7 | 70 | virtual void CorrectForDetectorEffects(); |
5c078806 | 71 | virtual void PrintOnTheScreen(); |
d66d46f7 | 72 | // 4.) Method GetOutputHistograms and method called within it: |
73 | virtual void GetOutputHistograms(TList *outputListHistos); | |
ecac11c2 | 74 | virtual void GetPointersForBaseHistograms(); |
d66d46f7 | 75 | virtual void GetPointersForCommonHistograms(); |
ecac11c2 | 76 | virtual void GetPointersForAllEventProfiles(); |
d66d46f7 | 77 | virtual void GetPointersForResultsHistograms(); |
78 | // 5.) Other methods: | |
79 | virtual void WriteHistograms(TString outputFileName); | |
80 | virtual void WriteHistograms(TDirectoryFile *outputFileName); | |
81 | // 6.) Setters and getters: | |
82 | void SetHistList(TList* const hl) {this->fHistList = hl;} | |
83 | TList* GetHistList() const {return this->fHistList;} | |
84 | void SetHistListName(const char *hln) {this->fHistListName->Append(*hln);}; | |
85 | TString *GetHistListName() const {return this->fHistListName;}; | |
86 | void SetAnalysisLabel(const char *al) {this->fAnalysisLabel->Append(*al);}; | |
87 | TString *GetAnalysisLabel() const {return this->fAnalysisLabel;}; | |
88 | void SetAnalysisSettings(TProfile* const as) {this->fAnalysisSettings = as;}; | |
ecac11c2 | 89 | TProfile* GetAnalysisSettings() const {return this->fAnalysisSettings;}; |
d66d46f7 | 90 | void SetCorrelatorInteger(Int_t const ci) {this->fCorrelatorInteger = ci;}; |
91 | Int_t GetCorrelatorInteger() const {return this->fCorrelatorInteger;}; | |
92 | void SetNoOfMultipicityBins(Int_t const nomb) {this->fNoOfMultipicityBins = nomb;}; | |
93 | Int_t GetNoOfMultipicityBins() const {return this->fNoOfMultipicityBins;}; | |
94 | void SetMultipicityBinWidth(Double_t const mbw) {this->fMultipicityBinWidth = mbw;}; | |
95 | Double_t GetMultipicityBinWidth() const {return this->fMultipicityBinWidth;}; | |
96 | void SetMinMultiplicity(Double_t const mm) {this->fMinMultiplicity = mm;}; | |
97 | Double_t GetMinMultiplicity() const {return this->fMinMultiplicity;}; | |
98 | void SetCorrectForDetectorEffects(Bool_t const cfde) {this->fCorrectForDetectorEffects = cfde;}; | |
99 | Bool_t GetCorrectForDetectorEffects() const {return this->fCorrectForDetectorEffects;}; | |
5c078806 | 100 | void SetPrintOnTheScreen(Bool_t const pots) {this->fPrintOnTheScreen = pots;}; |
101 | Bool_t GetPrintOnTheScreen() const {return this->fPrintOnTheScreen;}; | |
d66d46f7 | 102 | void SetCommonHists(AliFlowCommonHist* const ch) {this->fCommonHists = ch;}; |
103 | AliFlowCommonHist* GetCommonHists() const {return this->fCommonHists;}; | |
104 | void SetWeightsList(TList* const wl) {this->fWeightsList = (TList*)wl->Clone();} | |
105 | TList* GetWeightsList() const {return this->fWeightsList;} | |
106 | void SetUsePhiWeights(Bool_t const uPhiW) {this->fUsePhiWeights = uPhiW;}; | |
107 | Bool_t GetUsePhiWeights() const {return this->fUsePhiWeights;}; | |
108 | void SetUsePtWeights(Bool_t const uPtW) {this->fUsePtWeights = uPtW;}; | |
109 | Bool_t GetUsePtWeights() const {return this->fUsePtWeights;}; | |
110 | void SetUseEtaWeights(Bool_t const uEtaW) {this->fUseEtaWeights = uEtaW;}; | |
111 | Bool_t GetUseEtaWeights() const {return this->fUseEtaWeights;}; | |
112 | void SetUseParticleWeights(TProfile* const uPW) {this->fUseParticleWeights = uPW;}; | |
113 | TProfile* GetUseParticleWeights() const {return this->fUseParticleWeights;}; | |
114 | void SetPhiWeights(TH1F* const histPhiWeights) {this->fPhiWeights = histPhiWeights;}; | |
115 | TH1F* GetPhiWeights() const {return this->fPhiWeights;}; | |
116 | void SetPtWeights(TH1D* const histPtWeights) {this->fPtWeights = histPtWeights;}; | |
117 | TH1D* GetPtWeights() const {return this->fPtWeights;}; | |
118 | void SetEtaWeights(TH1D* const histEtaWeights) {this->fEtaWeights = histEtaWeights;}; | |
119 | TH1D* GetEtaWeights() const {return this->fEtaWeights;}; | |
ecac11c2 | 120 | void SetProfileList(TList* const plist) {this->fProfileList = plist;} |
121 | TList* GetProfileList() const {return this->fProfileList;} | |
d66d46f7 | 122 | void Set3pCorrelatorPro(TProfile* const s3pPro) {this->f3pCorrelatorPro = s3pPro;}; |
123 | TProfile* Get3pCorrelatorPro() const {return this->f3pCorrelatorPro;}; | |
ecac11c2 | 124 | void SetNonIsotropicTermsPro(TProfile* const nitPro) {this->fNonIsotropicTermsPro = nitPro;}; |
d66d46f7 | 125 | TProfile* GetNonIsotropicTermsPro() const {return this->fNonIsotropicTermsPro;}; |
ecac11c2 | 126 | void Set3pCorrelatorVsMPro(TProfile* const s3pVsMPro) {this->f3pCorrelatorVsMPro = s3pVsMPro;}; |
d66d46f7 | 127 | TProfile* Get3pCorrelatorVsMPro() const {return this->f3pCorrelatorVsMPro;}; |
ecac11c2 | 128 | void SetNonIsotropicTermsVsMPro(TProfile2D* const nitVsMPro) {this->fNonIsotropicTermsVsMPro = nitVsMPro;}; |
d66d46f7 | 129 | TProfile2D* GetNonIsotropicTermsVsMPro() const {return this->fNonIsotropicTermsVsMPro;}; |
130 | void SetResultsList(TList* const rlist) {this->fResultsList = rlist;} | |
5c078806 | 131 | TList* GetResultsList() const {return this->fResultsList;} |
132 | void Set3pCorrelatorHist(TH1D* const s3pHist) {this->f3pCorrelatorHist = s3pHist;}; | |
133 | TH1D* Get3pCorrelatorHist() const {return this->f3pCorrelatorHist;}; | |
d66d46f7 | 134 | void SetDetectorBiasHist(TH1D* const dbHist) {this->fDetectorBiasHist = dbHist;}; |
135 | TH1D* GetDetectorBiasHist() const {return this->fDetectorBiasHist;}; | |
d66d46f7 | 136 | void SetDetectorBiasVsMHist(TH1D* const dbVsMHist) {this->fDetectorBiasVsMHist = dbVsMHist;}; |
137 | TH1D* GetDetectorBiasVsMHist() const {return this->fDetectorBiasVsMHist;}; | |
5c078806 | 138 | void Set3pCorrelatorVsPtSumDiffPro(TProfile* const s3pcvpsd, Int_t const sd) {this->f3pCorrelatorVsPtSumDiffPro[sd] = s3pcvpsd;}; |
139 | TProfile* Get3pCorrelatorVsPtSumDiffPro(Int_t sd) const {return this->f3pCorrelatorVsPtSumDiffPro[sd];}; | |
140 | ||
d66d46f7 | 141 | private: |
142 | AliFlowAnalysisWithMixedHarmonics(const AliFlowAnalysisWithMixedHarmonics& afawQc); | |
143 | AliFlowAnalysisWithMixedHarmonics& operator=(const AliFlowAnalysisWithMixedHarmonics& afawQc); | |
144 | // 0.) Base: | |
145 | TList *fHistList; // base list to hold all output objects | |
146 | TString *fHistListName; // name of base list | |
147 | TString *fAnalysisLabel; // analysis label | |
148 | TProfile *fAnalysisSettings; // profile to hold analysis settings | |
149 | Int_t fCorrelatorInteger; // integer n in cos[n(2phi1-phi2-phi3)] | |
150 | Int_t fNoOfMultipicityBins; // number of multiplicity bins | |
151 | Double_t fMultipicityBinWidth; // width of multiplicity bin | |
152 | Double_t fMinMultiplicity; // minimal multiplicity | |
153 | Bool_t fCorrectForDetectorEffects; // correct 3-p correlator for detector effects | |
5c078806 | 154 | Bool_t fPrintOnTheScreen; // print or not the final results on the screen |
d66d46f7 | 155 | // 1.) Common: |
156 | AliFlowCommonHist *fCommonHists; // common control histograms (filled only with events with 3 or more tracks for 3-p correlators) | |
157 | Int_t fnBinsPhi; // number of phi bins | |
158 | Double_t fPhiMin; // minimum phi | |
159 | Double_t fPhiMax; // maximum phi | |
160 | Double_t fPhiBinWidth; // bin width for phi histograms | |
161 | Int_t fnBinsPt; // number of pt bins | |
162 | Double_t fPtMin; // minimum pt | |
163 | Double_t fPtMax; // maximum pt | |
164 | Double_t fPtBinWidth; // bin width for pt histograms | |
165 | Int_t fnBinsEta; // number of eta bins | |
166 | Double_t fEtaMin; // minimum eta | |
167 | Double_t fEtaMax; // maximum eta | |
168 | Double_t fEtaBinWidth; // bin width for eta histograms | |
169 | // 2a.) Particle weights: | |
170 | TList *fWeightsList; // list to hold all histograms with particle weights: fUseParticleWeights, fPhiWeights, fPtWeights and fEtaWeights | |
171 | Bool_t fUsePhiWeights; // use phi weights | |
172 | Bool_t fUsePtWeights; // use pt weights | |
173 | Bool_t fUseEtaWeights; // use eta weights | |
174 | TProfile *fUseParticleWeights; // profile with three bins to hold values of fUsePhiWeights, fUsePtWeights and fUseEtaWeights | |
175 | TH1F *fPhiWeights; // histogram holding phi weights | |
176 | TH1D *fPtWeights; // histogram holding phi weights | |
177 | TH1D *fEtaWeights; // histogram holding phi weights | |
178 | // 3.) Event-by-event quantities: | |
179 | TMatrixD *fReQnk; // fReQ[n][k] = Re[Q_{n,k}] = sum_{i=1}^{M} w_{i}^{k} cos(n*phi_{i}) | |
180 | TMatrixD *fImQnk; // fImQ[n][k] = Im[Q_{n,k}] = sum_{i=1}^{M} w_{i}^{k} sin(n*phi_{i}) | |
181 | TMatrixD *fSpk; // fS[p][k] = S_{p,k} = (sum_{i=1}^{M} w_{i}^{k})^{p+1} // note p+1 in the power to use 0th index in p in non-trivial way | |
182 | TH1D *f3pCorrelatorEBE; // 3-p correlator <cos[n(2phi1-phi2-phi3)]> for single event | |
183 | TH1D *fNonIsotropicTermsEBE; // correction terms to 3-p correlator <cos[n(2phi1-phi2-phi3)]> for single event | |
5c078806 | 184 | TProfile *fRePEBE[2]; // real part of p_n vs [(p1+p2)/2,|p1-p2|] |
185 | TProfile *fImPEBE[2]; // imaginary part of p_n vs [(p1+p2)/2,|p1-p2|] | |
d66d46f7 | 186 | // 4.) Profiles: |
ecac11c2 | 187 | TList *fProfileList; // list holding all all-event profiles |
5c078806 | 188 | TProfile *f3pCorrelatorPro; // 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> not corrected for detector effects |
189 | TProfile *fNonIsotropicTermsPro; // non-isotropic terms in the decomposition of 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> | |
190 | TProfile *f3pCorrelatorVsMPro; // 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> vs multiplicity | |
191 | TProfile2D *fNonIsotropicTermsVsMPro; // non-isotropic terms in the decomposition of <cos[n(2phi1-phi2-phi3)]> vs multiplicity | |
192 | TProfile *f3pCorrelatorVsPtSumDiffPro[2]; // differential 3-p correlator <<cos[n(2phi1-psi2-psi3)]>> vs [(p1+p2)/2,|p1-p2|] | |
d66d46f7 | 193 | // 5.) Final results: |
194 | TList *fResultsList; // list holding objects with final results | |
5c078806 | 195 | TH1D *f3pCorrelatorHist; // 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> corrected for detector effects |
d66d46f7 | 196 | TH1D *fDetectorBiasHist; // bias comming from detector inefficiencies to 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> (in %) |
d66d46f7 | 197 | TH1D *fDetectorBiasVsMHist; // bias comming from detector inefficiencies to 3-p correlator <<cos[n(2phi1-phi2-phi3)]>> (in %) versus multiplicity |
5c078806 | 198 | |
d66d46f7 | 199 | ClassDef(AliFlowAnalysisWithMixedHarmonics, 0); |
200 | ||
201 | }; | |
202 | ||
203 | //================================================================================================================ | |
204 | ||
205 | #endif | |
206 | ||
207 | ||
208 | ||
209 | ||
210 |