By default correction for non-uniform acceptance is not applied.

diff --git a/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.cxx b/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.cxx
index 293d5570f922205862e61e27e117e8907c50be36..ac599a6c4baa33b8f5f2472ca1c2783889172745 100644 (file)
--- a/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.cxx
+++ b/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.cxx
@@ -43,11 +43,13 @@ AliAnalysisTaskQCumulants::AliAnalysisTaskQCumulants(const char *name, Bool_t us
  fQC(NULL), 
  fListHistos(NULL),
  fHarmonic(2),  
- fApplyCorrectionForNUA(kTRUE), 
+ fApplyCorrectionForNUA(kFALSE), 
+ fApplyCorrectionForNUAVsM(kFALSE), 
+ fPropagateErrorAlsoFromNUATerms(kFALSE),
  fCalculate2DFlow(kFALSE),
  fStoreDistributions(kFALSE),
  fCalculateCumulantsVsM(kTRUE), 
- fMinimumBiasReferenceFlow(kFALSE), 
+ fMinimumBiasReferenceFlow(kTRUE), 
  fnBinsMult(10000),
  fMinMult(0.),  
  fMaxMult(10000.), 
@@ -84,6 +86,8 @@ AliAnalysisTaskQCumulants::AliAnalysisTaskQCumulants():
  fListHistos(NULL),
  fHarmonic(0),  
  fApplyCorrectionForNUA(kFALSE), 
+ fApplyCorrectionForNUAVsM(kFALSE), 
+ fPropagateErrorAlsoFromNUATerms(kFALSE),
  fCalculate2DFlow(kFALSE),
  fStoreDistributions(kFALSE),
  fCalculateCumulantsVsM(kFALSE),  
@@ -115,6 +119,8 @@ void AliAnalysisTaskQCumulants::UserCreateOutputObjects()
  // Common:
  fQC->SetHarmonic(fHarmonic);
  fQC->SetApplyCorrectionForNUA(fApplyCorrectionForNUA);
+ fQC->SetApplyCorrectionForNUAVsM(fApplyCorrectionForNUAVsM);
+ fQC->SetPropagateErrorAlsoFromNUATerms(fPropagateErrorAlsoFromNUATerms);
  fQC->SetCalculate2DFlow(fCalculate2DFlow);
  fQC->SetStoreDistributions(fStoreDistributions);
  fQC->SetCalculateCumulantsVsM(fCalculateCumulantsVsM);

diff --git a/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.h b/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.h
index 3598e8704f9fc83b95797e68be00daed8e1e1237..c87935426fd58127cbb098808ac7a06a191703a0 100644 (file)
--- a/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.h
+++ b/PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.h
@@ -42,6 +42,10 @@ class AliAnalysisTaskQCumulants : public AliAnalysisTaskSE{
   Int_t GetHarmonic() const {return this->fHarmonic;};
   void SetApplyCorrectionForNUA(Bool_t const applyCorrectionForNUA) {this->fApplyCorrectionForNUA = applyCorrectionForNUA;};
   Bool_t GetApplyCorrectionForNUA() const {return this->fApplyCorrectionForNUA;};
+  void SetApplyCorrectionForNUAVsM(Bool_t const applyCorrectionForNUAVsM) {this->fApplyCorrectionForNUAVsM = applyCorrectionForNUAVsM;};
+  Bool_t GetApplyCorrectionForNUAVsM() const {return this->fApplyCorrectionForNUAVsM;};      
+  void SetPropagateErrorAlsoFromNUATerms(Bool_t const pefc) {this->fPropagateErrorAlsoFromNUATerms = pefc;};
+  Bool_t GetPropagateErrorAlsoFromNUATerms() const {return this->fPropagateErrorAlsoFromNUATerms;};  
   void SetCalculate2DFlow(Bool_t const calculate2DFlow) {this->fCalculate2DFlow = calculate2DFlow;};
   Bool_t GetCalculate2DFlow() const {return this->fCalculate2DFlow;};
   void SetStoreDistributions(Bool_t const storeDistributions) {this->fStoreDistributions = storeDistributions;};
@@ -75,12 +79,14 @@ class AliAnalysisTaskQCumulants : public AliAnalysisTaskSE{
   AliFlowAnalysisWithQCumulants *fQC; // Q-cumulant object
   TList *fListHistos;                 // collection of output 
   // Common:
-  Int_t fHarmonic;                    // harmonic  
-  Bool_t fApplyCorrectionForNUA;      // apply correction for non-uniform acceptance 
-  Bool_t fCalculate2DFlow;            // calculate differential flow in (pt,eta) (Remark: this is very expensive in terms of CPU time)
-  Bool_t fStoreDistributions;         // store or not distributions of correlations
-  Bool_t fCalculateCumulantsVsM;      // calculate cumulants versus multiplicity  
-  Bool_t fMinimumBiasReferenceFlow;   // store as reference flow in AliFlowCommonHistResults the minimum bias result (kFALSE by default)     
+  Int_t fHarmonic;                        // harmonic  
+  Bool_t fApplyCorrectionForNUA;          // apply correction for non-uniform acceptance 
+  Bool_t fApplyCorrectionForNUAVsM;       // apply correction for non-uniform acceptance versus M    
+  Bool_t fPropagateErrorAlsoFromNUATerms; // propagate error by taking into account also non-isotrpic terms  
+  Bool_t fCalculate2DFlow;                // calculate differential flow in (pt,eta) (Remark: this is very expensive in terms of CPU time)
+  Bool_t fStoreDistributions;             // store or not distributions of correlations
+  Bool_t fCalculateCumulantsVsM;          // calculate cumulants versus multiplicity  
+  Bool_t fMinimumBiasReferenceFlow;       // store as reference flow in AliFlowCommonHistResults the minimum bias result (kFALSE by default)     
   // Multiparticle correlations vs multiplicity:
   Int_t fnBinsMult;                   // number of multiplicity bins for flow analysis versus multiplicity  
   Double_t fMinMult;                  // minimal multiplicity for flow analysis versus multiplicity