Skip to content

Commit 0dae45f

Browse files
jinhyunniJinhyun Park
authored andcommitted
[Xic0Xicp] Xicp added to selector and tree creator
1 parent 6558bb2 commit 0dae45f

File tree

3 files changed

+1163
-99
lines changed

3 files changed

+1163
-99
lines changed

PWGHF/Core/SelectorCuts.h

Lines changed: 75 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -846,23 +846,24 @@ static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Daug
846846
namespace hf_cuts_xic0_xicp_to_hadronic
847847
{
848848
// cuts for xic0
849-
849+
namespace xic0
850+
{
850851
static constexpr int NBinsPt = 10;
851852
static constexpr int NCutVars = 10;
852853
// default values for the pT bin edges
853854
// offset by 1 from the bin numbers in cuts array
854855
constexpr double BinsPt[NBinsPt+1] = {
855-
0.,
856-
1.,
857-
2.,
858-
3.,
859-
4.,
860-
5.,
861-
6.,
862-
8.,
863-
12.,
864-
24.,
865-
36.
856+
0.,
857+
1.,
858+
2.,
859+
3.,
860+
4.,
861+
5.,
862+
6.,
863+
8.,
864+
12.,
865+
24.,
866+
36.
866867
};
867868
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt+NBinsPt+1};
868869

@@ -903,6 +904,68 @@ static const std::vector<std::string> labelsCutVar = {
903904
"cosine pointing angle XY",
904905
"max impParXY Xi",
905906
"max impParXY Pi"};
907+
}// namspace hf_cuts_xic0_xicp_to_hadronic::xic0
908+
909+
namespace xicp
910+
{
911+
static constexpr int NBinsPt= 10;
912+
static constexpr int NCutVars= 12;
913+
// default values for the pT bin edges (can be used to configure histogram axis)
914+
// offset by 1 from the bin numbers in cuts array
915+
constexpr double BinsPt[NBinsPt+ 1] = {
916+
0.,
917+
1.,
918+
2.,
919+
3.,
920+
4.,
921+
5.,
922+
6.,
923+
8.,
924+
12.,
925+
24.,
926+
36.};
927+
auto vecBinsPt= std::vector<double>{BinsPt, BinsPt+ NBinsPt+ 1};
928+
929+
// default values for the cuts m ptXi ptPi0 ptPi1 chi2PCA dL dLXY cosp cospXY impParXY Xi Pi0 Pi1
930+
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 */
931+
{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 */
932+
{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 */
933+
{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 */
934+
{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 */
935+
{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 */
936+
{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 */
937+
{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 */
938+
{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 */
939+
{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 */
940+
// row labels
941+
static const std::vector<std::string> labelsPt= {
942+
"pT bin 0",
943+
"pT bin 1",
944+
"pT bin 2",
945+
"pT bin 3",
946+
"pT bin 4",
947+
"pT bin 5",
948+
"pT bin 6",
949+
"pT bin 7",
950+
"pT bin 8",
951+
"pT bin 9"};
952+
953+
// column labels
954+
static const std::vector<std::string> labelsCutVar= {
955+
"m",
956+
"pT Xi",
957+
"pT Pi0",
958+
"pT Pi1",
959+
"chi2PCA",
960+
"max decay length",
961+
"max decay length XY",
962+
"cos pointing angle",
963+
"cos pointing angle XY",
964+
"max impParXY Xi",
965+
"max impParXY Pi0",
966+
"max impParXY Pi1"}; // end Xicp settings
967+
968+
}// namspace hf_cuts_xic0_xicp_to_hadronic::xicp
906969

907970
}// namespace hf_cuts_xic0_xicp_to_hadronic
908971

0 commit comments

Comments
 (0)