Systematic errors for Ds PbPb 2011 (Gian Michele)
authorprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 17 Jul 2012 22:42:55 +0000 (22:42 +0000)
committerprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 17 Jul 2012 22:42:55 +0000 (22:42 +0000)
PWGHF/vertexingHF/AliHFSystErr.cxx
PWGHF/vertexingHF/AliHFSystErr.h

index 47765b6..64106fe 100644 (file)
@@ -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;
     
@@ -632,6 +636,60 @@ void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
 }
 
 //--------------------------------------------------------------------------
+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() {
   // 
   // D+->Kpipi syst errors. Responsible: ??
index 17d107e..4c9fea6 100644 (file)
@@ -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();