AliPhysics  2c6b7ad (2c6b7ad)
AliAnalysisTaskQCumulants.h
Go to the documentation of this file.
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  * analysis task for Q-cumulants *
9  * *
10  * authors: Naomi van der Kolk *
11  * (kolk@nikhef.nl) *
12  * Raimond Snellings *
13  * (snelling@nikhef.nl) *
14  * Ante Bilandzic *
15  * (anteb@nikhef.nl) *
16  * ***********************************/
17 
18 #ifndef ALIANALYSISTASKQCUMULANTS_H
19 #define ALIANALYSISTASKQCUMULANTS_H
20 
21 #include "AliAnalysisTaskSE.h"
22 #include "AliFlowCommonConstants.h"
23 
24 class TString;
25 class TList;
26 class AliFlowEventSimple;
28 
29 //================================================================================================================
30 
32  public:
34  AliAnalysisTaskQCumulants(const char *name, Bool_t useParticleWeights=kFALSE);
36 
37  virtual void UserCreateOutputObjects();
38  virtual void UserExec(Option_t *option);
39  virtual void Terminate(Option_t *);
40 
41  // Common:
42  void SetBookOnlyBasicCCH(Bool_t const bobcch) {this->fBookOnlyBasicCCH = bobcch;};
43  Bool_t GetBookOnlyBasicCCH() const {return this->fBookOnlyBasicCCH;};
48  void SetHarmonic(Int_t const harmonic) {this->fHarmonic = harmonic;};
49  Int_t GetHarmonic() const {return this->fHarmonic;};
50  void SetApplyCorrectionForNUA(Bool_t const applyCorrectionForNUA) {this->fApplyCorrectionForNUA = applyCorrectionForNUA;};
52  void SetApplyCorrectionForNUAVsM(Bool_t const applyCorrectionForNUAVsM) {this->fApplyCorrectionForNUAVsM = applyCorrectionForNUAVsM;};
54  void SetPropagateErrorAlsoFromNIT(Bool_t const peafNIT) {this->fPropagateErrorAlsoFromNIT = peafNIT;};
56  void SetCalculateDiffFlow(Bool_t const calculateDiffFlow) {this->fCalculateDiffFlow = calculateDiffFlow;};
58  void SetCalculate2DDiffFlow(Bool_t const calculate2DDiffFlow) {this->fCalculate2DDiffFlow = calculate2DDiffFlow;};
60  void SetCalculateDiffFlowVsEta(Bool_t const cdfve) {this->fCalculateDiffFlowVsEta = cdfve;};
62  void SetStoreDistributions(Bool_t const storeDistributions) {this->fStoreDistributions = storeDistributions;};
64  void SetCalculateCumulantsVsM(Bool_t const ccvm) {this->fCalculateCumulantsVsM = ccvm;};
82  void SetExactNoRPs(Int_t const enr) {this->fExactNoRPs = enr;};
83  Int_t GetExactNoRPs() const {return this->fExactNoRPs;};
84  void SetUse2DHistograms(Bool_t const u2dh){this->fUse2DHistograms = u2dh;if(u2dh){this->fStoreControlHistograms = kTRUE;}};
85  Bool_t GetUse2DHistograms() const {return this->fUse2DHistograms;};
88  void SetUseQvectorTerms(Bool_t const uqvt){this->fUseQvectorTerms = uqvt;if(uqvt){this->fStoreControlHistograms = kTRUE;}};
89  Bool_t GetUseQvectorTerms() const {return this->fUseQvectorTerms;};
90 
91  // Multiparticle correlations vs multiplicity:
92  void SetnBinsMult(Int_t const nbm) {this->fnBinsMult = nbm;};
93  Int_t GetnBinsMult() const {return this->fnBinsMult;};
94  void SetMinMult(Double_t const minm) {this->fMinMult = minm;};
95  Double_t GetMinMult() const {return this->fMinMult;};
96  void SetMaxMult(Double_t const maxm) {this->fMaxMult = maxm;};
97  Double_t GetMaxMult() const {return this->fMaxMult;};
98  // Particle weights:
99  void SetUsePhiWeights(Bool_t const uPhiW) {this->fUsePhiWeights = uPhiW;};
100  Bool_t GetUsePhiWeights() const {return this->fUsePhiWeights;};
101  void SetUsePtWeights(Bool_t const uPtW) {this->fUsePtWeights = uPtW;};
102  Bool_t GetUsePtWeights() const {return this->fUsePtWeights;};
103  void SetUseEtaWeights(Bool_t const uEtaW) {this->fUseEtaWeights = uEtaW;};
104  Bool_t GetUseEtaWeights() const {return this->fUseEtaWeights;};
105  void SetUseTrackWeights(Bool_t const uTrackW) {this->fUseTrackWeights = uTrackW;};
106  Bool_t GetUseTrackWeights() const {return this->fUseTrackWeights;};
107  // Event weights:
108  void SetMultiplicityWeight(const char *multiplicityWeight) {*this->fMultiplicityWeight = multiplicityWeight;};
110  // # of bins for correlation axis in fDistributions[4], fCorrelation2468VsMult[4] and fCorrelationProduct2468VsMult[1]
113  // Boundaries for distributions of correlations:
114  void SetMinValueOfCorrelation(Int_t const ci, Double_t const minValue) {this->fMinValueOfCorrelation[ci] = minValue;};
116  void SetMaxValueOfCorrelation(Int_t const ci, Double_t const maxValue) {this->fMaxValueOfCorrelation[ci] = maxValue;};
118  // min and max values of correlation products:
119  void SetMinValueOfCorrelationProduct(Int_t const cpi, Double_t const minValue) {this->fMinValueOfCorrelationProduct[cpi] = minValue;};
121  void SetMaxValueOfCorrelationProduct(Int_t const cpi, Double_t const maxValue) {this->fMaxValueOfCorrelationProduct[cpi] = maxValue;};
123  // min and max values of QvectorTerms:
124  void SetMinValueOfQvectorTerms(Int_t const qvti, Double_t const minValue) {this->fMinValueOfQvectorTerms[qvti] = minValue;};
126  void SetMaxValueOfQvectorTerms(Int_t const qvti, Double_t const maxValue) {this->fMaxValueOfQvectorTerms[qvti] = maxValue;};
128  // bootstrap:
129  void SetUseBootstrap(Bool_t const ub) {this->fUseBootstrap = ub;};
130  Bool_t GetUseBootstrap() const {return this->fUseBootstrap;};
131  void SetUseBootstrapVsM(Bool_t const ubVsM) {this->fUseBootstrapVsM = ubVsM;};
132  Bool_t GetUseBootstrapVsM() const {return this->fUseBootstrapVsM;};
133  void SetnSubsamples(Int_t const ns) {this->fnSubsamples = ns;};
134  Int_t GetnSubsamples() const {return this->fnSubsamples;};
135 
136  private:
139 
140  AliFlowEventSimple *fEvent; // the input event
141  AliFlowAnalysisWithQCumulants *fQC; // Q-cumulant object
142  TList *fListHistos; // collection of output
143  // Common:
144  Bool_t fBookOnlyBasicCCH; // book only basis common control histrograms (by default book them all)
145  Bool_t fFillMultipleControlHistograms; // fill separately control histos for events with >= 2, 4, 6 and 8 particles
146  Int_t fMaxCommonResultsHistogram; // can be [2468], e.g. if set to 2, AliFlowCommonHistResults[468]thOrderQC won't be booked
147  Int_t fHarmonic; // harmonic
148  Bool_t fApplyCorrectionForNUA; // apply correction for non-uniform acceptance
149  Bool_t fApplyCorrectionForNUAVsM; // apply correction for non-uniform acceptance versus M
150  Bool_t fPropagateErrorAlsoFromNIT; // propagate error by taking into account also non-isotrpic terms
151  Bool_t fCalculateDiffFlow; // calculate differential flow in pt or eta
152  Bool_t fCalculate2DDiffFlow; // calculate differential flow in (pt,eta) (Remark: this is very expensive in terms of CPU time)
153  Bool_t fCalculateDiffFlowVsEta; // if you set kFALSE only differential flow vs pt is calculated
154  Bool_t fStoreDistributions; // store or not distributions of correlations
155  Bool_t fCalculateCumulantsVsM; // calculate cumulants versus multiplicity
156  Bool_t fCalculateAllCorrelationsVsM; // calculate all correlations versus multiplicity
157  Bool_t fCalculateMixedHarmonics; // calculate all mixed harmonics correlations
158  Bool_t fCalculateMixedHarmonicsVsM; // calculate all mixed harmonics correlations versus multiplicity
159  Bool_t fStoreControlHistograms; // store or not control histograms
160  Bool_t fMinimumBiasReferenceFlow; // store as reference flow in AliFlowCommonHistResults the minimum bias result (kFALSE by default)
161  Bool_t fForgetAboutCovariances; // when propagating error forget about the covariances
162  Bool_t fStorePhiDistributionForOneEvent; // store phi distribution for one event to illustrate flow
163  Double_t fPhiDistributionForOneEventSettings[4]; // [v_min,v_max,refMult_min,refMult_max]
164  Int_t fExactNoRPs; // when shuffled, select only this number of RPs for the analysis
165  Bool_t fUse2DHistograms; // use TH2D instead of TProfile to improve numerical stability in reference flow calculation
166  Bool_t fFillProfilesVsMUsingWeights; // if the width of multiplicity bin is 1, weights are not needed
167  Bool_t fUseQvectorTerms; // use TH2D with separate Q-vector terms instead of TProfile to improve numerical stability in reference flow calculation
168  // Multiparticle correlations vs multiplicity:
169  Int_t fnBinsMult; // number of multiplicity bins for flow analysis versus multiplicity
170  Double_t fMinMult; // minimal multiplicity for flow analysis versus multiplicity
171  Double_t fMaxMult; // maximal multiplicity for flow analysis versus multiplicity
172  // Particle weights:
173  Bool_t fUseParticleWeights; // use any particle weights
174  Bool_t fUsePhiWeights; // use phi weights
175  Bool_t fUsePtWeights; // use pt weights
176  Bool_t fUseEtaWeights; // use eta weights
177  Bool_t fUseTrackWeights; // use track weights (e.g. VZERO sector weights)
178  TList *fWeightsList; // list with weights
179  // Event weights:
180  TString *fMultiplicityWeight; // event-by-event weights for multiparticle correlations ("combinations","unit" or "multiplicity")
181  AliFlowCommonConstants::ERefMultSource fMultiplicityIs; // by default "#RPs", other supported options are "RefMultFromESD" = ref. mult. from ESD, and "#POIs"
182  Int_t fnBinsForCorrelations; // # of bins for correlation axis in fDistributions[4], fCorrelation2468VsMult[4] and fCorrelationProduct2468VsMult[1]
183  // Boundaries for distributions of correlations:
184  Double_t fMinValueOfCorrelation[4]; // min values of <2>, <4>, <6> and <8>
185  Double_t fMaxValueOfCorrelation[4]; // max values of <2>, <4>, <6> and <8>
186  Double_t fMinValueOfCorrelationProduct[1]; // min values of <2><4>, <2><6>, <2><8>, <4><6> etc. TBI add the other ones when needed first time
187  Double_t fMaxValueOfCorrelationProduct[1]; // max values of <2><4>, <2><6>, <2><8>, <4><6> etc. TBI add the other ones when needed first time
188  Double_t fMinValueOfQvectorTerms[4]; // min value of Q-vector terms
189  Double_t fMaxValueOfQvectorTerms[4]; // max value of Q-vector terms
190  // Bootstrap:
191  Bool_t fUseBootstrap; // use bootstrap to estimate statistical spread
192  Bool_t fUseBootstrapVsM; // use bootstrap to estimate statistical spread for results vs M
193  Int_t fnSubsamples; // number of subsamples (SS), by default 10
194 
195  ClassDef(AliAnalysisTaskQCumulants, 2);
196 };
197 
198 //================================================================================================================
199 
200 #endif
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
void SetUseBootstrapVsM(Bool_t const ubVsM)
Bool_t GetCalculateAllCorrelationsVsM() const
void SetApplyCorrectionForNUA(Bool_t const applyCorrectionForNUA)
double Double_t
Definition: External.C:58
void SetMinMult(Double_t const minm)
void SetCalculateMixedHarmonicsVsM(Bool_t const cmhvm)
void SetMultiplicityWeight(const char *multiplicityWeight)
void SetUseBootstrap(Bool_t const ub)
void SetForgetAboutCovariances(Bool_t const fac)
void SetCalculateDiffFlow(Bool_t const calculateDiffFlow)
Double_t GetMinValueOfQvectorTerms(Int_t qvti) const
void SetStorePhiDistributionForOneEvent(Bool_t const spdfoe)
void SetCalculateDiffFlowVsEta(Bool_t const cdfve)
void SetMinimumBiasReferenceFlow(Bool_t const mmrf)
AliFlowCommonConstants::ERefMultSource fMultiplicityIs
void SetCalculateMixedHarmonics(Bool_t const cmh)
void SetUseTrackWeights(Bool_t const uTrackW)
void SetMaxCommonResultsHistogram(Int_t const mcrh)
virtual void UserExec(Option_t *option)
Bool_t GetStorePhiDistributionForOneEvent() const
void SetPhiDistributionForOneEventSettings(Double_t const pdfoes, Int_t const i)
void SetUse2DHistograms(Bool_t const u2dh)
void SetMaxValueOfCorrelation(Int_t const ci, Double_t const maxValue)
int Int_t
Definition: External.C:63
void SetMinValueOfQvectorTerms(Int_t const qvti, Double_t const minValue)
Double_t GetMinValueOfCorrelationProduct(Int_t cpi) const
void SetMultiplicityIs(AliFlowCommonConstants::ERefMultSource mi)
void SetMinValueOfCorrelation(Int_t const ci, Double_t const minValue)
void SetUsePhiWeights(Bool_t const uPhiW)
void SetnBinsMult(Int_t const nbm)
AliAnalysisTaskQCumulants & operator=(const AliAnalysisTaskQCumulants &aatqc)
void SetMaxValueOfCorrelationProduct(Int_t const cpi, Double_t const maxValue)
void SetUseQvectorTerms(Bool_t const uqvt)
void SetCalculateCumulantsVsM(Bool_t const ccvm)
void SetnBinsForCorrelations(Int_t const nb)
void SetUseEtaWeights(Bool_t const uEtaW)
void SetFillMultipleControlHistograms(Bool_t const fmch)
void SetBookOnlyBasicCCH(Bool_t const bobcch)
void SetPropagateErrorAlsoFromNIT(Bool_t const peafNIT)
void SetMaxValueOfQvectorTerms(Int_t const qvti, Double_t const maxValue)
void SetMinValueOfCorrelationProduct(Int_t const cpi, Double_t const minValue)
AliFlowAnalysisWithQCumulants * fQC
void SetCalculate2DDiffFlow(Bool_t const calculate2DDiffFlow)
Double_t GetMaxValueOfCorrelationProduct(Int_t cpi) const
Double_t GetMinValueOfCorrelation(Int_t ci) const
virtual void Terminate(Option_t *)
void SetCalculateAllCorrelationsVsM(Bool_t const cacvm)
const char Option_t
Definition: External.C:48
void SetHarmonic(Int_t const harmonic)
void SetApplyCorrectionForNUAVsM(Bool_t const applyCorrectionForNUAVsM)
Double_t GetPhiDistributionForOneEventSettings(Int_t const i) const
bool Bool_t
Definition: External.C:53
void SetMaxMult(Double_t const maxm)
Double_t GetMaxValueOfQvectorTerms(Int_t qvti) const
Bool_t GetFillProfilesVsMUsingWeights() const
void SetStoreControlHistograms(Bool_t const sch)
Bool_t GetFillMultipleControlHistograms() const
void SetFillProfilesVsMUsingWeights(Bool_t const fpvmuw)
Double_t GetMaxValueOfCorrelation(Int_t ci) const
void SetStoreDistributions(Bool_t const storeDistributions)
void SetUsePtWeights(Bool_t const uPtW)
void SetExactNoRPs(Int_t const enr)