PWGHF/vertexingHF/AliHFSystErr.cxx

   1
   2 /**************************************************************************
   3  * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
   4  *                                                                        *
   5  * Author: The ALICE Off-line Project.                                    *
   6  * Contributors are mentioned in the code where appropriate.              *
   7  *                                                                        *
   8  * Permission to use, copy, modify and distribute this software and its   *
   9  * documentation strictly for non-commercial purposes is hereby granted   *
  10  * without fee, provided that the above copyright notice appears in all   *
  11  * copies and that both the copyright notice and this permission notice   *
  12  * appear in the supporting documentation. The authors make no claims     *
  13  * about the suitability of this software for any purpose. It is          *
  14  * provided "as is" without express or implied warranty.                  *
  15  **************************************************************************/
  16
  17 /* $Id$ */
  18
  19 /////////////////////////////////////////////////////////////
  20 //
  21 // Class to handle systematic errors for charm hadrons
  22 //
  23 // Usage:
  24 // AliHFSystEff syst;           // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
  25 // syst.SetRunNumber(YEAR);     // YEAR = two last numbers of the year (is 10 for 2010)
  26 // syst.SetCollisionType(TYPE);  // TYPE =  0 is pp, 1 is PbPb
  27 // syst.SetCentrality(CENT);     // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
  28 // syst.Init(DECAY);             // DECAY = 1 for D0, 2, for D+, 3 for D*
  29 // syst.DrawErrors(); // to see a plot of the error contributions
  30 // syst.GetTotalSystErr(pt); // to get the total err at pt
  31 //
  32 // Author: A.Dainese, andrea.dainese@pd.infn.it
  33 /////////////////////////////////////////////////////////////
  34
  35 #include <TStyle.h>
  36 #include <TGraphAsymmErrors.h>
  37 #include <TMath.h>
  38 #include <TCanvas.h>
  39 #include <TH2F.h>
  40 #include <TLegend.h>
  41 #include <TColor.h>
  42
  43 #include "AliLog.h"
  44 #include "AliHFSystErr.h"
  45
  46
  47 ClassImp(AliHFSystErr)
  48
  49 //--------------------------------------------------------------------------
  50 AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
  51 TNamed(name,title),
  52 fNorm(0),
  53 fRawYield(0),
  54 fTrackingEff(0),
  55 fBR(0),
  56 fCutsEff(0),
  57 fPIDEff(0),
  58 fMCPtShape(0),
  59 fPartAntipart(0),
  60 fRunNumber(10),
  61 fCollisionType(0),
  62 fCentralityClass("0100"),
  63 fRapidityRange("0101"),
  64 fIsLowEnergy(false),
  65 fIsCentScan(false),
  66 fIsRapidityScan(false)
  67 {
  68   //
  69   // Default Constructor
  70   //
  71 }
  72
  73 //--------------------------------------------------------------------------
  74 AliHFSystErr::~AliHFSystErr() {
  75   //
  76   // Default Destructor
  77   //
  78   /*
  79
  80   if(fNorm)         { delete fNorm; fNorm=0; }
  81   if(fRawYield)     { delete fRawYield; fRawYield=0; }
  82   if(fTrackingEff)  { delete fTrackingEff; fTrackingEff=0; }
  83   if(fBR)           { delete fBR; fBR=0; }
  84   if(fCutsEff)      { delete fCutsEff; fCutsEff=0; }
  85   if(fPIDEff)       { delete fPIDEff; fPIDEff=0; }
  86   if(fMCPtShape)    { delete fMCPtShape; fMCPtShape=0; }
  87   if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
  88   */
  89 }
  90
  91 //--------------------------------------------------------------------------
  92 void AliHFSystErr::Init(Int_t decay){
  93   //
  94   // Variables/histos initialization
  95   //
  96
  97   if ((fRunNumber>11) && fIsLowEnergy==false) {
  98     AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
  99   }
 100
 101   switch(decay) {
 102   case 1: // D0->Kpi
 103     if (fCollisionType==0) {
 104       if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
 105       else InitD0toKpi2010pp();
 106     }
 107     else if (fCollisionType==1) {
 108       if (fRunNumber == 10){
 109         if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
 110         else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
 111         else if (fCentralityClass=="020")  InitD0toKpi2010PbPb020();
 112         else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
 113         else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
 114         else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
 115         else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
 116         else AliFatal("Not yet implemented");
 117       }
 118       else if (fRunNumber == 11){
 119         if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
 120         else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
 121         else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
 122         else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
 123     else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
 124         else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
 125         else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
 126         else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
 127         else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
 128         else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
 129         else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
 130         else AliFatal("Not yet implemented");
 131       }
 132     }
 133     else if (fCollisionType==2) {
 134       if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
 135
 136       if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
 137       if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
 138       if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
 139       if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
 140
 141       if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
 142       if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
 143       if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
 144       if (fCentralityClass=="60100ZNA") InitD0toKpi2013pPb60100ZNA();
 145
 146       if (fIsRapidityScan) {
 147         if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
 148         if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
 149         if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
 150         if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
 151         if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
 152       }
 153     }
 154     else AliFatal("Not yet implemented");
 155     break;
 156
 157   case 2: // D+->Kpipi
 158     if (fCollisionType==0) {
 159       if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
 160       else InitDplustoKpipi2010pp();
 161     }
 162     else if (fCollisionType==1) {
 163       if (fRunNumber == 10){
 164         if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
 165         else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
 166         else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
 167         else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
 168         else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
 169         else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
 170         else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
 171         else AliFatal("Not yet implemented");
 172       }
 173       if(fRunNumber == 11){
 174         if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
 175         else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
 176         else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
 177         else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
 178         else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
 179         else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
 180         else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
 181         else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
 182     else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
 183         else AliFatal("Not yet implemented");
 184       }
 185     }
 186     else if (fCollisionType==2) {
 187       if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
 188       if (fIsRapidityScan) {
 189         if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
 190         if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
 191         if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
 192         if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
 193         if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
 194       }
 195     }
 196     else AliFatal("Not yet implemented");
 197     break;
 198   case 3: // D*->D0pi
 199     if (fCollisionType==0) {
 200       if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
 201       else InitDstartoD0pi2010pp();
 202     }
 203     else if (fCollisionType==1) {
 204       if (fRunNumber == 10){
 205         if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
 206         else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
 207         else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
 208         else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
 209         else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
 210         else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
 211         else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
 212         else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
 213         else AliFatal("Not yet implemented");
 214       }
 215       if (fRunNumber == 11){
 216         if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
 217         else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
 218         else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
 219         else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
 220         else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
 221         else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
 222         else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
 223         else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
 224     else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
 225         else AliFatal("Not yet implemented");
 226       }
 227     }
 228     else if (fCollisionType==2) {
 229       if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
 230
 231         if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
 232         if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
 233         if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
 234         if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
 235
 236         if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
 237         if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
 238         if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
 239         if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
 240
 241       if (fIsRapidityScan) {
 242         if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
 243         if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
 244         if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
 245         if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
 246         if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
 247       }
 248     }
 249     else AliFatal("Not yet implemented");
 250     break;
 251   case 4: // D+s->KKpi
 252     if (fCollisionType==0) InitDstoKKpi2010pp();
 253     else if (fCollisionType==1) {
 254       if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
 255       else AliFatal("Not yet implemented");
 256     }
 257     else if (fCollisionType==2) {
 258       if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
 259     }
 260     else AliFatal("Not yet implemented");
 261     break;
 262   case 5: // Lc->pKpi
 263     if (fCollisionType==0) InitLctopKpi2010pp();
 264     else AliFatal("Not yet implemented");
 265     break;
 266   case 6: // Lc->pK0S
 267     if (fCollisionType==0) InitLctopK0S2010pp();
 268     else AliFatal("Not yet implemented");
 269     break;
 270
 271   default:
 272     printf("Invalid decay type: %d\n",decay);
 273     break;
 274   }
 275
 276 }
 277
 278 //--------------------------------------------------------------------------
 279 void AliHFSystErr::InitD0toKpi2010pp() {
 280   //
 281   // D0->Kpi syst errors. Responsible: A. Rossi
 282   //   2010 pp sample
 283   //
 284
 285   // Normalization
 286   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 287   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 288
 289   // Branching ratio
 290   fBR = new TH1F("fBR","fBR",24,0,24);
 291   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 292
 293   // Tracking efficiency
 294   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 295   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
 296
 297   // Raw yield extraction
 298   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 299   fRawYield->SetBinContent(1,1);
 300   fRawYield->SetBinContent(2,0.22);
 301   fRawYield->SetBinContent(3,0.1);
 302   for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
 303   for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
 304   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 305   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
 306
 307   // Cuts efficiency (from cuts variation)
 308   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 309   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 310
 311   // PID efficiency (from PID/noPID)
 312   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 313   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
 314   fPIDEff->SetBinContent(2,0.05); // 5%
 315
 316   // MC dN/dpt
 317   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 318   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 319   fMCPtShape->SetBinContent(1,0.03);
 320   fMCPtShape->SetBinContent(2,0.03);
 321
 322   // particle-antiparticle
 323   //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
 324   //  fPartAntipart->SetBinContent(1,1);
 325   //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
 326
 327   return;
 328 }
 329 //--------------------------------------------------------------------------
 330 void AliHFSystErr::InitD0toKpi2010PbPb020() {
 331   //
 332   // D0->Kpi syst errors. Responsible: A. Rossi
 333   //   2010 PbPb sample, 0-20 CC
 334   //
 335
 336   // Normalization
 337   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 338   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 339
 340   // Branching ratio
 341   fBR = new TH1F("fBR","fBR",20,0,20);
 342   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 343
 344   // Tracking efficiency
 345   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 346   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 347
 348   // Raw yield extraction
 349   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 350   fRawYield->SetBinContent(1,0);
 351   fRawYield->SetBinContent(2,0);
 352   fRawYield->SetBinContent(3,0.08);
 353   for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
 354   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 355
 356   // Cuts efficiency (from cuts variation)
 357   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 358   fCutsEff->SetBinContent(1,0.);
 359   fCutsEff->SetBinContent(2,0.);
 360   fCutsEff->SetBinContent(3,0.13);
 361   fCutsEff->SetBinContent(4,0.11);
 362   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 363   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 364
 365   // PID efficiency (from PID/noPID)
 366   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 367   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 368
 369   // MC dN/dpt
 370   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 371   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 372   fMCPtShape->SetBinContent(3,0.04);
 373   fMCPtShape->SetBinContent(4,0.02);
 374   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 375
 376 //   // particle-antiparticle
 377 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 378 //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 379 //   fPartAntipart->SetBinContent(3,0.10);
 380 //   fPartAntipart->SetBinContent(4,0.10);
 381 //   fPartAntipart->SetBinContent(7,0.10);
 382 //   fPartAntipart->SetBinContent(8,0.10);
 383
 384   return;
 385 }
 386 //--------------------------------------------------------------------------
 387 //--------------------------------------------------------------------------
 388 void AliHFSystErr::InitD0toKpi2011PbPb07half() {
 389   //
 390   // D0->Kpi syst errors. Responsible: A. Rossi
 391   //   2011 PbPb sample, 0-7.5 CC
 392   //
 393
 394   // Normalization
 395   fNorm = new TH1F("fNorm","fNorm",36,0,36);
 396   for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
 397
 398   // Branching ratio
 399   fBR = new TH1F("fBR","fBR",36,0,36);
 400   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 401
 402   // Tracking efficiency
 403   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
 404   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 405   for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 406
 407   // Raw yield extraction
 408   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 409   fRawYield->SetBinContent(1,0);
 410   fRawYield->SetBinContent(2,0.1);
 411   fRawYield->SetBinContent(3,0.05);
 412   for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 413   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 414   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
 415   for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 416
 417   // Cuts efficiency (from cuts variation)
 418   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 419   fCutsEff->SetBinContent(1,0.);
 420   fCutsEff->SetBinContent(2,0.15);
 421   fCutsEff->SetBinContent(3,0.13);
 422   fCutsEff->SetBinContent(4,0.11);
 423   fCutsEff->SetBinContent(5,0.08);
 424   for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
 425   for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 426
 427   // PID efficiency (from PID/noPID)
 428   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
 429   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
 430   for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 431
 432   // MC dN/dpt
 433   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
 434   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
 435   fMCPtShape->SetBinContent(2,0.06);
 436   fMCPtShape->SetBinContent(3,0.04);
 437   fMCPtShape->SetBinContent(4,0.02);
 438   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 439   for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
 440   for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 441
 442 //   // particle-antiparticle
 443 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
 444 //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 445 //   fPartAntipart->SetBinContent(3,0.10);
 446 //   fPartAntipart->SetBinContent(4,0.10);
 447 //   fPartAntipart->SetBinContent(7,0.10);
 448 //   fPartAntipart->SetBinContent(8,0.10);
 449
 450   return;
 451 }
 452
 453 //--------------------------------------------------------------------------
 454 void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
 455   //
 456   // D0->Kpi syst errors. Responsible: D. Caffarri
 457   //   2011 PbPb sample, 30-50 CC InPlane
 458   //
 459   InitD0toKpi2011PbPb07half();
 460   // Raw yield extraction
 461   // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 462   fRawYield->SetBinContent(1,0);
 463   fRawYield->SetBinContent(2,0);
 464   fRawYield->SetBinContent(3,0.05);
 465   fRawYield->SetBinContent(4,0.08);
 466   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
 467   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 468   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
 469   for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 470
 471   // Cuts efficiency (from cuts variation)
 472   //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 473   fCutsEff->SetBinContent(1,0.);
 474   fCutsEff->SetBinContent(2,0.0);
 475   fCutsEff->SetBinContent(3,0.10);
 476   fCutsEff->SetBinContent(4,0.10);
 477   fCutsEff->SetBinContent(5,0.10);
 478   fCutsEff->SetBinContent(6,0.10);
 479   for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
 480   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
 481   for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 482 }
 483
 484 //--------------------------------------------------------------------------
 485 void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
 486   //
 487   // D0->Kpi syst errors. Responsible: D. Caffarri
 488   //   2011 PbPb sample, 30-50 CC OutOfPlane
 489   //
 490   InitD0toKpi2011PbPb07half();
 491   // Raw yield extraction
 492   //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 493   fRawYield->SetBinContent(1,0);
 494   fRawYield->SetBinContent(2,0.);
 495   fRawYield->SetBinContent(3,0.05);
 496   for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
 497   for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
 498   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
 499   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
 500   for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 501
 502   // Cuts efficiency (from cuts variation)
 503   // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 504   fCutsEff->SetBinContent(1,0.);
 505   fCutsEff->SetBinContent(2,0.);
 506   fCutsEff->SetBinContent(3,0.10);
 507   fCutsEff->SetBinContent(4,0.10);
 508   fCutsEff->SetBinContent(5,0.10);
 509   fCutsEff->SetBinContent(6,0.10);
 510   for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
 511   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
 512   for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 513 }
 514
 515 //--------------------------------------------------------------------------
 516 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
 517   //
 518   // D0->Kpi syst errors. Responsible: A. Rossi
 519   //   2010 PbPb sample, 40-80 CC
 520   //
 521
 522   // Normalization
 523   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 524   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 525
 526   // Branching ratio
 527   fBR = new TH1F("fBR","fBR",20,0,20);
 528   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 529
 530   // Tracking efficiency
 531   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 532   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 533
 534
 535   // Raw yield extraction
 536   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 537   fRawYield->SetBinContent(1,0);
 538   fRawYield->SetBinContent(2,0);
 539   for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
 540   //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
 541
 542   // Cuts efficiency (from cuts variation)
 543   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 544   fCutsEff->SetBinContent(1,0.);
 545   fCutsEff->SetBinContent(2,0.);
 546   fCutsEff->SetBinContent(3,0.13);
 547   fCutsEff->SetBinContent(4,0.11);
 548   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 549   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 550
 551   // PID efficiency (from PID/noPID)
 552   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 553 //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
 554 //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 555   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 556
 557   // MC dN/dpt
 558   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 559   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 560   //  fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
 561   //  fMCPtShape->SetBinContent(4,0.02);
 562   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 563
 564   //   // particle-antiparticle
 565   //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 566   //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 567
 568   return;
 569 }
 570
 571 //--------------------------------------------------------------------------
 572 void AliHFSystErr::InitD0toKpi2011PbPb010() {
 573     //
 574     // D0->Kpi syst errors. Responsible: A. Rossi
 575     //   2011 PbPb sample, 0-10 CC
 576     //
 577
 578   AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
 579
 580     // Normalization
 581     fNorm = new TH1F("fNorm","fNorm",24,0,24);
 582     for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
 583
 584     // Branching ratio
 585     fBR = new TH1F("fBR","fBR",24,0,24);
 586     for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 587
 588     // Tracking efficiency
 589     fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 590     for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 591
 592
 593     // Raw yield extraction
 594     fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 595     fRawYield->SetBinContent(1,0);
 596     fRawYield->SetBinContent(2,0.15);
 597     fRawYield->SetBinContent(3,0.09);
 598     fRawYield->SetBinContent(4,0.06);
 599     fRawYield->SetBinContent(5,0.06);
 600     fRawYield->SetBinContent(6,0.06);
 601     for (Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.06);
 602     for (Int_t i=9; i<=12; i++)fRawYield->SetBinContent(i,0.07);
 603     for (Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.11);
 604     for (Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.30);
 605
 606     fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 607     fCutsEff->SetBinContent(1,0.);
 608     fCutsEff->SetBinContent(2,0.15);
 609     fCutsEff->SetBinContent(3,0.12);
 610     fCutsEff->SetBinContent(4,0.05);
 611     fCutsEff->SetBinContent(5,0.05);
 612     fCutsEff->SetBinContent(6,0.05);
 613     for (Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i,0.05);
 614     for (Int_t i=9; i<=12; i++)fCutsEff->SetBinContent(i,0.05);
 615     for (Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i,0.05);
 616     for (Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i,0.05);
 617
 618
 619     // PID efficiency (from PID/noPID)
 620     fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 621     //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
 622     //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 623     for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 624
 625     // MC dN/dpt
 626     fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 627     fMCPtShape->SetBinContent(1,0.);
 628     fMCPtShape->SetBinContent(2,0.15);
 629     fMCPtShape->SetBinContent(3,0.12);
 630     fMCPtShape->SetBinContent(4,0.02);
 631     for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
 632
 633     //   // particle-antiparticle
 634     //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 635     //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 636
 637     return;
 638 }
 639
 640 //--------------------------------------------------------------------------
 641 void AliHFSystErr::InitD0toKpi2011PbPb3050() {
 642     //
 643     // D0->Kpi syst errors. Responsible: D. Caffarri
 644     //   2011 PbPb sample, 30-50 CC
 645     //
 646
 647     // Normalization
 648     fNorm = new TH1F("fNorm","fNorm",20,0,20);
 649     for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 650
 651     // Branching ratio
 652     fBR = new TH1F("fBR","fBR",20,0,20);
 653     for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 654
 655     // Tracking efficiency
 656     fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 657     for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 658
 659
 660     // Raw yield extraction
 661     fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 662     fRawYield->SetBinContent(1,0);
 663     fRawYield->SetBinContent(2,0.102);
 664     fRawYield->SetBinContent(3,0.071);
 665     fRawYield->SetBinContent(4,0.052);
 666     fRawYield->SetBinContent(5,0.063);
 667     fRawYield->SetBinContent(6,0.071);
 668     fRawYield->SetBinContent(7,0.071);
 669     fRawYield->SetBinContent(8,0.071);
 670     fRawYield->SetBinContent(9,0.071);
 671     fRawYield->SetBinContent(10,0.071);
 672     fRawYield->SetBinContent(11,0.071);
 673     fRawYield->SetBinContent(12,0.071);
 674
 675     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
 676     for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
 677
 678     fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 679     fCutsEff->SetBinContent(1,0.);
 680     fCutsEff->SetBinContent(2,0.10);
 681     fCutsEff->SetBinContent(3,0.10);
 682     fCutsEff->SetBinContent(4,0.10);
 683     fCutsEff->SetBinContent(5,0.10);
 684     fCutsEff->SetBinContent(6,0.10);
 685     for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
 686     for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
 687
 688
 689     // PID efficiency (from PID/noPID)
 690     fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 691     //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
 692     //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 693     for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 694
 695     // MC dN/dpt
 696     fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
 697     for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 698     fMCPtShape->SetBinContent(2,0.06);
 699     fMCPtShape->SetBinContent(3,0.04);
 700     fMCPtShape->SetBinContent(4,0.02);
 701     for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 702     for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
 703     //   // particle-antiparticle
 704     //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 705     //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 706
 707     return;
 708 }
 709
 710
 711 //--------------------------------------------------------------------------
 712 void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 713   //
 714   // D0->Kpi syst errors. Low energy run
 715   //   2011 2.76 TeV pp sample
 716   //
 717   AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
 718
 719   // Normalization
 720   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 721   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
 722
 723   // Branching ratio
 724   fBR = new TH1F("fBR","fBR",20,0,20);
 725   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 726
 727   // Tracking efficiency
 728   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 729   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
 730
 731   // Raw yield extraction
 732   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 733   fRawYield->SetBinContent(1,1);
 734   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
 735
 736   // Cuts efficiency (from cuts variation)
 737   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 738   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 739   fCutsEff->SetBinContent(2,0.20);
 740   for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 741
 742
 743   // PID efficiency (from PID/noPID)
 744   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 745   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
 746   //  fPIDEff->SetBinContent(2,0.20);
 747   for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
 748
 749   // MC dN/dpt
 750   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 751   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 752   fMCPtShape->SetBinContent(1,0.03);
 753   fMCPtShape->SetBinContent(2,0.03);
 754
 755 //   // particle-antiparticle
 756 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 757 //   fPartAntipart->SetBinContent(1,1);
 758 //   fPartAntipart->SetBinContent(2,1);
 759 //   for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
 760 //   for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
 761
 762   return;
 763 }
 764 //_________________________________________________________________________
 765 void AliHFSystErr::InitD0toKpi2013pPb0100(){
 766   //
 767   // D0->Kpi syst errors. p-Pb data sample
 768   //
 769   AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
 770
 771   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 772   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
 773
 774   // Branching ratio
 775   fBR = new TH1F("fBR","fBR",24,0,24);
 776   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010
 777
 778   // Tracking efficiency
 779   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 780   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
 781
 782   // Raw yield extraction
 783   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 784   fRawYield->SetBinContent(1,0.0);
 785   fRawYield->SetBinContent(2,0.08);
 786   for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
 787   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
 788   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
 789
 790   // Cuts efficiency (from cuts variation)
 791   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 792   fCutsEff->SetBinContent(1,0.0);
 793   fCutsEff->SetBinContent(2,0.08);
 794   for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
 795
 796   // PID efficiency (from PID/noPID)
 797   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 798   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
 799
 800   // MC dN/dpt
 801   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 802   fMCPtShape->SetBinContent(1,0.0);
 803   fMCPtShape->SetBinContent(2,0.02);
 804   for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
 805
 806   // particle-antiparticle
 807   /*
 808   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 809   for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 810   */
 811   return;
 812
 813
 814 }
 815
 816
 817 //--------------------------------------------------------------------------
 818 void AliHFSystErr::InitDplustoKpipi2010pp() {
 819   //
 820   // D+->Kpipi syst errors. Responsible: R. Bala
 821   //  2010 pp sample
 822   //
 823
 824
 825 // Normalization
 826   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 827   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 828
 829   // Branching ratio
 830   fBR = new TH1F("fBR","fBR",24,0,24);
 831   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 832
 833   // Tracking efficiency
 834   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 835   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 836
 837
 838   // Raw yield extraction
 839   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 840   fRawYield->SetBinContent(1,1);
 841   fRawYield->SetBinContent(2,0.25);
 842   fRawYield->SetBinContent(3,0.25);
 843   fRawYield->SetBinContent(4,0.20);
 844   fRawYield->SetBinContent(5,0.09);
 845   fRawYield->SetBinContent(6,0.09);
 846   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 847   for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
 848
 849   // Cuts efficiency (from cuts variation)
 850   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 851   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 852
 853   // PID efficiency (from PID/noPID)
 854   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 855   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 856   fPIDEff->SetBinContent(1,0.15); // 15%
 857   fPIDEff->SetBinContent(2,0.15); // 15%
 858   fPIDEff->SetBinContent(3,0.15); // 15%
 859   fPIDEff->SetBinContent(4,0.15); // 15%
 860   for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
 861
 862   // MC dN/dpt
 863   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 864   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 865   fMCPtShape->SetBinContent(1,0.03);
 866   fMCPtShape->SetBinContent(2,0.03);
 867
 868
 869   // particle-antiparticle
 870   /*
 871   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 872   fPartAntipart->SetBinContent(1,1);
 873   fPartAntipart->SetBinContent(2,1);
 874   fPartAntipart->SetBinContent(3,0.12);
 875   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 876   */
 877   return;
 878 }
 879
 880 //--------------------------------------------------------------------------
 881 void AliHFSystErr::InitDstoKKpi2010pp() {
 882   //
 883   // D+s->KKpi syst errors. Responsible: G.M. Innocenti
 884   //  2010 pp sample
 885   //
 886
 887
 888 // Normalization
 889   fNorm = new TH1F("fNorm","fNorm",12,0,12);
 890   for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
 891
 892   // Branching ratio
 893   fBR = new TH1F("fBR","fBR",12,0,12);
 894   for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
 895
 896   // Tracking efficiency
 897   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
 898   for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 899
 900
 901   // Raw yield extraction
 902   fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
 903   fRawYield->SetBinContent(1,1);
 904   fRawYield->SetBinContent(2,1);
 905   fRawYield->SetBinContent(3,0.20);
 906   fRawYield->SetBinContent(4,0.20);
 907   fRawYield->SetBinContent(5,0.15);
 908   fRawYield->SetBinContent(6,0.15);
 909   fRawYield->SetBinContent(7,0.15);
 910   fRawYield->SetBinContent(8,0.15);
 911   fRawYield->SetBinContent(9,0.20);
 912   fRawYield->SetBinContent(10,0.20);
 913   fRawYield->SetBinContent(11,0.20);
 914   fRawYield->SetBinContent(12,0.20);
 915
 916   // Cuts efficiency (from cuts variation)
 917   fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
 918   for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
 919
 920   // PID efficiency (from PID/noPID)
 921   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
 922   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
 923
 924   // MC dN/dpt
 925   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
 926   for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
 927   for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
 928   for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
 929   for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
 930   for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
 931
 932
 933   // particle-antiparticle
 934   /*
 935   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
 936   fPartAntipart->SetBinContent(1,1);
 937   fPartAntipart->SetBinContent(2,1);
 938   fPartAntipart->SetBinContent(3,0.12);
 939   for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 940   */
 941   return;
 942 }
 943
 944
 945 //--------------------------------------------------------------------------
 946 void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
 947   //
 948   // D+->Kpipi syst errors. Responsible: E. Bruna
 949   //  2011 PbPb sample, 0-7.5% CC
 950   //
 951
 952  // Normalization
 953   fNorm = new TH1F("fNorm","fNorm",36,0,36);
 954   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
 955
 956   // Branching ratio
 957   fBR = new TH1F("fBR","fBR",36,0,36);
 958   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 959
 960   // Tracking efficiency
 961   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
 962   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 963
 964   // Raw yield extraction
 965   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 966   for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
 967   fRawYield->SetBinContent(4,0.30);
 968   fRawYield->SetBinContent(5,0.20);
 969   fRawYield->SetBinContent(6,0.20);
 970   for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
 971   for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
 972   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
 973   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
 974   for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
 975
 976
 977   // Cuts efficiency (from cuts variation)
 978   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 979   for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
 980   for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
 981
 982
 983   // PID efficiency (from PID/noPID)
 984   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
 985   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 986
 987   // MC dN/dpt  (24/7/2012)
 988   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
 989   for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 990   for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
 991   //  for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 992
 993
 994   // particle-antiparticle
 995   /*
 996   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 997   fPartAntipart->SetBinContent(1,1);
 998   fPartAntipart->SetBinContent(2,1);
 999   fPartAntipart->SetBinContent(3,0.12);
1000   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1001   */
1002
1003   return;
1004 }
1005 //_________________________________________________________________
1006 void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
1007   //
1008   // D+->Kpipi syst errors. Responsible: E. Bruna, C. Bedda
1009   //  2011 PbPb sample, 0-10% CC
1010   //
1011
1012  // Normalization
1013   fNorm = new TH1F("fNorm","fNorm",36,0,36);
1014   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1015
1016   // Branching ratio
1017   fBR = new TH1F("fBR","fBR",36,0,36);
1018   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1019
1020   // Tracking efficiency
1021   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1022   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1023
1024   // Raw yield extraction
1025   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1026   for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i,0.);  //5 to 10%
1027   fRawYield->SetBinContent(4,0.10);
1028   fRawYield->SetBinContent(5,0.08);
1029   fRawYield->SetBinContent(6,0.08);
1030   for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.08);
1031   for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1032   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.08);
1033   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
1034   for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.08);
1035
1036
1037   // Cuts efficiency (from cuts variation)
1038   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1039   for(Int_t i=1; i<=36; i++) fCutsEff->SetBinContent(i,0.10);
1040
1041
1042   // PID efficiency (from PID/noPID)
1043   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1044   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1045
1046   // MC dN/dpt
1047   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1048   for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.);
1049   fMCPtShape->SetBinContent(4,0.06);
1050   fMCPtShape->SetBinContent(5,0.02);
1051   for (Int_t i=6; i<=36; i++) fMCPtShape->SetBinContent(i,0.01);
1052   //  for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1053
1054
1055   // particle-antiparticle
1056   /*
1057   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1058   fPartAntipart->SetBinContent(1,1);
1059   fPartAntipart->SetBinContent(2,1);
1060   fPartAntipart->SetBinContent(3,0.12);
1061   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1062   */
1063
1064   return;
1065 }
1066
1067 //___________________________________________________________________
1068 void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
1069     //
1070     // D+->Kpipi syst errors. Responsible: R. Bala
1071     //  2011 PbPb sample, 30-50% CC
1072     //
1073     AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1074     // Normalization
1075     fNorm = new TH1F("fNorm","fNorm",16,0,16);
1076     for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1077
1078     // Branching ratio
1079     fBR = new TH1F("fBR","fBR",16,0,16);
1080     for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1081
1082     // Tracking efficiency
1083     fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
1084     for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1085
1086     // Raw yield extraction
1087     fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
1088     for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09);  //5 to 10%
1089     for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
1090
1091
1092     // Cuts efficiency (from cuts variation)
1093     fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
1094     for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
1095     for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
1096
1097     // PID efficiency (from PID/noPID)
1098     fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
1099     for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1100
1101     // MC dN/dpt  (24/7/2012)
1102     fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
1103     fMCPtShape->SetBinContent(3,0.10);
1104     fMCPtShape->SetBinContent(4,0.04);
1105     for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
1106
1107
1108     // particle-antiparticle
1109     /*
1110      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1111      fPartAntipart->SetBinContent(1,1);
1112      fPartAntipart->SetBinContent(2,1);
1113      fPartAntipart->SetBinContent(3,0.12);
1114      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1115      */
1116
1117     return;
1118 }
1119
1120
1121 //--------------------------------------------------------------------------
1122 void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
1123   //
1124   // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
1125   //  2011 PbPb sample, 0-7.5% CC
1126   //
1127
1128  // Normalization
1129   fNorm = new TH1F("fNorm","fNorm",12,0,12);
1130   for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1131
1132   // Branching ratio
1133   fBR = new TH1F("fBR","fBR",12,0,12);
1134   for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1135
1136   // Tracking efficiency
1137   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1138   for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1139
1140   // Raw yield extraction
1141   fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1142   for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
1143   for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
1144
1145
1146
1147   // Cuts efficiency (from cuts variation)
1148   fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1149   for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1150
1151   // PID efficiency (from PID/noPID)
1152   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1153   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
1154
1155    // MC dN/dpt
1156   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1157   for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
1158   for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
1159   for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
1160   for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
1161   for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
1162
1163   // particle-antiparticle
1164   /*
1165   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
1166   fPartAntipart->SetBinContent(1,1);
1167   fPartAntipart->SetBinContent(2,1);
1168   fPartAntipart->SetBinContent(3,0.12);
1169   for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1170   */
1171
1172   return;
1173 }
1174
1175 //_________________________________________________________________________
1176 void AliHFSystErr::InitDstoKKpi2013pPb0100(){
1177   //
1178   // Dstar syst
1179   AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1180
1181   //Normalization
1182   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1183   for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
1184
1185  // Branching ratio
1186   fBR = new TH1F("fBR","fBR",12,0,12);
1187   for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1188
1189   // Tracking efficiency
1190   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1191   for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1192
1193   // Raw yield extraction
1194   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1195   for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
1196   for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
1197   for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
1198   for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1199
1200   // Cuts efficiency (from cuts variation)
1201   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1202   for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
1203
1204   // PID efficiency (from PID/noPID)
1205   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1206   for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
1207   for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
1208
1209   // MC dN/dpt
1210   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1211   for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
1212
1213   // particle-antiparticle
1214   /*
1215   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1216   fPartAntipart->SetBinContent(1,1);
1217   fPartAntipart->SetBinContent(2,1);
1218   fPartAntipart->SetBinContent(3,0.12);
1219   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1220   */
1221   return;
1222 }
1223
1224
1225 //--------------------------------------------------------------------------
1226 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
1227   //
1228   // D+->Kpipi syst errors. Responsible: ??
1229   //  2010 PbPb sample, 0-20 CC
1230   //
1231
1232  // Normalization
1233   fNorm = new TH1F("fNorm","fNorm",20,0,20);
1234   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1235
1236   // Branching ratio
1237   fBR = new TH1F("fBR","fBR",20,0,20);
1238   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1239
1240   // Tracking efficiency
1241   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1242   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1243
1244   // Raw yield extraction
1245   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1246   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
1247   // fRawYield->SetBinContent(5,0.23);
1248   //fRawYield->SetBinContent(6,0.23);
1249   fRawYield->SetBinContent(7,0.20);
1250   fRawYield->SetBinContent(8,0.20);
1251   fRawYield->SetBinContent(9,0.15);
1252   fRawYield->SetBinContent(10,0.15);
1253   fRawYield->SetBinContent(11,0.15);
1254   fRawYield->SetBinContent(12,0.15);
1255
1256   // Cuts efficiency (from cuts variation)
1257   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1258   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1259
1260   // PID efficiency (from PID/noPID)
1261   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1262   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1263
1264   // MC dN/dpt  (2/2/2012)
1265   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1266   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
1267   for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1268   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1269   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1270
1271
1272   // particle-antiparticle
1273   /*
1274   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1275   fPartAntipart->SetBinContent(1,1);
1276   fPartAntipart->SetBinContent(2,1);
1277   fPartAntipart->SetBinContent(3,0.12);
1278   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1279   */
1280
1281   return;
1282 }
1283
1284 //--------------------------------------------------------------------------
1285 void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
1286   //
1287   // D+->Kpipi syst errors. Responsible: ??
1288   //  2010 PbPb sample, 40-80 CC
1289   //
1290
1291
1292  // Normalization
1293   fNorm = new TH1F("fNorm","fNorm",20,0,20);
1294   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1295
1296   // Branching ratio
1297   fBR = new TH1F("fBR","fBR",20,0,20);
1298   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
1299
1300   // Tracking efficiency
1301   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1302   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1303
1304
1305   // Raw yield extraction
1306   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1307   fRawYield->SetBinContent(1,1);
1308   fRawYield->SetBinContent(2,1);
1309   fRawYield->SetBinContent(3,1);
1310   fRawYield->SetBinContent(4,0.15);
1311   fRawYield->SetBinContent(5,0.05);
1312   fRawYield->SetBinContent(6,0.05);
1313   fRawYield->SetBinContent(7,0.15);
1314   fRawYield->SetBinContent(8,0.15);
1315   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
1316   for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1);  //5 to 10%
1317
1318   // Cuts efficiency (from cuts variation)
1319   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1320   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1321
1322   // PID efficiency (from PID/noPID)
1323   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1324   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1325   fPIDEff->SetBinContent(3,0.13); // 13%
1326
1327
1328   // MC dN/dpt  (2/2/2012)
1329   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1330   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
1331   for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1332   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1333   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1334
1335
1336   // particle-antiparticle
1337   /*
1338   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1339   fPartAntipart->SetBinContent(1,1);
1340   fPartAntipart->SetBinContent(2,1);
1341   fPartAntipart->SetBinContent(3,0.12);
1342   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1343   */
1344   return;
1345 }
1346
1347 //--------------------------------------------------------------------------
1348 void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
1349
1350   //
1351   // D+->Kpipi syst errors. Responsible: R. Bala
1352   //  2011 2.76 TeV pp sample
1353   //
1354   AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1355
1356   // Normalization
1357   fNorm = new TH1F("fNorm","fNorm",20,0,20);
1358   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
1359
1360   // Branching ratio
1361   fBR = new TH1F("fBR","fBR",20,0,20);
1362   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1363
1364   // Tracking efficiency
1365   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1366   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
1367
1368   // Raw yield extraction
1369   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1370   fRawYield->SetBinContent(1,1);
1371   fRawYield->SetBinContent(2,1);
1372   for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10);  //5 to 10%
1373   fRawYield->SetBinContent(7,0.15);
1374   fRawYield->SetBinContent(8,0.15);
1375   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055);  //5 to 10%
1376
1377   // Cuts efficiency (from cuts variation)
1378   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1379   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1380
1381   // PID efficiency (from PID/noPID)
1382   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1383   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1384   fPIDEff->SetBinContent(3,0.10); // 13%
1385   fPIDEff->SetBinContent(4,0.10); // 13%
1386
1387   // MC dN/dpt  (copied from D0 : will update later)
1388   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1389   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1390   fMCPtShape->SetBinContent(1,0.03);
1391   fMCPtShape->SetBinContent(2,0.03);
1392
1393   return;
1394 }
1395
1396 //_________________________________________________________________________
1397 void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
1398   //
1399   // Dplus syst
1400   AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1401
1402   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1403   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
1404
1405   fBR = new TH1F("fBR","fBR",20,0,24);
1406   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1407
1408   // Tracking efficiency
1409   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1410   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1411
1412   // Raw yield extraction
1413   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1414   fRawYield->SetBinContent(1,0.);
1415   fRawYield->SetBinContent(2,0.1);
1416   fRawYield->SetBinContent(3,0.08);
1417   fRawYield->SetBinContent(4,0.05);
1418   fRawYield->SetBinContent(5,0.05);
1419   fRawYield->SetBinContent(6,0.05);
1420   fRawYield->SetBinContent(7,0.05);
1421   fRawYield->SetBinContent(8,0.05);
1422   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1423   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1424   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
1425
1426   // Cuts efficiency (from cuts variation)
1427   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1428   fCutsEff->SetBinContent(1,0.);
1429   fCutsEff->SetBinContent(2,0.1);
1430   fCutsEff->SetBinContent(3,0.1);
1431   fCutsEff->SetBinContent(4,0.06);
1432   fCutsEff->SetBinContent(5,0.06);
1433   fCutsEff->SetBinContent(6,0.06);
1434   fCutsEff->SetBinContent(7,0.05);
1435   fCutsEff->SetBinContent(8,0.05);
1436   for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
1437   for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1438   for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1439
1440   // PID efficiency (from PID/noPID)
1441   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1442   fPIDEff->SetBinContent(1,1);
1443   fPIDEff->SetBinContent(2,0.1);
1444   for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
1445
1446   // MC dN/dpt
1447   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1448   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
1449
1450   // particle-antiparticle
1451   /*
1452   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1453   fPartAntipart->SetBinContent(1,1);
1454   fPartAntipart->SetBinContent(2,1);
1455   fPartAntipart->SetBinContent(3,0.12);
1456   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1457   */
1458   return;
1459
1460
1461 }
1462
1463 //--------------------------------------------------------------------------
1464 void AliHFSystErr::InitDstartoD0pi2010pp() {
1465   //
1466   // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
1467   //  2010 pp sample
1468   //
1469
1470  // Normalization
1471   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1472   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1473
1474   // Branching ratio
1475   fBR = new TH1F("fBR","fBR",24,0,24);
1476   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1477
1478   // Tracking efficiency
1479   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1480   fTrackingEff->SetBinContent(1,1.0);
1481   fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
1482   fTrackingEff->SetBinContent(3,0.12);
1483   fTrackingEff->SetBinContent(3,0.12);
1484   for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
1485
1486
1487   // Raw yield extraction
1488   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1489   fRawYield->SetBinContent(1,1.0);
1490   fRawYield->SetBinContent(2,0.10);
1491   fRawYield->SetBinContent(3,0.04);
1492   fRawYield->SetBinContent(4,0.03);
1493   fRawYield->SetBinContent(5,0.03);
1494   fRawYield->SetBinContent(6,0.05);
1495   fRawYield->SetBinContent(7,0.05);
1496   fRawYield->SetBinContent(8,0.05);
1497   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);  //4%
1498   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);  //4%
1499   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
1500
1501   // Cuts efficiency (from cuts variation)
1502   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1503   fCutsEff->SetBinContent(2,0.22);
1504   for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1505
1506   // PID efficiency (from PID/noPID)
1507   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1508   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1509
1510
1511   // MC dN/dpt  (copied from D0 : will update later)
1512   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1513   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1514   fMCPtShape->SetBinContent(1,0.03);
1515   fMCPtShape->SetBinContent(2,0.03);
1516
1517   return;
1518
1519
1520 }
1521 //--------------------------------------------------------------------------
1522 void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
1523
1524   //
1525   // D+->Kpipi syst errors. Responsible: A. Grelli
1526   //  2011 2.76 TeV pp sample
1527   //
1528   AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1529
1530 // Normalization
1531   fNorm = new TH1F("fNorm","fNorm",20,0,20);
1532   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
1533
1534   // Branching ratio
1535   fBR = new TH1F("fBR","fBR",20,0,20);
1536   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1537
1538   // Tracking efficiency
1539   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1540   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
1541
1542   // Raw yield extraction
1543   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1544   fRawYield->SetBinContent(1,1);
1545   fRawYield->SetBinContent(2,1);
1546   fRawYield->SetBinContent(3,0.14);
1547   fRawYield->SetBinContent(4,0.14);
1548   fRawYield->SetBinContent(5,0.12);
1549   fRawYield->SetBinContent(6,0.12);
1550   fRawYield->SetBinContent(7,0.06);
1551   fRawYield->SetBinContent(8,0.06);
1552   fRawYield->SetBinContent(9,0.08);
1553   fRawYield->SetBinContent(10,0.08);
1554   fRawYield->SetBinContent(11,0.08);
1555   fRawYield->SetBinContent(12,0.08);
1556   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
1557
1558   // Cuts efficiency (from cuts variation)
1559   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1560   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
1561   fCutsEff->SetBinContent(3,0.15);
1562   fCutsEff->SetBinContent(4,0.15);
1563   fCutsEff->SetBinContent(5,0.15);
1564   fCutsEff->SetBinContent(6,0.15);
1565   fCutsEff->SetBinContent(7,0.10);
1566   fCutsEff->SetBinContent(8,0.10);
1567   fCutsEff->SetBinContent(9,0.10);
1568   fCutsEff->SetBinContent(10,0.10);
1569   fCutsEff->SetBinContent(11,0.10);
1570   fCutsEff->SetBinContent(12,0.10);
1571
1572   // PID efficiency (from PID/noPID)
1573   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1574   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
1575
1576   // MC dN/dpt
1577   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1578   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1579   fMCPtShape->SetBinContent(1,0.03);
1580   fMCPtShape->SetBinContent(2,0.03);
1581
1582   return;
1583 }
1584
1585 //------------------------------------------------------------------------
1586 void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
1587   //
1588   // D*+->D0pi syst errors. Responsible: A. Grelli
1589   //  2010 PbPb sample, 0-20 CC
1590   //
1591
1592   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1593
1594  // Normalization
1595   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1596   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1597
1598   // Branching ratio
1599   fBR = new TH1F("fBR","fBR",24,0,24);
1600   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1601
1602   // Tracking efficiency
1603   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1604   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1605
1606
1607   // Raw yield extraction
1608   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1609   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
1610   fRawYield->SetBinContent(3,0.2);
1611   fRawYield->SetBinContent(4,0.2);
1612   fRawYield->SetBinContent(5,0.2);
1613   fRawYield->SetBinContent(6,0.2);
1614
1615   // Cuts efficiency (from cuts variation)
1616   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1617   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1618   fCutsEff->SetBinContent(4,0.15);
1619   fCutsEff->SetBinContent(5,0.15);
1620   fCutsEff->SetBinContent(6,0.15);
1621
1622   // PID efficiency (from PID/noPID)
1623   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1624   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1625
1626
1627   // MC dN/dpt  (from study on D* pt shape)
1628   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1629   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
1630   fMCPtShape->SetBinContent(4,0.025);
1631   fMCPtShape->SetBinContent(5,0.025);
1632   fMCPtShape->SetBinContent(6,0.025);
1633   fMCPtShape->SetBinContent(7,0.04);
1634   fMCPtShape->SetBinContent(8,0.04);
1635   fMCPtShape->SetBinContent(9,0.03);
1636   fMCPtShape->SetBinContent(10,0.03);
1637   fMCPtShape->SetBinContent(11,0.03);
1638   fMCPtShape->SetBinContent(12,0.03);
1639
1640   return;
1641
1642 }
1643 // ----------------------------- 2011 ----------------------------------
1644 void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
1645   //
1646   // D*+->D0pi syst errors. Responsible: A. Grelli
1647   //  2011 PbPb sample, 0-7.5 CC
1648   //
1649
1650   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1651
1652
1653  // Normalization
1654   fNorm = new TH1F("fNorm","fNorm",36,0,36);
1655   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1656
1657   // Branching ratio
1658   fBR = new TH1F("fBR","fBR",36,0,36);
1659   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1660
1661   // Tracking efficiency
1662   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1663   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1664
1665
1666   // Raw yield extraction
1667   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1668   for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05);  //4%
1669   fRawYield->SetBinContent(4,0.2);
1670   fRawYield->SetBinContent(5,0.10);
1671   fRawYield->SetBinContent(6,0.10);
1672   fRawYield->SetBinContent(7,0.08);
1673   for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15);  //4%
1674
1675   // Cuts efficiency (from cuts variation)
1676   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1677   for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1678
1679   // PID efficiency (from PID/noPID)
1680   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1681   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1682   fPIDEff->SetBinContent(4,0.09);
1683
1684   // MC dN/dpt  (from study on D* pt shape)
1685   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1686   for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
1687   fMCPtShape->SetBinContent(4,0.015);
1688   fMCPtShape->SetBinContent(5,0.015);
1689   fMCPtShape->SetBinContent(6,0.015);
1690   fMCPtShape->SetBinContent(7,0.02);
1691   fMCPtShape->SetBinContent(8,0.02);
1692   fMCPtShape->SetBinContent(9,0.03);
1693   fMCPtShape->SetBinContent(10,0.03);
1694   fMCPtShape->SetBinContent(11,0.03);
1695   fMCPtShape->SetBinContent(12,0.03);
1696
1697
1698
1699   return;
1700
1701 }
1702 //-------------------------------------------------------------------------
1703 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1704   //
1705   // D*+->D0pi syst errors. Responsible: A. Grelli
1706   //  2010 PbPb sample, 20-40 CC
1707   //
1708
1709   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1710
1711  // Normalization
1712   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1713   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1714
1715   // Branching ratio
1716   fBR = new TH1F("fBR","fBR",24,0,24);
1717   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1718
1719   // Tracking efficiency
1720   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1721   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1722
1723
1724   // Raw yield extraction
1725   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1726   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15);  //4%
1727
1728   // Cuts efficiency (from cuts variation)
1729   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1730   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1731
1732   // PID efficiency (from PID/noPID)
1733   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1734   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1735
1736
1737   // MC dN/dpt  (copied from D0 : will update later)
1738   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1739   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1740   fMCPtShape->SetBinContent(1,0.03);
1741   fMCPtShape->SetBinContent(2,0.03);
1742
1743   return;
1744
1745 }
1746
1747 //--------------------------------------------------------------------------
1748 void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1749   //
1750   // D*+->D0pi syst errors. Responsible: A. Grelli
1751   //  2010 PbPb sample, 40-80 CC
1752   //
1753
1754   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1755
1756  // Normalization
1757   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1758   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1759
1760   // Branching ratio
1761   fBR = new TH1F("fBR","fBR",24,0,24);
1762   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1763
1764   // Tracking efficiency
1765   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1766   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1767
1768
1769   // Raw yield extraction
1770   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1771   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
1772   fRawYield->SetBinContent(1,1);
1773   fRawYield->SetBinContent(2,0.15);
1774   fRawYield->SetBinContent(3,0.15);
1775   fRawYield->SetBinContent(4,0.15);
1776   fRawYield->SetBinContent(5,0.15);
1777   fRawYield->SetBinContent(6,0.10);
1778   fRawYield->SetBinContent(7,0.10);
1779   fRawYield->SetBinContent(8,0.10);
1780   fRawYield->SetBinContent(9,0.11);
1781   fRawYield->SetBinContent(10,0.11);
1782   fRawYield->SetBinContent(11,0.11);
1783   fRawYield->SetBinContent(12,0.11);
1784   fRawYield->SetBinContent(13,0.08);
1785   fRawYield->SetBinContent(14,0.08);
1786   fRawYield->SetBinContent(15,0.08);
1787   fRawYield->SetBinContent(16,0.08);
1788
1789
1790   // Cuts efficiency (from cuts variation)
1791   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1792   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1793
1794   // PID efficiency (from PID/noPID)
1795   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1796   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1797
1798
1799   // MC dN/dpt  (copied from D0 : will update later)
1800   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1801   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1802   fMCPtShape->SetBinContent(2,0.05);
1803   fMCPtShape->SetBinContent(3,0.05);
1804   fMCPtShape->SetBinContent(4,0.05);
1805   fMCPtShape->SetBinContent(5,0.04);
1806   fMCPtShape->SetBinContent(6,0.02);
1807   fMCPtShape->SetBinContent(7,0.04);
1808   fMCPtShape->SetBinContent(8,0.04);
1809
1810   return;
1811
1812 }
1813
1814 // ----------------------------- 2011 ----------------------------------
1815 void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
1816   //
1817   // D*+->D0pi syst errors. Responsible: A. Grelli
1818   //  2011 PbPb sample, 0-10% CC
1819   //
1820
1821   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
1822
1823
1824  // Normalization
1825   fNorm = new TH1F("fNorm","fNorm",36,0,36);
1826   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1827
1828   // Branching ratio
1829   fBR = new TH1F("fBR","fBR",36,0,36);
1830   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1831
1832   // Tracking efficiency
1833   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1834   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1835
1836
1837   // Raw yield extraction
1838   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1839   for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
1840   fRawYield->SetBinContent(4,0.10);
1841   fRawYield->SetBinContent(5,0.08);
1842   fRawYield->SetBinContent(6,0.05);
1843   for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
1844   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1845   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
1846   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
1847   for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10);  //4%
1848
1849   // Cuts efficiency (from cuts variation)
1850   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1851   for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
1852   fCutsEff->SetBinContent(4,0.15);
1853   for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1854
1855   // PID efficiency (from PID/noPID)
1856   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1857   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1858
1859   // MC dN/dpt  (from study on D* pt shape)
1860   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1861   for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
1862   fMCPtShape->SetBinContent(4,0.06);
1863   fMCPtShape->SetBinContent(5,0.03);
1864   for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
1865
1866   return;
1867
1868 }
1869
1870 //----------------------------------------------------------------------
1871 void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
1872     //
1873     // D*+->D0pi syst errors. Responsible: A. Grelli
1874     //  2011 PbPb sample, 30-50 CC
1875     //
1876
1877     AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
1878
1879     // Normalization
1880     fNorm = new TH1F("fNorm","fNorm",24,0,24);
1881     for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1882
1883     // Branching ratio
1884     fBR = new TH1F("fBR","fBR",24,0,24);
1885     for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1886
1887     // Tracking efficiency
1888     fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1889     for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1890
1891
1892     // Raw yield extraction
1893     fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1894     for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
1895     fRawYield->SetBinContent(1,0.12);
1896     fRawYield->SetBinContent(2,0.12);
1897     fRawYield->SetBinContent(3,0.12);
1898     fRawYield->SetBinContent(4,0.06);
1899     fRawYield->SetBinContent(5,0.06);
1900     fRawYield->SetBinContent(6,0.05);
1901     fRawYield->SetBinContent(7,0.08);
1902     fRawYield->SetBinContent(8,0.05);
1903     fRawYield->SetBinContent(9,0.05);
1904     fRawYield->SetBinContent(10,0.05);
1905     fRawYield->SetBinContent(11,0.05);
1906     fRawYield->SetBinContent(12,0.05);
1907     fRawYield->SetBinContent(13,0.05);
1908     fRawYield->SetBinContent(14,0.05);
1909     fRawYield->SetBinContent(15,0.05);
1910     fRawYield->SetBinContent(16,0.05);
1911
1912
1913     // Cuts efficiency (from cuts variation)
1914     fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1915     for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
1916     fCutsEff->SetBinContent(1,0.16);
1917     fCutsEff->SetBinContent(2,0.16);
1918     fCutsEff->SetBinContent(3,0.16);
1919     fCutsEff->SetBinContent(4,0.10);
1920     fCutsEff->SetBinContent(5,0.10);
1921     fCutsEff->SetBinContent(6,0.10);
1922     // PID efficiency (from PID/noPID)
1923     fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1924     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1925
1926
1927     // MC dN/dpt  (copied from D0 : will update later)
1928     fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1929     for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1930     fMCPtShape->SetBinContent(2,0.11);
1931     fMCPtShape->SetBinContent(3,0.11);
1932     fMCPtShape->SetBinContent(4,0.03);
1933     fMCPtShape->SetBinContent(5,0.02);
1934     fMCPtShape->SetBinContent(6,0.01);
1935     fMCPtShape->SetBinContent(7,0.01);
1936     fMCPtShape->SetBinContent(8,0.01);
1937
1938     return;
1939 }
1940 //_________________________________________________________________________
1941 void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
1942   // Dstar syst in pPb 2013 MB
1943   //
1944   AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
1945
1946   //Normalization
1947   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1948   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
1949
1950   // Branching ratio
1951   fBR = new TH1F("fBR","fBR",24,0,24);
1952   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1953
1954   // Tracking efficiency
1955   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1956   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1957
1958   // Raw yield extraction
1959   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1960   fRawYield->SetBinContent(1,1.0);
1961   fRawYield->SetBinContent(2,0.08);
1962   fRawYield->SetBinContent(3,0.02);
1963   fRawYield->SetBinContent(4,0.02);
1964   fRawYield->SetBinContent(5,0.05);
1965   fRawYield->SetBinContent(6,0.02);
1966   fRawYield->SetBinContent(7,0.04);
1967   fRawYield->SetBinContent(8,0.04);
1968   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1969   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
1970   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1971
1972   // Cuts efficiency (from cuts variation)
1973   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1974   fCutsEff->SetBinContent(1,0.);
1975   fCutsEff->SetBinContent(2,0.10);
1976   for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1977
1978   // PID efficiency (from PID/noPID)
1979   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1980   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
1981
1982   // MC dN/dpt
1983   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1984   for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
1985   for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
1986   for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1987
1988   // particle-antiparticle
1989   /*
1990   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1991   fPartAntipart->SetBinContent(1,1);
1992   fPartAntipart->SetBinContent(2,1);
1993   fPartAntipart->SetBinContent(3,0.12);
1994   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
1995   */
1996
1997     return;
1998
1999 }
2000
2001 //--------------------------------------------------------------------------
2002 void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
2003   // define errors for RAA vs. centrality
2004   InitD0toKpi2010PbPb020();
2005   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2006   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
2007   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2008 }
2009 //--------------------------------------------------------------------------
2010 void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
2011   // define errors for RAA vs. centrality
2012   InitD0toKpi2010PbPb020();
2013   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2014   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);
2015   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2016 }
2017 //--------------------------------------------------------------------------
2018 void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
2019   // define errors for RAA vs. centrality
2020   InitD0toKpi2010PbPb4080();
2021   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2022   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
2023   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2024 }
2025 //--------------------------------------------------------------------------
2026 void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
2027    // define errors for RAA vs. centrality
2028   InitD0toKpi2010PbPb4080();
2029   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2030   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);
2031   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2032 }
2033 //--------------------------------------------------------------------------
2034 void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
2035    // define errors for RAA vs. centrality
2036   InitD0toKpi2010PbPb4080();
2037   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2038   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
2039   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2040 }
2041 //--------------------------------------------------------------------------
2042 void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
2043   // define errors for RAA vs. centrality
2044   InitDplustoKpipi2010PbPb020();
2045   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
2046   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
2047
2048 }
2049 //--------------------------------------------------------------------------
2050 void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
2051   // define errors for RAA vs. centrality
2052   InitDplustoKpipi2010PbPb020();
2053   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
2054   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2055 }
2056 //--------------------------------------------------------------------------
2057 void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
2058   // define errors for RAA vs. centrality
2059   InitDplustoKpipi2010PbPb020();
2060   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2061   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
2062 }
2063 //--------------------------------------------------------------------------
2064 void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
2065   // define errors for RAA vs. centrality
2066   InitDplustoKpipi2010PbPb4080();
2067   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2068   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2069 }
2070 //--------------------------------------------------------------------------
2071 void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
2072   // define errors for RAA vs. centrality
2073   InitDplustoKpipi2010PbPb4080();
2074   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
2075   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
2076 }
2077
2078 //--------------------------------------------------------------------------
2079 void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
2080   // define errors for RAA vs. centrality
2081   InitDstartoD0pi2010PbPb020();
2082   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
2083   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2084 }
2085 //--------------------------------------------------------------------------
2086 void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
2087   // define errors for RAA vs. centrality
2088   InitDstartoD0pi2010PbPb020();
2089   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2090   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2091 }
2092 //--------------------------------------------------------------------------
2093 void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
2094   // define errors for RAA vs. centrality
2095   InitDstartoD0pi2010PbPb2040();
2096   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2097   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2098 }
2099 //--------------------------------------------------------------------------
2100 void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
2101   // define errors for RAA vs. centrality
2102   InitDstartoD0pi2010PbPb4080();
2103   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2104   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2105 }
2106 //--------------------------------------------------------------------------
2107 void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
2108   // define errors for RAA vs. centrality
2109   InitDstartoD0pi2010PbPb4080();
2110   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2111   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2112 }
2113
2114 //--------------------------------------------------------------------------
2115 void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
2116   // define errors for RAA vs. centrality 2011
2117   InitD0toKpi2011PbPb07half();
2118
2119 }
2120 //--------------------------------------------------------------------------
2121 void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
2122   // define errors for RAA vs. centrality 2011
2123   InitD0toKpi2011PbPb07half();
2124
2125 }
2126 //--------------------------------------------------------------------------
2127 void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
2128   // define errors for RAA vs. centrality 2011
2129   InitD0toKpi2011PbPb07half();
2130
2131 }
2132 //--------------------------------------------------------------------------
2133 void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
2134   // define errors for RAA vs. centrality 2011
2135   InitD0toKpi2011PbPb07half();
2136
2137 }
2138 //--------------------------------------------------------------------------
2139 void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
2140   // define errors for RAA vs. centrality 2011
2141   InitD0toKpi2011PbPb07half();
2142
2143 }
2144 //--------------------------------------------------------------------------
2145 void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
2146   // define errors for RAA vs. centrality 2011
2147   InitD0toKpi2011PbPb07half();
2148
2149 }
2150
2151 //--------------------------------------------------------------------------
2152 void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
2153   // define errors for RAA vs. centrality 2011
2154   InitDplustoKpipi2011PbPb07half();
2155
2156   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
2157   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2158   for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
2159
2160   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2161   for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2162   for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2163
2164   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2165   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2166   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
2167 }
2168
2169 //--------------------------------------------------------------------------
2170 void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
2171   // define errors for RAA vs. centrality 2011
2172   InitDplustoKpipi2011PbPb07half();
2173
2174   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2175   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2176   for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
2177
2178   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2179   for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2180   for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2181
2182   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2183   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2184   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
2185 }
2186
2187 //--------------------------------------------------------------------------
2188 void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
2189   // define errors for RAA vs. centrality 2011
2190   InitDplustoKpipi2011PbPb07half();
2191
2192   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2193   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2194   for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
2195
2196   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2197   for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2198   for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2199
2200   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2201   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2202   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2203 }
2204
2205 //--------------------------------------------------------------------------
2206 void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
2207   // define errors for RAA vs. centrality 2011
2208   InitDplustoKpipi2011PbPb07half();
2209
2210   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
2211   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2212   for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
2213
2214   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2215   for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2216   for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2217
2218   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2219   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2220   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2221 }
2222
2223 //--------------------------------------------------------------------------
2224 void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
2225   // define errors for RAA vs. centrality 2011
2226   InitDplustoKpipi2011PbPb07half();
2227
2228   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2229   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2230   for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
2231
2232   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2233   for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2234   for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2235
2236   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2237   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2238   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2239 }
2240
2241 //--------------------------------------------------------------------------
2242 void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
2243   // define errors for RAA vs. centrality 2011
2244   InitDplustoKpipi2010PbPb4080();
2245
2246   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2247   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2248   for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
2249
2250   for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2251   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2252   for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2253   for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2254 }
2255
2256 //--------------------------------------------------------------------------
2257 void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
2258   // define errors for RAA vs. centrality 2011
2259   InitDstartoD0pi2011PbPb07half();
2260
2261 }
2262 //--------------------------------------------------------------------------
2263 void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
2264   // define errors for RAA vs. centrality 2011
2265   InitDstartoD0pi2011PbPb07half();
2266
2267 }
2268 //--------------------------------------------------------------------------
2269 void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
2270   // define errors for RAA vs. centrality 2011
2271   InitDstartoD0pi2011PbPb07half();
2272
2273 }
2274 //--------------------------------------------------------------------------
2275 void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
2276   // define errors for RAA vs. centrality 2011
2277   InitDstartoD0pi2011PbPb07half();
2278
2279 }
2280 //--------------------------------------------------------------------------
2281 void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
2282   // define errors for RAA vs. centrality 2011
2283   InitDstartoD0pi2011PbPb07half();
2284
2285 }
2286 //--------------------------------------------------------------------------
2287 void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
2288   // define errors for RAA vs. centrality 2011
2289   InitDstartoD0pi2011PbPb07half();
2290
2291 }
2292
2293 //_________________________________________________________________________
2294 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
2295   InitD0toKpi2013pPb0100();
2296
2297   AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2298
2299   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2300   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
2301   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2302
2303   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2304
2305   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2306   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2307   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2308
2309   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2310   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2311   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2312
2313
2314 };
2315 //_________________________________________________________________________
2316 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
2317   InitD0toKpi2013pPb0100();
2318
2319    AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2320
2321   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2322   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2323   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2324
2325   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2326
2327    for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2328    for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2329    for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2330
2331   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2332   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2333   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2334 };
2335 //_________________________________________________________________________
2336 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
2337   InitD0toKpi2013pPb0100();
2338
2339   AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2340
2341   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2342   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
2343   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2344
2345   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2346
2347   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2348   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2349   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2350
2351   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2352   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2353   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2354
2355 };
2356 //_________________________________________________________________________
2357 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
2358   InitD0toKpi2013pPb0100();
2359
2360   AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2361
2362   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2363   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2364   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2365
2366   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2367
2368   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2369   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2370   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2371
2372   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2373   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2374   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2375
2376 };
2377 //_________________________________________________________________________
2378 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
2379   InitD0toKpi2013pPb0100();
2380
2381   AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2382
2383    for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2384    for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2385    for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2386
2387   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2388
2389   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2390   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2391   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2392
2393   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2394   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2395   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2396
2397 };
2398
2399 //_________________________________________________________________________
2400 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
2401   InitDplustoKpipi2013pPb0100();
2402
2403   AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2404
2405   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2406   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2407   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2408
2409   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2410
2411   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2412   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2413   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2414
2415   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2416   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
2417   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
2418
2419 };
2420 //_________________________________________________________________________
2421 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
2422   InitDplustoKpipi2013pPb0100();
2423   AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2424
2425    for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2426   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2427   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2428
2429   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2430
2431   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2432   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2433   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2434
2435   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2436   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2437   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2438
2439 };
2440 //_________________________________________________________________________
2441 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
2442   InitDplustoKpipi2013pPb0100();
2443   AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2444
2445   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2446   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2447   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2448
2449   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2450
2451   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2452   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2453   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2454
2455   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2456   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2457   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2458
2459 };
2460 //_________________________________________________________________________
2461 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
2462   InitDplustoKpipi2013pPb0100();
2463   AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2464
2465   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2466   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2467   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2468
2469   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2470
2471   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2472   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2473   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2474
2475   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2476   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2477   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2478
2479 };
2480 //_________________________________________________________________________
2481 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
2482   InitDplustoKpipi2013pPb0100();
2483
2484   AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2485
2486    for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2487    for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2488    for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2489
2490   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2491
2492   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2493   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2494   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2495
2496   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2497   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2498   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2499
2500 };
2501
2502
2503 //_________________________________________________________________________
2504 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
2505   InitDstartoD0pi2013pPb0100();
2506
2507   AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2508
2509   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2510   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2511   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2512
2513   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2514
2515   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2516   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2517   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2518
2519   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2520   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2521   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2522
2523 };
2524 //_________________________________________________________________________
2525 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
2526   InitDstartoD0pi2013pPb0100();
2527
2528   AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2529
2530   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2531   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2532   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2533
2534   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2535
2536   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2537   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2538   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2539
2540   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2541   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2542   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2543
2544 };
2545 //_________________________________________________________________________
2546 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
2547 InitDstartoD0pi2013pPb0100();
2548
2549    AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2550
2551   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2552   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2553   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2554
2555   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2556
2557   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2558   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2559   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2560
2561   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2562   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2563   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2564
2565 };
2566 //_________________________________________________________________________
2567 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
2568   InitDstartoD0pi2013pPb0100();
2569
2570   AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2571
2572   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2573   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2574   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2575
2576   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2577
2578   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2579   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2580   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2581
2582   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2583   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2584   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2585 };
2586 //_________________________________________________________________________
2587 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
2588   InitDstartoD0pi2013pPb0100();
2589
2590   AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2591
2592   for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2593   for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
2594   for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2595
2596   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2597
2598   for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2599   for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2600   for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2601
2602   for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2603   for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2604   for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2605
2606 };
2607
2608 //_________________________________________________________________________
2609 void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
2610     InitD0toKpi2013pPb0100();
2611
2612     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2613
2614
2615     // Raw yield extraction
2616     fRawYield->SetBinContent(1,0.0);
2617     fRawYield->SetBinContent(2,0.085);
2618     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2619     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2620     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2621     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2622     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2623     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2624
2625     // Cuts efficiency (from cuts variation)
2626     fCutsEff->SetBinContent(1,0.0);
2627     fCutsEff->SetBinContent(2,0.10);
2628     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2629
2630     // PID efficiency (from PID/noPID)
2631     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2632
2633
2634 }
2635
2636 //_________________________________________________________________________
2637 void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
2638     InitD0toKpi2013pPb0100();
2639
2640     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2641
2642
2643     // Raw yield extraction
2644     fRawYield->SetBinContent(1,0.0);
2645     fRawYield->SetBinContent(2,0.085);
2646     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2647     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2648     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2649     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2650     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2651     //for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2652
2653
2654     // Cuts efficiency (from cuts variation)
2655     fCutsEff->SetBinContent(1,0.0);
2656     fCutsEff->SetBinContent(2,0.10);
2657     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2658
2659
2660     // PID efficiency (from PID/noPID)
2661     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2662
2663 }
2664
2665 //_________________________________________________________________________
2666 void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
2667     InitD0toKpi2013pPb0100();
2668
2669     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2670
2671
2672     // Raw yield extraction
2673     fRawYield->SetBinContent(1,0.0);
2674     fRawYield->SetBinContent(2,0.085);
2675     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2676     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2677     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2678     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2679     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2680    // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2681
2682     // Cuts efficiency (from cuts variation)
2683     fCutsEff->SetBinContent(1,0.0);
2684     fCutsEff->SetBinContent(2,0.10);
2685     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2686
2687
2688     // PID efficiency (from PID/noPID)
2689     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2690
2691 }
2692
2693 //_________________________________________________________________________
2694 void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
2695     InitD0toKpi2013pPb0100();
2696
2697     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2698
2699     // Raw yield extraction
2700     fRawYield->SetBinContent(1,0.0);
2701     fRawYield->SetBinContent(2,0.085);
2702     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2703     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2704     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2705     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2706     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2707
2708
2709     // Cuts efficiency (from cuts variation)
2710     fCutsEff->SetBinContent(1,0.0);
2711     fCutsEff->SetBinContent(2,0.10);
2712     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2713
2714
2715     // PID efficiency (from PID/noPID)
2716     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2717
2718 }
2719
2720 //_________________________________________________________________________
2721 void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
2722     InitD0toKpi2013pPb0100();
2723
2724     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2725
2726     // Raw yield extraction
2727     fRawYield->SetBinContent(1,0.0);
2728     fRawYield->SetBinContent(2,0.085);
2729     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2730     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2731     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2732     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2733     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2734     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2735
2736     // Cuts efficiency (from cuts variation)
2737     fCutsEff->SetBinContent(1,0.0);
2738     fCutsEff->SetBinContent(2,0.10);
2739     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2740
2741
2742     // PID efficiency (from PID/noPID)
2743     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2744
2745 }
2746
2747 //_________________________________________________________________________
2748 void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
2749     InitD0toKpi2013pPb0100();
2750
2751     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
2752
2753  // Raw yield extraction
2754     fRawYield->SetBinContent(1,0.0);
2755     fRawYield->SetBinContent(2,0.085);
2756     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2757     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2758     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2759     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2760     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2761     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2762
2763
2764     // Cuts efficiency (from cuts variation)
2765     fCutsEff->SetBinContent(1,0.0);
2766     fCutsEff->SetBinContent(2,0.10);
2767     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2768
2769     // PID efficiency (from PID/noPID)
2770     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2771
2772 }
2773
2774 //_________________________________________________________________________
2775 void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
2776     InitD0toKpi2013pPb0100();
2777
2778     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
2779
2780  // Raw yield extraction
2781     fRawYield->SetBinContent(1,0.0);
2782     fRawYield->SetBinContent(2,0.085);
2783     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2784     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2785     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2786     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2787     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2788     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2789
2790     // Cuts efficiency (from cuts variation)
2791     fCutsEff->SetBinContent(1,0.0);
2792     fCutsEff->SetBinContent(2,0.10);
2793     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2794
2795     // PID efficiency (from PID/noPID)
2796     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2797
2798 }
2799
2800 //_________________________________________________________________________
2801 void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
2802     InitD0toKpi2013pPb0100();
2803
2804     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
2805
2806  // Raw yield extraction
2807     fRawYield->SetBinContent(1,0.0);
2808     fRawYield->SetBinContent(2,0.085);
2809     for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2810     for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2811     for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2812     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2813     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2814    // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2815
2816     // Cuts efficiency (from cuts variation)
2817     fCutsEff->SetBinContent(1,0.0);
2818     fCutsEff->SetBinContent(2,0.10);
2819     for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2820
2821     // PID efficiency (from PID/noPID)
2822     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2823
2824 }
2825
2826 //_________________________________________________________________________
2827 void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
2828     InitDstartoD0pi2013pPb0100();
2829
2830     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2831
2832     // Raw yield extraction
2833     fRawYield->SetBinContent(1,0.0);
2834     fRawYield->SetBinContent(2,0.10);
2835     fRawYield->SetBinContent(3,0.06);
2836     fRawYield->SetBinContent(4,0.06);
2837     fRawYield->SetBinContent(5,0.05);
2838     fRawYield->SetBinContent(6,0.05);
2839     fRawYield->SetBinContent(7,0.05);
2840     fRawYield->SetBinContent(8,0.05);
2841     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2842     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2843     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2844
2845     // Cuts efficiency (from cuts variation)
2846     fCutsEff->SetBinContent(1,0.0);
2847     fCutsEff->SetBinContent(2,0.10);
2848     fCutsEff->SetBinContent(3,0.06);
2849     fCutsEff->SetBinContent(4,0.06);
2850     fCutsEff->SetBinContent(5,0.06);
2851     fCutsEff->SetBinContent(6,0.06);
2852     fCutsEff->SetBinContent(7,0.06);
2853     fCutsEff->SetBinContent(8,0.06);
2854     for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2855     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2856
2857     // PID efficiency (from PID/noPID)
2858     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2859
2860
2861 }
2862
2863 //_________________________________________________________________________
2864 void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
2865     InitDstartoD0pi2013pPb0100();
2866
2867     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2868
2869     // Raw yield extraction
2870     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2871     fRawYield->SetBinContent(1,0.0);
2872     fRawYield->SetBinContent(2,0.10);
2873     fRawYield->SetBinContent(3,0.06);
2874     fRawYield->SetBinContent(4,0.06);
2875     fRawYield->SetBinContent(5,0.05);
2876     fRawYield->SetBinContent(6,0.05);
2877     fRawYield->SetBinContent(7,0.05);
2878     fRawYield->SetBinContent(8,0.05);
2879     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2880     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2881     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
2882
2883     // Cuts efficiency (from cuts variation)
2884     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2885     fCutsEff->SetBinContent(1,0.);
2886     fCutsEff->SetBinContent(2,0.10);
2887     fCutsEff->SetBinContent(3,0.06);
2888     fCutsEff->SetBinContent(4,0.06);
2889     fCutsEff->SetBinContent(5,0.06);
2890     fCutsEff->SetBinContent(6,0.06);
2891     fCutsEff->SetBinContent(7,0.06);
2892     fCutsEff->SetBinContent(8,0.06);
2893     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2894     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2895     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2896
2897     // PID efficiency (from PID/noPID)
2898     //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2899     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2900
2901     // MC dN/dpt
2902     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2903      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2904      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2905      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2906
2907     // particle-antiparticle
2908     /*
2909      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2910      fPartAntipart->SetBinContent(1,1);
2911      fPartAntipart->SetBinContent(2,1);
2912      fPartAntipart->SetBinContent(3,0.12);
2913      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
2914     */
2915
2916 }
2917
2918 //_________________________________________________________________________
2919 void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
2920     InitDstartoD0pi2013pPb0100();
2921
2922     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2923
2924     // Raw yield extraction
2925     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2926     fRawYield->SetBinContent(1,0.0);
2927     fRawYield->SetBinContent(2,0.0);
2928     fRawYield->SetBinContent(3,0.06);
2929     fRawYield->SetBinContent(4,0.06);
2930     fRawYield->SetBinContent(5,0.05);
2931     fRawYield->SetBinContent(6,0.05);
2932     fRawYield->SetBinContent(7,0.05);
2933     fRawYield->SetBinContent(8,0.05);
2934     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2935     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2936     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2937
2938     // Cuts efficiency (from cuts variation)
2939     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2940     fCutsEff->SetBinContent(1,0.);
2941     fCutsEff->SetBinContent(2,0.);
2942     fCutsEff->SetBinContent(3,0.06);
2943     fCutsEff->SetBinContent(4,0.06);
2944     fCutsEff->SetBinContent(5,0.06);
2945     fCutsEff->SetBinContent(6,0.06);
2946     fCutsEff->SetBinContent(7,0.06);
2947     fCutsEff->SetBinContent(8,0.06);
2948     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2949     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2950     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2951
2952     // PID efficiency (from PID/noPID)
2953     //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2954     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2955
2956     // MC dN/dpt
2957     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2958      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2959      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2960      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2961
2962     // particle-antiparticle
2963     /*
2964      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2965      fPartAntipart->SetBinContent(1,1);
2966      fPartAntipart->SetBinContent(2,1);
2967      fPartAntipart->SetBinContent(3,0.12);
2968      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
2969      */
2970
2971
2972 }
2973
2974 //_________________________________________________________________________
2975 void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
2976     InitDstartoD0pi2013pPb0100();
2977
2978     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2979
2980     // Raw yield extraction
2981     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2982     fRawYield->SetBinContent(1,0.0);
2983     fRawYield->SetBinContent(2,0.0);
2984     fRawYield->SetBinContent(3,0.06);
2985     fRawYield->SetBinContent(4,0.06);
2986     fRawYield->SetBinContent(5,0.05);
2987     fRawYield->SetBinContent(6,0.05);
2988     fRawYield->SetBinContent(7,0.05);
2989     fRawYield->SetBinContent(8,0.05);
2990     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2991     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2992     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2993
2994     // Cuts efficiency (from cuts variation)
2995     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2996     fCutsEff->SetBinContent(1,0.);
2997     fCutsEff->SetBinContent(2,0.);
2998     fCutsEff->SetBinContent(3,0.06);
2999     fCutsEff->SetBinContent(4,0.06);
3000     fCutsEff->SetBinContent(5,0.06);
3001     fCutsEff->SetBinContent(6,0.06);
3002     fCutsEff->SetBinContent(7,0.06);
3003     fCutsEff->SetBinContent(8,0.06);
3004     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3005     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3006     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3007
3008     // PID efficiency (from PID/noPID)
3009     //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3010     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3011
3012     // MC dN/dpt
3013     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3014      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3015      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3016      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3017
3018     // particle-antiparticle
3019     /*
3020      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3021      fPartAntipart->SetBinContent(1,1);
3022      fPartAntipart->SetBinContent(2,1);
3023      fPartAntipart->SetBinContent(3,0.12);
3024      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
3025      */
3026
3027
3028 }
3029
3030 //_________________________________________________________________________
3031 void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
3032     InitDstartoD0pi2013pPb0100();
3033
3034     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3035
3036     // Raw yield extraction
3037     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3038     fRawYield->SetBinContent(1,0.);
3039     fRawYield->SetBinContent(2,0.10);
3040     fRawYield->SetBinContent(3,0.06);
3041     fRawYield->SetBinContent(4,0.06);
3042     fRawYield->SetBinContent(5,0.05);
3043     fRawYield->SetBinContent(6,0.05);
3044     fRawYield->SetBinContent(7,0.05);
3045     fRawYield->SetBinContent(8,0.05);
3046     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3047     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3048     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3049
3050     // Cuts efficiency (from cuts variation)
3051     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3052     fCutsEff->SetBinContent(1,0.);
3053     fCutsEff->SetBinContent(2,0.10);
3054     fCutsEff->SetBinContent(3,0.06);
3055     fCutsEff->SetBinContent(4,0.06);
3056     fCutsEff->SetBinContent(5,0.06);
3057     fCutsEff->SetBinContent(6,0.06);
3058     fCutsEff->SetBinContent(7,0.06);
3059     fCutsEff->SetBinContent(8,0.06);
3060     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3061     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3062     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3063
3064     // PID efficiency (from PID/noPID)
3065     // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3066     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
3067
3068     // MC dN/dpt
3069     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3070      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3071      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3072      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3073
3074     // particle-antiparticle
3075     /*
3076      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3077      fPartAntipart->SetBinContent(1,1);
3078      fPartAntipart->SetBinContent(2,1);
3079      fPartAntipart->SetBinContent(3,0.12);
3080      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
3081      */
3082
3083
3084 }
3085
3086 //_________________________________________________________________________
3087 void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
3088     InitDstartoD0pi2013pPb0100();
3089
3090     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3091
3092     // Raw yield extraction
3093     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3094     fRawYield->SetBinContent(1,0.0);
3095     fRawYield->SetBinContent(2,0.0);
3096     fRawYield->SetBinContent(3,0.06);
3097     fRawYield->SetBinContent(4,0.06);
3098     fRawYield->SetBinContent(5,0.05);
3099     fRawYield->SetBinContent(6,0.05);
3100     fRawYield->SetBinContent(7,0.05);
3101     fRawYield->SetBinContent(8,0.05);
3102     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3103     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3104     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3105
3106     // Cuts efficiency (from cuts variation)
3107     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3108     fCutsEff->SetBinContent(1,0.);
3109     fCutsEff->SetBinContent(2,0.);
3110     fCutsEff->SetBinContent(3,0.06);
3111     fCutsEff->SetBinContent(4,0.06);
3112     fCutsEff->SetBinContent(5,0.06);
3113     fCutsEff->SetBinContent(6,0.06);
3114     fCutsEff->SetBinContent(7,0.06);
3115     fCutsEff->SetBinContent(8,0.06);
3116     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3117     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3118     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3119
3120     // PID efficiency (from PID/noPID)
3121     //  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3122     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3123
3124     // MC dN/dpt
3125     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3126      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3127      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3128      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3129
3130     // particle-antiparticle
3131     /*
3132      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3133      fPartAntipart->SetBinContent(1,1);
3134      fPartAntipart->SetBinContent(2,1);
3135      fPartAntipart->SetBinContent(3,0.12);
3136      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
3137      */
3138
3139
3140 }
3141
3142 //_________________________________________________________________________
3143 void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
3144     InitDstartoD0pi2013pPb0100();
3145
3146     AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3147
3148     // Raw yield extraction
3149     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3150     fRawYield->SetBinContent(1,0.0);
3151     fRawYield->SetBinContent(2,0.0);
3152     fRawYield->SetBinContent(3,0.06);
3153     fRawYield->SetBinContent(4,0.06);
3154     fRawYield->SetBinContent(5,0.05);
3155     fRawYield->SetBinContent(6,0.05);
3156     fRawYield->SetBinContent(7,0.05);
3157     fRawYield->SetBinContent(8,0.05);
3158     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3159     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3160     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3161
3162     // Cuts efficiency (from cuts variation)
3163     // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3164     fCutsEff->SetBinContent(1,0.);
3165     fCutsEff->SetBinContent(2,0.);
3166     fCutsEff->SetBinContent(3,0.06);
3167     fCutsEff->SetBinContent(4,0.06);
3168     fCutsEff->SetBinContent(5,0.06);
3169     fCutsEff->SetBinContent(6,0.06);
3170     fCutsEff->SetBinContent(7,0.06);
3171     fCutsEff->SetBinContent(8,0.06);
3172     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3173     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3174     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3175
3176     // PID efficiency (from PID/noPID)
3177     // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3178     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3179
3180     // MC dN/dpt
3181     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3182      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3183      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3184      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3185
3186     // particle-antiparticle
3187     /*
3188      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3189      fPartAntipart->SetBinContent(1,1);
3190      fPartAntipart->SetBinContent(2,1);
3191      fPartAntipart->SetBinContent(3,0.12);
3192      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
3193      */
3194
3195
3196
3197 }
3198
3199 //_________________________________________________________________________
3200 void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
3201     InitDstartoD0pi2013pPb0100();
3202
3203     AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3204
3205     // Raw yield extraction
3206     //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3207     fRawYield->SetBinContent(1,0.);
3208     fRawYield->SetBinContent(2,0.);
3209     fRawYield->SetBinContent(3,0.06);
3210     fRawYield->SetBinContent(4,0.06);
3211     fRawYield->SetBinContent(5,0.05);
3212     fRawYield->SetBinContent(6,0.05);
3213     fRawYield->SetBinContent(7,0.05);
3214     fRawYield->SetBinContent(8,0.05);
3215     for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3216     for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3217     for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3218
3219     // Cuts efficiency (from cuts variation)
3220     //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3221     fCutsEff->SetBinContent(1,0.);
3222     fCutsEff->SetBinContent(2,0.);
3223     fCutsEff->SetBinContent(3,0.06);
3224     fCutsEff->SetBinContent(4,0.06);
3225     fCutsEff->SetBinContent(5,0.06);
3226     fCutsEff->SetBinContent(6,0.06);
3227     fCutsEff->SetBinContent(7,0.06);
3228     fCutsEff->SetBinContent(8,0.06);
3229     for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3230     for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3231     for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3232
3233     // PID efficiency (from PID/noPID)
3234     // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3235     for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3236
3237     // MC dN/dpt
3238     /*    fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3239      for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3240      for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3241      for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3242
3243     // particle-antiparticle
3244     /*
3245      fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3246      fPartAntipart->SetBinContent(1,1);
3247      fPartAntipart->SetBinContent(2,1);
3248      fPartAntipart->SetBinContent(3,0.12);
3249      for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
3250      */
3251
3252 }
3253
3254
3255
3256
3257 //--------------------------------------------------------------------------
3258 void AliHFSystErr::InitLctopKpi2010pp() {
3259   //
3260   // Lc->pKpi syst errors. Responsible: R. Romita
3261   //  2010 pp sample
3262   //
3263
3264   // Normalization
3265   fNorm = new TH1F("fNorm","fNorm",4,2,6);
3266   for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3267
3268   // Tracking efficiency
3269   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
3270   for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
3271
3272   // Raw yield extraction
3273   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3274   fRawYield->SetBinContent(1,0.20);
3275   fRawYield->SetBinContent(2,0.15);
3276   fRawYield->SetBinContent(3,0.10);
3277   fRawYield->SetBinContent(4,0.15);
3278
3279   fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
3280   fCutsEff->SetBinContent(1,0.50);
3281   fCutsEff->SetBinContent(2,0.07);
3282   fCutsEff->SetBinContent(3,0.30);
3283   fCutsEff->SetBinContent(4,0.40);
3284
3285   // PID efficiency (from PID/noPID)
3286   fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
3287   fPIDEff->SetBinContent(1,0.13); // 15%
3288   fPIDEff->SetBinContent(2,0.30); // 15%
3289   fPIDEff->SetBinContent(3,0.5); // 15%
3290   fPIDEff->SetBinContent(4,0.30); // 15%
3291
3292   // MC dN/dpt
3293   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
3294   for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
3295   for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
3296   for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
3297   for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
3298   for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
3299 }
3300
3301 //--------------------------------------------------------------------------
3302 void AliHFSystErr::InitLctopK0S2010pp() {
3303   //
3304   // Lc->pK0S syst errors. Responsible: A. De Caro
3305   //  2010 pp sample
3306   //  Only tool preparation: no physical mean for these values!
3307   //
3308
3309   // Normalization
3310   fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
3311   for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and // only a try
3312
3313   // Tracking efficiency
3314   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
3315   for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1); // 10% (4% per track) // only a try
3316
3317   // Raw yield extraction
3318   fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
3319   for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.); // only a try
3320   for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1); // only a try
3321
3322   fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
3323   for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1); // only a try
3324
3325   // PID efficiency (from PID/noPID)
3326   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
3327   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // only a try
3328
3329   // MC dN/dpt
3330   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
3331   for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.); // only a try
3332   for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1); // only a try
3333
3334 }
3335
3336 //--------------------------------------------------------------------------
3337 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
3338   //
3339   // Get error
3340   //
3341
3342   Int_t bin=fCutsEff->FindBin(pt);
3343
3344   return fCutsEff->GetBinContent(bin);
3345 }
3346 //--------------------------------------------------------------------------
3347 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
3348   //
3349   // Get error
3350   //
3351
3352   Int_t bin=fMCPtShape->FindBin(pt);
3353
3354   return fMCPtShape->GetBinContent(bin);
3355 }
3356 //--------------------------------------------------------------------------
3357 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
3358   //
3359   // Get error
3360   //
3361
3362   Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3363
3364   return TMath::Sqrt(err);
3365 }
3366 //--------------------------------------------------------------------------
3367 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
3368   //
3369   // Get error
3370   //
3371
3372   Int_t bin=fPIDEff->FindBin(pt);
3373
3374   return fPIDEff->GetBinContent(bin);
3375 }
3376 //--------------------------------------------------------------------------
3377 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
3378   //
3379   // Get error
3380   //
3381
3382   Int_t bin=fTrackingEff->FindBin(pt);
3383
3384   return fTrackingEff->GetBinContent(bin);
3385 }
3386 //--------------------------------------------------------------------------
3387 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
3388   //
3389   // Get error
3390   //
3391
3392   Int_t bin=fRawYield->FindBin(pt);
3393
3394   return fRawYield->GetBinContent(bin);
3395 }
3396 //--------------------------------------------------------------------------
3397 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
3398   //
3399   // Get error
3400   //
3401
3402   Int_t bin=fPartAntipart->FindBin(pt);
3403
3404   return fPartAntipart->GetBinContent(bin);
3405 }
3406 //--------------------------------------------------------------------------
3407 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
3408   //
3409   // Get total syst error (except norm. error)
3410   //
3411
3412   Double_t err=0.;
3413
3414   if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
3415   if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
3416   //  if(fBR) err += GetBRErr()*GetBRErr();
3417   if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
3418   if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
3419   if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3420   if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
3421
3422   err += feeddownErr*feeddownErr;
3423
3424   return TMath::Sqrt(err);
3425 }
3426 //---------------------------------------------------------------------------
3427 void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
3428   //
3429   // Draw errors
3430   //
3431   gStyle->SetOptStat(0);
3432
3433   TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
3434   cSystErr->Range(0.20,-0.5,18.4,0.34);
3435   cSystErr->SetRightMargin(0.318);
3436   cSystErr->SetFillColor(0);
3437
3438   TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
3439   hFrame->SetAxisRange(1.,35.9,"X");
3440   hFrame->SetAxisRange(-0.5,0.5,"Y");
3441   hFrame->Draw();
3442
3443   TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
3444   leg->SetTextSize(0.03601695);
3445   leg->SetFillStyle(0);
3446   leg->SetBorderSize(0);
3447
3448   TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
3449   Int_t nbins = fNorm->GetNbinsX();
3450   TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
3451   for(Int_t i=1;i<=36;i++) {
3452     Double_t pt = hTotErr->GetBinCenter(i);
3453     Double_t ptwidth = hTotErr->GetBinWidth(i);
3454
3455     if(grErrFeeddown) {
3456       Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
3457       Double_t toterryl=0., toterryh=0.;
3458       for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
3459         grErrFeeddown->GetPoint(j,x,y);
3460         errxh = grErrFeeddown->GetErrorXhigh(j);
3461         errxl = grErrFeeddown->GetErrorXlow(j);
3462         if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
3463           erryh = grErrFeeddown->GetErrorYhigh(j);
3464           erryl = grErrFeeddown->GetErrorYlow(j);
3465         }
3466       }
3467       if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
3468       else toterryl = GetTotalSystErr(pt);
3469       if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
3470       else toterryh = GetTotalSystErr(pt);
3471
3472       hTotErr->SetBinContent(i,toterryh);
3473       gTotErr->SetPoint(i,pt,0.);
3474       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
3475     }
3476     else {
3477       hTotErr->SetBinContent(i,GetTotalSystErr(pt));
3478       gTotErr->SetPoint(i,pt,0.);
3479       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
3480     }
3481
3482   }
3483   gTotErr->SetLineColor(kBlack);
3484   gTotErr->SetFillColor(kRed);
3485   gTotErr->SetFillStyle(3002);
3486   gTotErr->Draw("2");
3487   leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
3488 //   hTotErr->SetLineColor(1);
3489 //   hTotErr->SetLineWidth(3);
3490 //   hTotErr->Draw("same");
3491 //   leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
3492
3493
3494   fNorm->SetFillColor(1);
3495   fNorm->SetFillStyle(3002);
3496   //fNorm->Draw("same");
3497   //TH1F *hNormRefl = ReflectHisto(fNorm);
3498   //hNormRefl->Draw("same");
3499   Double_t norm = fNorm->GetBinContent(1)*100;
3500   leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
3501
3502   if(grErrFeeddown) {
3503     grErrFeeddown->SetFillColor(kTeal-8);
3504     grErrFeeddown->SetFillStyle(3001);
3505     grErrFeeddown->Draw("2");
3506     leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
3507   }
3508   if(fTrackingEff) {
3509     fTrackingEff->SetFillColor(4);
3510     fTrackingEff->SetFillStyle(3006);
3511     fTrackingEff->Draw("same");
3512     TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
3513     hTrackingEffRefl->Draw("same");
3514     leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
3515   }
3516   if(fBR) {
3517     fBR->SetFillColor(6);
3518     fBR->SetFillStyle(3005);
3519     //fBR->SetFillStyle(3020);
3520     fBR->Draw("same");
3521     TH1F *hBRRefl = ReflectHisto(fBR);
3522     hBRRefl->Draw("same");
3523     leg->AddEntry(fBR,"Branching ratio","f");
3524   }
3525   if(fRawYield) {
3526     Int_t ci;   // for color index setting
3527     ci = TColor::GetColor("#00cc00");
3528     fRawYield->SetLineColor(ci);
3529     //    fRawYield->SetLineColor(3);
3530     fRawYield->SetLineWidth(3);
3531     fRawYield->Draw("same");
3532     TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
3533     hRawYieldRefl->Draw("same");
3534     leg->AddEntry(fRawYield,"Yield extraction","l");
3535   }
3536   if(fCutsEff) {
3537     fCutsEff->SetLineColor(4);
3538     fCutsEff->SetLineWidth(3);
3539     fCutsEff->Draw("same");
3540     TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
3541     hCutsEffRefl->Draw("same");
3542     leg->AddEntry(fCutsEff,"Cuts efficiency","l");
3543   }
3544   if(fPIDEff) {
3545     fPIDEff->SetLineColor(7);
3546     fPIDEff->SetLineWidth(3);
3547     fPIDEff->Draw("same");
3548     TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
3549     hPIDEffRefl->Draw("same");
3550     leg->AddEntry(fPIDEff,"PID efficiency","l");
3551   }
3552   if(fMCPtShape) {
3553     Int_t ci = TColor::GetColor("#9933ff");
3554     fMCPtShape->SetLineColor(ci);
3555     //    fMCPtShape->SetLineColor(8);
3556     fMCPtShape->SetLineWidth(3);
3557     fMCPtShape->Draw("same");
3558     TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
3559     hMCPtShapeRefl->Draw("same");
3560     leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
3561   }
3562   if(fPartAntipart) {
3563     Int_t ci = TColor::GetColor("#ff6600");
3564     fPartAntipart->SetLineColor(ci);
3565     //    fPartAntipart->SetLineColor(9);
3566     fPartAntipart->SetLineWidth(3);
3567     fPartAntipart->Draw("same");
3568     TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
3569     hPartAntipartRefl->Draw("same");
3570     leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
3571   }
3572
3573
3574   leg->Draw();
3575
3576   cSystErr->SaveAs("RelativeSystematics.eps");
3577
3578   return;
3579 }
3580 //-------------------------------------------------------------------------
3581 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
3582   //
3583   // Clones and reflects histogram
3584   //
3585   TH1F *hout=(TH1F*)hin->Clone("hout");
3586   hout->Scale(-1.);
3587
3588   return hout;
3589 }
3590
3591
3592