Fix in Dstar syst. errors + updates in macro (Zaida)

diff --git a/PWG3/vertexingHF/AliHFSystErr.cxx b/PWG3/vertexingHF/AliHFSystErr.cxx
index 33a6c526ac7d62ef36474eebf0a1cc924747807f..ce5f3f4557bf886189058d544d489625e760a2d2 100644 (file)
--- a/PWG3/vertexingHF/AliHFSystErr.cxx
+++ b/PWG3/vertexingHF/AliHFSystErr.cxx
@@ -715,6 +715,8 @@ void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
   fCutsEff->SetBinContent(8,0.12);
   fCutsEff->SetBinContent(9,0.2);
   fCutsEff->SetBinContent(10,0.2);
+  fCutsEff->SetBinContent(11,0.2);
+  fCutsEff->SetBinContent(12,0.2);
 
   // PID efficiency (from PID/noPID)
   fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);

diff --git a/PWG3/vertexingHF/macros/CombineFeedDownMCSubtractionMethodsUncertainties.C b/PWG3/vertexingHF/macros/CombineFeedDownMCSubtractionMethodsUncertainties.C
index bd8dcace0503d7f849f10d31774fb50ad293797d..ccfd5dc8970e0bf8a598e85a8b4d8cb370d041fc 100644 (file)
--- a/PWG3/vertexingHF/macros/CombineFeedDownMCSubtractionMethodsUncertainties.C
+++ b/PWG3/vertexingHF/macros/CombineFeedDownMCSubtractionMethodsUncertainties.C
@@ -35,13 +35,13 @@
 //
 //_________________________________________________________________________________________
 
-enum centrality{ kpp, k010, k020, k2040, k4060, k6080, k4080, k80100 };
+enum centrality{ kpp7, kpp276, k010, k020, k2040, k4060, k6080, k4080, k80100 };
 
 void CombineFeedDownMCSubtractionMethodsUncertainties(const char *fcfilename="HFPtSpectrum_D0Kpi_method1_221110_newnorm.root",
                                                      const char *nbfilename="HFPtSpectrum_D0Kpi_method2_221110_newnorm.root",
                                                      const char *outfilename="HFPtSpectrum_D0Kpi_combinedFD.root",
                                                      const char *thfilename="D0DplusDstarPredictions_y05.root",
-                                                     Int_t decay=1, Int_t centrality=kpp)
+                                                     Int_t decay=1, Int_t centrality=kpp7)
 {
   
   // 
@@ -108,13 +108,15 @@ void CombineFeedDownMCSubtractionMethodsUncertainties(const char *fcfilename="HF
   // Call the systematics uncertainty class for a given decay
   //   will help to compute the systematical unc. (but FD) 
   AliHFSystErr systematics;
-  if( centrality!=kpp)  {
+  if( centrality==kpp276 ) {
+    systematics.SetIsLowEnergy(true);
+  } else if( centrality!=kpp7 )  {
     systematics.SetCollisionType(1);
     if ( centrality == k020 ) {
-      systematics.SetCentrality(020);
+      systematics.SetCentrality("020");
     }
     else if ( centrality == k4080 ) {
-      systematics.SetCentrality(4080);
+      systematics.SetCentrality("4080");
     }
     else { 
       cout << " Systematics not yet implemented " << endl;
@@ -140,6 +142,7 @@ void CombineFeedDownMCSubtractionMethodsUncertainties(const char *fcfilename="HF
     valFcErrstat = histoSigmaCorrFc->GetBinError(ibin);
     Double_t value =0., ptt=0.;
     gSigmaCorrConservativeFc->GetPoint(ibin,ptt,value);
+    if (value<=0.) continue;
     if ( TMath::Abs(valFc-value)>0.1 || TMath::Abs(pt-ptt)>0.1 ) 
       cout << "Hey you ! There might be a problem with the fc input file, please, have a look !" << endl;
     valFcErrx = gSigmaCorrFc->GetErrorXlow(ibin);