PWGHF/vertexingHF/AliHFSystErr.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 /////////////////////////////////////////////////////////////
  19 //
  20 // Class to handle systematic errors for charm hadrons
  21 //
  22 // Usage:
  23 // AliHFSystEff syst;           // DECAY = 1 for D0, 2, for D+, 3 for D*
  24 // syst.SetRunNumber(YEAR);     // YEAR = two last numbers of the year (is 10 for 2010)
  25 // syst.SetCollisionType(TYPE);  // TYPE =  0 is pp, 1 is PbPb
  26 // syst.SetCentrality(CENT);     // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
  27 // syst.Init(DECAY);             // DECAY = 1 for D0, 2, for D+, 3 for D*
  28 // syst.DrawErrors(); // to see a plot of the error contributions
  29 // syst.GetTotalSystErr(pt); // to get the total err at pt
  30 //
  31 // Author: A.Dainese, andrea.dainese@pd.infn.it
  32 /////////////////////////////////////////////////////////////
  33
  34 #include <TStyle.h>
  35 #include <TGraphAsymmErrors.h>
  36 #include <TMath.h>
  37 #include <TCanvas.h>
  38 #include <TH2F.h>
  39 #include <TLegend.h>
  40 #include <TColor.h>
  41
  42 #include "AliLog.h"
  43 #include "AliHFSystErr.h"
  44
  45
  46 ClassImp(AliHFSystErr)
  47
  48 //--------------------------------------------------------------------------
  49 AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
  50 TNamed(name,title),
  51 fNorm(0),
  52 fRawYield(0),
  53 fTrackingEff(0),
  54 fBR(0),
  55 fCutsEff(0),
  56 fPIDEff(0),
  57 fMCPtShape(0),
  58 fPartAntipart(0),
  59 fRunNumber(10),
  60 fCollisionType(0),
  61 fCentralityClass("0100"),
  62 fIsLowEnergy(false),
  63 fIsCentScan(false)
  64 {
  65   //
  66   // Default Constructor
  67   //
  68 }
  69
  70 //--------------------------------------------------------------------------
  71 AliHFSystErr::~AliHFSystErr() {
  72   //
  73   // Default Destructor
  74   //
  75   /*
  76
  77   if(fNorm)         { delete fNorm; fNorm=0; }
  78   if(fRawYield)     { delete fRawYield; fRawYield=0; }
  79   if(fTrackingEff)  { delete fTrackingEff; fTrackingEff=0; }
  80   if(fBR)           { delete fBR; fBR=0; }
  81   if(fCutsEff)      { delete fCutsEff; fCutsEff=0; }
  82   if(fPIDEff)       { delete fPIDEff; fPIDEff=0; }
  83   if(fMCPtShape)    { delete fMCPtShape; fMCPtShape=0; }
  84   if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
  85   */
  86 }
  87
  88 //--------------------------------------------------------------------------
  89 void AliHFSystErr::Init(Int_t decay){
  90   //
  91   // Variables/histos initialization
  92   //
  93
  94   if (fRunNumber!=10 && fIsLowEnergy==false) {
  95     AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
  96   }
  97
  98   switch(decay) {
  99   case 1: // D0->Kpi
 100     if (fCollisionType==0) {
 101       if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
 102       else InitD0toKpi2010pp();
 103     } else if (fCollisionType==1) {
 104       if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
 105       else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
 106       else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
 107       else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
 108       else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
 109       else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
 110       else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
 111       else AliFatal("Not yet implemented");
 112     }
 113     //    else if (fCollisionType==2) InitD0toKpi2010ppLowEn();
 114     break;
 115   case 2: // D+->Kpipi
 116     if (fCollisionType==0) {
 117       if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
 118       else InitDplustoKpipi2010pp();
 119     } else if (fCollisionType==1) {
 120       if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
 121       else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
 122       else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
 123       else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
 124       else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
 125       else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
 126       else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
 127       else AliFatal("Not yet implemented");
 128     }
 129     break;
 130   case 3: // D*->D0pi
 131     if (fCollisionType==0) {
 132       if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
 133       else InitDstartoD0pi2010pp();
 134     }else if (fCollisionType==1) {
 135       if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
 136       else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
 137       else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
 138       else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
 139       else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
 140       else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
 141       else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
 142       else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
 143       else AliFatal("Not yet implemented");
 144     }
 145     break;
 146   case 4: // D+s->KKpi
 147     if (fCollisionType==0) InitDstoKKpi2010pp();
 148     else AliFatal("Not yet implemented");
 149     break;
 150
 151   default:
 152     printf("Invalid decay type: %d\n",decay);
 153     break;
 154   }
 155
 156 }
 157
 158 //--------------------------------------------------------------------------
 159 void AliHFSystErr::InitD0toKpi2010pp() {
 160   //
 161   // D0->Kpi syst errors. Responsible: A. Rossi
 162   //   2010 pp sample
 163   //
 164
 165   // Normalization
 166   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 167   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 168
 169   // Branching ratio
 170   fBR = new TH1F("fBR","fBR",24,0,24);
 171   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 172
 173   // Tracking efficiency
 174   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 175   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
 176
 177   // Raw yield extraction
 178   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 179   fRawYield->SetBinContent(1,1);
 180   fRawYield->SetBinContent(2,0.22);
 181   fRawYield->SetBinContent(3,0.1);
 182   for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
 183   for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
 184   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 185   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
 186
 187   // Cuts efficiency (from cuts variation)
 188   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 189   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 190
 191   // PID efficiency (from PID/noPID)
 192   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 193   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
 194   fPIDEff->SetBinContent(2,0.05); // 5%
 195
 196   // MC dN/dpt
 197   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 198   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 199   fMCPtShape->SetBinContent(1,0.03);
 200   fMCPtShape->SetBinContent(2,0.03);
 201
 202   // particle-antiparticle
 203   //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
 204   //  fPartAntipart->SetBinContent(1,1);
 205   //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
 206
 207   return;
 208 }
 209 //--------------------------------------------------------------------------
 210 void AliHFSystErr::InitD0toKpi2010PbPb020() {
 211   //
 212   // D0->Kpi syst errors. Responsible: A. Rossi
 213   //   2010 PbPb sample, 0-20 CC
 214   //
 215
 216   // Normalization
 217   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 218   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 219
 220   // Branching ratio
 221   fBR = new TH1F("fBR","fBR",20,0,20);
 222   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 223
 224   // Tracking efficiency
 225   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 226   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 227
 228   // Raw yield extraction
 229   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 230   fRawYield->SetBinContent(1,0);
 231   fRawYield->SetBinContent(2,0);
 232   fRawYield->SetBinContent(3,0.08);
 233   for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
 234   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 235
 236   // Cuts efficiency (from cuts variation)
 237   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 238   fCutsEff->SetBinContent(1,0.);
 239   fCutsEff->SetBinContent(2,0.);
 240   fCutsEff->SetBinContent(3,0.13);
 241   fCutsEff->SetBinContent(4,0.11);
 242   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 243   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 244
 245   // PID efficiency (from PID/noPID)
 246   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 247   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 248
 249   // MC dN/dpt
 250   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 251   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 252   fMCPtShape->SetBinContent(3,0.04);
 253   fMCPtShape->SetBinContent(4,0.02);
 254   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 255
 256 //   // particle-antiparticle
 257 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 258 //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 259 //   fPartAntipart->SetBinContent(3,0.10);
 260 //   fPartAntipart->SetBinContent(4,0.10);
 261 //   fPartAntipart->SetBinContent(7,0.10);
 262 //   fPartAntipart->SetBinContent(8,0.10);
 263
 264   return;
 265 }
 266 //--------------------------------------------------------------------------
 267 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
 268   //
 269   // D0->Kpi syst errors. Responsible: A. Rossi
 270   //   2010 PbPb sample, 40-80 CC
 271   //
 272
 273   // Normalization
 274   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 275   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 276
 277   // Branching ratio
 278   fBR = new TH1F("fBR","fBR",20,0,20);
 279   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 280
 281   // Tracking efficiency
 282   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 283   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 284
 285
 286   // Raw yield extraction
 287   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 288   fRawYield->SetBinContent(1,0);
 289   fRawYield->SetBinContent(2,0);
 290   for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
 291   //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
 292
 293   // Cuts efficiency (from cuts variation)
 294   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 295   fCutsEff->SetBinContent(1,0.);
 296   fCutsEff->SetBinContent(2,0.);
 297   fCutsEff->SetBinContent(3,0.13);
 298   fCutsEff->SetBinContent(4,0.11);
 299   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 300   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 301
 302   // PID efficiency (from PID/noPID)
 303   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 304 //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
 305 //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 306   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 307
 308   // MC dN/dpt
 309   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 310   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 311   //  fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
 312   //  fMCPtShape->SetBinContent(4,0.02);
 313   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 314
 315   //   // particle-antiparticle
 316   //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 317   //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 318
 319   return;
 320 }
 321
 322 //--------------------------------------------------------------------------
 323 void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 324   //
 325   // D0->Kpi syst errors. Low energy run
 326   //   2011 2.76 TeV pp sample
 327   //
 328   AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
 329
 330   // Normalization
 331   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 332   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
 333
 334   // Branching ratio
 335   fBR = new TH1F("fBR","fBR",20,0,20);
 336   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 337
 338   // Tracking efficiency
 339   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 340   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
 341
 342   // Raw yield extraction
 343   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 344   fRawYield->SetBinContent(1,1);
 345   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
 346
 347   // Cuts efficiency (from cuts variation)
 348   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 349   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 350   fCutsEff->SetBinContent(2,0.20);
 351   for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 352
 353
 354   // PID efficiency (from PID/noPID)
 355   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 356   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
 357   //  fPIDEff->SetBinContent(2,0.20);
 358   for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
 359
 360   // MC dN/dpt
 361   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 362   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 363   fMCPtShape->SetBinContent(1,0.03);
 364   fMCPtShape->SetBinContent(2,0.03);
 365
 366 //   // particle-antiparticle
 367 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 368 //   fPartAntipart->SetBinContent(1,1);
 369 //   fPartAntipart->SetBinContent(2,1);
 370 //   for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
 371 //   for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
 372
 373   return;
 374 }
 375
 376 //--------------------------------------------------------------------------
 377 void AliHFSystErr::InitDplustoKpipi2010pp() {
 378   //
 379   // D+->Kpipi syst errors. Responsible: R. Bala
 380   //  2010 pp sample
 381   //
 382
 383
 384 // Normalization
 385   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 386   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 387
 388   // Branching ratio
 389   fBR = new TH1F("fBR","fBR",24,0,24);
 390   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 391
 392   // Tracking efficiency
 393   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 394   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 395
 396
 397   // Raw yield extraction
 398   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 399   fRawYield->SetBinContent(1,1);
 400   fRawYield->SetBinContent(2,0.25);
 401   fRawYield->SetBinContent(3,0.25);
 402   fRawYield->SetBinContent(4,0.20);
 403   fRawYield->SetBinContent(5,0.09);
 404   fRawYield->SetBinContent(6,0.09);
 405   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 406   for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
 407
 408   // Cuts efficiency (from cuts variation)
 409   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 410   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 411
 412   // PID efficiency (from PID/noPID)
 413   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 414   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 415   fPIDEff->SetBinContent(1,0.15); // 15%
 416   fPIDEff->SetBinContent(2,0.15); // 15%
 417   fPIDEff->SetBinContent(3,0.15); // 15%
 418   fPIDEff->SetBinContent(4,0.15); // 15%
 419   for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
 420
 421   // MC dN/dpt
 422   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 423   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 424   fMCPtShape->SetBinContent(1,0.03);
 425   fMCPtShape->SetBinContent(2,0.03);
 426
 427
 428   // particle-antiparticle
 429   /*
 430   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 431   fPartAntipart->SetBinContent(1,1);
 432   fPartAntipart->SetBinContent(2,1);
 433   fPartAntipart->SetBinContent(3,0.12);
 434   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 435   */
 436   return;
 437 }
 438
 439 //--------------------------------------------------------------------------
 440 void AliHFSystErr::InitDstoKKpi2010pp() {
 441   //
 442   // D+s->KKpi syst errors. Responsible: G.M. Innocenti
 443   //  2010 pp sample
 444   //
 445
 446
 447 // Normalization
 448   fNorm = new TH1F("fNorm","fNorm",12,0,12);
 449   for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 7% error on sigmaV0and
 450
 451   // Branching ratio
 452   fBR = new TH1F("fBR","fBR",12,0,12);
 453   for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.05); // 5% PDG2010
 454
 455   // Tracking efficiency
 456   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
 457   for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.03); // 3% (1% per track)
 458
 459
 460   // Raw yield extraction
 461   fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
 462   fRawYield->SetBinContent(1,1);
 463   fRawYield->SetBinContent(2,1);
 464   fRawYield->SetBinContent(3,0.20);
 465   fRawYield->SetBinContent(4,0.20);
 466   fRawYield->SetBinContent(5,0.15);
 467   fRawYield->SetBinContent(6,0.15);
 468   fRawYield->SetBinContent(7,0.15);
 469   fRawYield->SetBinContent(8,0.15);
 470   fRawYield->SetBinContent(9,0.20);
 471   fRawYield->SetBinContent(10,0.20);
 472   fRawYield->SetBinContent(11,0.20);
 473   fRawYield->SetBinContent(12,0.20);
 474
 475   // Cuts efficiency (from cuts variation)
 476   fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
 477   for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
 478
 479   // PID efficiency (from PID/noPID)
 480   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
 481   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
 482
 483   // MC dN/dpt  (copied from D0 : will update later)
 484   //fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
 485   //for(Int_t i=1;i<=12;i++) fMCPtShape->SetBinContent(i,(Float_t)i*0.006);
 486
 487
 488   // particle-antiparticle
 489   /*
 490   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
 491   fPartAntipart->SetBinContent(1,1);
 492   fPartAntipart->SetBinContent(2,1);
 493   fPartAntipart->SetBinContent(3,0.12);
 494   for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 495   */
 496   return;
 497 }
 498
 499
 500 //--------------------------------------------------------------------------
 501 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
 502   //
 503   // D+->Kpipi syst errors. Responsible: ??
 504   //  2010 PbPb sample, 0-20 CC
 505   //
 506
 507  // Normalization
 508   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 509   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 510
 511   // Branching ratio
 512   fBR = new TH1F("fBR","fBR",20,0,20);
 513   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 514
 515   // Tracking efficiency
 516   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 517   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 518
 519   // Raw yield extraction
 520   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 521   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
 522   // fRawYield->SetBinContent(5,0.23);
 523   //fRawYield->SetBinContent(6,0.23);
 524   fRawYield->SetBinContent(7,0.20);
 525   fRawYield->SetBinContent(8,0.20);
 526   fRawYield->SetBinContent(9,0.15);
 527   fRawYield->SetBinContent(10,0.15);
 528   fRawYield->SetBinContent(11,0.15);
 529   fRawYield->SetBinContent(12,0.15);
 530
 531   // Cuts efficiency (from cuts variation)
 532   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 533   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 534
 535   // PID efficiency (from PID/noPID)
 536   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 537   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 538
 539   // MC dN/dpt  (2/2/2012)
 540   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 541   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
 542   for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
 543   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 544   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 545
 546
 547   // particle-antiparticle
 548   /*
 549   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 550   fPartAntipart->SetBinContent(1,1);
 551   fPartAntipart->SetBinContent(2,1);
 552   fPartAntipart->SetBinContent(3,0.12);
 553   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 554   */
 555
 556   return;
 557 }
 558
 559 //--------------------------------------------------------------------------
 560 void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
 561   //
 562   // D+->Kpipi syst errors. Responsible: ??
 563   //  2010 PbPb sample, 40-80 CC
 564   //
 565
 566
 567  // Normalization
 568   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 569   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 570
 571   // Branching ratio
 572   fBR = new TH1F("fBR","fBR",20,0,20);
 573   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
 574
 575   // Tracking efficiency
 576   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 577   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 578
 579
 580   // Raw yield extraction
 581   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 582   fRawYield->SetBinContent(1,1);
 583   fRawYield->SetBinContent(2,1);
 584   fRawYield->SetBinContent(3,1);
 585   fRawYield->SetBinContent(4,0.15);
 586   fRawYield->SetBinContent(5,0.05);
 587   fRawYield->SetBinContent(6,0.05);
 588   fRawYield->SetBinContent(7,0.15);
 589   fRawYield->SetBinContent(8,0.15);
 590   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
 591   for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1);  //5 to 10%
 592
 593   // Cuts efficiency (from cuts variation)
 594   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 595   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 596
 597   // PID efficiency (from PID/noPID)
 598   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 599   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 600   fPIDEff->SetBinContent(3,0.13); // 13%
 601
 602
 603   // MC dN/dpt  (2/2/2012)
 604   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 605   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
 606   for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
 607   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
 608   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
 609
 610
 611   // particle-antiparticle
 612   /*
 613   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 614   fPartAntipart->SetBinContent(1,1);
 615   fPartAntipart->SetBinContent(2,1);
 616   fPartAntipart->SetBinContent(3,0.12);
 617   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 618   */
 619   return;
 620 }
 621
 622 //--------------------------------------------------------------------------
 623 void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
 624
 625   //
 626   // D+->Kpipi syst errors. Responsible: R. Bala
 627   //  2011 2.76 TeV pp sample
 628   //
 629   AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
 630
 631   // Normalization
 632   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 633   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
 634
 635   // Branching ratio
 636   fBR = new TH1F("fBR","fBR",20,0,20);
 637   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 638
 639   // Tracking efficiency
 640   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 641   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
 642
 643   // Raw yield extraction
 644   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 645   fRawYield->SetBinContent(1,1);
 646   fRawYield->SetBinContent(2,1);
 647   for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10);  //5 to 10%
 648   fRawYield->SetBinContent(7,0.15);
 649   fRawYield->SetBinContent(8,0.15);
 650   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055);  //5 to 10%
 651
 652   // Cuts efficiency (from cuts variation)
 653   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 654   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 655
 656   // PID efficiency (from PID/noPID)
 657   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 658   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 659   fPIDEff->SetBinContent(3,0.10); // 13%
 660   fPIDEff->SetBinContent(4,0.10); // 13%
 661
 662   // MC dN/dpt  (copied from D0 : will update later)
 663   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 664   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 665   fMCPtShape->SetBinContent(1,0.03);
 666   fMCPtShape->SetBinContent(2,0.03);
 667
 668   return;
 669 }
 670
 671 //--------------------------------------------------------------------------
 672 void AliHFSystErr::InitDstartoD0pi2010pp() {
 673   //
 674   // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
 675   //  2010 pp sample
 676   //
 677
 678  // Normalization
 679   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 680   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 681
 682   // Branching ratio
 683   fBR = new TH1F("fBR","fBR",24,0,24);
 684   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 685
 686   // Tracking efficiency
 687   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 688   fTrackingEff->SetBinContent(1,1.0);
 689   fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
 690   fTrackingEff->SetBinContent(3,0.12);
 691   fTrackingEff->SetBinContent(3,0.12);
 692   for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 693
 694
 695   // Raw yield extraction
 696   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 697   fRawYield->SetBinContent(1,1.0);
 698   fRawYield->SetBinContent(2,0.10);
 699   fRawYield->SetBinContent(3,0.04);
 700   fRawYield->SetBinContent(4,0.03);
 701   fRawYield->SetBinContent(5,0.03);
 702   fRawYield->SetBinContent(6,0.05);
 703   fRawYield->SetBinContent(7,0.05);
 704   fRawYield->SetBinContent(8,0.05);
 705   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);  //4%
 706   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);  //4%
 707   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
 708
 709   // Cuts efficiency (from cuts variation)
 710   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 711   fCutsEff->SetBinContent(2,0.22);
 712   for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 713
 714   // PID efficiency (from PID/noPID)
 715   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 716   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
 717
 718
 719   // MC dN/dpt  (copied from D0 : will update later)
 720   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 721   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 722   fMCPtShape->SetBinContent(1,0.03);
 723   fMCPtShape->SetBinContent(2,0.03);
 724
 725   return;
 726
 727
 728 }
 729 //--------------------------------------------------------------------------
 730 void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 731
 732   //
 733   // D+->Kpipi syst errors. Responsible: A. Grelli
 734   //  2011 2.76 TeV pp sample
 735   //
 736   AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
 737
 738 // Normalization
 739   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 740   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
 741
 742   // Branching ratio
 743   fBR = new TH1F("fBR","fBR",20,0,20);
 744   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 745
 746   // Tracking efficiency
 747   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 748   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
 749
 750   // Raw yield extraction
 751   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 752   fRawYield->SetBinContent(1,1);
 753   fRawYield->SetBinContent(2,1);
 754   fRawYield->SetBinContent(3,0.14);
 755   fRawYield->SetBinContent(4,0.14);
 756   fRawYield->SetBinContent(5,0.12);
 757   fRawYield->SetBinContent(6,0.12);
 758   fRawYield->SetBinContent(7,0.06);
 759   fRawYield->SetBinContent(8,0.06);
 760   fRawYield->SetBinContent(9,0.08);
 761   fRawYield->SetBinContent(10,0.08);
 762   fRawYield->SetBinContent(11,0.08);
 763   fRawYield->SetBinContent(12,0.08);
 764   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
 765
 766   // Cuts efficiency (from cuts variation)
 767   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 768   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
 769   fCutsEff->SetBinContent(3,0.15);
 770   fCutsEff->SetBinContent(4,0.15);
 771   fCutsEff->SetBinContent(5,0.15);
 772   fCutsEff->SetBinContent(6,0.15);
 773   fCutsEff->SetBinContent(7,0.10);
 774   fCutsEff->SetBinContent(8,0.10);
 775   fCutsEff->SetBinContent(9,0.10);
 776   fCutsEff->SetBinContent(10,0.10);
 777   fCutsEff->SetBinContent(11,0.10);
 778   fCutsEff->SetBinContent(12,0.10);
 779
 780   // PID efficiency (from PID/noPID)
 781   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 782   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
 783
 784   // MC dN/dpt
 785   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 786   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 787   fMCPtShape->SetBinContent(1,0.03);
 788   fMCPtShape->SetBinContent(2,0.03);
 789
 790   return;
 791 }
 792
 793 //------------------------------------------------------------------------
 794 void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
 795   //
 796   // D*+->D0pi syst errors. Responsible: A. Grelli
 797   //  2010 PbPb sample, 0-20 CC
 798   //
 799
 800   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
 801
 802  // Normalization
 803   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 804   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 805
 806   // Branching ratio
 807   fBR = new TH1F("fBR","fBR",24,0,24);
 808   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 809
 810   // Tracking efficiency
 811   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
 812   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 813
 814
 815   // Raw yield extraction
 816   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 817   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
 818   fRawYield->SetBinContent(3,0.2);
 819   fRawYield->SetBinContent(4,0.2);
 820   fRawYield->SetBinContent(5,0.2);
 821   fRawYield->SetBinContent(6,0.2);
 822
 823   // Cuts efficiency (from cuts variation)
 824   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 825   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 826   fCutsEff->SetBinContent(4,0.15);
 827   fCutsEff->SetBinContent(5,0.15);
 828   fCutsEff->SetBinContent(6,0.15);
 829
 830   // PID efficiency (from PID/noPID)
 831   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 832   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
 833
 834
 835   // MC dN/dpt  (from study on D* pt shape)
 836   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 837   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
 838   fMCPtShape->SetBinContent(4,0.025);
 839   fMCPtShape->SetBinContent(5,0.025);
 840   fMCPtShape->SetBinContent(6,0.025);
 841   fMCPtShape->SetBinContent(7,0.04);
 842   fMCPtShape->SetBinContent(8,0.04);
 843   fMCPtShape->SetBinContent(9,0.03);
 844   fMCPtShape->SetBinContent(10,0.03);
 845   fMCPtShape->SetBinContent(11,0.03);
 846   fMCPtShape->SetBinContent(12,0.03);
 847
 848
 849
 850   return;
 851
 852 }
 853
 854 //-------------------------------------------------------------------------
 855 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
 856   //
 857   // D*+->D0pi syst errors. Responsible: A. Grelli
 858   //  2010 PbPb sample, 20-40 CC
 859   //
 860
 861   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
 862
 863  // Normalization
 864   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 865   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
 866
 867   // Branching ratio
 868   fBR = new TH1F("fBR","fBR",24,0,24);
 869   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 870
 871   // Tracking efficiency
 872   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
 873   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 874
 875
 876   // Raw yield extraction
 877   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 878   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15);  //4%
 879
 880   // Cuts efficiency (from cuts variation)
 881   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 882   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 883
 884   // PID efficiency (from PID/noPID)
 885   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 886   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
 887
 888
 889   // MC dN/dpt  (copied from D0 : will update later)
 890   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 891   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
 892   fMCPtShape->SetBinContent(1,0.03);
 893   fMCPtShape->SetBinContent(2,0.03);
 894
 895   return;
 896
 897 }
 898
 899 //--------------------------------------------------------------------------
 900 void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 901   //
 902   // D*+->D0pi syst errors. Responsible: A. Grelli
 903   //  2010 PbPb sample, 40-80 CC
 904   //
 905
 906   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
 907
 908  // Normalization
 909   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 910   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 911
 912   // Branching ratio
 913   fBR = new TH1F("fBR","fBR",24,0,24);
 914   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 915
 916   // Tracking efficiency
 917   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
 918   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 919
 920
 921   // Raw yield extraction
 922   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 923   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
 924   fRawYield->SetBinContent(1,1);
 925   fRawYield->SetBinContent(2,0.15);
 926   fRawYield->SetBinContent(3,0.15);
 927   fRawYield->SetBinContent(4,0.15);
 928   fRawYield->SetBinContent(5,0.15);
 929   fRawYield->SetBinContent(6,0.10);
 930   fRawYield->SetBinContent(7,0.10);
 931   fRawYield->SetBinContent(8,0.10);
 932   fRawYield->SetBinContent(9,0.11);
 933   fRawYield->SetBinContent(10,0.11);
 934   fRawYield->SetBinContent(11,0.11);
 935   fRawYield->SetBinContent(12,0.11);
 936   fRawYield->SetBinContent(13,0.08);
 937   fRawYield->SetBinContent(14,0.08);
 938   fRawYield->SetBinContent(15,0.08);
 939   fRawYield->SetBinContent(16,0.08);
 940
 941
 942   // Cuts efficiency (from cuts variation)
 943   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 944   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 945
 946   // PID efficiency (from PID/noPID)
 947   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 948   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
 949
 950
 951   // MC dN/dpt  (copied from D0 : will update later)
 952   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 953   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
 954   fMCPtShape->SetBinContent(2,0.05);
 955   fMCPtShape->SetBinContent(3,0.05);
 956   fMCPtShape->SetBinContent(4,0.05);
 957   fMCPtShape->SetBinContent(5,0.04);
 958   fMCPtShape->SetBinContent(6,0.02);
 959   fMCPtShape->SetBinContent(7,0.04);
 960   fMCPtShape->SetBinContent(8,0.04);
 961
 962   return;
 963
 964 }
 965
 966 //--------------------------------------------------------------------------
 967 void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
 968   // define errors for RAA vs. centrality
 969   InitD0toKpi2010PbPb020();
 970   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 971   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
 972   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 973 }
 974 //--------------------------------------------------------------------------
 975 void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
 976   // define errors for RAA vs. centrality
 977   InitD0toKpi2010PbPb020();
 978   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 979   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);
 980   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 981 }
 982 //--------------------------------------------------------------------------
 983 void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
 984   // define errors for RAA vs. centrality
 985   InitD0toKpi2010PbPb4080();
 986   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 987   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
 988   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 989 }
 990 //--------------------------------------------------------------------------
 991 void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
 992    // define errors for RAA vs. centrality
 993   InitD0toKpi2010PbPb4080();
 994   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
 995   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);
 996   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 997 }
 998 //--------------------------------------------------------------------------
 999 void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
1000    // define errors for RAA vs. centrality
1001   InitD0toKpi2010PbPb4080();
1002   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1003   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
1004   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1005 }
1006 //--------------------------------------------------------------------------
1007 void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
1008   // define errors for RAA vs. centrality
1009   InitDplustoKpipi2010PbPb020();
1010   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
1011   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
1012
1013 }
1014 //--------------------------------------------------------------------------
1015 void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
1016   // define errors for RAA vs. centrality
1017   InitDplustoKpipi2010PbPb020();
1018   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
1019   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1020 }
1021 //--------------------------------------------------------------------------
1022 void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
1023   // define errors for RAA vs. centrality
1024   InitDplustoKpipi2010PbPb020();
1025   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1026   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
1027 }
1028 //--------------------------------------------------------------------------
1029 void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
1030   // define errors for RAA vs. centrality
1031   InitDplustoKpipi2010PbPb4080();
1032   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1033   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1034 }
1035 //--------------------------------------------------------------------------
1036 void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
1037   // define errors for RAA vs. centrality
1038   InitDplustoKpipi2010PbPb4080();
1039   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
1040   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
1041 }
1042
1043 //--------------------------------------------------------------------------
1044 void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
1045   // define errors for RAA vs. centrality
1046   InitDstartoD0pi2010PbPb020();
1047   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
1048   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
1049 }
1050 //--------------------------------------------------------------------------
1051 void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
1052   // define errors for RAA vs. centrality
1053   InitDstartoD0pi2010PbPb020();
1054   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1055   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
1056 }
1057 //--------------------------------------------------------------------------
1058 void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
1059   // define errors for RAA vs. centrality
1060   InitDstartoD0pi2010PbPb2040();
1061   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1062   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1063 }
1064 //--------------------------------------------------------------------------
1065 void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
1066   // define errors for RAA vs. centrality
1067   InitDstartoD0pi2010PbPb4080();
1068   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1069   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
1070 }
1071 //--------------------------------------------------------------------------
1072 void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
1073   // define errors for RAA vs. centrality
1074   InitDstartoD0pi2010PbPb4080();
1075   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1076   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
1077 }
1078
1079
1080 //--------------------------------------------------------------------------
1081 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
1082   //
1083   // Get error
1084   //
1085
1086   Int_t bin=fCutsEff->FindBin(pt);
1087
1088   return fCutsEff->GetBinContent(bin);
1089 }
1090 //--------------------------------------------------------------------------
1091 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
1092   //
1093   // Get error
1094   //
1095
1096   Int_t bin=fMCPtShape->FindBin(pt);
1097
1098   return fMCPtShape->GetBinContent(bin);
1099 }
1100 //--------------------------------------------------------------------------
1101 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
1102   //
1103   // Get error
1104   //
1105
1106   Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
1107
1108   return TMath::Sqrt(err);
1109 }
1110 //--------------------------------------------------------------------------
1111 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
1112   //
1113   // Get error
1114   //
1115
1116   Int_t bin=fPIDEff->FindBin(pt);
1117
1118   return fPIDEff->GetBinContent(bin);
1119 }
1120 //--------------------------------------------------------------------------
1121 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
1122   //
1123   // Get error
1124   //
1125
1126   Int_t bin=fTrackingEff->FindBin(pt);
1127
1128   return fTrackingEff->GetBinContent(bin);
1129 }
1130 //--------------------------------------------------------------------------
1131 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
1132   //
1133   // Get error
1134   //
1135
1136   Int_t bin=fRawYield->FindBin(pt);
1137
1138   return fRawYield->GetBinContent(bin);
1139 }
1140 //--------------------------------------------------------------------------
1141 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
1142   //
1143   // Get error
1144   //
1145
1146   Int_t bin=fPartAntipart->FindBin(pt);
1147
1148   return fPartAntipart->GetBinContent(bin);
1149 }
1150 //--------------------------------------------------------------------------
1151 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
1152   //
1153   // Get total syst error (except norm. error)
1154   //
1155
1156   Double_t err=0.;
1157
1158   if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
1159   if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
1160   //  if(fBR) err += GetBRErr()*GetBRErr();
1161   if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
1162   if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
1163   if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
1164   if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
1165
1166   err += feeddownErr*feeddownErr;
1167
1168   return TMath::Sqrt(err);
1169 }
1170 //---------------------------------------------------------------------------
1171 void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
1172   //
1173   // Draw errors
1174   //
1175   gStyle->SetOptStat(0);
1176
1177   TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
1178   cSystErr->Range(0.20,-0.5,18.4,0.34);
1179   cSystErr->SetRightMargin(0.318);
1180   cSystErr->SetFillColor(0);
1181
1182   TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",30,0,30,100,-1,+1);
1183   hFrame->SetAxisRange(1.,24.9,"X");
1184   hFrame->SetAxisRange(-0.5,0.5,"Y");
1185   hFrame->Draw();
1186
1187   TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
1188   leg->SetTextSize(0.03601695);
1189   leg->SetFillStyle(0);
1190   leg->SetBorderSize(0);
1191
1192   TH1F *hTotErr=new TH1F("hTotErr","",24,0,24);
1193   Int_t nbins = fNorm->GetNbinsX();
1194   TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
1195   for(Int_t i=1;i<=24;i++) {
1196     Double_t pt = hTotErr->GetBinCenter(i);
1197     Double_t ptwidth = hTotErr->GetBinWidth(i);
1198
1199     if(grErrFeeddown) {
1200       Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
1201       Double_t toterryl=0., toterryh=0.;
1202       for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
1203         grErrFeeddown->GetPoint(j,x,y);
1204         errxh = grErrFeeddown->GetErrorXhigh(j);
1205         errxl = grErrFeeddown->GetErrorXlow(j);
1206         if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
1207           erryh = grErrFeeddown->GetErrorYhigh(j);
1208           erryl = grErrFeeddown->GetErrorYlow(j);
1209         }
1210       }
1211       if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
1212       else toterryl = GetTotalSystErr(pt);
1213       if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
1214       else toterryh = GetTotalSystErr(pt);
1215
1216       hTotErr->SetBinContent(i,toterryh);
1217       gTotErr->SetPoint(i,pt,0.);
1218       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
1219     }
1220     else {
1221       hTotErr->SetBinContent(i,GetTotalSystErr(pt));
1222       gTotErr->SetPoint(i,pt,0.);
1223       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
1224     }
1225
1226   }
1227   gTotErr->SetLineColor(kBlack);
1228   gTotErr->SetFillColor(kRed);
1229   gTotErr->SetFillStyle(3002);
1230   gTotErr->Draw("2");
1231   leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
1232 //   hTotErr->SetLineColor(1);
1233 //   hTotErr->SetLineWidth(3);
1234 //   hTotErr->Draw("same");
1235 //   leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
1236
1237
1238   fNorm->SetFillColor(1);
1239   fNorm->SetFillStyle(3002);
1240   //fNorm->Draw("same");
1241   //TH1F *hNormRefl = ReflectHisto(fNorm);
1242   //hNormRefl->Draw("same");
1243   Double_t norm = fNorm->GetBinContent(1)*100;
1244   leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
1245
1246   if(grErrFeeddown) {
1247     grErrFeeddown->SetFillColor(kTeal-8);
1248     grErrFeeddown->SetFillStyle(3001);
1249     grErrFeeddown->Draw("2");
1250     leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
1251   }
1252   if(fTrackingEff) {
1253     fTrackingEff->SetFillColor(4);
1254     fTrackingEff->SetFillStyle(3006);
1255     fTrackingEff->Draw("same");
1256     TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
1257     hTrackingEffRefl->Draw("same");
1258     leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
1259   }
1260   if(fBR) {
1261     fBR->SetFillColor(6);
1262     fBR->SetFillStyle(3005);
1263     //fBR->SetFillStyle(3020);
1264     fBR->Draw("same");
1265     TH1F *hBRRefl = ReflectHisto(fBR);
1266     hBRRefl->Draw("same");
1267     leg->AddEntry(fBR,"Branching ratio","f");
1268   }
1269   if(fRawYield) {
1270     Int_t ci;   // for color index setting
1271     ci = TColor::GetColor("#00cc00");
1272     fRawYield->SetLineColor(ci);
1273     //    fRawYield->SetLineColor(3);
1274     fRawYield->SetLineWidth(3);
1275     fRawYield->Draw("same");
1276     TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
1277     hRawYieldRefl->Draw("same");
1278     leg->AddEntry(fRawYield,"Yield extraction","l");
1279   }
1280   if(fCutsEff) {
1281     fCutsEff->SetLineColor(4);
1282     fCutsEff->SetLineWidth(3);
1283     fCutsEff->Draw("same");
1284     TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
1285     hCutsEffRefl->Draw("same");
1286     leg->AddEntry(fCutsEff,"Cuts efficiency","l");
1287   }
1288   if(fPIDEff) {
1289     fPIDEff->SetLineColor(7);
1290     fPIDEff->SetLineWidth(3);
1291     fPIDEff->Draw("same");
1292     TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
1293     hPIDEffRefl->Draw("same");
1294     leg->AddEntry(fPIDEff,"PID efficiency","l");
1295   }
1296   if(fMCPtShape) {
1297     Int_t ci = TColor::GetColor("#9933ff");
1298     fMCPtShape->SetLineColor(ci);
1299     //    fMCPtShape->SetLineColor(8);
1300     fMCPtShape->SetLineWidth(3);
1301     fMCPtShape->Draw("same");
1302     TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
1303     hMCPtShapeRefl->Draw("same");
1304     leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
1305   }
1306   if(fPartAntipart) {
1307     Int_t ci = TColor::GetColor("#ff6600");
1308     fPartAntipart->SetLineColor(ci);
1309     //    fPartAntipart->SetLineColor(9);
1310     fPartAntipart->SetLineWidth(3);
1311     fPartAntipart->Draw("same");
1312     TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
1313     hPartAntipartRefl->Draw("same");
1314     leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
1315   }
1316
1317
1318   leg->Draw();
1319
1320   cSystErr->SaveAs("RelativeSystematics.eps");
1321
1322   return;
1323 }
1324 //-------------------------------------------------------------------------
1325 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
1326   //
1327   // Clones and reflects histogram
1328   //
1329   TH1F *hout=(TH1F*)hin->Clone("hout");
1330   hout->Scale(-1.);
1331
1332   return hout;
1333 }