[Xic0Xicp] Xicp added to selector and tree creator

Author: jinhyunni

PWGHF/Core/SelectorCuts.h

@@ -803,23 +803,24 @@ static const std::vector<std::string> labelsCutVar = {"m", "pT Xi", "pT Pi0", "p
 namespace hf_cuts_xic0_xicp_to_hadronic
 {
 // cuts for xic0
-
+namespace xic0
+{
 static constexpr int NBinsPt = 10;
 static constexpr int NCutVars = 10;
 // default values for the pT bin edges
 // offset by 1 from the bin numbers in cuts array
 constexpr double BinsPt[NBinsPt+1] = {
-	0.,
-	1.,
-	2.,
-	3.,
-	4.,
-	5.,
-	6.,
-	8.,
-	12.,
-	24.,
-	36.
+  0.,
+  1.,
+  2.,
+  3.,
+  4.,
+  5.,
+  6.,
+  8.,
+  12.,
+  24.,
+  36.
 };
 auto vecBinsPt = std::vector<double>{BinsPt, BinsPt+NBinsPt+1};
 
@@ -860,6 +861,68 @@ static const std::vector<std::string> labelsCutVar = {
 	"cosine pointing angle XY",
 	"max impParXY Xi",
 	"max impParXY Pi"};
+}// namspace hf_cuts_xic0_xicp_to_hadronic::xic0
+
+namespace xicp
+{
+static constexpr int NBinsPt= 10;
+static constexpr int NCutVars= 12;
+// default values for the pT bin edges (can be used to configure histogram axis)
+// offset by 1 from the bin numbers in cuts array
+constexpr double BinsPt[NBinsPt+ 1] = {
+  0.,
+  1.,
+  2.,
+  3.,
+  4.,
+  5.,
+  6.,
+  8.,
+  12.,
+  24.,
+  36.};
+auto vecBinsPt= std::vector<double>{BinsPt, BinsPt+ NBinsPt+ 1};
+
+// default values for the cuts               m   ptXi ptPi0 ptPi1 chi2PCA dL dLXY cosp cospXY impParXY Xi Pi0 Pi1
+constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 0  < pT < 1  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 1  < pT < 2  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 2  < pT < 3  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 3  < pT < 4  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 4  < pT < 5  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 5  < pT < 6  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 6  < pT < 8  */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 8  < pT < 10 */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1},  /* 12 < pT < 24 */
+                                            {0.4, 0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1, 0.1}}; /* 24 < pT < 36 */
+// row labels
+static const std::vector<std::string> labelsPt= {
+  "pT bin 0",
+  "pT bin 1",
+  "pT bin 2",
+  "pT bin 3",
+  "pT bin 4",
+  "pT bin 5",
+  "pT bin 6",
+  "pT bin 7",
+  "pT bin 8",
+  "pT bin 9"};
+
+// column labels
+static const std::vector<std::string> labelsCutVar= {
+  "m", 
+  "pT Xi", 
+  "pT Pi0", 
+  "pT Pi1", 
+  "chi2PCA", 
+  "max decay length", 
+  "max decay length XY", 
+  "cos pointing angle", 
+  "cos pointing angle XY", 
+  "max impParXY Xi", 
+  "max impParXY Pi0", 
+  "max impParXY Pi1"}; // end Xicp settings
+
+}// namspace hf_cuts_xic0_xicp_to_hadronic::xicp
 
 }// namespace hf_cuts_xic0_xicp_to_hadronic