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