Systematic errors for Ds PbPb 2011 (Gian Michele)

diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 47765b693e682ed5b4d87b2eaf45e77981584d6d..64106fe7dc9f32e77db8ec3726aa868a14365a55 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -149,6 +149,10 @@ void AliHFSystErr::Init(Int_t decay){
     break;
   case 4: // D+s->KKpi
     if (fCollisionType==0) InitDstoKKpi2010pp();
+    else if (fCollisionType==1) {
+      if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
+      else AliFatal("Not yet implemented");
+    }
     else AliFatal("Not yet implemented");
     break;
     
@@ -631,6 +635,60 @@ void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
   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.05); // 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,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.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() {
   // 

diff --git a/PWGHF/vertexingHF/AliHFSystErr.h b/PWGHF/vertexingHF/AliHFSystErr.h
index 17d107e375d88885864514cd32daee3cf500af60..4c9fea66333cd179232da025f27cd3b144b06e12 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.h
+++ b/PWGHF/vertexingHF/AliHFSystErr.h
@@ -87,6 +87,7 @@ class AliHFSystErr : public TNamed
   void InitDstartoD0pi2010PbPb2040CentScan();
   void InitDstartoD0pi2010PbPb4060CentScan();
   void InitDstartoD0pi2010PbPb6080CentScan();
+  
 
  private:
 
@@ -97,6 +98,7 @@ class AliHFSystErr : public TNamed
   void InitDplustoKpipi2010pp();
   void InitDstartoD0pi2010pp();
   void InitDstoKKpi2010pp();
+  void InitDstoKKpi2011PbPb07half();
 
   void InitD0toKpi2010PbPb020();
   void InitD0toKpi2011PbPb07half();