From 66247513bd9a7600bfa5cbf42d9521d2a846d976 Mon Sep 17 00:00:00 2001
From: dainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 18 Oct 2011 13:04:01 +0000
Subject: [PATCH] Added Ds pp 2010 (Gian Michele)

---
 PWG3/vertexingHF/AliHFSystErr.cxx | 66 +++++++++++++++++++++++++++++++
 PWG3/vertexingHF/AliHFSystErr.h   |  1 +
 2 files changed, 67 insertions(+)

diff --git a/PWG3/vertexingHF/AliHFSystErr.cxx b/PWG3/vertexingHF/AliHFSystErr.cxx
index b24d4af9534..778d81a35aa 100644
--- a/PWG3/vertexingHF/AliHFSystErr.cxx
+++ b/PWG3/vertexingHF/AliHFSystErr.cxx
@@ -130,6 +130,11 @@ void AliHFSystErr::Init(Int_t decay){
       if (fCentralityClass!="4080" && fCentralityClass!="2040" && fCentralityClass!="020")  AliError("Not yet implemented");
     }
     break;
+  case 4: // D+s->KKpi
+    if (fCollisionType==0) InitDstoKKpi2010pp();
+    else AliError("Not yet implemented");
+    break;
+    
   default:
     printf("Invalid decay type: %d\n",decay);
     break;
@@ -438,6 +443,67 @@ void AliHFSystErr::InitDplustoKpipi2010pp() {
   return;
 }
  
+//--------------------------------------------------------------------------
+void AliHFSystErr::InitDstoKKpi2010pp() {
+  // 
+  // D+s->KKpi syst errors. Responsible: G.M. Innocenti
+  //  2010 pp sample
+  //
+
+
+// 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
+
+  // 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
+
+  // 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)
+
+
+  // Raw yield extraction
+  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
+  fRawYield->SetBinContent(1,1);
+  fRawYield->SetBinContent(2,1);
+  fRawYield->SetBinContent(3,0.20);
+  fRawYield->SetBinContent(4,0.20);
+  fRawYield->SetBinContent(5,0.15);
+  fRawYield->SetBinContent(6,0.15);
+  fRawYield->SetBinContent(7,0.15);
+  fRawYield->SetBinContent(8,0.15);
+  fRawYield->SetBinContent(9,0.20);
+  fRawYield->SetBinContent(10,0.20);
+  fRawYield->SetBinContent(11,0.20);
+  fRawYield->SetBinContent(12,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.15); // 15%
+
+  // 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.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);
+
+
+  // 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/PWG3/vertexingHF/AliHFSystErr.h b/PWG3/vertexingHF/AliHFSystErr.h
index a5e16208df8..ad7640dd04b 100644
--- a/PWG3/vertexingHF/AliHFSystErr.h
+++ b/PWG3/vertexingHF/AliHFSystErr.h
@@ -74,6 +74,7 @@ class AliHFSystErr : public TNamed
   void InitD0toKpi2010pp();
   void InitDplustoKpipi2010pp();
   void InitDstartoD0pi2010pp();
+  void InitDstoKKpi2010pp();
 
   void InitD0toKpi2010PbPb020();
   void InitDplustoKpipi2010PbPb020();
-- 
2.43.5