AliHFSystErr.cxx

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

//
// Class to handle systematic errors for charm hadrons
//
// Usage:
// AliHFSystEff syst;           // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
// syst.SetRunNumber(YEAR);     // YEAR = two last numbers of the year (is 10 for 2010)
// syst.SetCollisionType(TYPE);  // TYPE =  0 is pp, 1 is PbPb
// syst.SetCentrality(CENT);     // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
// syst.Init(DECAY);             // DECAY = 1 for D0, 2, for D+, 3 for D* 
// syst.DrawErrors(); // to see a plot of the error contributions
// syst.GetTotalSystErr(pt); // to get the total err at pt 
//
// Author: A.Dainese, andrea.dainese@pd.infn.it

#include <TStyle.h>
#include <TGraphAsymmErrors.h>
#include <TMath.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <TLegend.h>
#include <TColor.h>

#include "AliLog.h"
#include "AliHFSystErr.h"

ClassImp(AliHFSystErr);

//--------------------------------------------------------------------------
AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) : 
  TNamed(name,title),
  fNorm(0),
  fRawYield(0),
  fTrackingEff(0),
  fBR(0),
  fCutsEff(0),
  fPIDEff(0),
  fMCPtShape(0),
  fPartAntipart(0),
  fRunNumber(10),
  fCollisionType(0),
  fCentralityClass("0100"),
  fRapidityRange("0101"),
  fIsLowEnergy(false),
  fIsLowPtAnalysis(false),
  fIsPass4Analysis(false),
  fIs5TeVAnalysis(false),
  fIsBDTAnalysis(false),
  fIsCentScan(false),
  fIsRapidityScan(false)
{
  //
  //
}

//--------------------------------------------------------------------------
AliHFSystErr::~AliHFSystErr() {
  //  
  //
  /*

    if(fNorm)         { delete fNorm; fNorm=0; }
    if(fRawYield)     { delete fRawYield; fRawYield=0; }
    if(fTrackingEff)  { delete fTrackingEff; fTrackingEff=0; }
    if(fBR)           { delete fBR; fBR=0; }
    if(fCutsEff)      { delete fCutsEff; fCutsEff=0; }
    if(fPIDEff)       { delete fPIDEff; fPIDEff=0; }
    if(fMCPtShape)    { delete fMCPtShape; fMCPtShape=0; }
    if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
  */
}

//--------------------------------------------------------------------------
void AliHFSystErr::Init(Int_t decay){
  //
  //

  //  if ((fRunNumber>11) && fIsLowEnergy==false) {
  //    AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
  //  }

  switch(decay) {
  case 1: // D0->Kpi
    if (fCollisionType==0) {
      if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
      else if(fRunNumber == 10){
  if(fIsLowPtAnalysis) InitD0toKpi2010ppLowPtAn();
  else if(fIsPass4Analysis) InitD0toKpi2010ppPass4();
  else if(fIs5TeVAnalysis) InitD0toKpi2015pp5TeV();
  else InitD0toKpi2010pp();
      } else AliFatal("Not yet implemented");
    } 
    else if (fCollisionType==1) {
      if (fRunNumber == 10){
  if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
  else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
  else if (fCentralityClass=="020")  InitD0toKpi2010PbPb020();
  else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
  else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
  else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
  else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
  else AliFatal("Not yet implemented");
      }
      else if (fRunNumber == 11){
  if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
  else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
  else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
  else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
  else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
  else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
  else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
  else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
  else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
  else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
  else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
  else AliFatal("Not yet implemented");
      }
      else if (fRunNumber == 15){
  if (fCentralityClass=="010") InitD0toKpi2015PbPb010();
  else if (fCentralityClass=="3050") InitD0toKpi2015PbPb3050();
  else if (fCentralityClass=="6080") InitD0toKpi2015PbPb6080();
  else AliFatal("Not yet implemented");
      }
    } 
    else if (fCollisionType==2) { 
      if (fCentralityClass=="0100"){
  if(fIsLowPtAnalysis){
    if(fRunNumber==16 || fRunNumber==2016) InitD0toKpi2016pPb0100LowPtAn();
    else InitD0toKpi2013pPb0100LowPtAn();
  }else{
    if(fRunNumber==16 || fRunNumber==2016) InitD0toKpi2016pPb0100();
    else InitD0toKpi2013pPb0100();
  }
      }
      if(fRunNumber==16 || fRunNumber==2016){
      if (fCentralityClass=="010ZNA") InitD0toKpi2016pPb010ZNA();
      else if(fCentralityClass=="60100ZNA") InitD0toKpi2016pPb60100ZNA();
      }else{
      if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
      if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
      if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
      if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
        
      if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
      if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
      if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
      if (fCentralityClass=="60100ZNA")InitD0toKpi2013pPb60100ZNA();
         
      if (fCentralityClass=="020CL1") InitD0toKpi2013pPb020CL1();
      if (fCentralityClass=="2040CL1") InitD0toKpi2013pPb2040CL1();
      if (fCentralityClass=="4060CL1") InitD0toKpi2013pPb4060CL1();
      if (fCentralityClass=="60100CL1") InitD0toKpi2013pPb60100CL1();
      
      if (fIsRapidityScan) {
  if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
  if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
  if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
  if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
  if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
      }
    }
   } 
    else AliFatal("Not yet implemented");
    break;
    
  case 2: // D+->Kpipi
    if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
    if (fCollisionType==0) {
      if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
      else if(fIs5TeVAnalysis) InitDplustoKpipi2015pp5TeV();
      else if(fRunNumber == 10){
  if(fIsPass4Analysis) InitDplustoKpipi2010ppPass4();
  else InitDplustoKpipi2010pp();
      } else if(fRunNumber == 12){
  InitDplustoKpipi2012pp();
      } else if(fRunNumber == 16){
  InitDplustoKpipi2016pp13TeV();
      } else AliFatal("Not yet implemented");
    }
    else if (fCollisionType==1) {
      if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
      if (fRunNumber == 10){
  if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
  else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
  else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
  else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
  else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
  else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
  else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
  else AliFatal("Not yet implemented");
      }
      if(fRunNumber == 11){
  if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
  else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
  else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
  else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
  else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
  else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
  else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
  else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
  else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
  else AliFatal("Not yet implemented");
      }
      if(fRunNumber == 15){
  if (fCentralityClass=="010") InitDplustoKpipi2015PbPb010();
  else if (fCentralityClass=="3050") InitDplustoKpipi2015PbPb3050();
  else if (fCentralityClass=="6080") InitDplustoKpipi2015PbPb6080();
  else AliFatal("Not yet implemented");
      }
    } 
    else if (fCollisionType==2) { 
        if(fRunNumber==16 || fRunNumber==2016) {
          if (fCentralityClass=="0100") InitDplustoKpipi2016pPb0100();
          
          if (fCentralityClass=="010ZNA") InitDplustoKpipi2016pPb010ZNA();
          if (fCentralityClass=="60100ZNA") InitDplustoKpipi2016pPb60100ZNA();
          
          if (fCentralityClass=="140trkl") InitDplustoKpipi2016pPb140trkl();
          if (fCentralityClass=="4070trkl") InitDplustoKpipi2016pPb4070trkl();
          if (fCentralityClass=="70200trkl") InitDplustoKpipi2016pPb70200trkl();
        }
        else {
      if (fCentralityClass=="0100"){InitDplustoKpipi2013pPb0100();}

      if (fCentralityClass=="020V0A") InitDplustoKpipi2013pPb020V0A();
      if (fCentralityClass=="2040V0A") InitDplustoKpipi2013pPb2040V0A();
      if (fCentralityClass=="4060V0A") InitDplustoKpipi2013pPb4060V0A();
      if (fCentralityClass=="60100V0A") InitDplustoKpipi2013pPb60100V0A();
      
      if (fCentralityClass=="020ZNA") InitDplustoKpipi2013pPb020ZNA();
      if (fCentralityClass=="2040ZNA") InitDplustoKpipi2013pPb2040ZNA();
      if (fCentralityClass=="4060ZNA") InitDplustoKpipi2013pPb4060ZNA();
      if (fCentralityClass=="60100ZNA") InitDplustoKpipi2013pPb60100ZNA();
      
      
      if (fCentralityClass=="020CL1") InitDplustoKpipi2013pPb020CL1();
      if (fCentralityClass=="2040CL1") InitDplustoKpipi2013pPb2040CL1();
      if (fCentralityClass=="4060CL1") InitDplustoKpipi2013pPb4060CL1();
      if (fCentralityClass=="60100CL1") InitDplustoKpipi2013pPb60100CL1();
      
      if (fIsRapidityScan) {
  if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
  if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
  if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
  if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
  if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
      }
    }
    }
    else AliFatal("Not yet implemented");
    break;
  case 3: // D*->D0pi
    if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
    if (fCollisionType==0) {
      if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
      else if(fRunNumber == 10 || fRunNumber==2010){
  if(fIsPass4Analysis) InitDstartoD0pi2010ppPass4();
  else InitDstartoD0pi2010pp();
      } else if(fRunNumber == 12 || fRunNumber==2012){
  InitDstartoD0pi2012pp();
      } else AliFatal("Not yet implemented");
    }
    else if (fCollisionType==1) {
      if (fRunNumber == 10  || fRunNumber==2010){
  if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
  else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
  else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
  else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
  else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
  else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
  else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
  else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
  else AliFatal("Not yet implemented");
      }
      if (fRunNumber == 11 || fRunNumber==2011){
  if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
  else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
  else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
  else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
  else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
  else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
  else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
  else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
  else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
  else AliFatal("Not yet implemented");
      }
      else if (fRunNumber == 15 || fRunNumber==2015){
  if (fCentralityClass=="010") InitDstartoD0pi2015PbPb010();
  else if (fCentralityClass=="3050") InitDstartoD0pi2015PbPb3050();
  else if (fCentralityClass=="6080") InitDstartoD0pi2015PbPb6080();
  else AliFatal("Not yet implemented");
      }
    }
    else if (fCollisionType==2) { 
      if (fRunNumber == 16 || fRunNumber==2016){
  if (fCentralityClass=="0100")InitDstartoD0pi2016pPb0100();    
  else if (fCentralityClass=="010ZNA")InitDstartoD0pi2016pPb010ZNA();
  else if (fCentralityClass=="60100ZNA")InitDstartoD0pi2016pPb60100ZNA();
  else AliFatal("Not yet implemented");
  

      }
      else if (fRunNumber == 13 || fRunNumber==2013){
  
  if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
        if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
  if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();  
        if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
  
        if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
        if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
        if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
        if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
      
        if (fCentralityClass=="020CL1") InitDstartoD0pi2013pPb020CL1();
        if (fCentralityClass=="2040CL1") InitDstartoD0pi2013pPb2040CL1();
        if (fCentralityClass=="4060CL1") InitDstartoD0pi2013pPb4060CL1();
        if (fCentralityClass=="60100CL1") InitDstartoD0pi2013pPb60100CL1();
      
        if (fIsRapidityScan) {
          if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
          if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
          if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
          if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
          if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
        } 
    }
} 
    else AliFatal("Not yet implemented");
    break;
  case 4: // D+s->KKpi
    if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
    if (fCollisionType==0) {
      if(fIsPass4Analysis) InitDstoKKpi2010ppPass4();
      else if (fRunNumber==16 || fRunNumber==2016) InitDstoKKpi2016pp13TeV();
      else InitDstoKKpi2010pp();
    }
    else if (fCollisionType==1) {
      if (fRunNumber == 15){
  if (fCentralityClass=="010") InitDstoKKpi2015PbPb010();
  else if (fCentralityClass=="3050") InitDstoKKpi2015PbPb3050();
  else if (fCentralityClass=="6080") InitDstoKKpi2015PbPb6080();
  else AliFatal("Not yet implemented");
      }
      else{
  if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
  else if (fCentralityClass=="010") InitDstoKKpi2011PbPb010();
  else if (fCentralityClass=="2050") InitDstoKKpi2011PbPb2050();
  else AliFatal("Not yet implemented");
      }
    }
    else if (fCollisionType==2) {
      if(fRunNumber==13 || fRunNumber==2013) {
  if (fCentralityClass=="0100")           InitDstoKKpi2013pPb0100();
      }
      if(fRunNumber==16 || fRunNumber==2016) {
  if (fCentralityClass=="0100")           InitDstoKKpi2016pPb0100();
  else if (fCentralityClass=="140trkl")   InitDstoKKpi2016pPb140trkl();
  else if (fCentralityClass=="4070trkl")  InitDstoKKpi2016pPb4070trkl();
  else if (fCentralityClass=="70200trkl") InitDstoKKpi2016pPb70200trkl();
      }
    }
    else AliFatal("Not yet implemented");
    break;
  case 5: // Lc->pKpi
    if (fCollisionType==0) {
      if (fIsBDTAnalysis) InitLctopKpi2010ppBDT();
      else                InitLctopKpi2010pp();
    }
    else if (fCollisionType==2) {
      if (fIsBDTAnalysis) InitLctopKpi2013pPbBDT();
      else                InitLctopKpi2013pPb();
    }
    else AliFatal("Not yet implemented");
    break;
  case 6: // Lc->pK0S
    if (fCollisionType==0) InitLctopK0S2010pp();
    else if (fCollisionType==2) {
      if (fIsBDTAnalysis) InitLctopK0S2013pPbBDT();
      else                InitLctopK0S2013pPb();
    }
    else AliFatal("Not yet implemented");
    break;
    
  default:
    printf("Invalid decay type: %d\n",decay);
    break;
  }

}
  
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010pp() {
  // 
  // D0->Kpi syst errors. Responsible: A. Rossi
  //   2010 pp sample
  //

  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010pp");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.22);
  fRawYield->SetBinContent(3,0.1);
  for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
  fPIDEff->SetBinContent(2,0.05); // 5%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  fPartAntipart->SetBinContent(1,1); 
  //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
  
  return;
}
//-------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010ppPass4(){
  // 
  // D0->Kpi syst errors. Responsible: C. Terrevoli
  //  2010 pp sample - pass4
  //
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppPass4");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  fNorm->SetBinContent(1,0.0); // bin 0.-1.
  for(Int_t i=2;i<=36;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  fBR->SetBinContent(1,0.0); // 1.2% PDG2010
  for(Int_t i=2;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
  
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  fTrackingEff->SetBinContent(1,0.0); //  
  fTrackingEff->SetBinContent(2,0.03); //  
  fTrackingEff->SetBinContent(3,0.035); // 
  for(Int_t i=4;i<=5;i++) fTrackingEff->SetBinContent(i,0.04); //  
  for(Int_t i=6;i<=9;i++) fTrackingEff->SetBinContent(i,0.045); //  
  for(Int_t i=10;i<=14;i++) fTrackingEff->SetBinContent(i,0.04); //  
  for(Int_t i=15;i<=36;i++) fTrackingEff->SetBinContent(i,0.035); //  

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0.04);
  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=13;i<=36;i++) fRawYield->SetBinContent(i,0.05);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.0);
  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.05);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); 

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(2,0.03);

  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  fPartAntipart->SetBinContent(1,1); 
  //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);

  return;
}
//-------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2015pp5TeV(){
  // 
  // D0->Kpi syst errors. Responsible: C. Terrevoli
  //  2015 pp sample - 5TeV
  //
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2015pp5TeV");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  fNorm->SetBinContent(1,0.0); // bin 0.-1.
  for(Int_t i=2;i<=36;i++) fNorm->SetBinContent(i,0.023); 

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  fBR->SetBinContent(1,0.0); // 1.2% PDG2010
  for(Int_t i=2;i<=36;i++) fBR->SetBinContent(i,0.01); //

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  fTrackingEff->SetBinContent(1,0.0); //  0-1
  fTrackingEff->SetBinContent(2,0.045); //  1-2
  fTrackingEff->SetBinContent(3,0.04); // 2-3
  fTrackingEff->SetBinContent(4,0.035); // 3-4
  for(Int_t i=5;i<=10;i++) fTrackingEff->SetBinContent(i,0.03); // 
  for(Int_t i=11;i<=12;i++) fTrackingEff->SetBinContent(i,0.035); // 
  for(Int_t i=13;i<=16;i++) fTrackingEff->SetBinContent(i,0.045); // 
  for(Int_t i=17;i<=36;i++) fTrackingEff->SetBinContent(i,0.05); //  

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0.03);
  for(Int_t i=2;i<=7;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=12;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.06);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.0);
  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.05);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(2,0.03);

  // particle-antiparticle                                              
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  fPartAntipart->SetBinContent(1,1);                                
  //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);    

  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb020() {
  // 
  // D0->Kpi syst errors. Responsible: A. Rossi
  //   2010 PbPb sample, 0-20 CC
  //

  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb020");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0);
  fRawYield->SetBinContent(3,0.08);
  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.);
  fCutsEff->SetBinContent(3,0.13);
  fCutsEff->SetBinContent(4,0.11);  
  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
  for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(3,0.04);
  fMCPtShape->SetBinContent(4,0.02);
  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03); 

  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
  //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
  //   fPartAntipart->SetBinContent(3,0.10);
  //   fPartAntipart->SetBinContent(4,0.10);
  //   fPartAntipart->SetBinContent(7,0.10);
  //   fPartAntipart->SetBinContent(8,0.10);
  
  return;
}
//--------------------------------------------------------------------------
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb07half() {
  // 
  // D0->Kpi syst errors. Responsible: A. Rossi
  //   2011 PbPb sample, 0-7.5 CC
  //

  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb07half");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0.1);
  fRawYield->SetBinContent(3,0.05);
  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.15);
  fCutsEff->SetBinContent(3,0.13);
  fCutsEff->SetBinContent(4,0.11);  
  fCutsEff->SetBinContent(5,0.08); 
  for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
  for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(2,0.06);
  fMCPtShape->SetBinContent(3,0.04);
  fMCPtShape->SetBinContent(4,0.02);
  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03); 
  for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE

  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
  //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
  //   fPartAntipart->SetBinContent(3,0.10);
  //   fPartAntipart->SetBinContent(4,0.10);
  //   fPartAntipart->SetBinContent(7,0.10);
  //   fPartAntipart->SetBinContent(8,0.10);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
  //
  // D0->Kpi syst errors. Responsible: D. Caffarri
  //   2011 PbPb sample, 30-50 CC InPlane
  //
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050InPlane");

  InitD0toKpi2011PbPb07half();
  // Raw yield extraction
  // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0);
  fRawYield->SetBinContent(3,0.05);
  fRawYield->SetBinContent(4,0.08);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE

  // Cuts efficiency (from cuts variation)
  //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.0);
  fCutsEff->SetBinContent(3,0.10);
  fCutsEff->SetBinContent(4,0.10);
  fCutsEff->SetBinContent(5,0.10);
  fCutsEff->SetBinContent(6,0.10);
  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
  //
  // D0->Kpi syst errors. Responsible: D. Caffarri
  //   2011 PbPb sample, 30-50 CC OutOfPlane
  //
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050OutOfPlane");
  InitD0toKpi2011PbPb07half();
  // Raw yield extraction
  //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0.);
  fRawYield->SetBinContent(3,0.05);
  for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE

  // Cuts efficiency (from cuts variation)
  // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.);
  fCutsEff->SetBinContent(3,0.10);
  fCutsEff->SetBinContent(4,0.10);
  fCutsEff->SetBinContent(5,0.10);
  fCutsEff->SetBinContent(6,0.10);
  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb4080() {
  //
  // D0->Kpi syst errors. Responsible: A. Rossi
  //   2010 PbPb sample, 40-80 CC
  //
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb4080");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0);
  for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.);
  fCutsEff->SetBinContent(3,0.13);
  fCutsEff->SetBinContent(4,0.11);  
  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
  for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
  //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
  for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  //  fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
  //  fMCPtShape->SetBinContent(4,0.02);
  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03); 
  
  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
  //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb010() {
  //
  // D0->Kpi syst errors. Responsible: A. Rossi
  //   2011 PbPb sample, 0-10 CC
  //

  AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb010");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
    
  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
    
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
    
    
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0.15);
  fRawYield->SetBinContent(3,0.09);
  fRawYield->SetBinContent(4,0.06);
  fRawYield->SetBinContent(5,0.06);
  fRawYield->SetBinContent(6,0.06);
  for (Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.06);
  for (Int_t i=9; i<=12; i++)fRawYield->SetBinContent(i,0.07);
  for (Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.11);
  for (Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.30);
    
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.15);
  fCutsEff->SetBinContent(3,0.12);
  fCutsEff->SetBinContent(4,0.05);
  fCutsEff->SetBinContent(5,0.05);
  fCutsEff->SetBinContent(6,0.05);
  for (Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i,0.05);
  for (Int_t i=9; i<=12; i++)fCutsEff->SetBinContent(i,0.05);
  for (Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i,0.05);
  for (Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i,0.05);

    
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
  //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
  for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
    
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  fMCPtShape->SetBinContent(1,0.);
  fMCPtShape->SetBinContent(2,0.15);
  fMCPtShape->SetBinContent(3,0.12);
  fMCPtShape->SetBinContent(4,0.02);
  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
  
  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
  //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
    
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb3050() {
  //
  // D0->Kpi syst errors. Responsible: D. Caffarri
  //   2011 PbPb sample, 30-50 CC
  //
    
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
    
  // Branching ratio
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
    
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
    
    
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,0);
  fRawYield->SetBinContent(2,0.102);
  fRawYield->SetBinContent(3,0.071);
  fRawYield->SetBinContent(4,0.052);
  fRawYield->SetBinContent(5,0.063);
  fRawYield->SetBinContent(6,0.071);
  fRawYield->SetBinContent(7,0.071);
  fRawYield->SetBinContent(8,0.071);
  fRawYield->SetBinContent(9,0.071);
  fRawYield->SetBinContent(10,0.071);
  fRawYield->SetBinContent(11,0.071);
  fRawYield->SetBinContent(12,0.071);
    
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
  for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
    
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.10);
  fCutsEff->SetBinContent(3,0.10);
  fCutsEff->SetBinContent(4,0.10);
  fCutsEff->SetBinContent(5,0.10);
  fCutsEff->SetBinContent(6,0.10);
  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);

    
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
  //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
  for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
    
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(2,0.06);
  fMCPtShape->SetBinContent(3,0.04);
  fMCPtShape->SetBinContent(4,0.02);
  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
  //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
    
  return;
}


//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010ppLowEn() {
  // 
  // D0->Kpi syst errors. Low energy run
  //   2011 2.76 TeV pp sample
  //
  AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppLowEn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); //  error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,1);
  for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10% 
  fCutsEff->SetBinContent(2,0.20);
  for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10% 


  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
  //  fPIDEff->SetBinContent(2,0.20);
  for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  //   // particle-antiparticle
  //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
  //   fPartAntipart->SetBinContent(1,1);
  //   fPartAntipart->SetBinContent(2,1);
  //   for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
  //   for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
  
  return;
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010ppLowPtAn() {
  // 
  // D0->Kpi syst errors. Responsible: 
  //   2010 pp sample, analysis without topological cuts
  //

  AliInfo(" Settings for D0 --> K pi, pp collisions at 7 TeV, analysis without topological cuts"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppLowPtAn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88) 

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
  fRawYield->SetBinContent(1,0.14);
  fRawYield->SetBinContent(2,0.14);
  fRawYield->SetBinContent(3,0.10);
  fRawYield->SetBinContent(4,0.10);
  fRawYield->SetBinContent(5,0.10);
  fRawYield->SetBinContent(6,0.14);
  fRawYield->SetBinContent(7,0.18);
  fRawYield->SetBinContent(8,0.18);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.2);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.); 

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
  fPIDEff->SetBinContent(1,0.05); // 5%
  fPIDEff->SetBinContent(2,0.05); // 5%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);

  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
  
  return;
}
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100(){
  //
  // D0->Kpi syst errors. p-Pb data sample
  //
  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100");

  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); // 
  
  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,0.0);
  fRawYield->SetBinContent(2,0.08);
  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(1,0.0);
  fCutsEff->SetBinContent(2,0.08);
  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  fMCPtShape->SetBinContent(1,0.0);
  fMCPtShape->SetBinContent(2,0.02);
  for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;


}

//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100LowPtAn(){
  //
  // D0->Kpi syst errors. p-Pb data sample
  // analysis without topological cuts

  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, analysis without topological cuts"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100LowPtAn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88) 

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
  fRawYield->SetBinContent(1,0.15);
  fRawYield->SetBinContent(2,0.15);
  fRawYield->SetBinContent(3,0.10);
  fRawYield->SetBinContent(4,0.10);
  fRawYield->SetBinContent(5,0.10);
  fRawYield->SetBinContent(6,0.10);
  fRawYield->SetBinContent(7,0.10);
  fRawYield->SetBinContent(8,0.10);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.11);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.); 

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); 

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);

  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);

  return;


}
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2016pPb0100(){
  //
  // D0->Kpi syst errors. p-Pb data sample
  //
  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV - 2016 run"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); // 

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88) 


  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); // 
  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); // 
  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); // 
  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); // 

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0.0);
  fRawYield->SetBinContent(2,0.05);
  for(Int_t i=3;i<=17;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=24;i<=36;i++) fRawYield->SetBinContent(i,0.04);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.0);
  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  fMCPtShape->SetBinContent(1,0.0);
  fMCPtShape->SetBinContent(2,0.00);
  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;


}
//_______________________________________________________________________
void AliHFSystErr::InitD0toKpi2016pPb010ZNA(){
  //
  // D0->Kpi syst errors. p-Pb data sample
  //
  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, 010 ZNA - 2016 run");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); // 

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88) 


  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); // 
  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); // 
  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); // 
  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); // 

  // Raw yield extraction //aggiorna
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0.0);
  fRawYield->SetBinContent(2,0.08);
  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.035);
  for(Int_t i=16;i<=24;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.0);
  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  fMCPtShape->SetBinContent(1,0.0);
  fMCPtShape->SetBinContent(2,0.00);
  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}
//________________________________________________________________________
void AliHFSystErr::InitD0toKpi2016pPb60100ZNA(){
  //
  // D0->Kpi syst errors. p-Pb data sample
  //
  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV 60100% ZNA - 2016 run");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); // 

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88) 


  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); // 
  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); // 
  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); // 
  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); // 

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,0.0);
  fRawYield->SetBinContent(2,0.06);
  for(Int_t i=3;i<=10;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=11;i<=16;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.06);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,0.0);
  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);


  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  fMCPtShape->SetBinContent(1,0.0);
  fMCPtShape->SetBinContent(2,0.00);
  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2016pPb0100LowPtAn(){
  //
  // D0->Kpi syst errors. p-Pb 2016 data sample
  // analysis without topological cuts

  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV 2016, analysis without topological cuts"); 
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100LowPtAn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.037); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88) 

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); // 
  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); // 
  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); // 

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
  fRawYield->SetBinContent(1,0.06);
  fRawYield->SetBinContent(2,0.06);
  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.05);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.); 

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); 

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);

  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);

  return;


}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010pp() {
  // 
  // D+->Kpipi syst errors. Responsible: R. Bala
  //  2010 pp sample
  //

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010pp");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.25);
  fRawYield->SetBinContent(3,0.25);
  fRawYield->SetBinContent(4,0.20);
  fRawYield->SetBinContent(5,0.09);
  fRawYield->SetBinContent(6,0.09);
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);  
  for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);  
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
  fPIDEff->SetBinContent(1,0.15); // 15%
  fPIDEff->SetBinContent(2,0.15); // 15%
  fPIDEff->SetBinContent(3,0.15); // 15%
  fPIDEff->SetBinContent(4,0.15); // 15%
  for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}
 
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010ppPass4(){
  // 
  // D+->Kpipi syst errors
  //  2010 pp sample - pass4
  //


  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010ppPass4");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  fTrackingEff->SetBinContent(1,0.04);
  fTrackingEff->SetBinContent(2,0.045);
  fTrackingEff->SetBinContent(3,0.045);
  fTrackingEff->SetBinContent(4,0.05);
  fTrackingEff->SetBinContent(5,0.05);
  fTrackingEff->SetBinContent(6,0.055);
  for(Int_t i=7; i<=16; i++) fTrackingEff->SetBinContent(i,0.06);
  for(Int_t i=17; i<=36; i++) fTrackingEff->SetBinContent(i,0.055);


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.06);
  fRawYield->SetBinContent(3,0.06);
  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.03);  
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);  
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);  
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.20);  
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i,0.1); 
  for(Int_t i=4;i<=36;i++) fCutsEff->SetBinContent(i,0.05); 

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.); // 0%

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(2,0.03);
  for(Int_t i=8;i<=24;i++) fMCPtShape->SetBinContent(i,0.02);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2012pp(){
  //
  // D+->Kpipi syst errors
  //  2012 pp sample
  //

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2012pp");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // 4% error on sigmaV0and

  // Branching ratio
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); //   PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.07);
  fRawYield->SetBinContent(3,0.05);
  for(Int_t i=4;i<=36;i++) fRawYield->SetBinContent(i,0.04);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,1.);
  fCutsEff->SetBinContent(2,0.07);
  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.05); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.02);
  for(Int_t i=4;i<=5;i++) fMCPtShape->SetBinContent(i,0.01);

  return;
}




//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2015pp5TeV(){
    //
    // D+->Kpipi syst errors
    //  2015 pp sample
    //Responsible R. Bala
 AliInfo(" Settings for D+ --> K pi pi, pp collisions for 5.02 TeV");
 SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2015pp");
    // Normalization
    fNorm = new TH1F("fNorm","fNorm",24,0,24);
    for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.023); // 4% error on sigmaV0and

    // Branching ratio
    fBR = new TH1F("fBR","fBR",24,0,24);
    for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); //   PDG2015

    // Tracking efficiency
    fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
    fTrackingEff->SetBinContent(2,0.06);
    fTrackingEff->SetBinContent(3,0.06);
    fTrackingEff->SetBinContent(4,0.06);
    fTrackingEff->SetBinContent(5,0.06);
    fTrackingEff->SetBinContent(6,0.05);
    fTrackingEff->SetBinContent(7,0.05);
    fTrackingEff->SetBinContent(8,0.05);
    fTrackingEff->SetBinContent(9,0.04);
    fTrackingEff->SetBinContent(10,0.04);
    for(Int_t i=11;i<=16;i++) fTrackingEff->SetBinContent(i,0.05); 
    for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); 
    
    // Raw yield extraction
    fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
    fRawYield->SetBinContent(2,0.08);
    fRawYield->SetBinContent(3,0.04);
    fRawYield->SetBinContent(4,0.02);
    fRawYield->SetBinContent(5,0.02);
    fRawYield->SetBinContent(6,0.02);
    fRawYield->SetBinContent(7,0.04);
    fRawYield->SetBinContent(8,0.04);
    fRawYield->SetBinContent(9,0.03);
    fRawYield->SetBinContent(10,0.03);
    fRawYield->SetBinContent(11,0.06);
    fRawYield->SetBinContent(12,0.06);


    for(Int_t i=13;i<=24;i++) fRawYield->SetBinContent(i,0.04);

    // Cuts efficiency (from cuts variation)
    fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
    fCutsEff->SetBinContent(2,0.07);
    fCutsEff->SetBinContent(3,0.02);
    fCutsEff->SetBinContent(4,0.02);
    fCutsEff->SetBinContent(5,0.02);
    fCutsEff->SetBinContent(6,0.02);
    fCutsEff->SetBinContent(7,0.03);
    fCutsEff->SetBinContent(8,0.03);
    fCutsEff->SetBinContent(9,0.04);
    fCutsEff->SetBinContent(10,0.04);
   
    for(Int_t i=11;i<=16;i++) fCutsEff->SetBinContent(i,0.03); 
    for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.04); 

    // PID efficiency (from PID/noPID)
    fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
    for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00); // 0%

    // MC dN/dpt
    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
    for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
       return;
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pp13TeV(){
  //
  // D+->Kpipi syst errors
  //  2016 pp sample
  //

 AliInfo(" Settings for D+ --> K pi pi, pp collisions for 13 TeV");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pp13TeV");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // 4% error on sigmaV0and

  // Branching ratio
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); //   PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  //  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)
  fTrackingEff->SetBinContent(2,0.055);
  fTrackingEff->SetBinContent(3,0.055);
  for(Int_t i=4;i<=6;i++) fTrackingEff->SetBinContent(i,0.06); // 
  for(Int_t i=7;i<=24;i++) fTrackingEff->SetBinContent(i,0.07); // 
  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.075); //
  


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.09);
  fRawYield->SetBinContent(3,0.03);
  for(Int_t i=4;i<=10;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=11;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.09);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(1,1.);
  fCutsEff->SetBinContent(2,0.07);
  fCutsEff->SetBinContent(3,0.06);
  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.04); // 10%
 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.04); // 10%
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.04);

  return;
}


//--------------------------------------------------------------------------
void AliHFSystErr::InitDstoKKpi2010pp() {
  // 
  // D+s->KKpi syst errors. Responsible: G.M. Innocenti
  //  2010 pp sample
  //


  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2010pp");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",12,0,12);
  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,1);
  fRawYield->SetBinContent(3,0.20);
  fRawYield->SetBinContent(4,0.20);
  fRawYield->SetBinContent(5,0.15);
  fRawYield->SetBinContent(6,0.15);
  fRawYield->SetBinContent(7,0.15);
  fRawYield->SetBinContent(8,0.15);
  fRawYield->SetBinContent(9,0.20);
  fRawYield->SetBinContent(10,0.20);
  fRawYield->SetBinContent(11,0.20);
  fRawYield->SetBinContent(12,0.20);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}
   
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstoKKpi2010ppPass4() {
  // 
  // D+s->KKpi syst errors. Responsible: A. Barbano
  // 2010 pp sample -- pass4
  //


  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2010ppPass4");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",12,0,12);
  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.035); // 0.08/2.27 PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
  for(Int_t i=1; i<=2; i++) fTrackingEff->SetBinContent(i,0.03);
  for(Int_t i=3; i<=4; i++) fTrackingEff->SetBinContent(i,0.05);
  for(Int_t i=5; i<6; i++) fTrackingEff->SetBinContent(i,0.055);
  for(Int_t i=6; i<=12; i++) fTrackingEff->SetBinContent(i,0.06);

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.05); //5%
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.07); // 7%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);

  return;
}
 
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
  // 
  // D+->Kpipi syst errors. Responsible: E. Bruna
  //  2011 PbPb sample, 0-7.5% CC
  //

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb07half");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
  fRawYield->SetBinContent(4,0.30);
  fRawYield->SetBinContent(5,0.20);
  fRawYield->SetBinContent(6,0.20);
  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);


  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
  for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);


  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%

  // MC dN/dpt  (24/7/2012)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
  for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
  //  for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */

  return;
}
//_________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
  // 
  // D+->Kpipi syst errors. Responsible: E. Bruna, C. Bedda 
  //  2011 PbPb sample, 0-10% CC
  //

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb010");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i,0.);  //5 to 10%
  fRawYield->SetBinContent(4,0.10);
  fRawYield->SetBinContent(5,0.08);
  fRawYield->SetBinContent(6,0.08);
  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.08);


  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  for(Int_t i=1; i<=36; i++) fCutsEff->SetBinContent(i,0.10);


  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.); 
  fMCPtShape->SetBinContent(4,0.06);
  fMCPtShape->SetBinContent(5,0.02);
  for (Int_t i=6; i<=36; i++) fMCPtShape->SetBinContent(i,0.01);
  //  for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */

  return;
}

//___________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
  //
  // D+->Kpipi syst errors. Responsible: R. Bala
  //  2011 PbPb sample, 30-50% CC
  //
  AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb3050");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",16,0,16);
  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
    
  // Branching ratio
  fBR = new TH1F("fBR","fBR",16,0,16);
  for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
    
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
    
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09);  //5 to 10%
  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
    

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
  for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
  for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
    
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05); // 5%
    
  // MC dN/dpt  (24/7/2012)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
  fMCPtShape->SetBinContent(3,0.10);
  fMCPtShape->SetBinContent(4,0.04);
  for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
    
    
  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
    
  return;
}


//--------------------------------------------------------------------------
void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
  // 
  // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
  //  2011 PbPb sample, 0-7.5% CC
  //
  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 0-7.5 CC");
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb07half");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",12,0,12);
  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
  for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30); 
  for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
 


  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */

  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstoKKpi2011PbPb010() {
  // 
  // D+s->Kpipi syst errors. Responsible: A.Barbano
  //  2011 PbPb sample, 0-10% CC
  //

  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 0-10 CC");
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb010");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",12,0,12);
  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.08); // 8%


  // Cuts efficiency (from cuts variation + phi mass selection cuts)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%  

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 10%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
  for(Int_t i=1; i<=6; i++) fMCPtShape->SetBinContent(i,0.02);  //2%
  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.01);  //1%
  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%

  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstoKKpi2011PbPb2050() {
  // 
  // D+s->Kpipi syst errors. Responsible: A.Barbano
  //  2011 PbPb sample, 20-50% CC
  //

  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 20-50 CC");
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb2050");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",12,0,12);
  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.051); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.08); //8%


  // Cuts efficiency (from cuts variation + phi selection cuts)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%

  // MC dN/dpt 
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
  for(Int_t i=1; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%

  return;
}

//_________________________________________________________________________
void AliHFSystErr::InitDstoKKpi2013pPb0100(){
  //
  // Dstar syst
  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2013pPb0100");
 
  //Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); // 
  
  // Branching ratio 
  fBR = new TH1F("fBR","fBR",12,0,12);
  for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10); 
  for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05); 

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}

//_________________________________________________________________________
void AliHFSystErr::InitDstoKKpi2016pPb0100(){
  // Ds systematics in pPb 2016 MB
  // Responsible: J. Hamon

  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb0100");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037); // Same as D0 low pt?

  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08

  // Tracking efficiency
  //   From C. Terrevoli's studies:  3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;  i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
  for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fRawYield->SetBinContent(i, 0.03); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fRawYield->SetBinContent(i, 0.01); // [4-6]
  for(Int_t i=7;  i<=8;  i++) fRawYield->SetBinContent(i, 0.01); // [6-8]
  for(Int_t i=9;  i<=12; i++) fRawYield->SetBinContent(i, 0.01); // [8-12]
  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.03); // [12-16]
  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]

  // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
  //   ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
  for(Int_t i=7;  i<=8;  i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
  for(Int_t i=9;  i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]

  // PID efficiency
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
  for(Int_t i=5;  i<=6;  i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
  for(Int_t i=7;  i<=8;  i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
  for(Int_t i=9;  i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.00); // [12-16] - ConservativePID
  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.00); // [16-24] - ConservativePID

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.00); // No assigned systematics

  return;
}
//_________________________________________________________________________
 void AliHFSystErr::InitDstoKKpi2016pPb140trkl(){
   // Ds systematics in pPb 2016, 1-40 tracklets
   // Responsible: A. Barbano
        
   AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [1,40) trkls");
   SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb140trkls");
        
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
   for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
        
   // Branching ratio
   fBR = new TH1F("fBR","fBR",24,0,24);
   for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
        
   // Tracking efficiency
   //   From C. Terrevoli's studies:  3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
   for(Int_t i=1;  i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
   for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
        
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fRawYield->SetBinContent(i, 0.03); // [6-8]
   for(Int_t i=9;  i<=12; i++) fRawYield->SetBinContent(i, 0.03); // [8-12]
   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.09); // [12-16]
   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
        
   // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
   // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
   for(Int_t i=9;  i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
   for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
   for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
        
   // PID efficiency
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
   for(Int_t i=5;  i<=6;  i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
   for(Int_t i=7;  i<=8;  i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
   for(Int_t i=9;  i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
   for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
   for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
        
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
   for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
     
   return;
 }
    
 //_________________________________________________________________________
 void AliHFSystErr::InitDstoKKpi2016pPb4070trkl(){
   // Ds systematics in pPb 2016, 40-70 tracklets
   // Responsible: A. Barbano
        
   AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [40,70) trkls");
   SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb4070trkls");
        
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
   for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
        
   // Branching ratio
   fBR = new TH1F("fBR","fBR",24,0,24);
   for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
        
   // Tracking efficiency
   //   From C. Terrevoli's studies:  3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
   for(Int_t i=1;  i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
   for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
        
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fRawYield->SetBinContent(i, 0.03); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fRawYield->SetBinContent(i, 0.03); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fRawYield->SetBinContent(i, 0.04); // [6-8]
   for(Int_t i=9;  i<=12; i++) fRawYield->SetBinContent(i, 0.05); // [8-12]
   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.12); // [12-16]
   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
        
   // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
   // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
   for(Int_t i=9;  i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
   for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
   for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
        
   // PID efficiency
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
   for(Int_t i=5;  i<=6;  i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
   for(Int_t i=7;  i<=8;  i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
   for(Int_t i=9;  i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
   for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
   for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
        
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
   for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
     
   return;
 }
    
 //_________________________________________________________________________
 void AliHFSystErr::InitDstoKKpi2016pPb70200trkl(){
   // Ds systematics in pPb 2016, 70-200 tracklets
   // Responsible: A. Barbano
        
   AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [70,200) trkls");
   SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb70200trkls");
        
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
   for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
        
   // Branching ratio
   fBR = new TH1F("fBR","fBR",24,0,24);
   for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
        
   // Tracking efficiency
   //   From C. Terrevoli's studies:  3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
   for(Int_t i=1;  i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
   for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
        
   // Raw yield extraction 
   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fRawYield->SetBinContent(i, 0.04); // [6-8]
   for(Int_t i=9;  i<=12; i++) fRawYield->SetBinContent(i, 0.04); // [8-12]
   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.10); // [12-16]
   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
        
   // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
   // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
   for(Int_t i=5;  i<=6;  i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
   for(Int_t i=7;  i<=8;  i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
   for(Int_t i=9;  i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
   for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
   for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
        
   // PID efficiency
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
   for(Int_t i=1;  i<=2;  i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
   for(Int_t i=3;  i<=4;  i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
   for(Int_t i=5;  i<=6;  i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
   for(Int_t i=7;  i<=8;  i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
   for(Int_t i=9;  i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
   for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
   for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
        
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
   for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
        
   return;
 }

//_________________________________________________________________________
void AliHFSystErr::InitDstoKKpi2016pp13TeV(){
  // Ds systematics in MB pp 13 TeV (2016)
  // Responsible: J. Hamon

  AliInfo(" Settings for Ds --> K K pi, MB pp collisions at 13 TeV - 2016 runs (LHC16k,l)");
  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pp13TeV");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.05); // from V0

  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fTrackingEff->SetBinContent(i, 1.000); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fTrackingEff->SetBinContent(i, 0.055); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fTrackingEff->SetBinContent(i, 0.060); // [4-6]
  for(Int_t i=7;  i<=12; i++) fTrackingEff->SetBinContent(i, 0.070); // [6-8] + [8-12]
  for(Int_t i=12; i<=16; i++) fTrackingEff->SetBinContent(i, 0.075); // [12-16]
  for(Int_t i=17; i<=24; i++) fTrackingEff->SetBinContent(i, 1.000); // [16-24]

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
  for(Int_t i=7;  i<=8;  i++) fRawYield->SetBinContent(i, 0.03); // [6-8]
  for(Int_t i=9;  i<=12; i++) fRawYield->SetBinContent(i, 0.04); // [8-12]
  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.08); // [12-16]
  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 1.00); // [16-24]

  // Cuts efficiency
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fCutsEff->SetBinContent(i, 0.18); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
  for(Int_t i=7;  i<=8;  i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
  for(Int_t i=9;  i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.08); // [12-16]
  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 1.00); // [16-24]

  // PID efficiency
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;  i<=2;  i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
  for(Int_t i=3;  i<=4;  i++) fPIDEff->SetBinContent(i, 0.00); // [2-4]
  for(Int_t i=5;  i<=6;  i++) fPIDEff->SetBinContent(i, 0.00); // [4-6]
  for(Int_t i=7;  i<=8;  i++) fPIDEff->SetBinContent(i, 0.00); // [6-8]
  for(Int_t i=9;  i<=12; i++) fPIDEff->SetBinContent(i, 0.00); // [8-12]
  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.00); // [12-16]
  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 1.00); // [16-24]

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.00); // No assigned systematics
  for(Int_t i=3; i<=4;  i++) fMCPtShape->SetBinContent(i, 0.01); // [2-4] from toy MC pt shape

  return;
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
  // 
  // D+->Kpipi syst errors. Responsible: ??
  //  2010 PbPb sample, 0-20 CC
  //

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb020");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
  // fRawYield->SetBinContent(5,0.23);
  //fRawYield->SetBinContent(6,0.23);
  fRawYield->SetBinContent(7,0.20);
  fRawYield->SetBinContent(8,0.20);
  fRawYield->SetBinContent(9,0.15);
  fRawYield->SetBinContent(10,0.15);
  fRawYield->SetBinContent(11,0.15);
  fRawYield->SetBinContent(12,0.15);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%

  // MC dN/dpt  (2/2/2012)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
  for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */

  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
  // 
  // D+->Kpipi syst errors. Responsible: ??
  //  2010 PbPb sample, 40-80 CC
  //
  

  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb4080");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% 

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,1);
  fRawYield->SetBinContent(3,1);
  fRawYield->SetBinContent(4,0.15);
  fRawYield->SetBinContent(5,0.05);
  fRawYield->SetBinContent(6,0.05);
  fRawYield->SetBinContent(7,0.15);
  fRawYield->SetBinContent(8,0.15);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
  for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1);  //5 to 10%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
  fPIDEff->SetBinContent(3,0.13); // 13%
 

  // MC dN/dpt  (2/2/2012)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
  for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);


  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {

  // 
  // D+->Kpipi syst errors. Responsible: R. Bala
  //  2011 2.76 TeV pp sample
  //
  AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010ppLowEn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,1);
  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10);  //5 to 10%
  fRawYield->SetBinContent(7,0.15);
  fRawYield->SetBinContent(8,0.15); 
  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055);  //5 to 10%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
  fPIDEff->SetBinContent(3,0.10); // 13%
  fPIDEff->SetBinContent(4,0.10); // 13%
 
  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  return;
}

//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
  //
  // Dplus syst 
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100");

  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); // 
  
  fBR = new TH1F("fBR","fBR",20,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,0.);
  fRawYield->SetBinContent(2,0.1);
  fRawYield->SetBinContent(3,0.08);
  fRawYield->SetBinContent(4,0.05);
  fRawYield->SetBinContent(5,0.05);
  fRawYield->SetBinContent(6,0.05);
  fRawYield->SetBinContent(7,0.05);
  fRawYield->SetBinContent(8,0.05);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.1);
  fCutsEff->SetBinContent(3,0.1);
  fCutsEff->SetBinContent(4,0.06);
  fCutsEff->SetBinContent(5,0.06);
  fCutsEff->SetBinContent(6,0.06);
  fCutsEff->SetBinContent(7,0.05);
  fCutsEff->SetBinContent(8,0.05);
  for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  fPIDEff->SetBinContent(1,1);
  fPIDEff->SetBinContent(2,0.1);
  for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;


}

//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2016pPb0100(){
  //
  // Dplus syst 
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV - 2016 run"); 
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb0100");

  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);

  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(2,0.03);
  fRawYield->SetBinContent(3,0.01);
  fRawYield->SetBinContent(4,0.01);
  fRawYield->SetBinContent(5,0.01);
  for(Int_t i=6;i<=12;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.04);
  

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  fCutsEff->SetBinContent(2,0.12);
  fCutsEff->SetBinContent(3,0.07);
  for(Int_t i=4;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
  for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.05);
  

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); //

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  return;

}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pPb010ZNA() {
  //
  // D+->Kpipi syst errors. Responsible: F. Grosa
  //  2016 pPb sample, 0-10% ZNA
  //
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 0-10% ZNA - 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb010ZNA");
  
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); //  0.24/9.46 PDG2016
  
  // Tracking efficiency (from MB)
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
  //for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(2,0.05);
  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=7;i<=16;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(2,0.07);
  fCutsEff->SetBinContent(3,0.04);
  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.03);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); //
  
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=7;i++) fMCPtShape->SetBinContent(i,0.01);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pPb60100ZNA() {
  //
  // D+->Kpipi syst errors. Responsible: F. Grosa
  //  2016 pPb sample, 60-100% ZNA
  //
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 60-100% ZNA - 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb60100ZNA");
  
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); //  0.24/9.46 PDG2016
  
  // Tracking efficiency (from MB)
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
  //for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(2,0.08);
  for(Int_t i=3;i<=7;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(2,0.07);
  fCutsEff->SetBinContent(3,0.04);
  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.03);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); //
  
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  fMCPtShape->SetBinContent(2,0.03);
  fMCPtShape->SetBinContent(3,0.02);
  for(Int_t i=4;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
  
  return;
}
  
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pPb140trkl() {
  
  //
  // D+->Kpipi syst errors. Responsible: F. Grosa
  //  2016 pPb sample, 1-40 tracklets
  //
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
  
  fNorm = new TH1F("fNorm","fNorm",16,0,16);
  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",16,0,16);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); //  0.24/9.46 PDG2016
  
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
  
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pPb4070trkl() {
  
  //
  // D+->Kpipi syst errors. Responsible: F. Grosa
  //  2016 pPb sample, 40-70 tracklets
  //
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
  
  fNorm = new TH1F("fNorm","fNorm",16,0,16);
  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",16,0,16);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); //  0.24/9.46 PDG2016
  
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
  
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2016pPb70200trkl() {
  
  //
  // D+->Kpipi syst errors. Responsible: F. Grosa
  //  2016 pPb sample, 70-200 tracklets
  //
  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
  
  fNorm = new TH1F("fNorm","fNorm",16,0,16);
  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
  
  fBR = new TH1F("fBR","fBR",16,0,16);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); //  0.24/9.46 PDG2016
  
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
  
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.02);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.04);
  
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
  
  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
  
  return;
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010pp() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
  //  2010 pp sample
  //

  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010pp");
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  fTrackingEff->SetBinContent(1,1.0);
  fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
  fTrackingEff->SetBinContent(3,0.12);
  fTrackingEff->SetBinContent(3,0.12);
  for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1.0);
  fRawYield->SetBinContent(2,0.10);
  fRawYield->SetBinContent(3,0.04);
  fRawYield->SetBinContent(4,0.03);
  fRawYield->SetBinContent(5,0.03);
  fRawYield->SetBinContent(6,0.05);
  fRawYield->SetBinContent(7,0.05);
  fRawYield->SetBinContent(8,0.05);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);  //4%
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);  //4%
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(2,0.22);
  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
 

  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  return;


}
void AliHFSystErr::InitDstartoD0pi2010ppPass4(){
  // 
  // D*+->Kpipi syst errors. Responsible: A. Grelli
  //  2010 pp sample - pass4
  //
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010ppPass4");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  fTrackingEff->SetBinContent(1,1.0);
  fTrackingEff->SetBinContent(2,0.06); //
  fTrackingEff->SetBinContent(3,0.06);
  fTrackingEff->SetBinContent(3,0.06);
  for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); //


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1.0);
  fRawYield->SetBinContent(2,0.07);
  fRawYield->SetBinContent(3,0.02);
  fRawYield->SetBinContent(4,0.02);
  fRawYield->SetBinContent(5,0.01);
  fRawYield->SetBinContent(6,0.01);
  fRawYield->SetBinContent(7,0.01);
  fRawYield->SetBinContent(8,0.01);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.02);  //4%
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.02);  //4%
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.02);  //4%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(2,0.07);
  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.0);  //4%
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);  //4%

  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0); // 3%
 

  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.02);


  // particle-antiparticle
  //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
  //  fPartAntipart->SetBinContent(1,1); 
  //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);

  return;
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {

  // 
  // D+->Kpipi syst errors. Responsible: A. Grelli
  //  2011 2.76 TeV pp sample
  //
  AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010ppLowEn");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",20,0,20);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",20,0,20);
  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,1);
  fRawYield->SetBinContent(3,0.14);
  fRawYield->SetBinContent(4,0.14);
  fRawYield->SetBinContent(5,0.12);
  fRawYield->SetBinContent(6,0.12);
  fRawYield->SetBinContent(7,0.06);
  fRawYield->SetBinContent(8,0.06);
  fRawYield->SetBinContent(9,0.08);
  fRawYield->SetBinContent(10,0.08);
  fRawYield->SetBinContent(11,0.08);
  fRawYield->SetBinContent(12,0.08);
  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);  
  fCutsEff->SetBinContent(3,0.15);
  fCutsEff->SetBinContent(4,0.15);
  fCutsEff->SetBinContent(5,0.15);
  fCutsEff->SetBinContent(6,0.15);
  fCutsEff->SetBinContent(7,0.10);
  fCutsEff->SetBinContent(8,0.10);
  fCutsEff->SetBinContent(9,0.10);
  fCutsEff->SetBinContent(10,0.10);
  fCutsEff->SetBinContent(11,0.10);
  fCutsEff->SetBinContent(12,0.10);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%

  // MC dN/dpt
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  return;
}


//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2012pp() {
  //
  // D*+->D0pi syst errors.
  //  2012 pp sample
  //

  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2012pp");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and

  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2015

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1.0);
  fRawYield->SetBinContent(2,0.08);
  fRawYield->SetBinContent(3,0.06);
  for(Int_t i=4;i<=8;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=9;i<=24;i++) fRawYield->SetBinContent(i,0.02);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(1,1.0);
  for(Int_t i=2;i<=3;i++) fCutsEff->SetBinContent(i,0.10);
  for(Int_t i=4;i<=7;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=8;i<=24;i++) fCutsEff->SetBinContent(i,0.01);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.02);
  for(Int_t i=4;i<=5;i++) fMCPtShape->SetBinContent(i,0.01);

  return;
}


//------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2010 PbPb sample, 0-20 CC
  //

  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb020");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
  fRawYield->SetBinContent(3,0.2);
  fRawYield->SetBinContent(4,0.2);
  fRawYield->SetBinContent(5,0.2);
  fRawYield->SetBinContent(6,0.2);
 
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
  fCutsEff->SetBinContent(4,0.15);
  fCutsEff->SetBinContent(5,0.15);
  fCutsEff->SetBinContent(6,0.15);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
 

  // MC dN/dpt  (from study on D* pt shape)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
  fMCPtShape->SetBinContent(4,0.025);
  fMCPtShape->SetBinContent(5,0.025);
  fMCPtShape->SetBinContent(6,0.025);
  fMCPtShape->SetBinContent(7,0.04);
  fMCPtShape->SetBinContent(8,0.04);
  fMCPtShape->SetBinContent(9,0.03);
  fMCPtShape->SetBinContent(10,0.03);
  fMCPtShape->SetBinContent(11,0.03);
  fMCPtShape->SetBinContent(12,0.03);
  
  return;

}
// ----------------------------- 2011 ----------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2011 PbPb sample, 0-7.5 CC
  //

  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb07half");


  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05);  //4%
  fRawYield->SetBinContent(4,0.2);
  fRawYield->SetBinContent(5,0.10);
  fRawYield->SetBinContent(6,0.10);
  fRawYield->SetBinContent(7,0.08);
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15);  //4%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
  fPIDEff->SetBinContent(4,0.09);

  // MC dN/dpt  (from study on D* pt shape)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
  fMCPtShape->SetBinContent(4,0.015);
  fMCPtShape->SetBinContent(5,0.015);
  fMCPtShape->SetBinContent(6,0.015);
  fMCPtShape->SetBinContent(7,0.02);
  fMCPtShape->SetBinContent(8,0.02);
  fMCPtShape->SetBinContent(9,0.03);
  fMCPtShape->SetBinContent(10,0.03);
  fMCPtShape->SetBinContent(11,0.03);
  fMCPtShape->SetBinContent(12,0.03);
  
  

  return;

}
//-------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2010 PbPb sample, 20-40 CC
  //

  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb2040");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15);  //4%
 
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
 

  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
  fMCPtShape->SetBinContent(1,0.03);
  fMCPtShape->SetBinContent(2,0.03);

  return;

}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2010 PbPb sample, 40-80 CC
  //

  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb4080");

  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
  fRawYield->SetBinContent(1,1);
  fRawYield->SetBinContent(2,0.15);
  fRawYield->SetBinContent(3,0.15);
  fRawYield->SetBinContent(4,0.15);
  fRawYield->SetBinContent(5,0.15);
  fRawYield->SetBinContent(6,0.10);
  fRawYield->SetBinContent(7,0.10);
  fRawYield->SetBinContent(8,0.10);
  fRawYield->SetBinContent(9,0.11);
  fRawYield->SetBinContent(10,0.11);
  fRawYield->SetBinContent(11,0.11);
  fRawYield->SetBinContent(12,0.11);
  fRawYield->SetBinContent(13,0.08);
  fRawYield->SetBinContent(14,0.08);
  fRawYield->SetBinContent(15,0.08);
  fRawYield->SetBinContent(16,0.08);


  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
 

  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(2,0.05);
  fMCPtShape->SetBinContent(3,0.05);
  fMCPtShape->SetBinContent(4,0.05);
  fMCPtShape->SetBinContent(5,0.04);
  fMCPtShape->SetBinContent(6,0.02);
  fMCPtShape->SetBinContent(7,0.04);
  fMCPtShape->SetBinContent(8,0.04);
 
  return;

}

// ----------------------------- 2011 ----------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
  // 
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2011 PbPb sample, 0-10% CC
  //

  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality "); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb010");


  // Normalization
  fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm

  // Branching ratio 
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track


  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
  fRawYield->SetBinContent(4,0.10);
  fRawYield->SetBinContent(5,0.08);
  fRawYield->SetBinContent(6,0.05);
  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10);  //4%

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
  fCutsEff->SetBinContent(4,0.15);
  for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%

  // MC dN/dpt  (from study on D* pt shape)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
  fMCPtShape->SetBinContent(4,0.06);
  fMCPtShape->SetBinContent(5,0.03);
  for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
 
  return;

}

//----------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
  //
  // D*+->D0pi syst errors. Responsible: A. Grelli
  //  2011 PbPb sample, 30-50 CC
  //
    
  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb3050");
   
  // Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
    
  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
    
  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
    
    
  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
  fRawYield->SetBinContent(1,0.12);
  fRawYield->SetBinContent(2,0.12);
  fRawYield->SetBinContent(3,0.12);
  fRawYield->SetBinContent(4,0.06);
  fRawYield->SetBinContent(5,0.06);
  fRawYield->SetBinContent(6,0.05);
  fRawYield->SetBinContent(7,0.08);
  fRawYield->SetBinContent(8,0.05);
  fRawYield->SetBinContent(9,0.05);
  fRawYield->SetBinContent(10,0.05);
  fRawYield->SetBinContent(11,0.05);
  fRawYield->SetBinContent(12,0.05);
  fRawYield->SetBinContent(13,0.05);
  fRawYield->SetBinContent(14,0.05);
  fRawYield->SetBinContent(15,0.05);
  fRawYield->SetBinContent(16,0.05);
    
    
  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
  fCutsEff->SetBinContent(1,0.16);
  fCutsEff->SetBinContent(2,0.16);
  fCutsEff->SetBinContent(3,0.16);
  fCutsEff->SetBinContent(4,0.10);
  fCutsEff->SetBinContent(5,0.10);
  fCutsEff->SetBinContent(6,0.10);
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
    
    
  // MC dN/dpt  (copied from D0 : will update later)
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
  fMCPtShape->SetBinContent(2,0.11);
  fMCPtShape->SetBinContent(3,0.11);
  fMCPtShape->SetBinContent(4,0.03);
  fMCPtShape->SetBinContent(5,0.02);
  fMCPtShape->SetBinContent(6,0.01);
  fMCPtShape->SetBinContent(7,0.01);
  fMCPtShape->SetBinContent(8,0.01);
    
  return;
}
//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
  // Dstar syst in pPb 2013 MB
  //
  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100");

  //Normalization
  fNorm = new TH1F("fNorm","fNorm",24,0,24);
  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
  
  // Branching ratio
  fBR = new TH1F("fBR","fBR",24,0,24);
  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
  fRawYield->SetBinContent(1,1.0);
  fRawYield->SetBinContent(2,0.08);
  fRawYield->SetBinContent(3,0.02);
  fRawYield->SetBinContent(4,0.02);
  fRawYield->SetBinContent(5,0.05);
  fRawYield->SetBinContent(6,0.02);
  fRawYield->SetBinContent(7,0.04);
  fRawYield->SetBinContent(8,0.04);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05); 
  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
  fCutsEff->SetBinContent(1,0.);
  fCutsEff->SetBinContent(2,0.10);
  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
  for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
  for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);

  // particle-antiparticle
  /*
    fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
    fPartAntipart->SetBinContent(1,1);
    fPartAntipart->SetBinContent(2,1);
    fPartAntipart->SetBinContent(3,0.12);
    for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
  */
  
  return;

}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2016pPb0100(){
  // Dstar syst in pPb 2016 MB
  // Responsible: A.M. Veen

  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV 2016"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2016pPb0100");

  //Normalization
fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.037); // Same as D0 low pt?
  
  // Branching ratio
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2016 (Dstar->d0+d0->kPipi uncertainty)

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36); //Dstar meson: 3.7 in 1-10, 4 in 10-24, 4.5 in 24-36
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);//new version
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);//
  for(Int_t i=25;i<=36; i++) fTrackingEff->SetBinContent(i,0.045);//
  

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(1,1.0);
  fRawYield->SetBinContent(2,0.06);
  fRawYield->SetBinContent(3,0.04);
  fRawYield->SetBinContent(4,0.03);
  fRawYield->SetBinContent(5,0.03);
  fRawYield->SetBinContent(6,0.02);
  fRawYield->SetBinContent(7,0.02);
  fRawYield->SetBinContent(8,0.02);
  for(Int_t i=9;i<=24;i++) fRawYield->SetBinContent(i,0.02); 
  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.04);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  // For now uncertainty added as if we use average with full improver as uncertainty
  fCutsEff->SetBinContent(2,0.04);// Original values are higher then effect improver in most bins
  fCutsEff->SetBinContent(3,0.02);
  for(Int_t i=4;i<=16;i++) fCutsEff->SetBinContent(i,0.02);// Very conservative first values rounded up
  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.03);// a little higher due to variations in d0 only
  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.02);

  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.02); // PID for 2 sigma in TPC

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0); //No systematic assigned final


  
  return;

}



//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2016pPb010ZNA(){
  // Dstar syst in pPb 2016 010ZNA
  // Responsible: C. Bedda

  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV 2016"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2016pPb010ZNA");

  //Normalization
fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.037); 
  
  // Branching ratio
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2016 (Dstar->d0+d0->kPipi uncertainty)

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36); //Dstar meson: 3.7 in 1-10, 4 in 10-24, 4.5 in 24-36
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);//same as MB
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);//
  for(Int_t i=25;i<=36; i++) fTrackingEff->SetBinContent(i,0.045);//
  

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(2,0.14);
  fRawYield->SetBinContent(3,0.07);
  
  for(Int_t i=4;i<=15;i++) fRawYield->SetBinContent(i,0.03); 
  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.05);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);//same as MB
  
 fCutsEff->SetBinContent(2,0.04);

  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.02);
  
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.02); // PID for 2 sigma in TPC

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
  fMCPtShape->SetBinContent(2,0.01);
  for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0); //Systematics on multiplicity weights (MC pt shape negligible)
  for(Int_t i=24;i<=36;i++) fMCPtShape->SetBinContent(i,0.01); //Systematics on multiplicity weights (MC pt shape negligible)


  
  return;

}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2016pPb60100ZNA(){
  // Dstar syst in pPb 2016 60100ZNA
  // Responsible: C. Bedda

  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV 2016"); 
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2016pPb60100ZNA");

  //Normalization
fNorm = new TH1F("fNorm","fNorm",36,0,36);
  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.037); 
  
  // Branching ratio
  fBR = new TH1F("fBR","fBR",36,0,36);
  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2016 (Dstar->d0+d0->kPipi uncertainty)

  // Tracking efficiency
  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36); //Dstar meson: 3.7 in 1-10, 4 in 10-24, 4.5 in 24-36
  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);//same as MB analysis
  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);//
  for(Int_t i=25;i<=36; i++) fTrackingEff->SetBinContent(i,0.045);//
  

  // Raw yield extraction
  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
  fRawYield->SetBinContent(2,0.11);
  fRawYield->SetBinContent(3,0.05);
  fRawYield->SetBinContent(4,0.05);
 
  for(Int_t i=5;i<=24;i++) fRawYield->SetBinContent(i,0.03); 
  for(Int_t i=24;i<=36;i++) fRawYield->SetBinContent(i,0.05);

  // Cuts efficiency (from cuts variation)
  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
  // Same as  MB
 fCutsEff->SetBinContent(2,0.04);
 
  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.02);
  // PID efficiency (from PID/noPID)
  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.02); // PID for 2 sigma in TPC

  // MC dN/dpt  
  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
 //Systematics on multiplicity weights (MC pt shape negligible)
  fMCPtShape->SetBinContent(2,0.03);
  fMCPtShape->SetBinContent(3,0.02);
  for(Int_t i=4;i<=36;i++) fMCPtShape->SetBinContent(i,0.01); 
  

  
  return;

}


//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
  // define errors for RAA vs. centrality
  InitD0toKpi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb010CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);      
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
  // define errors for RAA vs. centrality
  InitD0toKpi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb1020CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);    
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
  // define errors for RAA vs. centrality
  InitD0toKpi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb2040CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
  // define errors for RAA vs. centrality
  InitD0toKpi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb4060CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);  
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
  // define errors for RAA vs. centrality
  InitD0toKpi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb6080CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
  // define errors for RAA vs. centrality
  InitDplustoKpipi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb010CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);

}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
  // define errors for RAA vs. centrality
  InitDplustoKpipi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb1020CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
  // define errors for RAA vs. centrality
  InitDplustoKpipi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb2040CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
  // define errors for RAA vs. centrality
  InitDplustoKpipi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb4060CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
  // define errors for RAA vs. centrality
  InitDplustoKpipi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb6080CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
  // define errors for RAA vs. centrality
  InitDstartoD0pi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb010CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16); 
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
  // define errors for RAA vs. centrality
  InitDstartoD0pi2010PbPb020();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb1020CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05); 
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
  // define errors for RAA vs. centrality
  InitDstartoD0pi2010PbPb2040();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb2040CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
  // define errors for RAA vs. centrality
  InitDstartoD0pi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb4060CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
  // define errors for RAA vs. centrality
  InitDstartoD0pi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb6080CentScan");
  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb010CentScan");
 
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb1020CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb2030CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3040CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb4050CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
  // define errors for RAA vs. centrality 2011
  InitD0toKpi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb5080CentScan");
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb010CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb1020CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb2030CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb3040CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb4050CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDplustoKpipi2010PbPb4080();
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb5080CentScan");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);

  for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
}

//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb010CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb1020CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb2030CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb3040CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb4050CentScan");
}
//--------------------------------------------------------------------------
void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
  // define errors for RAA vs. centrality 2011
  InitDstartoD0pi2011PbPb07half();
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb5080CentScan");
}

//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
  InitD0toKpi2013pPb0100();
  
  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100RapScan0804");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.058);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.058);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.028);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
  
  
};
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
  InitD0toKpi2013pPb0100();

  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100RapScan0401");
  
  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.042);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.058);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.058);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.028);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
};
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
  InitD0toKpi2013pPb0100();
  
  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100RapScan0101");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.104);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.104);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.104);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.028);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
  
};
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
  InitD0toKpi2013pPb0100();
  
  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100RapScan0104");
  
  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.042);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.058);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.058);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.028);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
  
};
//_________________________________________________________________________
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
  InitD0toKpi2013pPb0100();
  
  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100RapScan0408");
 
  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.058);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.058);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.028);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
  
};

//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
  InitDplustoKpipi2013pPb0100();

  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100RapScan0804");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);

};

//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
  InitDplustoKpipi2013pPb0100();

  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100RapScan0401");
  
  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);

  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);

};
//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
  InitDplustoKpipi2013pPb0100();

  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100RapScan0101");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);

  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);

};
//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
  InitDplustoKpipi2013pPb0100();

  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100RapScan0104");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
};
//_________________________________________________________________________
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
  InitDplustoKpipi2013pPb0100();

  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100RapScan0408");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);

  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
  
};


//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
  InitDstartoD0pi2013pPb0100();

  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100RapScan0804");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
};
//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
  InitDstartoD0pi2013pPb0100();

  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");  
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100RapScan0401");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);

  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);

};
//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
  InitDstartoD0pi2013pPb0100();

  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100RapScan0101");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);

  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%

  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);

};
//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
  InitDstartoD0pi2013pPb0100();

  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");  
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100RapScan0104");
 
  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);


};
//_________________________________________________________________________
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
  InitDstartoD0pi2013pPb0100();

  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");  
  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100RapScan0408");

  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
  
  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
  
  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
  
  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.117);
  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.054);
  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.064);
};

void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
  InitD0toKpi2013pPb0100();
  
  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb020V0A");
  
  // Raw yield extraction
  fRawYield->SetBinContent(1,0.0);
  fRawYield->SetBinContent(2,0.085);
  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
  for(