PWGHF/vertexingHF/AliHFSystErr.cxx

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