D* systematics in pPb (G.Luparello)

diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 080832be127743653c8d33f36071bbc3424a734c..72b8e9cc08ae7838a3296b9faa7aaef4482f0aef 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -1482,14 +1482,14 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
 
 //_________________________________________________________________________
 void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
+  // Dstar syst in pPb 2013 MB
+  //
 
-  // Dstar syst
-  
   //Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
   for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.00); // 
   
- // Branching ratio 
+  // Branching ratio
   fBR = new TH1F("fBR","fBR",24,0,24);
   for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
 
@@ -1499,21 +1499,36 @@ void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
-  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.0);  
-  
+  fRawYield->SetBinContent(1,1.0);
+  fRawYield->SetBinContent(2,0.10);
+  fRawYield->SetBinContent(3,0.05);
+  fRawYield->SetBinContent(4,0.02);
+  fRawYield->SetBinContent(5,0.02);
+  fRawYield->SetBinContent(6,0.04);
+  fRawYield->SetBinContent(7,0.04);
+  fRawYield->SetBinContent(8,0.04);
+  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); 
+  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
+
   // 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.0); //
+  fCutsEff->SetBinContent(1,0.);
+  fCutsEff->SetBinContent(2,0.10);
+  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
 
   // PID efficiency (from PID/noPID)
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
+  fPIDEff->SetBinContent(1,0.0);
   fPIDEff->SetBinContent(2,0.05);
-  for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.00); 
-  for(Int_t i=12;i<=24;i++) fPIDEff->SetBinContent(i,0.05); 
+  for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.00);
+  for(Int_t i=12;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
 
   // MC dN/dpt  
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
-  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
+  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.02);
+  for(Int_t i=4;i<=8;i++) fMCPtShape->SetBinContent(i,0.);
+  for(Int_t i=9;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
 
   // particle-antiparticle
   /*