From cede55f5b68d0ab291ef74cb2e2103edf5c644d0 Mon Sep 17 00:00:00 2001
From: prino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 6 Feb 2012 15:22:19 +0000
Subject: [PATCH] Add option for centrality scan (Zaida)

---
 PWGHF/vertexingHF/AliHFSystErr.cxx | 30 +++++++++++++++++++++++-------
 PWGHF/vertexingHF/AliHFSystErr.h   |  9 +++++++--
 2 files changed, 30 insertions(+), 9 deletions(-)

diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 064ca537f35..537cb38a418 100644
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -59,7 +59,8 @@ fPartAntipart(0),
 fRunNumber(10),
 fCollisionType(0),
 fCentralityClass("0100"),
-fIsLowEnergy(false)
+fIsLowEnergy(false),
+fIsCentScan(false)
 {
   //
   // Default Constructor
@@ -103,7 +104,12 @@ void AliHFSystErr::Init(Int_t decay){
       if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
       else InitD0toKpi2010pp();
     } else if (fCollisionType==1) {
-      if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
+      if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
+      else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
+      else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
+      else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
       else AliError("Not yet implemented");
     }
@@ -114,7 +120,12 @@ void AliHFSystErr::Init(Int_t decay){
       if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
       else InitDplustoKpipi2010pp();
     } else if (fCollisionType==1) {
-      if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
+      if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
+      else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
+      else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
+      else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
       else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
       else AliError("Not yet implemented");
     }
@@ -124,10 +135,15 @@ void AliHFSystErr::Init(Int_t decay){
       if(fIsLowEnergy)  InitDstartoD0pi2010ppLowEn();
       else InitDstartoD0pi2010pp();
     }else if (fCollisionType==1) {
-      if (fCentralityClass=="020")  InitDstartoD0pi2010PbPb020();
-      if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
-      if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
-      if (fCentralityClass!="4080" && fCentralityClass!="2040" && fCentralityClass!="020")  AliError("Not yet implemented");
+      if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
+      else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
+      else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
+      else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
+      else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
+      else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
+      else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
+      else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
+      else AliError("Not yet implemented");
     }
     break;
   case 4: // D+s->KKpi
diff --git a/PWGHF/vertexingHF/AliHFSystErr.h b/PWGHF/vertexingHF/AliHFSystErr.h
index d45f946d8c4..abe4b76c2b1 100644
--- a/PWGHF/vertexingHF/AliHFSystErr.h
+++ b/PWGHF/vertexingHF/AliHFSystErr.h
@@ -62,6 +62,10 @@ class AliHFSystErr : public TNamed
     fIsLowEnergy = flag; 
     AliInfo(" Settings for the low energy run");
   }
+  void SetIsPbPb2010EnergyScan(Bool_t flag) {
+    fIsCentScan = flag;
+    AliInfo(" Settings for the PbPb 2010 energy scan");
+  }
 
   // Function to initialize the variables/histograms
   void Init(Int_t decay);
@@ -107,7 +111,7 @@ class AliHFSystErr : public TNamed
   void InitDplustoKpipi2010ppLowEn();
   void InitDstartoD0pi2010ppLowEn();
 
- void InitDstartoD0pi2010PbPb2040();
+  void InitDstartoD0pi2010PbPb2040();
 
   TH1F* ReflectHisto(TH1F *hin) const;
 
@@ -125,8 +129,9 @@ class AliHFSystErr : public TNamed
   TString fCentralityClass;  // Centrality class
                            // MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
   Bool_t fIsLowEnergy;     // flag for the low energy (2.76TeV) run
+  Bool_t fIsCentScan;      // flag fot the PbPb centrality scan
 
-  ClassDef(AliHFSystErr,3);  // class for systematic errors of charm hadrons
+  ClassDef(AliHFSystErr,4);  // class for systematic errors of charm hadrons
 };
 
 #endif
-- 
2.39.3