AliAnalysisTaskSEDvsEventShapes.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSEDVSEVENTSHAPES_H
#define ALIANALYSISTASKSEDVSEVENTSHAPES_H

/* Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

//*************************************************************************
// Class AliAnalysisTaskSEDvsEventShapes
// AliAnalysisTaskSE for the D meson vs. Event shape analysis in different mutiplicity window
// Authors: Renu Bala, Manoj Bhanudas Jadhav
//*************************************************************************

#include <TROOT.h>
#include <TSystem.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <THnSparse.h>
#include <TArrayD.h>
#include <TFile.h>
#include <TRandom.h>
#include <TProfile.h>
#include "AliRDHFCutsDplustoKpipi.h"
#include "AliAnalysisTaskSE.h"
#include "AliAnalysisVertexingHF.h"
#include "AliNormalizationCounter.h"
#include "AliAODMCHeader.h"
#include "AliAODMCParticle.h"
#include "AliVertexingHFUtils.h"
#include "AliVEvent.h"


class AliAnalysisTaskSEDvsEventShapes : public AliAnalysisTaskSE
{
public:
    
    AliAnalysisTaskSEDvsEventShapes();
    AliAnalysisTaskSEDvsEventShapes(const char *name, Int_t pdgMeson, AliRDHFCuts* cuts, Bool_t switchPPb);
    virtual ~AliAnalysisTaskSEDvsEventShapes();
    
    
    void SetMassLimits(Double_t lowlimit, Double_t uplimit);
    void SetMassLimits(Int_t pdg, Double_t range);
    Double_t GetUpperMassLimit() const {return fUpmasslimit;}
    Double_t GetLowerMassLimit() const {return fLowmasslimit;}
    void SetNMassBins(Int_t nbins){fNMassBins=nbins;}
    Int_t GetNMassBins() const{return fNMassBins;}
    Bool_t GetSubtractTrackletsFromDaughters() const {return fSubtractTrackletsFromDau;}
    
    void SetImpactParameterBinning(Int_t nbins, Double_t dmin, Double_t dmax){
        fNImpParBins=nbins;
        fLowerImpPar=dmin;
        fHigherImpPar=dmax;
    }
    
    void SetReadMC(Bool_t readMC=kTRUE){fReadMC=readMC;}
    void SetMCOption(Int_t option=0){ fMCOption = option; }
    void SetIsPPbData(Bool_t flag=kTRUE){
        fisPPbData=flag;
    }
    void SetUseBit(Bool_t use=kTRUE){fUseBit=use;}
    void SetDoImpactParameterHistos(Bool_t doImp=kTRUE){fDoImpPar=doImp;}
    void SetKeepEstimatorCorrelationPlots(Bool_t use=kTRUE){fKeepCorrPlots=use;}
    
    void SetMultiplVsZProfileLHC10b(TProfile* hprof){
        if(fMultEstimatorAvg[0]) delete fMultEstimatorAvg[0];
        fMultEstimatorAvg[0]=new TProfile(*hprof);
    }
    void SetMultiplVsZProfileLHC10c(TProfile* hprof){
        if(fMultEstimatorAvg[1]) delete fMultEstimatorAvg[1];
        fMultEstimatorAvg[1]=new TProfile(*hprof);
    }
    void SetMultiplVsZProfileLHC10d(TProfile* hprof){
        if(fMultEstimatorAvg[2]) delete fMultEstimatorAvg[2];
        fMultEstimatorAvg[2]=new TProfile(*hprof);
    }
    void SetMultiplVsZProfileLHC10e(TProfile* hprof){
        if(fMultEstimatorAvg[3]) delete fMultEstimatorAvg[3];
        fMultEstimatorAvg[3]=new TProfile(*hprof);
    }
    
    void SetMultiplVsZProfileLHC13b(TProfile* hprof){
        if(fMultEstimatorAvg[0]) delete fMultEstimatorAvg[0];
        fMultEstimatorAvg[0]=new TProfile(*hprof);
    }
    void SetMultiplVsZProfileLHC13c(TProfile* hprof){
        if(fMultEstimatorAvg[1]) delete fMultEstimatorAvg[1];
        fMultEstimatorAvg[1]=new TProfile(*hprof);
    }
    
    void SetReferenceMultiplcity(Double_t rmu){fRefMult=rmu;}
    
    // Nch Ntrk weights on MC
    void UseMCNchWeight(Int_t flag) { fUseNchWeight = flag; }
    void SetHistoNchWeight(TH1F *h){
        if(fHistoMCNch) delete fHistoMCNch;
        fHistoMCNch = new TH1F(*h);
    }
    void SetMeasuredNchHisto(TH1F* h){
        if(fHistoMeasNch) delete fHistoMeasNch;
        fHistoMeasNch = new TH1F(*h);
    }
    
    void SetSubtractTrackletsFromDaughters(Bool_t opt){fSubtractTrackletsFromDau=opt;}
    Int_t CheckOrigin(TClonesArray* arrayMC, AliAODMCParticle *mcPartCandidate) const;
    
    // Flag to use the zvtx correction from ( 0= none, 1= usual d2h, 2=AliESDUtils for VZERO multiplicity)
    void SetUseVZEROParameterizedVertexCorr(Int_t flag) { fDoVZER0ParamVertexCorr=flag; }
    // Flag to fill THnSparse with MultUncorr and NoPid cases ( 0 = only Mult, 1 = Mult and multUncorr, 2 = NoPid and 3 is All)
    void SetFillSoSparseForMultUncorrNoPid(Int_t flag) { fFillSoSparseChecks=flag; }
    void SetEventShapeParameters(Double_t ptMin, Double_t ptMax, Double_t etaMin, Double_t etaMax, Int_t minMult, Double_t phiStepSizeDeg, Int_t filtbit1, Int_t filtbit2) { fptMin=ptMin; fptMax=ptMax; fetaMin=etaMin; fetaMax=etaMax; fminMult=minMult; fphiStepSizeDeg=phiStepSizeDeg; ffiltbit1=filtbit1; ffiltbit2=filtbit2;}
    void SetCalculationsForSphericity(Bool_t CalSpheri){fCalculateSphericity=CalSpheri;}
    
    Int_t GetUseVZEROParameterizedVertexCorr() { return fDoVZER0ParamVertexCorr; }
    
    enum { kNtrk10=0, kNtrk10to16=1, kVZERO=2, kNtrk03=3, kNtrk05=4, kVZEROA=5, kVZEROEq=6, kVZEROAEq=7 };
    void SetMultiplicityEstimator(Int_t value){ fMultiplicityEstimator=value; }
    Int_t GetMultiplicityEstimator(){ return fMultiplicityEstimator; }
    enum { kEta10=0, kEta10to16=1, kEtaVZERO=2, kEta03=3, kEta05=5, kEtaVZEROA=5 };
    void SetMCPrimariesEstimator(Int_t value){ fMCPrimariesEstimator=value; }
    Int_t GetMCPrimariesEstimator(){ return fMCPrimariesEstimator; }
    
    // Implementation of interface methods
    virtual void UserCreateOutputObjects();
    virtual void Init();
    virtual void LocalInit() {Init();}
    virtual void UserExec(Option_t *option);
    virtual void Terminate(Option_t *option);
    
private:
    
    AliAnalysisTaskSEDvsEventShapes(const AliAnalysisTaskSEDvsEventShapes &source);
    AliAnalysisTaskSEDvsEventShapes& operator=(const AliAnalysisTaskSEDvsEventShapes& source);
    
    TProfile* GetEstimatorHistogram(const AliVEvent *event);
    void CreateImpactParameterHistos();
    void CreateMeasuredNchHisto();
    void FillMCMassHistos(TClonesArray *arrayMC, Int_t labD, Int_t countMult,Double_t nchWeight);
    
    TList  *fOutput; 
    TList  *fListCuts; //list of cuts
    TList  *fOutputCounters; 
    TList  *fListProfiles; //list of profile histos for z-vtx correction
    
    TH1F *fHistNEvents;     
    
    TH2F* fHistNtrEta16vsNtrEta1EvSel; 
    TH2F* fHistNtrEta05vsNtrEta1EvSel; 
    TH2F* fHistNtrEta03vsNtrEta1EvSel; 
    TH2F* fHistNtrEtaV0AvsNtrEta1EvSel; 
    TH2F* fHistNtrEtaV0MvsNtrEta1EvSel; 
    TH2F* fHistNtrEtaV0AvsV0AEqEvSel;   
    TH2F* fHistNtrEtaV0MvsV0MEqEvSel;   
    TH2F* fHistNtrCorrEta1vsNtrRawEta1EvSel; 
    TH2F* fHistMultCorrvsMultRawEvSel;       
    TH2F* fHistNtrEta16vsNtrEta1EvWithCand; 
    TH2F* fHistNtrEta05vsNtrEta1EvWithCand; 
    TH2F* fHistNtrEta03vsNtrEta1EvWithCand; 
    TH2F* fHistNtrEtaV0AvsNtrEta1EvWithCand; 
    TH2F* fHistNtrEtaV0MvsNtrEta1EvWithCand; 
    TH2F* fHistNtrEtaV0AvsV0AEqEvWithCand;     
    TH2F* fHistNtrEtaV0MvsV0MEqEvWithCand;     
    TH2F* fHistNtrCorrEta1vsNtrRawEta1EvWithCand; 
    TH2F* fHistMultCorrvsMultRawEvWithCand;       
    TH2F* fHistNtrEta16vsNtrEta1EvWithD; 
    TH2F* fHistNtrEta05vsNtrEta1EvWithD; 
    TH2F* fHistNtrEta03vsNtrEta1EvWithD; 
    TH2F* fHistNtrEtaV0AvsNtrEta1EvWithD; 
    TH2F* fHistNtrEtaV0MvsNtrEta1EvWithD; 
    TH2F* fHistNtrEtaV0AvsV0AEqEvWithD;   
    TH2F* fHistNtrEtaV0MvsV0MEqEvWithD;   
    TH2F* fHistNtrCorrEta1vsNtrRawEta1EvWithD; 
    TH2F* fHistMultCorrvsMultRawEvWithD;       
    
    TH2F* fHistNtrVsZvtx; 
    TH2F* fHistNtrCorrVsZvtx; 
    TH2F* fHistNtrVsSo; 
    TH2F* fHistNtrCorrVsSo; 
    
    TH2F* fHistNtrVsNchMC; 
    TH2F* fHistNtrCorrVsNchMC; 
    TH2F* fHistNtrVsNchMCPrimary; 
    TH2F* fHistNtrCorrVsNchMCPrimary; 
    TH2F* fHistNtrVsNchMCPhysicalPrimary; 
    TH2F* fHistNtrCorrVsNchMCPhysicalPrimary; 
    TH1F* fHistGenPrimaryParticlesInelGt0; 
    TH3F* fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary; 
    
    TH1F* fHistNtrUnCorrPSSel; 
    TH1F* fHistNtrUnCorrPSTrigSel; 
    TH1F* fHistNtrUnCorrPSTrigPileUpSel; 
    TH1F* fHistNtrUnCorrPSTrigPileUpVtxSel; 
    TH1F* fHistNtrUnCorrPSTrigPileUpVtxContSel; 
    TH1F* fHistNtrUnCorrPSTrigPileUpVtxRangeSel; 
    TH1F* fHistNtrUnCorrPSTrigPileUpVtxRangeCentrSel; 
    TH1F* fHistNtrUnCorrEvSel; 
    TH1F* fHistNtrUnCorrEvWithCand; 
    TH1F* fHistNtrUnCorrEvWithD;
    TH1F* fHistNtrCorrPSSel; 
    TH1F* fHistNtrCorrEvSel; 
    TH1F* fHistNtrCorrEvWithCand; 
    TH1F* fHistNtrCorrEvWithD;
    TH3F *fPtVsMassVsMultMC;  
    
    THnSparseD *fSparseSpherocity;
    THnSparseD *fSparseSpherocitywithNoPid;
    
    THnSparseF *fHistMassPtImpPar[5];
    
    Double_t fUpmasslimit;  //upper inv mass limit for histos
    Double_t fLowmasslimit; //lower inv mass limit for histos
    Int_t   fNMassBins;    // nbins for invariant mass histos
    
    AliRDHFCuts *fRDCutsAnalysis; // Cuts for Analysis
    AliNormalizationCounter *fCounterC;           
    AliNormalizationCounter *fCounterU;           
    
    AliNormalizationCounter *fCounterCandidates;  
    
    
    Bool_t fDoImpPar;  //swicth for D impact parameter THnSparse
    Int_t  fNImpParBins;   // nunber of bins in impact parameter histos
    Double_t fLowerImpPar;  // lower limit in impact parameter (um)
    Double_t fHigherImpPar; // higher limit in impact parameter (um)
    
    Bool_t fReadMC;    //flag for access to MC
    Int_t  fMCOption;  // 0=keep all cand, 1=keep only signal, 2= keep only back
    Bool_t fisPPbData; // flag to run on pPb data (differen histogram bining)
    Bool_t fUseBit;    // flag to use bitmask
    Bool_t fSubtractTrackletsFromDau; // flag for subtracting D meson daughter contribution to N of tracklets
    Bool_t fKeepCorrPlots; // flag to look at the correlation of different estimators (eta ranges)
    Bool_t fCalculateSphericity;
    Int_t fUseNchWeight; // weight on the MC on the generated multiplicity (0->no weights, 1->Nch weights, 2->Ntrk weights)
    TH1F* fHistoMCNch;    // weight histogram for the MC on the generated multiplicity
    TH1F* fHistoMeasNch;  // weight histogram on the true measured multiplicity
    
    TProfile* fMultEstimatorAvg[4]; //TProfile with mult vs. Z per period
    Double_t fRefMult;   // refrence multiplcity (period b)
    Int_t fPdgMeson;   // pdg code of analyzed meson
    
    Int_t fMultiplicityEstimator; // Definition of the multiplicity estimator: kNtrk10=0, kNtrk10to16=1, kVZERO=2
    Int_t fMCPrimariesEstimator;  // Definition of the primaries estimator eta range: |eta|<1.0=0, -1.6<|eta|<1.0=1, VZEROrange=2
    
    Int_t fDoVZER0ParamVertexCorr; // Flag to use the zvtx correction from (0=none, 1=usual d2h, 2=AliESDUtils for VZERO multiplicity)
    
    Int_t fFillSoSparseChecks; // Flag to fill THnSparse with MultUncorr and NoPid cases ( 0 = only Mult, 1 = Mult and multUncorr, 2 = NoPid and 3 is All)
    Double_t fetaMin;
    Double_t fetaMax;
    Double_t fptMin;
    Double_t fptMax;
    Int_t fminMult;
    Int_t ffiltbit1;
    Int_t ffiltbit2;
    Double_t fphiStepSizeDeg;
    
    ClassDef(AliAnalysisTaskSEDvsEventShapes,1); // event shape task
};

#endif