AliHFSystErr.h

Go to the documentation of this file.

#ifndef ALIHFSYSTERR_H
#define ALIHFSYSTERR_H
/* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */ 

//***********************************************************
//***********************************************************

#include <TNamed.h>
#include <TString.h>
#include <TH1F.h>
#include "AliLog.h"
#include "TGraphAsymmErrors.h"


class AliHFSystErr : public TNamed 
{
 public:

  AliHFSystErr(const Char_t* name="HFSystErr", const Char_t* title="");
    
  virtual ~AliHFSystErr();
  
  void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const; 

  Double_t GetNormErr() const {return (fNorm ? fNorm->GetBinContent(0) : 0.);}
  Double_t GetBRErr() const {return (fBR ? fBR->GetBinContent(0) : 0.);}
  Double_t GetCutsEffErr(Double_t pt) const;
  Double_t GetMCPtShapeErr(Double_t pt) const;
  Double_t GetSeleEffErr(Double_t pt) const;
  Double_t GetPartAntipartErr(Double_t pt) const;
  Double_t GetPIDEffErr(Double_t pt) const;
  Double_t GetRawYieldErr(Double_t pt) const;
  Double_t GetTrackingEffErr(Double_t pt) const;
  Double_t GetTotalSystErr(Double_t pt,Double_t feeddownErr=0) const;

  void SetRunNumber(Int_t number) { 
    fRunNumber = number; 
    AliInfo(Form(" Settings for run year 20%2d",fRunNumber));
  }
  Int_t GetRunNumber() const { return fRunNumber; }
  void SetCollisionType(Int_t type) { 
    fCollisionType = type; 
    if (fCollisionType==0) { AliInfo(" Settings for p-p collisions"); }
    else if(fCollisionType==1) { AliInfo(" Settings for Pb-Pb collisions"); }
    else if(fCollisionType==2) { AliInfo(" Settings for p-Pb collisions"); }
  }
  Int_t GetCollisionType() const { return fCollisionType; }
  void SetCentrality(TString centrality) { 
    fCentralityClass = centrality; 
    AliInfo(Form(" Settings for centrality class %s",fCentralityClass.Data()));
  }
  void SetIsLowEnergy(Bool_t flag) { 
    fIsLowEnergy = flag; 
    if(flag) AliInfo(" Settings for the low energy run");
  }
  void SetIsLowPtAnalysis(Bool_t flag){
    fIsLowPtAnalysis = flag;
    if(flag) AliInfo("Settings for the low pt analysis");  
  }
  void SetIsPass4Analysis(Bool_t flag){
    fIsPass4Analysis = flag;
    if(flag) AliInfo("Settings for the pass4 analysis");  
  }
  void SetIs5TeVAnalysis(Bool_t flag){
    fIs5TeVAnalysis = flag;
    if(flag) AliInfo("Settings for the 5TeV analysis");  
  }
  void SetIsBDTAnalysis(Bool_t flag){
    fIsBDTAnalysis = flag;
    if(flag) AliInfo("Settings for the Lc BDT analysis");  
  }


  void SetIsPbPb2010EnergyScan(Bool_t flag) {
    fIsCentScan = flag;
    if(flag) AliInfo(" Settings for the PbPb 2010 energy scan");
  }

  void SetRapidity(TString rapidity) {
    fRapidityRange = rapidity;
    AliInfo(Form(" Settings for rapidity interval %s",fRapidityRange.Data()));
  }
  void SetIspPb2011RapidityScan(Bool_t flag){
    fIsRapidityScan = flag; 
    if(flag) AliInfo("Settings for the pPb vs y measurement");
  }

  void Init(Int_t decay);

  void InitD0toKpi2010PbPb010CentScan();
  void InitD0toKpi2010PbPb1020CentScan();
  void InitD0toKpi2010PbPb2040CentScan();
  void InitD0toKpi2010PbPb4060CentScan();
  void InitD0toKpi2010PbPb6080CentScan();

  void InitD0toKpi2011PbPb3050InPlane();
  void InitD0toKpi2011PbPb3050OutOfPlane();

  void InitDplustoKpipi2010PbPb010CentScan();
  void InitDplustoKpipi2010PbPb1020CentScan();
  void InitDplustoKpipi2010PbPb2040CentScan();
  void InitDplustoKpipi2010PbPb4060CentScan();
  void InitDplustoKpipi2010PbPb6080CentScan();

  void InitDstartoD0pi2010PbPb010CentScan();
  void InitDstartoD0pi2010PbPb1020CentScan();
  void InitDstartoD0pi2010PbPb2040CentScan();
  void InitDstartoD0pi2010PbPb4060CentScan();
  void InitDstartoD0pi2010PbPb6080CentScan();

  void InitD0toKpi2011PbPb010CentScan();
  void InitD0toKpi2011PbPb1020CentScan();
  void InitD0toKpi2011PbPb2030CentScan();
  void InitD0toKpi2011PbPb3040CentScan();
  void InitD0toKpi2011PbPb4050CentScan();
  void InitD0toKpi2010PbPb5080CentScan();

  void InitDplustoKpipi2011PbPb010CentScan();
  void InitDplustoKpipi2011PbPb1020CentScan();
  void InitDplustoKpipi2011PbPb2030CentScan();
  void InitDplustoKpipi2011PbPb3040CentScan();
  void InitDplustoKpipi2011PbPb4050CentScan();
  void InitDplustoKpipi2010PbPb5080CentScan();

  void InitDstartoD0pi2011PbPb010CentScan();
  void InitDstartoD0pi2011PbPb1020CentScan();
  void InitDstartoD0pi2011PbPb2030CentScan();
  void InitDstartoD0pi2011PbPb3040CentScan();
  void InitDstartoD0pi2011PbPb4050CentScan();
  void InitDstartoD0pi2010PbPb5080CentScan();

  void InitD0toKpi2013pPb0100RapScan0804();
  void InitD0toKpi2013pPb0100RapScan0401();
  void InitD0toKpi2013pPb0100RapScan0101();
  void InitD0toKpi2013pPb0100RapScan0104();
  void InitD0toKpi2013pPb0100RapScan0408();

  void InitDplustoKpipi2013pPb0100RapScan0804();
  void InitDplustoKpipi2013pPb0100RapScan0401();
  void InitDplustoKpipi2013pPb0100RapScan0101();
  void InitDplustoKpipi2013pPb0100RapScan0104();
  void InitDplustoKpipi2013pPb0100RapScan0408();

  void InitDstartoD0pi2013pPb0100RapScan0804();
  void InitDstartoD0pi2013pPb0100RapScan0401();
  void InitDstartoD0pi2013pPb0100RapScan0101();
  void InitDstartoD0pi2013pPb0100RapScan0104();
  void InitDstartoD0pi2013pPb0100RapScan0408();
    
  void InitD0toKpi2013pPb020V0A();
  void InitD0toKpi2013pPb2040V0A();
  void InitD0toKpi2013pPb4060V0A();
  void InitD0toKpi2013pPb60100V0A();
  
  void InitD0toKpi2013pPb020ZNA();
  void InitD0toKpi2013pPb2040ZNA();
  void InitD0toKpi2013pPb4060ZNA();
  void InitD0toKpi2013pPb60100ZNA();

  void InitD0toKpi2016pPb010ZNA();
  void InitD0toKpi2016pPb60100ZNA();
  
  void InitD0toKpi2013pPb020CL1();
  void InitD0toKpi2013pPb2040CL1();
  void InitD0toKpi2013pPb4060CL1();
  void InitD0toKpi2013pPb60100CL1();
  
  void InitDstartoD0pi2013pPb020V0A();
  void InitDstartoD0pi2013pPb2040V0A();
  void InitDstartoD0pi2013pPb4060V0A();
  void InitDstartoD0pi2013pPb60100V0A();
  
  void InitDstartoD0pi2013pPb020ZNA();
  void InitDstartoD0pi2013pPb2040ZNA();
  void InitDstartoD0pi2013pPb4060ZNA();
  void InitDstartoD0pi2013pPb60100ZNA();

   void InitDstartoD0pi2016pPb010ZNA();
   void InitDstartoD0pi2016pPb60100ZNA();
 
  void InitDstartoD0pi2013pPb020CL1();
  void InitDstartoD0pi2013pPb2040CL1();
  void InitDstartoD0pi2013pPb4060CL1();
  void InitDstartoD0pi2013pPb60100CL1();

  void InitDplustoKpipi2013pPb020V0A();
  void InitDplustoKpipi2013pPb2040V0A();
  void InitDplustoKpipi2013pPb4060V0A();
  void InitDplustoKpipi2013pPb60100V0A();
  
  void InitDplustoKpipi2013pPb020ZNA();
  void InitDplustoKpipi2013pPb2040ZNA();
  void InitDplustoKpipi2013pPb4060ZNA();
  void InitDplustoKpipi2013pPb60100ZNA();
  
  void InitDplustoKpipi2013pPb020CL1();
  void InitDplustoKpipi2013pPb2040CL1();
  void InitDplustoKpipi2013pPb4060CL1();
  void InitDplustoKpipi2013pPb60100CL1();
    
  void InitDplustoKpipi2016pPb140trkl();
  void InitDplustoKpipi2016pPb4070trkl();
  void InitDplustoKpipi2016pPb70200trkl();
  void InitDplustoKpipi2016pPb010ZNA();
  void InitDplustoKpipi2016pPb60100ZNA();

 private:

  AliHFSystErr(const AliHFSystErr& source);
  AliHFSystErr& operator=(const AliHFSystErr& source); 
 
  void InitD0toKpi2010pp();
  void InitD0toKpi2010ppLowEn();
  void InitD0toKpi2010ppLowPtAn();
  void InitD0toKpi2010ppPass4();
  void InitD0toKpi2015pp5TeV();
  void InitD0toKpi2011PbPb07half();
  void InitD0toKpi2010PbPb020();
  void InitD0toKpi2010PbPb4080();
  void InitD0toKpi2011PbPb3050(); 
  void InitD0toKpi2011PbPb010();
  void InitD0toKpi2013pPb0100();
  void InitD0toKpi2016pPb0100();
  void InitD0toKpi2013pPb0100LowPtAn();
  void InitD0toKpi2016pPb0100LowPtAn();

  void InitDplustoKpipi2010pp();
  void InitDplustoKpipi2010ppPass4();
  void InitDplustoKpipi2010ppLowEn();
  void InitDplustoKpipi2012pp();
  void InitDplustoKpipi2015pp5TeV();
  void InitDplustoKpipi2016pp13TeV();
  void InitDplustoKpipi2011PbPb07half();
  void InitDplustoKpipi2010PbPb020();
  void InitDplustoKpipi2010PbPb4080();
  void InitDplustoKpipi2011PbPb3050();
  void InitDplustoKpipi2011PbPb010();
  void InitDplustoKpipi2013pPb0100();
  void InitDplustoKpipi2016pPb0100();

  void InitDstartoD0pi2010pp();
  void InitDstartoD0pi2010ppLowEn();
  void InitDstartoD0pi2012pp();
  void InitDstartoD0pi2011PbPb07half();
  void InitDstartoD0pi2010PbPb020();
  void InitDstartoD0pi2010PbPb2040();
  void InitDstartoD0pi2010PbPb4080();
  void InitDstartoD0pi2011PbPb3050();
  void InitDstartoD0pi2011PbPb010();
  void InitDstartoD0pi2013pPb0100();
  void InitDstartoD0pi2016pPb0100();
  void InitDstartoD0pi2010ppPass4();
 
  void InitDstoKKpi2010pp();
  void InitDstoKKpi2010ppPass4();
  void InitDstoKKpi2011PbPb07half();
  void InitDstoKKpi2011PbPb010();
  void InitDstoKKpi2011PbPb2050();
  void InitDstoKKpi2013pPb0100();
  void InitDstoKKpi2016pPb0100();
  void InitDstoKKpi2016pPb140trkl();
  void InitDstoKKpi2016pPb4070trkl();
  void InitDstoKKpi2016pPb70200trkl();
  void InitDstoKKpi2016pp13TeV();

  void InitLctopKpi2010pp();
  void InitLctopKpi2010ppBDT();
  void InitLctopKpi2013pPb();
  void InitLctopKpi2013pPbBDT();

  void InitLctopK0S2010pp();
  void InitLctopK0S2013pPb();
  void InitLctopK0S2013pPbBDT();

  void InitD0toKpi2015PbPb010();
  void InitD0toKpi2015PbPb3050();
  void InitD0toKpi2015PbPb6080();
  
  void InitDplustoKpipi2015PbPb010();
  void InitDplustoKpipi2015PbPb3050();
  void InitDplustoKpipi2015PbPb6080();

  void InitDstoKKpi2015PbPb010();
  void InitDstoKKpi2015PbPb3050();
  void InitDstoKKpi2015PbPb6080();

  void InitDstartoD0pi2015PbPb010();
  void InitDstartoD0pi2015PbPb3050();
  void InitDstartoD0pi2015PbPb6080();


  TH1F* ReflectHisto(TH1F *hin) const;

  TH1F *fNorm;            
  TH1F *fRawYield;        
  TH1F *fTrackingEff;     
  TH1F *fBR;              
  TH1F *fCutsEff;         
  TH1F *fPIDEff;          
  TH1F *fMCPtShape;       
  TH1F *fPartAntipart;    

  Int_t fRunNumber;        
  Int_t fCollisionType;    
  TString fCentralityClass;  
  TString fRapidityRange;  

  Bool_t fIsLowEnergy;     
  Bool_t fIsLowPtAnalysis; 
  Bool_t fIsPass4Analysis; 
  Bool_t fIs5TeVAnalysis; 
  Bool_t fIsBDTAnalysis;   
  Bool_t fIsCentScan;      
  Bool_t fIsRapidityScan;  

  ClassDef(AliHFSystErr,9);  
};

#endif