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=="07half") InitD0toKpi2011PbPb07half();
 105       else if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
 106       else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
 107       else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
 108       else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
 109       else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
 110       else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
 111       else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
 112       else AliFatal("Not yet implemented");
 113     }
 114     //    else if (fCollisionType==2) InitD0toKpi2010ppLowEn();
 115     break;
 116   case 2: // D+->Kpipi
 117     if (fCollisionType==0) {
 118       if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
 119       else InitDplustoKpipi2010pp();
 120     } else if (fCollisionType==1) {
 121       if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
 122       else if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
 123       else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
 124       else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
 125       else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
 126       else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
 127       else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
 128       else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
 129       else AliFatal("Not yet implemented");
 130     }
 131     break;
 132   case 3: // D*->D0pi
 133     if (fCollisionType==0) {
 134       if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
 135       else InitDstartoD0pi2010pp();
 136     }else if (fCollisionType==1) {
 137       if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
 138       else if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
 139       else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
 140       else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
 141       else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
 142       else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
 143       else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
 144       else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
 145       else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
 146       else AliFatal("Not yet implemented");
 147     }
 148     break;
 149   case 4: // D+s->KKpi
 150     if (fCollisionType==0) InitDstoKKpi2010pp();
 151     else AliFatal("Not yet implemented");
 152     break;
 153
 154   default:
 155     printf("Invalid decay type: %d\n",decay);
 156     break;
 157   }
 158
 159 }
 160
 161 //--------------------------------------------------------------------------
 162 void AliHFSystErr::InitD0toKpi2010pp() {
 163   //
 164   // D0->Kpi syst errors. Responsible: A. Rossi
 165   //   2010 pp sample
 166   //
 167
 168   // Normalization
 169   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 170   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 171
 172   // Branching ratio
 173   fBR = new TH1F("fBR","fBR",24,0,24);
 174   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 175
 176   // Tracking efficiency
 177   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 178   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
 179
 180   // Raw yield extraction
 181   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 182   fRawYield->SetBinContent(1,1);
 183   fRawYield->SetBinContent(2,0.22);
 184   fRawYield->SetBinContent(3,0.1);
 185   for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
 186   for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
 187   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 188   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
 189
 190   // Cuts efficiency (from cuts variation)
 191   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 192   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 193
 194   // PID efficiency (from PID/noPID)
 195   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 196   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
 197   fPIDEff->SetBinContent(2,0.05); // 5%
 198
 199   // MC dN/dpt
 200   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 201   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 202   fMCPtShape->SetBinContent(1,0.03);
 203   fMCPtShape->SetBinContent(2,0.03);
 204
 205   // particle-antiparticle
 206   //  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
 207   //  fPartAntipart->SetBinContent(1,1);
 208   //  for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
 209
 210   return;
 211 }
 212 //--------------------------------------------------------------------------
 213 void AliHFSystErr::InitD0toKpi2010PbPb020() {
 214   //
 215   // D0->Kpi syst errors. Responsible: A. Rossi
 216   //   2010 PbPb sample, 0-20 CC
 217   //
 218
 219   // Normalization
 220   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 221   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 222
 223   // Branching ratio
 224   fBR = new TH1F("fBR","fBR",20,0,20);
 225   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 226
 227   // Tracking efficiency
 228   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 229   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 230
 231   // Raw yield extraction
 232   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 233   fRawYield->SetBinContent(1,0);
 234   fRawYield->SetBinContent(2,0);
 235   fRawYield->SetBinContent(3,0.08);
 236   for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
 237   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
 238
 239   // Cuts efficiency (from cuts variation)
 240   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 241   fCutsEff->SetBinContent(1,0.);
 242   fCutsEff->SetBinContent(2,0.);
 243   fCutsEff->SetBinContent(3,0.13);
 244   fCutsEff->SetBinContent(4,0.11);
 245   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 246   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 247
 248   // PID efficiency (from PID/noPID)
 249   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 250   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 251
 252   // MC dN/dpt
 253   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 254   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 255   fMCPtShape->SetBinContent(3,0.04);
 256   fMCPtShape->SetBinContent(4,0.02);
 257   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 258
 259 //   // particle-antiparticle
 260 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 261 //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 262 //   fPartAntipart->SetBinContent(3,0.10);
 263 //   fPartAntipart->SetBinContent(4,0.10);
 264 //   fPartAntipart->SetBinContent(7,0.10);
 265 //   fPartAntipart->SetBinContent(8,0.10);
 266
 267   return;
 268 }
 269 //--------------------------------------------------------------------------
 270 //--------------------------------------------------------------------------
 271 void AliHFSystErr::InitD0toKpi2011PbPb07half() {
 272   //
 273   // D0->Kpi syst errors. Responsible: A. Rossi
 274   //   2011 PbPb sample, 0-7.5 CC
 275   //
 276
 277   // Normalization
 278   fNorm = new TH1F("fNorm","fNorm",36,0,36);
 279   for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 280
 281   // Branching ratio
 282   fBR = new TH1F("fBR","fBR",36,0,36);
 283   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 284
 285   // Tracking efficiency
 286   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
 287   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 288   for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 289
 290   // Raw yield extraction
 291   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 292   fRawYield->SetBinContent(1,0);
 293   fRawYield->SetBinContent(2,0.1);
 294   fRawYield->SetBinContent(3,0.05);
 295   for(Int_t i=4;i<=16;i++) fRawYield->SetBinContent(i,0.05);
 296   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
 297   for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
 298
 299   // Cuts efficiency (from cuts variation)
 300   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 301   fCutsEff->SetBinContent(1,0.);
 302   fCutsEff->SetBinContent(2,0.15);
 303   fCutsEff->SetBinContent(3,0.13);
 304   fCutsEff->SetBinContent(4,0.11);
 305   fCutsEff->SetBinContent(5,0.08);
 306   for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
 307   for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 308
 309   // PID efficiency (from PID/noPID)
 310   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
 311   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
 312   for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 313
 314   // MC dN/dpt
 315   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
 316   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
 317   fMCPtShape->SetBinContent(2,0.06);
 318   fMCPtShape->SetBinContent(3,0.04);
 319   fMCPtShape->SetBinContent(4,0.02);
 320   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 321   for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
 322   for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
 323
 324 //   // particle-antiparticle
 325 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
 326 //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 327 //   fPartAntipart->SetBinContent(3,0.10);
 328 //   fPartAntipart->SetBinContent(4,0.10);
 329 //   fPartAntipart->SetBinContent(7,0.10);
 330 //   fPartAntipart->SetBinContent(8,0.10);
 331
 332   return;
 333 }
 334 //--------------------------------------------------------------------------
 335 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
 336   //
 337   // D0->Kpi syst errors. Responsible: A. Rossi
 338   //   2010 PbPb sample, 40-80 CC
 339   //
 340
 341   // Normalization
 342   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 343   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 344
 345   // Branching ratio
 346   fBR = new TH1F("fBR","fBR",20,0,20);
 347   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 348
 349   // Tracking efficiency
 350   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 351   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 352
 353
 354   // Raw yield extraction
 355   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 356   fRawYield->SetBinContent(1,0);
 357   fRawYield->SetBinContent(2,0);
 358   for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
 359   //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
 360
 361   // Cuts efficiency (from cuts variation)
 362   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 363   fCutsEff->SetBinContent(1,0.);
 364   fCutsEff->SetBinContent(2,0.);
 365   fCutsEff->SetBinContent(3,0.13);
 366   fCutsEff->SetBinContent(4,0.11);
 367   for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
 368   for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
 369
 370   // PID efficiency (from PID/noPID)
 371   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 372 //   for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
 373 //   for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 374   for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
 375
 376   // MC dN/dpt
 377   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 378   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 379   //  fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
 380   //  fMCPtShape->SetBinContent(4,0.02);
 381   for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
 382
 383   //   // particle-antiparticle
 384   //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 385   //   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
 386
 387   return;
 388 }
 389
 390 //--------------------------------------------------------------------------
 391 void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 392   //
 393   // D0->Kpi syst errors. Low energy run
 394   //   2011 2.76 TeV pp sample
 395   //
 396   AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
 397
 398   // Normalization
 399   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 400   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
 401
 402   // Branching ratio
 403   fBR = new TH1F("fBR","fBR",20,0,20);
 404   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
 405
 406   // Tracking efficiency
 407   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 408   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
 409
 410   // Raw yield extraction
 411   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 412   fRawYield->SetBinContent(1,1);
 413   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
 414
 415   // Cuts efficiency (from cuts variation)
 416   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 417   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 418   fCutsEff->SetBinContent(2,0.20);
 419   for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 420
 421
 422   // PID efficiency (from PID/noPID)
 423   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 424   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
 425   //  fPIDEff->SetBinContent(2,0.20);
 426   for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
 427
 428   // MC dN/dpt
 429   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 430   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 431   fMCPtShape->SetBinContent(1,0.03);
 432   fMCPtShape->SetBinContent(2,0.03);
 433
 434 //   // particle-antiparticle
 435 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 436 //   fPartAntipart->SetBinContent(1,1);
 437 //   fPartAntipart->SetBinContent(2,1);
 438 //   for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
 439 //   for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
 440
 441   return;
 442 }
 443
 444 //--------------------------------------------------------------------------
 445 void AliHFSystErr::InitDplustoKpipi2010pp() {
 446   //
 447   // D+->Kpipi syst errors. Responsible: R. Bala
 448   //  2010 pp sample
 449   //
 450
 451
 452 // Normalization
 453   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 454   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 455
 456   // Branching ratio
 457   fBR = new TH1F("fBR","fBR",24,0,24);
 458   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 459
 460   // Tracking efficiency
 461   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 462   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 463
 464
 465   // Raw yield extraction
 466   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 467   fRawYield->SetBinContent(1,1);
 468   fRawYield->SetBinContent(2,0.25);
 469   fRawYield->SetBinContent(3,0.25);
 470   fRawYield->SetBinContent(4,0.20);
 471   fRawYield->SetBinContent(5,0.09);
 472   fRawYield->SetBinContent(6,0.09);
 473   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
 474   for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
 475
 476   // Cuts efficiency (from cuts variation)
 477   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 478   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 479
 480   // PID efficiency (from PID/noPID)
 481   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 482   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 483   fPIDEff->SetBinContent(1,0.15); // 15%
 484   fPIDEff->SetBinContent(2,0.15); // 15%
 485   fPIDEff->SetBinContent(3,0.15); // 15%
 486   fPIDEff->SetBinContent(4,0.15); // 15%
 487   for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
 488
 489   // MC dN/dpt
 490   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 491   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 492   fMCPtShape->SetBinContent(1,0.03);
 493   fMCPtShape->SetBinContent(2,0.03);
 494
 495
 496   // particle-antiparticle
 497   /*
 498   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 499   fPartAntipart->SetBinContent(1,1);
 500   fPartAntipart->SetBinContent(2,1);
 501   fPartAntipart->SetBinContent(3,0.12);
 502   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 503   */
 504   return;
 505 }
 506
 507 //--------------------------------------------------------------------------
 508 void AliHFSystErr::InitDstoKKpi2010pp() {
 509   //
 510   // D+s->KKpi syst errors. Responsible: G.M. Innocenti
 511   //  2010 pp sample
 512   //
 513
 514
 515 // Normalization
 516   fNorm = new TH1F("fNorm","fNorm",12,0,12);
 517   for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
 518
 519   // Branching ratio
 520   fBR = new TH1F("fBR","fBR",12,0,12);
 521   for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
 522
 523   // Tracking efficiency
 524   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
 525   for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 526
 527
 528   // Raw yield extraction
 529   fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
 530   fRawYield->SetBinContent(1,1);
 531   fRawYield->SetBinContent(2,1);
 532   fRawYield->SetBinContent(3,0.20);
 533   fRawYield->SetBinContent(4,0.20);
 534   fRawYield->SetBinContent(5,0.15);
 535   fRawYield->SetBinContent(6,0.15);
 536   fRawYield->SetBinContent(7,0.15);
 537   fRawYield->SetBinContent(8,0.15);
 538   fRawYield->SetBinContent(9,0.20);
 539   fRawYield->SetBinContent(10,0.20);
 540   fRawYield->SetBinContent(11,0.20);
 541   fRawYield->SetBinContent(12,0.20);
 542
 543   // Cuts efficiency (from cuts variation)
 544   fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
 545   for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
 546
 547   // PID efficiency (from PID/noPID)
 548   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
 549   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
 550
 551   // MC dN/dpt
 552   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
 553   for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
 554   for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
 555   for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
 556   for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
 557   for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
 558
 559
 560   // particle-antiparticle
 561   /*
 562   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
 563   fPartAntipart->SetBinContent(1,1);
 564   fPartAntipart->SetBinContent(2,1);
 565   fPartAntipart->SetBinContent(3,0.12);
 566   for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 567   */
 568   return;
 569 }
 570
 571
 572 //--------------------------------------------------------------------------
 573 void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
 574   //
 575   // D+->Kpipi syst errors. Responsible: E. Bruna
 576   //  2011 PbPb sample, 0-7.5% CC
 577   //
 578
 579  // Normalization
 580   fNorm = new TH1F("fNorm","fNorm",36,0,36);
 581   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 582
 583   // Branching ratio
 584   fBR = new TH1F("fBR","fBR",36,0,36);
 585   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 586
 587   // Tracking efficiency
 588   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
 589   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 590
 591   // Raw yield extraction
 592   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
 593   for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
 594   fRawYield->SetBinContent(4,0.20);
 595   fRawYield->SetBinContent(5,0.15);
 596   fRawYield->SetBinContent(6,0.15);
 597   for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.12);
 598   for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.10);
 599   for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
 600   for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.07);
 601   for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
 602
 603
 604   // Cuts efficiency (from cuts variation)
 605   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
 606   for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.15); // 15%
 607
 608   // PID efficiency (from PID/noPID)
 609   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
 610   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 611
 612   // MC dN/dpt  (2/2/2012)
 613   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
 614   for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.);
 615   for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
 616   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 617   for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 618
 619
 620   // particle-antiparticle
 621   /*
 622   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 623   fPartAntipart->SetBinContent(1,1);
 624   fPartAntipart->SetBinContent(2,1);
 625   fPartAntipart->SetBinContent(3,0.12);
 626   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 627   */
 628
 629   return;
 630 }
 631
 632 //--------------------------------------------------------------------------
 633 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
 634   //
 635   // D+->Kpipi syst errors. Responsible: ??
 636   //  2010 PbPb sample, 0-20 CC
 637   //
 638
 639  // Normalization
 640   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 641   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 642
 643   // Branching ratio
 644   fBR = new TH1F("fBR","fBR",20,0,20);
 645   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 646
 647   // Tracking efficiency
 648   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 649   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 650
 651   // Raw yield extraction
 652   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 653   for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
 654   // fRawYield->SetBinContent(5,0.23);
 655   //fRawYield->SetBinContent(6,0.23);
 656   fRawYield->SetBinContent(7,0.20);
 657   fRawYield->SetBinContent(8,0.20);
 658   fRawYield->SetBinContent(9,0.15);
 659   fRawYield->SetBinContent(10,0.15);
 660   fRawYield->SetBinContent(11,0.15);
 661   fRawYield->SetBinContent(12,0.15);
 662
 663   // Cuts efficiency (from cuts variation)
 664   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 665   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 666
 667   // PID efficiency (from PID/noPID)
 668   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 669   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 670
 671   // MC dN/dpt  (2/2/2012)
 672   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 673   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
 674   for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
 675   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 676   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
 677
 678
 679   // particle-antiparticle
 680   /*
 681   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 682   fPartAntipart->SetBinContent(1,1);
 683   fPartAntipart->SetBinContent(2,1);
 684   fPartAntipart->SetBinContent(3,0.12);
 685   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 686   */
 687
 688   return;
 689 }
 690
 691 //--------------------------------------------------------------------------
 692 void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
 693   //
 694   // D+->Kpipi syst errors. Responsible: ??
 695   //  2010 PbPb sample, 40-80 CC
 696   //
 697
 698
 699  // Normalization
 700   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 701   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 702
 703   // Branching ratio
 704   fBR = new TH1F("fBR","fBR",20,0,20);
 705   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
 706
 707   // Tracking efficiency
 708   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 709   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 710
 711
 712   // Raw yield extraction
 713   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 714   fRawYield->SetBinContent(1,1);
 715   fRawYield->SetBinContent(2,1);
 716   fRawYield->SetBinContent(3,1);
 717   fRawYield->SetBinContent(4,0.15);
 718   fRawYield->SetBinContent(5,0.05);
 719   fRawYield->SetBinContent(6,0.05);
 720   fRawYield->SetBinContent(7,0.15);
 721   fRawYield->SetBinContent(8,0.15);
 722   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
 723   for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1);  //5 to 10%
 724
 725   // Cuts efficiency (from cuts variation)
 726   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 727   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 728
 729   // PID efficiency (from PID/noPID)
 730   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 731   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 732   fPIDEff->SetBinContent(3,0.13); // 13%
 733
 734
 735   // MC dN/dpt  (2/2/2012)
 736   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 737   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
 738   for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
 739   for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
 740   for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
 741
 742
 743   // particle-antiparticle
 744   /*
 745   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
 746   fPartAntipart->SetBinContent(1,1);
 747   fPartAntipart->SetBinContent(2,1);
 748   fPartAntipart->SetBinContent(3,0.12);
 749   for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
 750   */
 751   return;
 752 }
 753
 754 //--------------------------------------------------------------------------
 755 void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
 756
 757   //
 758   // D+->Kpipi syst errors. Responsible: R. Bala
 759   //  2011 2.76 TeV pp sample
 760   //
 761   AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
 762
 763   // Normalization
 764   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 765   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
 766
 767   // Branching ratio
 768   fBR = new TH1F("fBR","fBR",20,0,20);
 769   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
 770
 771   // Tracking efficiency
 772   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 773   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
 774
 775   // Raw yield extraction
 776   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 777   fRawYield->SetBinContent(1,1);
 778   fRawYield->SetBinContent(2,1);
 779   for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10);  //5 to 10%
 780   fRawYield->SetBinContent(7,0.15);
 781   fRawYield->SetBinContent(8,0.15);
 782   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055);  //5 to 10%
 783
 784   // Cuts efficiency (from cuts variation)
 785   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 786   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
 787
 788   // PID efficiency (from PID/noPID)
 789   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 790   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 791   fPIDEff->SetBinContent(3,0.10); // 13%
 792   fPIDEff->SetBinContent(4,0.10); // 13%
 793
 794   // MC dN/dpt  (copied from D0 : will update later)
 795   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 796   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 797   fMCPtShape->SetBinContent(1,0.03);
 798   fMCPtShape->SetBinContent(2,0.03);
 799
 800   return;
 801 }
 802
 803 //--------------------------------------------------------------------------
 804 void AliHFSystErr::InitDstartoD0pi2010pp() {
 805   //
 806   // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
 807   //  2010 pp sample
 808   //
 809
 810  // Normalization
 811   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 812   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 813
 814   // Branching ratio
 815   fBR = new TH1F("fBR","fBR",24,0,24);
 816   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 817
 818   // Tracking efficiency
 819   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
 820   fTrackingEff->SetBinContent(1,1.0);
 821   fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
 822   fTrackingEff->SetBinContent(3,0.12);
 823   fTrackingEff->SetBinContent(3,0.12);
 824   for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 825
 826
 827   // Raw yield extraction
 828   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 829   fRawYield->SetBinContent(1,1.0);
 830   fRawYield->SetBinContent(2,0.10);
 831   fRawYield->SetBinContent(3,0.04);
 832   fRawYield->SetBinContent(4,0.03);
 833   fRawYield->SetBinContent(5,0.03);
 834   fRawYield->SetBinContent(6,0.05);
 835   fRawYield->SetBinContent(7,0.05);
 836   fRawYield->SetBinContent(8,0.05);
 837   for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);  //4%
 838   for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);  //4%
 839   for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
 840
 841   // Cuts efficiency (from cuts variation)
 842   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 843   fCutsEff->SetBinContent(2,0.22);
 844   for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 845
 846   // PID efficiency (from PID/noPID)
 847   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 848   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
 849
 850
 851   // MC dN/dpt  (copied from D0 : will update later)
 852   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 853   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
 854   fMCPtShape->SetBinContent(1,0.03);
 855   fMCPtShape->SetBinContent(2,0.03);
 856
 857   return;
 858
 859
 860 }
 861 //--------------------------------------------------------------------------
 862 void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 863
 864   //
 865   // D+->Kpipi syst errors. Responsible: A. Grelli
 866   //  2011 2.76 TeV pp sample
 867   //
 868   AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
 869
 870 // Normalization
 871   fNorm = new TH1F("fNorm","fNorm",20,0,20);
 872   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
 873
 874   // Branching ratio
 875   fBR = new TH1F("fBR","fBR",20,0,20);
 876   for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 877
 878   // Tracking efficiency
 879   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
 880   for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
 881
 882   // Raw yield extraction
 883   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
 884   fRawYield->SetBinContent(1,1);
 885   fRawYield->SetBinContent(2,1);
 886   fRawYield->SetBinContent(3,0.14);
 887   fRawYield->SetBinContent(4,0.14);
 888   fRawYield->SetBinContent(5,0.12);
 889   fRawYield->SetBinContent(6,0.12);
 890   fRawYield->SetBinContent(7,0.06);
 891   fRawYield->SetBinContent(8,0.06);
 892   fRawYield->SetBinContent(9,0.08);
 893   fRawYield->SetBinContent(10,0.08);
 894   fRawYield->SetBinContent(11,0.08);
 895   fRawYield->SetBinContent(12,0.08);
 896   for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
 897
 898   // Cuts efficiency (from cuts variation)
 899   fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
 900   for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
 901   fCutsEff->SetBinContent(3,0.15);
 902   fCutsEff->SetBinContent(4,0.15);
 903   fCutsEff->SetBinContent(5,0.15);
 904   fCutsEff->SetBinContent(6,0.15);
 905   fCutsEff->SetBinContent(7,0.10);
 906   fCutsEff->SetBinContent(8,0.10);
 907   fCutsEff->SetBinContent(9,0.10);
 908   fCutsEff->SetBinContent(10,0.10);
 909   fCutsEff->SetBinContent(11,0.10);
 910   fCutsEff->SetBinContent(12,0.10);
 911
 912   // PID efficiency (from PID/noPID)
 913   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
 914   for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
 915
 916   // MC dN/dpt
 917   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
 918   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
 919   fMCPtShape->SetBinContent(1,0.03);
 920   fMCPtShape->SetBinContent(2,0.03);
 921
 922   return;
 923 }
 924
 925 //------------------------------------------------------------------------
 926 void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
 927   //
 928   // D*+->D0pi syst errors. Responsible: A. Grelli
 929   //  2010 PbPb sample, 0-20 CC
 930   //
 931
 932   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
 933
 934  // Normalization
 935   fNorm = new TH1F("fNorm","fNorm",24,0,24);
 936   for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 937
 938   // Branching ratio
 939   fBR = new TH1F("fBR","fBR",24,0,24);
 940   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 941
 942   // Tracking efficiency
 943   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
 944   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 945
 946
 947   // Raw yield extraction
 948   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
 949   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);  //4%
 950   fRawYield->SetBinContent(3,0.2);
 951   fRawYield->SetBinContent(4,0.2);
 952   fRawYield->SetBinContent(5,0.2);
 953   fRawYield->SetBinContent(6,0.2);
 954
 955   // Cuts efficiency (from cuts variation)
 956   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
 957   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 958   fCutsEff->SetBinContent(4,0.15);
 959   fCutsEff->SetBinContent(5,0.15);
 960   fCutsEff->SetBinContent(6,0.15);
 961
 962   // PID efficiency (from PID/noPID)
 963   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
 964   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
 965
 966
 967   // MC dN/dpt  (from study on D* pt shape)
 968   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
 969   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
 970   fMCPtShape->SetBinContent(4,0.025);
 971   fMCPtShape->SetBinContent(5,0.025);
 972   fMCPtShape->SetBinContent(6,0.025);
 973   fMCPtShape->SetBinContent(7,0.04);
 974   fMCPtShape->SetBinContent(8,0.04);
 975   fMCPtShape->SetBinContent(9,0.03);
 976   fMCPtShape->SetBinContent(10,0.03);
 977   fMCPtShape->SetBinContent(11,0.03);
 978   fMCPtShape->SetBinContent(12,0.03);
 979
 980   return;
 981
 982 }
 983 // ----------------------------- 2011 ----------------------------------
 984 void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
 985   //
 986   // D*+->D0pi syst errors. Responsible: A. Grelli
 987   //  2011 PbPb sample, 0-7.5 CC
 988   //
 989
 990   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
 991
 992
 993  // Normalization
 994   fNorm = new TH1F("fNorm","fNorm",36,0,36);
 995   for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 996
 997   // Branching ratio
 998   fBR = new TH1F("fBR","fBR",36,0,36);
 999   for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1000
1001   // Tracking efficiency
1002   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1003   for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1004
1005
1006   // Raw yield extraction
1007   fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1008   for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05);  //4%
1009   fRawYield->SetBinContent(3,0.2);
1010   fRawYield->SetBinContent(4,0.10);
1011   fRawYield->SetBinContent(5,0.10);
1012   fRawYield->SetBinContent(6,0.08);
1013   for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15);  //4%
1014
1015
1016   // Cuts efficiency (from cuts variation)
1017   fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1018   for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1019   fCutsEff->SetBinContent(4,0.1);
1020   fCutsEff->SetBinContent(5,0.05);
1021   fCutsEff->SetBinContent(6,0.1);
1022
1023   // PID efficiency (from PID/noPID)
1024   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1025   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1026   fPIDEff->SetBinContent(3,0.09);
1027
1028   // MC dN/dpt  (from study on D* pt shape)
1029   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1030   for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.045);
1031   fMCPtShape->SetBinContent(4,0.025);
1032   fMCPtShape->SetBinContent(5,0.025);
1033   fMCPtShape->SetBinContent(6,0.025);
1034   fMCPtShape->SetBinContent(7,0.04);
1035   fMCPtShape->SetBinContent(8,0.04);
1036   fMCPtShape->SetBinContent(9,0.03);
1037   fMCPtShape->SetBinContent(10,0.03);
1038   fMCPtShape->SetBinContent(11,0.03);
1039   fMCPtShape->SetBinContent(12,0.03);
1040
1041
1042
1043   return;
1044
1045 }
1046 //-------------------------------------------------------------------------
1047 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1048   //
1049   // D*+->D0pi syst errors. Responsible: A. Grelli
1050   //  2010 PbPb sample, 20-40 CC
1051   //
1052
1053   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1054
1055  // Normalization
1056   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1057   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1058
1059   // Branching ratio
1060   fBR = new TH1F("fBR","fBR",24,0,24);
1061   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1062
1063   // Tracking efficiency
1064   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1065   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1066
1067
1068   // Raw yield extraction
1069   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1070   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15);  //4%
1071
1072   // Cuts efficiency (from cuts variation)
1073   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1074   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1075
1076   // PID efficiency (from PID/noPID)
1077   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1078   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1079
1080
1081   // MC dN/dpt  (copied from D0 : will update later)
1082   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1083   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1084   fMCPtShape->SetBinContent(1,0.03);
1085   fMCPtShape->SetBinContent(2,0.03);
1086
1087   return;
1088
1089 }
1090
1091 //--------------------------------------------------------------------------
1092 void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1093   //
1094   // D*+->D0pi syst errors. Responsible: A. Grelli
1095   //  2010 PbPb sample, 40-80 CC
1096   //
1097
1098   AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1099
1100  // Normalization
1101   fNorm = new TH1F("fNorm","fNorm",24,0,24);
1102   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1103
1104   // Branching ratio
1105   fBR = new TH1F("fBR","fBR",24,0,24);
1106   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1107
1108   // Tracking efficiency
1109   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1110   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1111
1112
1113   // Raw yield extraction
1114   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1115   for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);  //4%
1116   fRawYield->SetBinContent(1,1);
1117   fRawYield->SetBinContent(2,0.15);
1118   fRawYield->SetBinContent(3,0.15);
1119   fRawYield->SetBinContent(4,0.15);
1120   fRawYield->SetBinContent(5,0.15);
1121   fRawYield->SetBinContent(6,0.10);
1122   fRawYield->SetBinContent(7,0.10);
1123   fRawYield->SetBinContent(8,0.10);
1124   fRawYield->SetBinContent(9,0.11);
1125   fRawYield->SetBinContent(10,0.11);
1126   fRawYield->SetBinContent(11,0.11);
1127   fRawYield->SetBinContent(12,0.11);
1128   fRawYield->SetBinContent(13,0.08);
1129   fRawYield->SetBinContent(14,0.08);
1130   fRawYield->SetBinContent(15,0.08);
1131   fRawYield->SetBinContent(16,0.08);
1132
1133
1134   // Cuts efficiency (from cuts variation)
1135   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1136   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1137
1138   // PID efficiency (from PID/noPID)
1139   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1140   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1141
1142
1143   // MC dN/dpt  (copied from D0 : will update later)
1144   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1145   for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1146   fMCPtShape->SetBinContent(2,0.05);
1147   fMCPtShape->SetBinContent(3,0.05);
1148   fMCPtShape->SetBinContent(4,0.05);
1149   fMCPtShape->SetBinContent(5,0.04);
1150   fMCPtShape->SetBinContent(6,0.02);
1151   fMCPtShape->SetBinContent(7,0.04);
1152   fMCPtShape->SetBinContent(8,0.04);
1153
1154   return;
1155
1156 }
1157
1158 //--------------------------------------------------------------------------
1159 void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
1160   // define errors for RAA vs. centrality
1161   InitD0toKpi2010PbPb020();
1162   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1163   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
1164   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1165 }
1166 //--------------------------------------------------------------------------
1167 void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
1168   // define errors for RAA vs. centrality
1169   InitD0toKpi2010PbPb020();
1170   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1171   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);
1172   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1173 }
1174 //--------------------------------------------------------------------------
1175 void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
1176   // define errors for RAA vs. centrality
1177   InitD0toKpi2010PbPb4080();
1178   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1179   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
1180   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1181 }
1182 //--------------------------------------------------------------------------
1183 void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
1184    // define errors for RAA vs. centrality
1185   InitD0toKpi2010PbPb4080();
1186   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
1187   for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);
1188   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1189 }
1190 //--------------------------------------------------------------------------
1191 void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
1192    // define errors for RAA vs. centrality
1193   InitD0toKpi2010PbPb4080();
1194   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1195   for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
1196   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1197 }
1198 //--------------------------------------------------------------------------
1199 void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
1200   // define errors for RAA vs. centrality
1201   InitDplustoKpipi2010PbPb020();
1202   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
1203   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
1204
1205 }
1206 //--------------------------------------------------------------------------
1207 void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
1208   // define errors for RAA vs. centrality
1209   InitDplustoKpipi2010PbPb020();
1210   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
1211   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1212 }
1213 //--------------------------------------------------------------------------
1214 void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
1215   // define errors for RAA vs. centrality
1216   InitDplustoKpipi2010PbPb020();
1217   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1218   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
1219 }
1220 //--------------------------------------------------------------------------
1221 void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
1222   // define errors for RAA vs. centrality
1223   InitDplustoKpipi2010PbPb4080();
1224   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
1225   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1226 }
1227 //--------------------------------------------------------------------------
1228 void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
1229   // define errors for RAA vs. centrality
1230   InitDplustoKpipi2010PbPb4080();
1231   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
1232   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
1233 }
1234
1235 //--------------------------------------------------------------------------
1236 void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
1237   // define errors for RAA vs. centrality
1238   InitDstartoD0pi2010PbPb020();
1239   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
1240   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
1241 }
1242 //--------------------------------------------------------------------------
1243 void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
1244   // define errors for RAA vs. centrality
1245   InitDstartoD0pi2010PbPb020();
1246   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1247   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
1248 }
1249 //--------------------------------------------------------------------------
1250 void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
1251   // define errors for RAA vs. centrality
1252   InitDstartoD0pi2010PbPb2040();
1253   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1254   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
1255 }
1256 //--------------------------------------------------------------------------
1257 void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
1258   // define errors for RAA vs. centrality
1259   InitDstartoD0pi2010PbPb4080();
1260   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1261   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
1262 }
1263 //--------------------------------------------------------------------------
1264 void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
1265   // define errors for RAA vs. centrality
1266   InitDstartoD0pi2010PbPb4080();
1267   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
1268   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
1269 }
1270
1271
1272 //--------------------------------------------------------------------------
1273 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
1274   //
1275   // Get error
1276   //
1277
1278   Int_t bin=fCutsEff->FindBin(pt);
1279
1280   return fCutsEff->GetBinContent(bin);
1281 }
1282 //--------------------------------------------------------------------------
1283 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
1284   //
1285   // Get error
1286   //
1287
1288   Int_t bin=fMCPtShape->FindBin(pt);
1289
1290   return fMCPtShape->GetBinContent(bin);
1291 }
1292 //--------------------------------------------------------------------------
1293 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
1294   //
1295   // Get error
1296   //
1297
1298   Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
1299
1300   return TMath::Sqrt(err);
1301 }
1302 //--------------------------------------------------------------------------
1303 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
1304   //
1305   // Get error
1306   //
1307
1308   Int_t bin=fPIDEff->FindBin(pt);
1309
1310   return fPIDEff->GetBinContent(bin);
1311 }
1312 //--------------------------------------------------------------------------
1313 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
1314   //
1315   // Get error
1316   //
1317
1318   Int_t bin=fTrackingEff->FindBin(pt);
1319
1320   return fTrackingEff->GetBinContent(bin);
1321 }
1322 //--------------------------------------------------------------------------
1323 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
1324   //
1325   // Get error
1326   //
1327
1328   Int_t bin=fRawYield->FindBin(pt);
1329
1330   return fRawYield->GetBinContent(bin);
1331 }
1332 //--------------------------------------------------------------------------
1333 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
1334   //
1335   // Get error
1336   //
1337
1338   Int_t bin=fPartAntipart->FindBin(pt);
1339
1340   return fPartAntipart->GetBinContent(bin);
1341 }
1342 //--------------------------------------------------------------------------
1343 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
1344   //
1345   // Get total syst error (except norm. error)
1346   //
1347
1348   Double_t err=0.;
1349
1350   if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
1351   if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
1352   //  if(fBR) err += GetBRErr()*GetBRErr();
1353   if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
1354   if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
1355   if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
1356   if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
1357
1358   err += feeddownErr*feeddownErr;
1359
1360   return TMath::Sqrt(err);
1361 }
1362 //---------------------------------------------------------------------------
1363 void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
1364   //
1365   // Draw errors
1366   //
1367   gStyle->SetOptStat(0);
1368
1369   TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
1370   cSystErr->Range(0.20,-0.5,18.4,0.34);
1371   cSystErr->SetRightMargin(0.318);
1372   cSystErr->SetFillColor(0);
1373
1374   TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
1375   hFrame->SetAxisRange(1.,35.9,"X");
1376   hFrame->SetAxisRange(-0.5,0.5,"Y");
1377   hFrame->Draw();
1378
1379   TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
1380   leg->SetTextSize(0.03601695);
1381   leg->SetFillStyle(0);
1382   leg->SetBorderSize(0);
1383
1384   TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
1385   Int_t nbins = fNorm->GetNbinsX();
1386   TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
1387   for(Int_t i=1;i<=36;i++) {
1388     Double_t pt = hTotErr->GetBinCenter(i);
1389     Double_t ptwidth = hTotErr->GetBinWidth(i);
1390
1391     if(grErrFeeddown) {
1392       Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
1393       Double_t toterryl=0., toterryh=0.;
1394       for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
1395         grErrFeeddown->GetPoint(j,x,y);
1396         errxh = grErrFeeddown->GetErrorXhigh(j);
1397         errxl = grErrFeeddown->GetErrorXlow(j);
1398         if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
1399           erryh = grErrFeeddown->GetErrorYhigh(j);
1400           erryl = grErrFeeddown->GetErrorYlow(j);
1401         }
1402       }
1403       if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
1404       else toterryl = GetTotalSystErr(pt);
1405       if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
1406       else toterryh = GetTotalSystErr(pt);
1407
1408       hTotErr->SetBinContent(i,toterryh);
1409       gTotErr->SetPoint(i,pt,0.);
1410       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
1411     }
1412     else {
1413       hTotErr->SetBinContent(i,GetTotalSystErr(pt));
1414       gTotErr->SetPoint(i,pt,0.);
1415       gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
1416     }
1417
1418   }
1419   gTotErr->SetLineColor(kBlack);
1420   gTotErr->SetFillColor(kRed);
1421   gTotErr->SetFillStyle(3002);
1422   gTotErr->Draw("2");
1423   leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
1424 //   hTotErr->SetLineColor(1);
1425 //   hTotErr->SetLineWidth(3);
1426 //   hTotErr->Draw("same");
1427 //   leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
1428
1429
1430   fNorm->SetFillColor(1);
1431   fNorm->SetFillStyle(3002);
1432   //fNorm->Draw("same");
1433   //TH1F *hNormRefl = ReflectHisto(fNorm);
1434   //hNormRefl->Draw("same");
1435   Double_t norm = fNorm->GetBinContent(1)*100;
1436   leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
1437
1438   if(grErrFeeddown) {
1439     grErrFeeddown->SetFillColor(kTeal-8);
1440     grErrFeeddown->SetFillStyle(3001);
1441     grErrFeeddown->Draw("2");
1442     leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
1443   }
1444   if(fTrackingEff) {
1445     fTrackingEff->SetFillColor(4);
1446     fTrackingEff->SetFillStyle(3006);
1447     fTrackingEff->Draw("same");
1448     TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
1449     hTrackingEffRefl->Draw("same");
1450     leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
1451   }
1452   if(fBR) {
1453     fBR->SetFillColor(6);
1454     fBR->SetFillStyle(3005);
1455     //fBR->SetFillStyle(3020);
1456     fBR->Draw("same");
1457     TH1F *hBRRefl = ReflectHisto(fBR);
1458     hBRRefl->Draw("same");
1459     leg->AddEntry(fBR,"Branching ratio","f");
1460   }
1461   if(fRawYield) {
1462     Int_t ci;   // for color index setting
1463     ci = TColor::GetColor("#00cc00");
1464     fRawYield->SetLineColor(ci);
1465     //    fRawYield->SetLineColor(3);
1466     fRawYield->SetLineWidth(3);
1467     fRawYield->Draw("same");
1468     TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
1469     hRawYieldRefl->Draw("same");
1470     leg->AddEntry(fRawYield,"Yield extraction","l");
1471   }
1472   if(fCutsEff) {
1473     fCutsEff->SetLineColor(4);
1474     fCutsEff->SetLineWidth(3);
1475     fCutsEff->Draw("same");
1476     TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
1477     hCutsEffRefl->Draw("same");
1478     leg->AddEntry(fCutsEff,"Cuts efficiency","l");
1479   }
1480   if(fPIDEff) {
1481     fPIDEff->SetLineColor(7);
1482     fPIDEff->SetLineWidth(3);
1483     fPIDEff->Draw("same");
1484     TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
1485     hPIDEffRefl->Draw("same");
1486     leg->AddEntry(fPIDEff,"PID efficiency","l");
1487   }
1488   if(fMCPtShape) {
1489     Int_t ci = TColor::GetColor("#9933ff");
1490     fMCPtShape->SetLineColor(ci);
1491     //    fMCPtShape->SetLineColor(8);
1492     fMCPtShape->SetLineWidth(3);
1493     fMCPtShape->Draw("same");
1494     TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
1495     hMCPtShapeRefl->Draw("same");
1496     leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
1497   }
1498   if(fPartAntipart) {
1499     Int_t ci = TColor::GetColor("#ff6600");
1500     fPartAntipart->SetLineColor(ci);
1501     //    fPartAntipart->SetLineColor(9);
1502     fPartAntipart->SetLineWidth(3);
1503     fPartAntipart->Draw("same");
1504     TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
1505     hPartAntipartRefl->Draw("same");
1506     leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
1507   }
1508
1509
1510   leg->Draw();
1511
1512   cSystErr->SaveAs("RelativeSystematics.eps");
1513
1514   return;
1515 }
1516 //-------------------------------------------------------------------------
1517 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
1518   //
1519   // Clones and reflects histogram
1520   //
1521   TH1F *hout=(TH1F*)hin->Clone("hout");
1522   hout->Scale(-1.);
1523
1524   return hout;
1525 }
1526