PWG2/FLOW/AliFlowTasks/AliAnalysisTaskQCumulants.h

   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
  23 class TString;
  24 class TList;
  25 class AliFlowEventSimple;
  26 class AliFlowAnalysisWithQCumulants;
  27
  28 //================================================================================================================
  29
  30 class AliAnalysisTaskQCumulants : public AliAnalysisTaskSE{
  31  public:
  32   AliAnalysisTaskQCumulants();
  33   AliAnalysisTaskQCumulants(const char *name, Bool_t useParticleWeights=kFALSE);
  34   virtual ~AliAnalysisTaskQCumulants(){};
  35
  36   virtual void UserCreateOutputObjects();
  37   virtual void UserExec(Option_t *option);
  38   virtual void Terminate(Option_t *);
  39
  40   // Common:
  41   void SetBookOnlyBasicCCH(Bool_t const bobcch) {this->fBookOnlyBasicCCH = bobcch;};
  42   Bool_t GetBookOnlyBasicCCH() const {return this->fBookOnlyBasicCCH;};
  43   void SetFillMultipleControlHistograms(Bool_t const fmch) {this->fFillMultipleControlHistograms = fmch;};
  44   Bool_t GetFillMultipleControlHistograms() const {return this->fFillMultipleControlHistograms;};
  45   void SetHarmonic(Int_t const harmonic) {this->fHarmonic = harmonic;};
  46   Int_t GetHarmonic() const {return this->fHarmonic;};
  47   void SetApplyCorrectionForNUA(Bool_t const applyCorrectionForNUA) {this->fApplyCorrectionForNUA = applyCorrectionForNUA;};
  48   Bool_t GetApplyCorrectionForNUA() const {return this->fApplyCorrectionForNUA;};
  49   void SetApplyCorrectionForNUAVsM(Bool_t const applyCorrectionForNUAVsM) {this->fApplyCorrectionForNUAVsM = applyCorrectionForNUAVsM;};
  50   Bool_t GetApplyCorrectionForNUAVsM() const {return this->fApplyCorrectionForNUAVsM;};
  51   void SetPropagateErrorAlsoFromNIT(Bool_t const peafNIT) {this->fPropagateErrorAlsoFromNIT = peafNIT;};
  52   Bool_t GetPropagateErrorAlsoFromNIT() const {return this->fPropagateErrorAlsoFromNIT;};
  53   void SetCalculateDiffFlow(Bool_t const calculateDiffFlow) {this->fCalculateDiffFlow = calculateDiffFlow;};
  54   Bool_t GetCalculateDiffFlow() const {return this->fCalculateDiffFlow;};
  55   void SetCalculate2DDiffFlow(Bool_t const calculate2DDiffFlow) {this->fCalculate2DDiffFlow = calculate2DDiffFlow;};
  56   Bool_t GetCalculate2DDiffFlow() const {return this->fCalculate2DDiffFlow;};
  57   void SetCalculateDiffFlowVsEta(Bool_t const cdfve) {this->fCalculateDiffFlowVsEta = cdfve;};
  58   Bool_t GetCalculateDiffFlowVsEta() const {return this->fCalculateDiffFlowVsEta;};
  59   void SetStoreDistributions(Bool_t const storeDistributions) {this->fStoreDistributions = storeDistributions;};
  60   Bool_t GetStoreDistributions() const {return this->fStoreDistributions;};
  61   void SetCalculateCumulantsVsM(Bool_t const ccvm) {this->fCalculateCumulantsVsM = ccvm;};
  62   Bool_t GetCalculateCumulantsVsM() const {return this->fCalculateCumulantsVsM;};
  63   void SetCalculateAllCorrelationsVsM(Bool_t const cacvm) {this->fCalculateAllCorrelationsVsM = cacvm;};
  64   Bool_t GetCalculateAllCorrelationsVsM() const {return this->fCalculateAllCorrelationsVsM;};
  65   void SetMinimumBiasReferenceFlow(Bool_t const mmrf) {this->fMinimumBiasReferenceFlow = mmrf;};
  66   Bool_t GetMinimumBiasReferenceFlow() const {return this->fMinimumBiasReferenceFlow;};
  67   void SetForgetAboutCovariances(Bool_t const fac) {this->fForgetAboutCovariances = fac;};
  68   Bool_t GetForgetAboutCovariances() const {return this->fForgetAboutCovariances;};
  69   void SetStorePhiDistributionForOneEvent(Bool_t const spdfoe) {this->fStorePhiDistributionForOneEvent = spdfoe;};
  70   Bool_t GetStorePhiDistributionForOneEvent() const {return this->fStorePhiDistributionForOneEvent;};
  71   void SetPhiDistributionForOneEventSettings(Double_t const pdfoes, Int_t const i) {this->fPhiDistributionForOneEventSettings[i] = pdfoes;};
  72   Double_t GetPhiDistributionForOneEventSettings(Int_t const i) const {return this->fPhiDistributionForOneEventSettings[i];};
  73
  74   // Multiparticle correlations vs multiplicity:
  75   void SetnBinsMult(Int_t const nbm) {this->fnBinsMult = nbm;};
  76   Int_t GetnBinsMult() const {return this->fnBinsMult;};
  77   void SetMinMult(Double_t const minm) {this->fMinMult = minm;};
  78   Double_t GetMinMult() const {return this->fMinMult;};
  79   void SetMaxMult(Double_t const maxm) {this->fMaxMult = maxm;};
  80   Double_t GetMaxMult() const {return this->fMaxMult;};
  81   // Particle weights:
  82   void SetUsePhiWeights(Bool_t const uPhiW) {this->fUsePhiWeights = uPhiW;};
  83   Bool_t GetUsePhiWeights() const {return this->fUsePhiWeights;};
  84   void SetUsePtWeights(Bool_t const uPtW) {this->fUsePtWeights = uPtW;};
  85   Bool_t GetUsePtWeights() const {return this->fUsePtWeights;};
  86   void SetUseEtaWeights(Bool_t const uEtaW) {this->fUseEtaWeights = uEtaW;};
  87   Bool_t GetUseEtaWeights() const {return this->fUseEtaWeights;};
  88   void SetUseTrackWeights(Bool_t const uTrackW) {this->fUseTrackWeights = uTrackW;};
  89   Bool_t GetUseTrackWeights() const {return this->fUseTrackWeights;};
  90   // Event weights:
  91   void SetMultiplicityWeight(const char *multiplicityWeight) {*this->fMultiplicityWeight = multiplicityWeight;};
  92
  93  private:
  94   AliAnalysisTaskQCumulants(const AliAnalysisTaskQCumulants& aatqc);
  95   AliAnalysisTaskQCumulants& operator=(const AliAnalysisTaskQCumulants& aatqc);
  96
  97   AliFlowEventSimple *fEvent;         // the input event
  98   AliFlowAnalysisWithQCumulants *fQC; // Q-cumulant object
  99   TList *fListHistos;                 // collection of output
 100   // Common:
 101   Bool_t fBookOnlyBasicCCH;              // book only basis common control histrograms (by default book them all)
 102   Bool_t fFillMultipleControlHistograms; // fill separately control histos for events with >= 2, 4, 6 and 8 particles
 103   Int_t fHarmonic;                       // harmonic
 104   Bool_t fApplyCorrectionForNUA;         // apply correction for non-uniform acceptance
 105   Bool_t fApplyCorrectionForNUAVsM;      // apply correction for non-uniform acceptance versus M
 106   Bool_t fPropagateErrorAlsoFromNIT;     // propagate error by taking into account also non-isotrpic terms
 107   Bool_t fCalculateDiffFlow;             // calculate differential flow in pt or eta
 108   Bool_t fCalculate2DDiffFlow;           // calculate differential flow in (pt,eta) (Remark: this is very expensive in terms of CPU time)
 109   Bool_t fCalculateDiffFlowVsEta;        // if you set kFALSE only differential flow vs pt is calculated
 110   Bool_t fStoreDistributions;            // store or not distributions of correlations
 111   Bool_t fCalculateCumulantsVsM;         // calculate cumulants versus multiplicity
 112   Bool_t fCalculateAllCorrelationsVsM;   // calculate all correlations versus multiplicity
 113   Bool_t fMinimumBiasReferenceFlow;      // store as reference flow in AliFlowCommonHistResults the minimum bias result (kFALSE by default)
 114   Bool_t fForgetAboutCovariances;        // when propagating error forget about the covariances
 115   Bool_t fStorePhiDistributionForOneEvent; // store phi distribution for one event to illustrate flow
 116   Double_t fPhiDistributionForOneEventSettings[4]; // [v_min,v_max,refMult_min,refMult_max]
 117   // Multiparticle correlations vs multiplicity:
 118   Int_t fnBinsMult;                   // number of multiplicity bins for flow analysis versus multiplicity
 119   Double_t fMinMult;                  // minimal multiplicity for flow analysis versus multiplicity
 120   Double_t fMaxMult;                  // maximal multiplicity for flow analysis versus multiplicity
 121   // Particle weights:
 122   Bool_t fUseParticleWeights;         // use any particle weights
 123   Bool_t fUsePhiWeights;              // use phi weights
 124   Bool_t fUsePtWeights;               // use pt weights
 125   Bool_t fUseEtaWeights;              // use eta weights
 126   Bool_t fUseTrackWeights;            // use track weights (e.g. VZERO sector weights)
 127   TList *fWeightsList;                // list with weights
 128   // Event weights:
 129   TString *fMultiplicityWeight;       // event-by-event weights for multiparticle correlations ("combinations","unit" or "multiplicity")
 130
 131   ClassDef(AliAnalysisTaskQCumulants, 1);
 132 };
 133
 134 //================================================================================================================
 135
 136 #endif
 137
 138
 139
 140
 141
 142
 143
 144
 145
 146
 147