PWGHF/vertexingHF/AliCFTaskVertexingHF.h

   1 #ifndef ALICFTASKVERTEXINGHF_H
   2 #define ALICFTASKVERTEXINGHF_H
   3 /**************************************************************************
   4  * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
   5  *                                                                        *
   6  * Author: The ALICE Off-line Project.                                    *
   7  * Contributors are mentioned in the code where appropriate.              *
   8  *                                                                        *
   9  * Permission to use, copy, modify and distribute this software and its   *
  10  * documentation strictly for non-commercial purposes is hereby granted   *
  11  * without fee, provided that the above copyright notice appears in all   *
  12  * copies and that both the copyright notice and this permission notice   *
  13  * appear in the supporting documentation. The authors make no claims     *
  14  * about the suitability of this software for any purpose. It is          *
  15  * provided "as is" without express or implied warranty.                  *
  16  **************************************************************************/
  17
  18 /* $Id$ */
  19
  20 //-----------------------------------------------------------------------
  21 // Class for HF corrections as a function of many variables and step
  22 // Author : C. Zampolli, CERN
  23 //                      D. Caffarri, Univ & INFN Padova caffarri@pd.infn.it
  24 // Base class for HF Unfolding - agrelli@uu.nl
  25 //-----------------------------------------------------------------------
  26
  27
  28 #include "AliAnalysisTaskSE.h"
  29 #include "AliCFVertexingHF2Prong.h"
  30 #include "AliCFVertexingHF3Prong.h"
  31 #include "AliCFVertexingHFLctoV0bachelor.h"
  32 #include "AliCFVertexingHF.h"
  33 #include <TH1F.h>
  34 #include <TProfile.h>
  35
  36 class TH1I;
  37 class TParticle ;
  38 class TFile ;
  39 class TClonesArray ;
  40 class AliCFManager;
  41 class AliAODRecoDecay;
  42 class AliAODRecoDecayHF2Prong;
  43 class AliAODMCParticle;
  44 class THnSparse;
  45 class TF1;
  46 class AliRDHFCuts;
  47 class AliCFVertexingHF2Prong;
  48 class AliCFVertexingHF3Prong;
  49
  50 class AliCFTaskVertexingHF: public AliAnalysisTaskSE {
  51 public:
  52
  53         enum {
  54                 kStepGeneratedLimAcc = 0,
  55                 kStepGenerated       = 1,
  56                 kStepAcceptance      = 2,
  57                 kStepVertex          = 3,
  58                 kStepRefit           = 4,
  59                 kStepReconstructed   = 5,
  60                 kStepRecoAcceptance  = 6,
  61                 kStepRecoITSClusters = 7,
  62                 kStepRecoPPR         = 8,
  63                 kStepRecoPID         = 9
  64         };
  65
  66         enum {
  67                 kSnail = 0,    // slow configuration, all variables
  68                 kCheetah = 1   // fast configuration, only a subset of variables
  69         };
  70
  71         enum {
  72           kAll = 0,   // all decays (resonant + non-resonant)
  73           kNonResonant = 1, // only non resonant
  74           kL1520 = 2,  // Lc --> L(1520) + p
  75           kKstar = 3,  // Lc --> K* + pi
  76           kDelta = 4   // Lc --> Delta + K
  77         };
  78
  79         enum { kNtrk10=0, kNtrk10to16=1, kVZERO=2 }; // multiplicity estimators
  80
  81         AliCFTaskVertexingHF();
  82         AliCFTaskVertexingHF(const Char_t* name, AliRDHFCuts* cuts, TF1* func = 0x0);
  83         AliCFTaskVertexingHF& operator= (const AliCFTaskVertexingHF& c);
  84         AliCFTaskVertexingHF(const AliCFTaskVertexingHF& c);
  85         virtual ~AliCFTaskVertexingHF();
  86
  87         // ANALYSIS FRAMEWORK STUFF to loop on data and fill output objects
  88         void     UserCreateOutputObjects();
  89         void     UserExec(Option_t *option);
  90         void     Init();
  91         void     LocalInit() {Init();}
  92         void     Terminate(Option_t *);
  93
  94         // UNFOLDING
  95         void     SetCorrelationMatrix(THnSparse* h) {fCorrelation=h;}
  96         void     SetAcceptanceUnf(Bool_t AcceptanceUnf) {fAcceptanceUnf = AcceptanceUnf;}
  97         Bool_t   GetAcceptanceUnf() const {return fAcceptanceUnf;}
  98
  99
 100         // CORRECTION FRAMEWORK RELATED FUNCTIONS
 101         void           SetCFManager(AliCFManager* io) {fCFManager = io;}   // global correction manager
 102         AliCFManager * GetCFManager()                 {return fCFManager;} // get corr manager
 103
 104         // Setters (and getters) for the config macro
 105         void    SetFillFromGenerated(Bool_t flag) {fFillFromGenerated = flag;}
 106         Bool_t  GetFillFromGenerated() const {return fFillFromGenerated;}
 107         void    SetDecayChannel (Int_t decayChannel) {fDecayChannel = decayChannel;}
 108         Int_t   GetDecayChannel () {return fDecayChannel;}
 109         void     SetUseWeight(Bool_t useWeight){fUseWeight=useWeight;}
 110         Bool_t   GetUseWeight() const {return fUseWeight;}
 111         Double_t GetWeight(Float_t pt);
 112         Double_t dNdptFit(Float_t pt, Double_t* par);
 113         Double_t GetPtWeightFromHistogram(Float_t pt);
 114
 115         void SetUseFlatPtWeight(Bool_t useWeight){fUseFlatPtWeight=useWeight; fUseWeight=useWeight;}
 116         Bool_t GetUseFlatPtWeight() const {return fUseFlatPtWeight;}
 117         void SetUseZWeight(Bool_t useWeight){fUseZWeight=useWeight;}
 118         Bool_t GetUseZWeight() const {return fUseZWeight;}
 119         Double_t GetZWeight(Float_t z, Int_t runnumber);
 120         Double_t DodzFit(Float_t z, Double_t* par);
 121
 122         void SetUseNchWeight(Bool_t useWeight){fUseNchWeight=useWeight;}
 123         Bool_t GetUseNchWeight() const {return fUseNchWeight;}
 124         void SetMCNchHisto(TH1F* h){
 125           if(fHistoMCNch) delete fHistoMCNch;
 126           fHistoMCNch=new TH1F(*h);
 127         }
 128         void CreateMeasuredNchHisto();
 129         void SetMeasuredNchHisto(TH1F* h){
 130           if(fHistoMeasNch) delete fHistoMeasNch;
 131           fHistoMeasNch=new TH1F(*h);
 132         }
 133         Double_t GetNchWeight(Int_t nch);
 134         void SetMultiplicityEstimator(Int_t value){ fMultiplicityEstimator=value; }
 135         Int_t GetMultiplicityEstimator(){ return fMultiplicityEstimator; }
 136         void SetIsPPData(Bool_t flag){ fIsPPData = flag; }
 137         void SetIsPPbData(Bool_t flag){ fIsPPbData = flag; }
 138
 139         void SetUseNchTrackletsWeight(Bool_t useWeight = kTRUE) { fUseNchWeight=useWeight; fUseTrackletsWeight=useWeight; fUseMultRatioAsWeight=useWeight; }
 140         Bool_t GetUseNchTrackletsWeight() const {return fUseTrackletsWeight;}
 141         void SetUseRatioMultiplicityDistributionsAsWeight(Bool_t flag=kTRUE){ fUseMultRatioAsWeight=flag; }
 142         Bool_t GetUseRatioMultiplicityDistributionsAsWeight() const {return fUseMultRatioAsWeight;}
 143
 144         void SetUseZvtxCorrectedNtrkEstimator(Bool_t flag) { fZvtxCorrectedNtrkEstimator=flag; }
 145         Bool_t GetUseZvtxCorrectedNtrkEstimator() { return fZvtxCorrectedNtrkEstimator; }
 146         void SetMultiplVsZProfileLHC10b(TProfile* hprof){
 147           if(fMultEstimatorAvg[0]) delete fMultEstimatorAvg[0];
 148           fMultEstimatorAvg[0]=new TProfile(*hprof);
 149         }
 150         void SetMultiplVsZProfileLHC10c(TProfile* hprof){
 151           if(fMultEstimatorAvg[1]) delete fMultEstimatorAvg[1];
 152           fMultEstimatorAvg[1]=new TProfile(*hprof);
 153         }
 154         void SetMultiplVsZProfileLHC10d(TProfile* hprof){
 155           if(fMultEstimatorAvg[2]) delete fMultEstimatorAvg[2];
 156           fMultEstimatorAvg[2]=new TProfile(*hprof);
 157         }
 158         void SetMultiplVsZProfileLHC10e(TProfile* hprof){
 159           if(fMultEstimatorAvg[3]) delete fMultEstimatorAvg[3];
 160           fMultEstimatorAvg[3]=new TProfile(*hprof);
 161         }
 162
 163         void SetMultiplVsZProfileLHC13b(TProfile* hprof){
 164           if(fMultEstimatorAvg[0]) delete fMultEstimatorAvg[0];
 165           fMultEstimatorAvg[0]=new TProfile(*hprof);
 166         }
 167         void SetMultiplVsZProfileLHC13c(TProfile* hprof){
 168           if(fMultEstimatorAvg[1]) delete fMultEstimatorAvg[1];
 169           fMultEstimatorAvg[1]=new TProfile(*hprof);
 170         }
 171
 172
 173         TProfile* GetEstimatorHistogram(const AliVEvent* event);
 174         void SetReferenceMultiplcity(Double_t rmu){fRefMult=rmu;}
 175
 176         void   SetDselection(UShort_t originDselection) {fOriginDselection=originDselection;}
 177         UShort_t GetDselection (){return fOriginDselection;}
 178         void SetSign(Char_t isSign) {fSign = isSign;}
 179         Char_t GetSign() {return fSign;}
 180
 181         void SetCentralitySelection(Bool_t centSelec = kTRUE) {fCentralitySelection = centSelec;}
 182         Bool_t GetCentralitySelection() {return fCentralitySelection;}
 183
 184         void SetFakeSelection(Int_t fakeSel = 0) {fFakeSelection=fakeSel;}
 185         Int_t GetFakeSelection(){return fFakeSelection;}
 186
 187         void SetRejectCandidateIfNotFromQuark(Bool_t opt){fRejectIfNoQuark=opt;}
 188         Bool_t GetRejectCandidateIfNotFromQuark(){return fRejectIfNoQuark;}
 189
 190         void SetUseMCVertex(Bool_t opt){fUseMCVertex=opt;}
 191         Bool_t GetUseMCVertex(){return fUseMCVertex;}
 192
 193
 194         void SetKeepDsViaPhi(){fDsOption=1;}
 195         void SetKeepDsViaK0star(){fDsOption=2;}
 196         void SetKeepAllDs(){fDsOption=3;}
 197         void SetCountAllDs(){fGenDsOption=AliCFVertexingHF3Prong::kCountAllDsKKpi;}
 198         void SetCountDsViaPhi(){fGenDsOption=AliCFVertexingHF3Prong::kCountPhipi;}
 199         void SetCountDsViaK0star(){fGenDsOption=AliCFVertexingHF3Prong::kCountK0stK;}
 200         void SetCountResonantDs(){fGenDsOption=AliCFVertexingHF3Prong::kCountResonant;}
 201         void SetCountNonResonantDs(){fGenDsOption=AliCFVertexingHF3Prong::kCountNonResonant;}
 202
 203         Bool_t ProcessDs(Int_t returnCodeDs) const;
 204
 205         void SetConfiguration(Int_t configuration) {(configuration == kSnail) ? Printf("Slow configuration chosen, all variables will be used!") : Printf("Fast configuration chosen, all variablesOnly pt, y, phi, ct, fake, z_vtx, centrality and multiplicity will be used!"); fConfiguration = configuration;}
 206         Int_t GetConfiguration() const {return fConfiguration;}
 207
 208         void SetWeightFunction(TF1* func) {fFuncWeight = func;}
 209         TF1* GetWeightFunction() const {return fFuncWeight;}
 210         void SetWeightHistogram(TH1F* histo) {
 211           if(fHistoPtWeight) delete fHistoPtWeight;
 212           fHistoPtWeight=new TH1F(*histo);
 213         }
 214         TH1F* GetWeightHistogram() const {return (TH1F*)fHistoPtWeight;}
 215
 216         void SetPtWeightsFromFONLL276overLHC12a17a();
 217         void SetPtWeightsFromDataPbPb276overLHC12a17a();
 218         void SetPtWeightsFromFONLL276overLHC12a17b();
 219         void SetPtWeightsFromFONLL276andBAMPSoverLHC12a17b();
 220         void SetPtWeightsFromFONLL276overLHC10f6a();
 221         void SetPtWeightsFromFONLL7overLHC10f6a();
 222         void SetPtWeightsFromFONLL7overLHC12a12();
 223         void SetPtWeightsFromFONLL7overLHC12a12bis();
 224         void SetPtWeightsFromFONLL7overLHC13e2fix();
 225         void SetPtWeightsFromFONLL5overLHC10f6a();
 226         void SetPtWeightsFromFONLL5overLHC13d3();
 227
 228         void SetResonantDecay(UInt_t resonantDecay) {fResonantDecay = resonantDecay;}
 229         UInt_t GetResonantDecay() const {return fResonantDecay;}
 230
 231         void SetKeepLctoK0Sp() {fLctoV0bachelorOption=1;}
 232         void SetKeepLctoLambdaBarpi() {fLctoV0bachelorOption=2;}
 233         void SetKeepLctoLambdapi() {fLctoV0bachelorOption=4;}
 234         void SetKeepLctoV0bachelor() {fLctoV0bachelorOption=7;}
 235
 236         void SetCountLctoK0Sp(){fGenLctoV0bachelorOption=AliCFVertexingHFLctoV0bachelor::kCountK0Sp;}
 237         void SetCountLctoLambdapi(){fGenLctoV0bachelorOption=AliCFVertexingHFLctoV0bachelor::kCountLambdapi;}
 238
 239         void SetUseSelectionBit(Bool_t flag) { fUseSelectionBit=flag; }
 240         Bool_t GetUseSelectionBit() const { return fUseSelectionBit; }
 241
 242         Bool_t ProcessLctoV0Bachelor(Int_t returnCodeDs) const;
 243
 244         void SetUseAdditionalCuts(Bool_t flag) { fDecayChannel == 22 ? fUseAdditionalCuts = flag : fUseAdditionalCuts = kFALSE;}
 245         Bool_t GetUseAdditionalCuts() const { return fUseAdditionalCuts; }
 246         void SetUseCutsForTMVA(Bool_t useCutsForTMVA) { fDecayChannel == 22 ? fUseCutsForTMVA = useCutsForTMVA : fUseAdditionalCuts = kFALSE;}
 247         Bool_t GetUseCutsForTMVA() const {return fUseCutsForTMVA;}
 248
 249         void SetUseCascadeTaskForLctoV0bachelor(Bool_t useCascadeTaskForLctoV0bachelor) {fUseCascadeTaskForLctoV0bachelor = useCascadeTaskForLctoV0bachelor;}
 250         Bool_t GetUseCascadeTaskForLctoV0bachelor() const {return fUseCascadeTaskForLctoV0bachelor;}
 251
 252         void SetCutOnMomConservation(Float_t cut) {fCutOnMomConservation = cut;}
 253         Bool_t GetCutOnMomConservation() const {return fCutOnMomConservation;}
 254
 255 protected:
 256         AliCFManager   *fCFManager;   //  pointer to the CF manager
 257         TH1I *fHistEventsProcessed;   //! simple histo for monitoring the number of events processed
 258         THnSparse* fCorrelation;      //  response matrix for unfolding
 259         TList  *fListProfiles; //list of profile histos for z-vtx correction
 260         Int_t fCountMC;               //  MC particle found
 261         Int_t fCountAcc;              //  MC particle found that satisfy acceptance cuts
 262         Int_t fCountVertex;       //  Reco particle found that satisfy vertex constrained
 263         Int_t fCountRefit;        //  Reco particle found that satisfy kTPCrefit and kITSrefit
 264         Int_t fCountReco;             //  Reco particle found that satisfy cuts
 265         Int_t fCountRecoAcc;          //  Reco particle found that satisfy cuts in requested acceptance
 266         Int_t fCountRecoITSClusters;  //  Reco particle found that satisfy cuts in n. of ITS clusters
 267         Int_t fCountRecoPPR;          //  Reco particle found that satisfy cuts in PPR
 268         Int_t fCountRecoPID;          //Reco PID step
 269         Int_t fEvents;                //  n. of events
 270         Int_t fDecayChannel;          // decay channel to configure the task
 271         Bool_t fFillFromGenerated;    //  flag to indicate whether data container should be filled with generated values also for reconstructed particles
 272         UShort_t fOriginDselection;      // flag to select D0 origins. 0 Only from charm 1 only from beauty 2 both from charm and beauty
 273         Bool_t fAcceptanceUnf;        //  flag for unfolding before or after cuts.
 274         AliRDHFCuts* fCuts;            // cuts
 275         Bool_t fUseWeight;             //flag to decide whether to use pt-weights != 1 when filling the container or not
 276         Double_t fWeight;              //weight used to fill the container
 277         Bool_t fUseFlatPtWeight;       // flag to decide to use a flat pt shape
 278         Bool_t fUseZWeight;           // flag to decide whether to use z-vtx weights != 1 when filling the container or not
 279         Bool_t fUseNchWeight;         // flag to decide whether to use Ncharged weights != 1 when filling the container or not
 280         Bool_t fUseTrackletsWeight;   // flag to decide whether to use Ncharged weights != 1 when filling the container or not
 281         Bool_t fUseMultRatioAsWeight; // flag to use directly the ratio of the distributions (fHistoMCNch) instead of computing it
 282         Int_t fNvar;                  // number of variables for the container
 283         TString fPartName;    // D meson name
 284         TString fDauNames;    // daughter in fin state
 285         Char_t fSign;                 // flag to decide wheter to keep D0 only (0), D0bar only (1), or both D0 and D0bar (2)
 286         Bool_t fCentralitySelection;  //flag to switch off the centrality selection
 287         Int_t  fFakeSelection;  //selection flag for fakes tracks
 288         Bool_t fRejectIfNoQuark;  // flag to remove events not geenrated with PYTHIA
 289         Bool_t fUseMCVertex;  // flag to use MC vertex (useful when runnign in pp)
 290         Int_t  fDsOption;     // Ds decay option (selection level)
 291         Int_t  fGenDsOption;     // Ds decay option (generation level)
 292         Int_t fConfiguration; // configuration (slow / fast) of the CF --> different variables will be allocated (all / reduced number)
 293         TF1* fFuncWeight;     // user-defined function to be used to calculate weights
 294         TH1F* fHistoPtWeight; // user-defined histogram to calculate the Pt weights
 295         TH1F* fHistoMeasNch;  // histogram with measured Nch distribution (pp 7 TeV)
 296         TH1F* fHistoMCNch;  // histogram with Nch distribution from MC production
 297         UInt_t fResonantDecay;  // resonant deacy channel to be used if the CF should be run on resonant channels only
 298         Int_t fLctoV0bachelorOption; // Lc->V0+bachelor decay option (selection level)
 299         Int_t fGenLctoV0bachelorOption; // Lc->V0+bachelor decay option (generation level)
 300         Bool_t fUseSelectionBit;     // flag to use selection bit
 301         UInt_t fPDGcode; // PDG code
 302
 303         Int_t fMultiplicityEstimator; // Definition of the multiplicity estimator: kNtrk10=0, kNtrk10to16=1, kVZERO=2
 304         TProfile* fMultEstimatorAvg[4]; // TProfile with mult vas. Z per period
 305         Double_t fRefMult;   // refrence multiplcity (period b)
 306         Bool_t fZvtxCorrectedNtrkEstimator; // flag to use the z-vtx corrected (if not use uncorrected) multiplicity estimator
 307         Bool_t fIsPPData; // flag for pp data (not checking centrality)
 308         Bool_t fIsPPbData; // flag for pPb data (used for multiplicity corrections)
 309         Bool_t fUseAdditionalCuts;  // flag to use additional cuts needed for Lc --> K0S + p, TMVA
 310         Bool_t fUseCutsForTMVA;     // flag to use additional cuts needed for Lc --> K0S + p, TMVA
 311                                     // these are the pre-selection cuts for the TMVA
 312         Bool_t fUseCascadeTaskForLctoV0bachelor;   // flag to define which task to use for Lc --> K0S+p
 313         Float_t fCutOnMomConservation; // cut on momentum conservation
 314
 315         ClassDef(AliCFTaskVertexingHF,24); // class for HF corrections as a function of many variables
 316 };
 317
 318 #endif