Fix in the SetLabelArray of CF for Lc->V0+bachelor (Annalisa)

[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / AliHFSystErr.cxx

diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 7c7f50721112b656f8f036aa970cf03ed5357ba3..81ce4302ccd1a136a55499504af94df67578c90c 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -1,3 +1,4 @@
+
 /**************************************************************************
  * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
@@ -59,7 +60,8 @@ fPartAntipart(0),
 fRunNumber(10),
 fCollisionType(0),
 fCentralityClass("0100"),
-fIsLowEnergy(false)
+fIsLowEnergy(false),
+fIsCentScan(false)
 {
   //
   // Default Constructor
@@ -91,10 +93,7 @@ void AliHFSystErr::Init(Int_t decay){
   //
 
   if (fRunNumber!=10 && fIsLowEnergy==false) { 
-    AliError("Only settings for 2010 and the low energy runs are implemented so far");
-  }
-  if (fCentralityClass!="020" && fCentralityClass!="4080" && fCentralityClass!="0100"){
-    AliError("Only settings for MB2010 are implemented so far");
+    AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
   }
 
   switch(decay) {
@@ -103,9 +102,15 @@ void AliHFSystErr::Init(Int_t decay){
       if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
       else InitD0toKpi2010pp();
     } else if (fCollisionType==1) {
-      if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
+      if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
+      else if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
+      else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
+      else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
+      else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
-      else AliError("Not yet implemented");
+      else AliFatal("Not yet implemented");
     }
     //    else if (fCollisionType==2) InitD0toKpi2010ppLowEn();
     break;
@@ -114,9 +119,15 @@ void AliHFSystErr::Init(Int_t decay){
       if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
       else InitDplustoKpipi2010pp();
     } else if (fCollisionType==1) {
-      if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
+      if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
+      else if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
+      else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
+      else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
+      else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
-      else AliError("Not yet implemented");
+      else AliFatal("Not yet implemented");
     }
     break;
   case 3: // D*->D0pi
@@ -124,15 +135,25 @@ void AliHFSystErr::Init(Int_t decay){
       if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
       else InitDstartoD0pi2010pp();
     }else if (fCollisionType==1) {
-      if (fCentralityClass=="020")  InitDstartoD0pi2010PbPb020();
-      if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
-      if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
-      if (fCentralityClass!="4080" && fCentralityClass!="2040" && fCentralityClass!="020")  AliError("Not yet implemented");
+      if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
+      else if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
+      else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
+      else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
+      else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
+      else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
+      else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
+      else AliFatal("Not yet implemented");
     }
     break;
   case 4: // D+s->KKpi
     if (fCollisionType==0) InitDstoKKpi2010pp();
-    else AliError("Not yet implemented");
+    else if (fCollisionType==1) {
+      if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
+      else AliFatal("Not yet implemented");
+    }
+    else AliFatal("Not yet implemented");
     break;
     
   default:
@@ -151,7 +172,7 @@ void AliHFSystErr::InitD0toKpi2010pp() {
 
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -202,7 +223,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb020() {
 
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",20,0,20);
-  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",20,0,20);
@@ -210,7 +231,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb020() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.08);//0.05+0.005*(Float_t)i);
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -238,6 +259,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb020() {
   for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
   fMCPtShape->SetBinContent(3,0.04);
   fMCPtShape->SetBinContent(4,0.02);
+  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03); 
 
 //   // particle-antiparticle
 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
@@ -250,6 +272,72 @@ void AliHFSystErr::InitD0toKpi2010PbPb020() {
   return;
 }
 //--------------------------------------------------------------------------
+//--------------------------------------------------------------------------
+void AliHFSystErr::InitD0toKpi2011PbPb07half() {
+  // 
+  // D0->Kpi syst errors. Responsible: A. Rossi
+  //   2011 PbPb sample, 0-7.5 CC
+  //
+
+  // Normalization
+  fNorm = new TH1F("fNorm","fNorm",36,0,36);
+  for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
+
+  // Branching ratio 
+  fBR = new TH1F("fBR","fBR",36,0,36);
+  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
+
+  // Tracking efficiency
+  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
+  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
+  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
+
+  // Raw yield extraction
+  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
+  fRawYield->SetBinContent(1,0);
+  fRawYield->SetBinContent(2,0.1);
+  fRawYield->SetBinContent(3,0.05);
+  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
+  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
+  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
+
+  // Cuts efficiency (from cuts variation)
+  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
+  fCutsEff->SetBinContent(1,0.);
+  fCutsEff->SetBinContent(2,0.15);
+  fCutsEff->SetBinContent(3,0.13);
+  fCutsEff->SetBinContent(4,0.11);  
+  fCutsEff->SetBinContent(5,0.08); 
+  for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
+
+  // PID efficiency (from PID/noPID)
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
+  for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
+
+  // MC dN/dpt
+  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
+  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
+  fMCPtShape->SetBinContent(2,0.06);
+  fMCPtShape->SetBinContent(3,0.04);
+  fMCPtShape->SetBinContent(4,0.02);
+  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03); 
+  for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
+  for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
+
+//   // particle-antiparticle
+//   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
+//   for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
+//   fPartAntipart->SetBinContent(3,0.10);
+//   fPartAntipart->SetBinContent(4,0.10);
+//   fPartAntipart->SetBinContent(7,0.10);
+//   fPartAntipart->SetBinContent(8,0.10);
+  
+  return;
+}
+//--------------------------------------------------------------------------
 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
   // 
   // D0->Kpi syst errors. Responsible: A. Rossi
@@ -258,7 +346,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb4080() {
 
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",20,0,20);
-  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",20,0,20);
@@ -266,7 +354,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb4080() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.08);//0.5*(0.05+0.005*(Float_t)i));
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
 
 
   // Raw yield extraction
@@ -323,7 +411,7 @@ void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.12); //10% (5% per track)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -345,9 +433,9 @@ void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
-//   fMCPtShape->SetBinContent(1,0.03);
-//   fMCPtShape->SetBinContent(2,0.03);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
+  fMCPtShape->SetBinContent(1,0.03);
+  fMCPtShape->SetBinContent(2,0.03);
 
 //   // particle-antiparticle
 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
@@ -369,7 +457,7 @@ void AliHFSystErr::InitDplustoKpipi2010pp() {
 
 // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -432,15 +520,15 @@ void AliHFSystErr::InitDstoKKpi2010pp() {
 
 // Normalization
   fNorm = new TH1F("fNorm","fNorm",12,0,12);
-  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.07); // 7% error on sigmaV0and
+  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",12,0,12);
-  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.05); // 5% PDG2010
+  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
-  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.03); // 3% (1% per track)
+  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
 
 
   // Raw yield extraction
@@ -466,9 +554,13 @@ void AliHFSystErr::InitDstoKKpi2010pp() {
   fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
   for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
 
-  // MC dN/dpt  (copied from D0 : will update later)
-  //fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
-  //for(Int_t i=1;i<=12;i++) fMCPtShape->SetBinContent(i,(Float_t)i*0.006);
+  // MC dN/dpt 
+  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
+  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
+  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
+  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
+  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
+  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
 
 
   // particle-antiparticle
@@ -483,6 +575,121 @@ void AliHFSystErr::InitDstoKKpi2010pp() {
 }
    
  
+//--------------------------------------------------------------------------
+void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
+  // 
+  // D+->Kpipi syst errors. Responsible: E. Bruna
+  //  2011 PbPb sample, 0-7.5% CC
+  //
+
+ // Normalization
+  fNorm = new TH1F("fNorm","fNorm",36,0,36);
+  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
+
+  // Branching ratio 
+  fBR = new TH1F("fBR","fBR",36,0,36);
+  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
+
+  // Tracking efficiency
+  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
+
+  // Raw yield extraction
+  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
+  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10);  //5 to 10%
+  fRawYield->SetBinContent(4,0.30);
+  fRawYield->SetBinContent(5,0.20);
+  fRawYield->SetBinContent(6,0.20);
+  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
+  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
+
+
+  // Cuts efficiency (from cuts variation)
+  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
+  for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
+  for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
+
+
+  // PID efficiency (from PID/noPID)
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
+
+  // MC dN/dpt  (24/7/2012)
+  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
+  for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
+  for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
+  //  for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
+
+
+  // particle-antiparticle
+  /*
+  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
+  fPartAntipart->SetBinContent(1,1);
+  fPartAntipart->SetBinContent(2,1);
+  fPartAntipart->SetBinContent(3,0.12);
+  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
+  */
+
+  return;
+}
+
+//--------------------------------------------------------------------------
+void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
+  // 
+  // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
+  //  2011 PbPb sample, 0-7.5% CC
+  //
+
+ // Normalization
+  fNorm = new TH1F("fNorm","fNorm",12,0,12);
+  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
+
+  // Branching ratio 
+  fBR = new TH1F("fBR","fBR",12,0,12);
+  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
+
+  // Tracking efficiency
+  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
+  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
+
+  // Raw yield extraction
+  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
+  for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30); 
+  for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
+ 
+
+
+  // Cuts efficiency (from cuts variation)
+  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
+  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
+
+  // PID efficiency (from PID/noPID)
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
+  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
+
+   // MC dN/dpt 
+  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
+  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
+  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
+  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
+  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
+  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
+
+  // particle-antiparticle
+  /*
+  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
+  fPartAntipart->SetBinContent(1,1);
+  fPartAntipart->SetBinContent(2,1);
+  fPartAntipart->SetBinContent(3,0.12);
+  for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);   //5 to 12%
+  */
+
+  return;
+}
+
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
   // 
@@ -492,7 +699,7 @@ void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",20,0,20);
-  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",20,0,20);
@@ -500,7 +707,7 @@ void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track) //i,0.075+0.005*(Float_t)i); 
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -552,7 +759,7 @@ void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",20,0,20);
-  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",20,0,20);
@@ -560,7 +767,7 @@ void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 
 
   // Raw yield extraction
@@ -624,7 +831,7 @@ void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.18); // 3% (1% per track)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -647,7 +854,7 @@ void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
  
   // MC dN/dpt  (copied from D0 : will update later)
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
   fMCPtShape->SetBinContent(1,0.03);
   fMCPtShape->SetBinContent(2,0.03);
 
@@ -663,7 +870,7 @@ void AliHFSystErr::InitDstartoD0pi2010pp() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -731,7 +938,7 @@ void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.18); //10% (to be checked!!)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -769,7 +976,7 @@ void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
   fMCPtShape->SetBinContent(1,0.03);
   fMCPtShape->SetBinContent(2,0.03);
 
@@ -787,7 +994,7 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -795,7 +1002,7 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
-  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)i,0.075+0.005*(Float_t)i); 
+  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 
 
   // Raw yield extraction
@@ -809,22 +1016,90 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
   // Cuts efficiency (from cuts variation)
   fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
   for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
+  fCutsEff->SetBinContent(4,0.15);
+  fCutsEff->SetBinContent(5,0.15);
+  fCutsEff->SetBinContent(6,0.15);
 
   // PID efficiency (from PID/noPID)
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
   for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
  
 
-  // MC dN/dpt  (copied from D0 : will update later)
+  // MC dN/dpt  (from study on D* pt shape)
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
-  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
-  fMCPtShape->SetBinContent(1,0.03);
-  fMCPtShape->SetBinContent(2,0.03);
-
+  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
+  fMCPtShape->SetBinContent(4,0.025);
+  fMCPtShape->SetBinContent(5,0.025);
+  fMCPtShape->SetBinContent(6,0.025);
+  fMCPtShape->SetBinContent(7,0.04);
+  fMCPtShape->SetBinContent(8,0.04);
+  fMCPtShape->SetBinContent(9,0.03);
+  fMCPtShape->SetBinContent(10,0.03);
+  fMCPtShape->SetBinContent(11,0.03);
+  fMCPtShape->SetBinContent(12,0.03);
+  
   return;
 
 }
+// ----------------------------- 2011 ----------------------------------
+void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
+  // 
+  // D*+->D0pi syst errors. Responsible: A. Grelli
+  //  2011 PbPb sample, 0-7.5 CC
+  //
+
+  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY"); 
+
+
+ // Normalization
+  fNorm = new TH1F("fNorm","fNorm",36,0,36);
+  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
+
+  // Branching ratio 
+  fBR = new TH1F("fBR","fBR",36,0,36);
+  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
+
+  // Tracking efficiency
+  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
+  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
+
+
+  // Raw yield extraction
+  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
+  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05);  //4%
+  fRawYield->SetBinContent(4,0.2);
+  fRawYield->SetBinContent(5,0.10);
+  fRawYield->SetBinContent(6,0.10);
+  fRawYield->SetBinContent(7,0.08);
+  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15);  //4%
+
+  // Cuts efficiency (from cuts variation)
+  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
 
+  // PID efficiency (from PID/noPID)
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
+  fPIDEff->SetBinContent(4,0.09);
+
+  // MC dN/dpt  (from study on D* pt shape)
+  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
+  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
+  fMCPtShape->SetBinContent(4,0.015);
+  fMCPtShape->SetBinContent(5,0.015);
+  fMCPtShape->SetBinContent(6,0.015);
+  fMCPtShape->SetBinContent(7,0.02);
+  fMCPtShape->SetBinContent(8,0.02);
+  fMCPtShape->SetBinContent(9,0.03);
+  fMCPtShape->SetBinContent(10,0.03);
+  fMCPtShape->SetBinContent(11,0.03);
+  fMCPtShape->SetBinContent(12,0.03);
+  
+  
+
+  return;
+
+}
 //-------------------------------------------------------------------------
 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
   // 
@@ -844,7 +1119,7 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
-  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 3% (1% per track)
+  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 
 
   // Raw yield extraction
@@ -881,7 +1156,7 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -889,7 +1164,7 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
-  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12);//0.5*(0.075+0.005*(Float_t)i));
+  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
 
 
   // Raw yield extraction
@@ -924,10 +1199,15 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 
   // MC dN/dpt  (copied from D0 : will update later)
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
-  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
-  fMCPtShape->SetBinContent(1,0.03);
-  fMCPtShape->SetBinContent(2,0.03);
-
+  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
+  fMCPtShape->SetBinContent(2,0.05);
+  fMCPtShape->SetBinContent(3,0.05);
+  fMCPtShape->SetBinContent(4,0.05);
+  fMCPtShape->SetBinContent(5,0.04);
+  fMCPtShape->SetBinContent(6,0.02);
+  fMCPtShape->SetBinContent(7,0.04);
+  fMCPtShape->SetBinContent(8,0.04);
+ 
   return;
 
 }
@@ -936,51 +1216,77 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb020();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);      
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb020();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);    
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
   // define errors for RAA vs. centrality
-  InitD0toKpi2010PbPb020();
+  InitD0toKpi2010PbPb4080();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
-  // define errors for RAA vs. centrality
+   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb4080();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);  
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
-  // define errors for RAA vs. centrality
+   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb4080();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
   // define errors for RAA vs. centrality
   InitDplustoKpipi2010PbPb020();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
+
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
   // define errors for RAA vs. centrality
   InitDplustoKpipi2010PbPb020();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
   // define errors for RAA vs. centrality
   InitDplustoKpipi2010PbPb020();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
   // define errors for RAA vs. centrality
   InitDplustoKpipi2010PbPb4080();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
   // define errors for RAA vs. centrality
   InitDplustoKpipi2010PbPb4080();
+  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
 }
 
 //--------------------------------------------------------------------------
@@ -988,30 +1294,35 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
   // define errors for RAA vs. centrality
   InitDstartoD0pi2010PbPb020();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16); 
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
   // define errors for RAA vs. centrality
   InitDstartoD0pi2010PbPb020();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05); 
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
   // define errors for RAA vs. centrality
   InitDstartoD0pi2010PbPb2040();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
   // define errors for RAA vs. centrality
   InitDstartoD0pi2010PbPb4080();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
 }
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
   // define errors for RAA vs. centrality
   InitDstartoD0pi2010PbPb4080();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); 
+  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
 }
 
 
@@ -1117,8 +1428,8 @@ void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
   cSystErr->SetRightMargin(0.318);
   cSystErr->SetFillColor(0);
 
-  TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} [GeV/c]; Relative Error",20,0,20,100,-1,+1);
-  hFrame->SetAxisRange(1.,15.9,"X");
+  TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
+  hFrame->SetAxisRange(1.,35.9,"X");
   hFrame->SetAxisRange(-0.5,0.5,"Y");
   hFrame->Draw();
 
@@ -1127,10 +1438,10 @@ void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
   leg->SetFillStyle(0);
   leg->SetBorderSize(0);
   
-  TH1F *hTotErr=new TH1F("hTotErr","",24,0,24);
+  TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
   Int_t nbins = fNorm->GetNbinsX();
   TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
-  for(Int_t i=1;i<=24;i++) {
+  for(Int_t i=1;i<=36;i++) {
     Double_t pt = hTotErr->GetBinCenter(i);
     Double_t ptwidth = hTotErr->GetBinWidth(i);
 
@@ -1178,7 +1489,8 @@ void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
   //fNorm->Draw("same");
   //TH1F *hNormRefl = ReflectHisto(fNorm);
   //hNormRefl->Draw("same");
-  leg->AddEntry(fNorm,"Normalization (10%)","");
+  Double_t norm = fNorm->GetBinContent(1)*100;
+  leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
 
   if(grErrFeeddown) {
     grErrFeeddown->SetFillColor(kTeal-8);
@@ -1268,3 +1580,6 @@ TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
 
   return hout;
 }
+
+
+