Add centrality classes
authorzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Nov 2012 18:37:56 +0000 (18:37 +0000)
committerzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Nov 2012 18:37:56 +0000 (18:37 +0000)
PWGHF/vertexingHF/macros/HFPtSpectrum.C
PWGHF/vertexingHF/macros/HFPtSpectrumRaa.C

index aa040c2..94f6bdc 100644 (file)
@@ -38,7 +38,7 @@
 //  9) Flag to decide if there is need to evaluate the dependence on the energy loss
 //
 
-enum centrality{ kpp7, kpp276, k07half, k010, k1020, k020, k2040, k3050, k4060, k6080, k4080, k80100 };
+enum centrality{ kpp7, kpp276, k07half, k010, k1020, k020, k2040, k2030, k3040, k4050, k3050, k5060, k4060, k6080, k4080, k80100 };
 enum BFDSubtrMethod { knone, kfc, kNb };
 enum RaavsEP {kPhiIntegrated, kInPlane, kOutOfPlane};
 
@@ -95,12 +95,20 @@ void HFPtSpectrum ( const char *mcfilename="FeedDownCorrectionMC.root",
     tab = 18.93; tabUnc = 0.74;
   } else if ( cc == k2040 ) {
     tab = 6.86; tabUnc = 0.28;
+  } else if ( cc == k2030 ) {
+    tab = 8.73769; tabUnc = 0.370219;
+  } else if ( cc == k3040 ) {
+    tab = 5.02755; tabUnc = 0.22099;
+  } else if ( cc == k4050 ) {
+    tab = 2.68327; tabUnc = 0.137073;
   } else if ( cc == k3050 ) {
     tab = 3.87011; tabUnc = 0.183847;
   } else if ( cc == k4060 ) {
     tab = 2.00;  tabUnc= 0.11;
   } else if ( cc == k4080 ) {
     tab = 1.20451; tabUnc = 0.071843;
+  } else if ( cc == k5060 ) {
+    tab = 1.32884; tabUnc = 0.0929536;
   } else if ( cc == k6080 ) {
     tab = 0.419; tabUnc = 0.033;
   } else if ( cc == k80100 ){
@@ -314,7 +322,7 @@ void HFPtSpectrum ( const char *mcfilename="FeedDownCorrectionMC.root",
     else if ( cc == k010 )  systematics->SetCentrality("010");
     else if ( cc == k1020 )  systematics->SetCentrality("1020");
     else if ( cc == k020 )  systematics->SetCentrality("020");
-    else if ( cc == k2040 ) {
+    else if ( cc == k2040 || cc == k2030 || cc == k3040 ) {
       systematics->SetCentrality("2040");
       systematics->SetIsPbPb2010EnergyScan(true);
     }
@@ -323,7 +331,7 @@ void HFPtSpectrum ( const char *mcfilename="FeedDownCorrectionMC.root",
       else if (isRaavsEP == kInPlane) systematics->SetCentrality("3050InPlane");
       else if (isRaavsEP == kOutOfPlane) systematics->SetCentrality("3050OutOfPlane");
     }
-    else if ( cc == k4060 )  systematics->SetCentrality("4060");
+    else if ( cc == k4060 || cc == k4050 || cc == k5060 )  systematics->SetCentrality("4060");
     else if ( cc == k6080 )  systematics->SetCentrality("6080");
     else if ( cc == k4080 ) systematics->SetCentrality("4080");
     else {
index 1ec4069..b88a70f 100644 (file)
@@ -47,7 +47,7 @@
 //
 ///////////////////////////////////////////////////////////////////////////////
 
-enum centrality{ kpp, k07half, k010, k1020, k020, k2040, k3050, k4060, k6080, k4080, k80100 };
+enum centrality{ kpp, k07half, k010, k1020, k020, k2040, k2030, k3040, k4050, k3050, k5060, k4060, k6080, k4080, k80100 };
 enum energy{ k276, k55 };
 enum BFDSubtrMethod { kfc, kNb };
 enum RaavsEP {kPhiIntegrated, kInPlane, kOutOfPlane};
@@ -119,19 +119,27 @@ void HFPtSpectrumRaa(const char *ppfile="HFPtSpectrum_D0Kpi_method2_rebinnedth_2
   } else if ( cc == k010 ) {
     Tab = 23.48; TabSyst = 0.97;
   } else if ( cc == k1020 ) {
-    Tab = 14.4318; TabSyst = 0.573289;
+    Tab = 14.4318; TabSyst = 0.5733;
   } else if ( cc == k020 ) {
     Tab = 18.93; TabSyst = 0.74;
   } else if ( cc == k2040 ) {
     Tab = 6.86; TabSyst = 0.28;
+  } else if ( cc == k2030 ) {
+    Tab = 8.73769; TabSyst = 0.370219;
+  } else if ( cc == k3040 ) {
+    Tab = 5.02755; TabSyst = 0.22099;
+  } else if ( cc == k4050 ) {
+    Tab = 2.68327; TabSyst = 0.137073;
   } else if ( cc == k3050 ) {
     Tab = 3.87011; TabSyst = 0.183847;
   } else if ( cc == k4060 ) {
-    Tab = 2.00; TabSyst = 0.11;
-  } else if ( cc == k6080 ) {
-    Tab = 0.419; TabSyst = 0.033;
+    Tab = 2.00;  TabSyst= 0.11;
   } else if ( cc == k4080 ) {
     Tab = 1.20451; TabSyst = 0.071843;
+  } else if ( cc == k5060 ) {
+    Tab = 1.32884; TabSyst = 0.0929536;
+  } else if ( cc == k6080 ) {
+    Tab = 0.419; TabSyst = 0.033;
   } else if ( cc == k80100 ){
     Tab = 0.0690; TabSyst = 0.0062;
   }
@@ -203,11 +211,11 @@ void HFPtSpectrumRaa(const char *ppfile="HFPtSpectrum_D0Kpi_method2_rebinnedth_2
   else if ( cc == k010 ) systematicsAB->SetCentrality("010");
   else if ( cc == k1020 ) systematicsAB->SetCentrality("1020");
   else if ( cc == k020 ) systematicsAB->SetCentrality("020");
-  else if ( cc == k2040 ) {
+  else if ( cc == k2040 || cc == k2030 || cc == k3040 ) {
     systematicsAB->SetCentrality("2040");
     systematicsAB->SetIsPbPb2010EnergyScan(true);
   }
-  else if ( cc == k4060 ) systematicsAB->SetCentrality("4060");
+  else if ( cc == k4060 || cc == k4050 || cc == k5060 ) systematicsAB->SetCentrality("4060");
   else if ( cc == k6080 ) systematicsAB->SetCentrality("6080");
   else if ( cc == k4080 ) systematicsAB->SetCentrality("4080");
   else if ( cc == k3050 ) {