D+ pPb systematics (R.Russo)
authorzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 1 Jul 2013 16:25:44 +0000 (16:25 +0000)
committerzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 1 Jul 2013 16:25:44 +0000 (16:25 +0000)
PWGHF/vertexingHF/AliHFSystErr.cxx

index 04c88e8..816ac37 100644 (file)
@@ -1101,24 +1101,46 @@ void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
   for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
-
+  
   // 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);
+  fRawYield->SetBinContent(2,1);
+  fRawYield->SetBinContent(3,0.06);
+  fRawYield->SetBinContent(4,0.12);
+  fRawYield->SetBinContent(5,0.05);
+  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.07);
+  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.17);
+  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.14);
+
   // 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,1);
+  fCutsEff->SetBinContent(2,1);
+  fCutsEff->SetBinContent(3,0.07);
+  fCutsEff->SetBinContent(4,0.04);
+  fCutsEff->SetBinContent(5,0.04);
+  fCutsEff->SetBinContent(6,0.08);
+  fCutsEff->SetBinContent(7,0.08);
+  fCutsEff->SetBinContent(8,0.08);
+  for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
+  for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
+  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.1); 
 
   // PID efficiency (from PID/noPID)
   fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
-  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); 
+  fPIDEff->SetBinContent(1,1);
+  fPIDEff->SetBinContent(2,1);
+  for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.0); //
+  for(Int_t i=12;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
 
-  // MC dN/dpt  
+  // 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<=12;i++) fMCPtShape->SetBinContent(i,0.02);
+  for(Int_t i=13;i<=24;i++) fMCPtShape->SetBinContent(i,0.04);
 
   // particle-antiparticle
   /*