]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGHF/vertexingHF/AliHFSystErr.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Updates syst err for pp low energy (Zaida) and D0 RAA vs centrality (Davide)
[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / AliHFSystErr.cxx
diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 537cb38a4182a819c68ed1afd129f9c43455d551..a96646578e02df566a392d6f12b46826d818f0ea 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -92,10 +92,7 @@ void AliHFSystErr::Init(Int_t decay){
   //
 
   if (fRunNumber!=10 && fIsLowEnergy==false) { 
-    AliError("Only settings for 2010 and the low energy runs are implemented so far");
-  }
-  if (fCentralityClass!="020" && fCentralityClass!="4080" && fCentralityClass!="0100"){
-    AliError("Only settings for MB2010 are implemented so far");
+    AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
   }
 
   switch(decay) {
@@ -111,7 +108,7 @@ void AliHFSystErr::Init(Int_t decay){
       else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
       else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
-      else AliError("Not yet implemented");
+      else AliFatal("Not yet implemented");
     }
     //    else if (fCollisionType==2) InitD0toKpi2010ppLowEn();
     break;
@@ -127,7 +124,7 @@ void AliHFSystErr::Init(Int_t decay){
       else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
       else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
-      else AliError("Not yet implemented");
+      else AliFatal("Not yet implemented");
     }
     break;
   case 3: // D*->D0pi
@@ -143,12 +140,12 @@ void AliHFSystErr::Init(Int_t decay){
       else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
       else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
-      else AliError("Not yet implemented");
+      else AliFatal("Not yet implemented");
     }
     break;
   case 4: // D+s->KKpi
     if (fCollisionType==0) InitDstoKKpi2010pp();
-    else AliError("Not yet implemented");
+    else AliFatal("Not yet implemented");
     break;
     
   default:
@@ -167,7 +164,7 @@ void AliHFSystErr::InitD0toKpi2010pp() {
 
   // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -340,7 +337,7 @@ void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.12); //10% (5% per track)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -362,9 +359,9 @@ void AliHFSystErr::InitD0toKpi2010ppLowEn() {
 
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
-//   fMCPtShape->SetBinContent(1,0.03);
-//   fMCPtShape->SetBinContent(2,0.03);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
+  fMCPtShape->SetBinContent(1,0.03);
+  fMCPtShape->SetBinContent(2,0.03);
 
 //   // particle-antiparticle
 //   fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
@@ -386,7 +383,7 @@ void AliHFSystErr::InitDplustoKpipi2010pp() {
 
 // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -449,7 +446,7 @@ void AliHFSystErr::InitDstoKKpi2010pp() {
 
 // 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
+  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 7% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",12,0,12);
@@ -641,7 +638,7 @@ void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.18); // 3% (1% per track)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -664,7 +661,7 @@ void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
  
   // MC dN/dpt  (copied from D0 : will update later)
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
   fMCPtShape->SetBinContent(1,0.03);
   fMCPtShape->SetBinContent(2,0.03);
 
@@ -680,7 +677,7 @@ void AliHFSystErr::InitDstartoD0pi2010pp() {
 
  // Normalization
   fNorm = new TH1F("fNorm","fNorm",24,0,24);
-  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.04); // 4% error on sigmaV0and
+  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
 
   // Branching ratio 
   fBR = new TH1F("fBR","fBR",24,0,24);
@@ -748,7 +745,7 @@ void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 
   // Tracking efficiency
   fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
-  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.18); //10% (to be checked!!)
+  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
 
   // Raw yield extraction
   fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
@@ -786,7 +783,7 @@ void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
 
   // MC dN/dpt
   fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
-  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
+  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
   fMCPtShape->SetBinContent(1,0.03);
   fMCPtShape->SetBinContent(2,0.03);
 
@@ -971,6 +968,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb020();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);      
   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
@@ -978,6 +976,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb020();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.17);    
   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
@@ -985,6 +984,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
   // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb4080();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.14);
   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
@@ -992,6 +992,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
    // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb4080();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=3;i<=5;i++)  fMCPtShape->SetBinContent(i,0.11);  
   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
@@ -999,6 +1000,7 @@ void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
    // define errors for RAA vs. centrality
   InitD0toKpi2010PbPb4080();
   for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
   for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08); 
 }
 //--------------------------------------------------------------------------
@@ -1177,8 +1179,8 @@ void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
   cSystErr->SetRightMargin(0.318);
   cSystErr->SetFillColor(0);
 
-  TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} [GeV/c]; Relative Error",20,0,20,100,-1,+1);
-  hFrame->SetAxisRange(1.,15.9,"X");
+  TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",30,0,30,100,-1,+1);
+  hFrame->SetAxisRange(1.,24.9,"X");
   hFrame->SetAxisRange(-0.5,0.5,"Y");
   hFrame->Draw();
 
@@ -1238,7 +1240,8 @@ void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
   //fNorm->Draw("same");
   //TH1F *hNormRefl = ReflectHisto(fNorm);
   //hNormRefl->Draw("same");
-  leg->AddEntry(fNorm,"Normalization (10%)","");
+  Double_t norm = fNorm->GetBinContent(1)*100;
+  leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
 
   if(grErrFeeddown) {
     grErrFeeddown->SetFillColor(kTeal-8);