Skip to content

Commit 16272dd

Browse files
jinhyunnijinhyunni
committed
[Xic0Omegac0 Update] First draft of updated Xic0Omegac0 DCAFitte & KFPartice integrated workflow
1 parent 763dd95 commit 16272dd

File tree

8 files changed

+2837
-2211
lines changed

8 files changed

+2837
-2211
lines changed

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 79 additions & 119 deletions
Original file line numberDiff line numberDiff line change
@@ -842,10 +842,10 @@ static const std::vector<std::string> labelsPt = {
842842
static const std::vector<std::string> labelsCutVar = {"m", "pT p", "pT K", "pT Pi", "chi2PCA", "decay length", "cos pointing angle", "decLengthXY", "normDecLXY", "ct", "impParXY"};
843843
} // namespace hf_cuts_xic_to_p_k_pi
844844

845-
namespace hf_cuts_xic_to_xi_pi_pi
845+
namespace hf_cuts_to_xi_pi
846846
{
847847
static constexpr int NBinsPt = 13;
848-
static constexpr int NCutVars = 13;
848+
static constexpr int NCutVars = 26;
849849
// default values for the pT bin edges (can be used to configure histogram axis)
850850
// offset by 1 from the bin numbers in cuts array
851851
constexpr double BinsPt[NBinsPt + 1] = {
@@ -865,20 +865,22 @@ constexpr double BinsPt[NBinsPt + 1] = {
865865
20.};
866866
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
867867

868-
// default values for the cuts m Y Eta EtaPi EtaXi pT Pi0 Pi1 Sum chi2SV dL dLXY invMass Xi-Pi pairs
869-
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
870-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
871-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
872-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
873-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
874-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
875-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
876-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
877-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
878-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
879-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
880-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
881-
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
868+
// default values for the cuts
869+
constexpr double Cuts[NBinsPt][NCutVars] = {
870+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
871+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
872+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
873+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
874+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
875+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
876+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
877+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
878+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
879+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
880+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
881+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
882+
{2.0, 3.1, 0.01, 0.01, 0.2, 0.15, 1.0, 1.0, 2.0, 0.06, 0.06, 0.04, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.0, 10.0, 0.97, 0.97, 0.6, 1.2, 1.0, 0.8},
883+
};
882884

883885
// row labels
884886
static const std::vector<std::string> labelsPt = {
@@ -897,105 +899,76 @@ static const std::vector<std::string> labelsPt = {
897899
"pT bin 12"};
898900

899901
// column labels
900-
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Pi from XicPlus", "eta Xi Daughters", "pT Pi0", "pT Pi1", "pT Pi0 + Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
901-
} // namespace hf_cuts_xic_to_xi_pi_pi
902-
903-
namespace hf_cuts_xic0_xicp_to_hadronic
904-
{
905-
// cuts for xic0
906-
namespace xic0
907-
{
908-
static constexpr int NBinsPt = 10;
909-
static constexpr int NCutVars = 10;
910-
// default values for the pT bin edges
911-
// offset by 1 from the bin numbers in cuts array
912-
constexpr double BinsPt[NBinsPt+1] = {
913-
0.,
914-
1.,
915-
2.,
916-
3.,
917-
4.,
918-
5.,
919-
6.,
920-
8.,
921-
12.,
922-
24.,
923-
36.
902+
static const std::vector<std::string> labelsCutVar = {
903+
"mCharmBaryonMin",
904+
"mCharmBaryonMax",
905+
"massWindowCascade",
906+
"massWindowV0", // Invariant mass
907+
"ptCharmBachelorMin",
908+
"ptBachelorMin", // pt of daughter
909+
"dcaCascDauMax",
910+
"dcaV0DauMax",
911+
"dcaCharmBaryonDauMax", // dca between daughters
912+
"dcaxyV0PosDauToPvMin",
913+
"dcaxyV0NegDauToPvMin",
914+
"dcaxyBachToPvMin", // DCAXY(impact parameter)
915+
"impactParXYCharmBachelorMin",
916+
"impactParXYCharmBachelorMax",
917+
"impactParZCharmBachelorMin",
918+
"impactParZCharmBachelorMax", // DCAXY(impact parameter)
919+
"impactParXYCascMin",
920+
"impactParXYCascMax",
921+
"impactParZCascMin",
922+
"impactParZCascMax", // DCAXY(impact parameter)
923+
"cosPaCascMin",
924+
"cosPaV0Min", // CosPa
925+
"radiusCascMin",
926+
"radiusV0Min", // radius
927+
"etaTrackLFDauMax",
928+
"etaTrackCharmBachMax" // eta cut to daughter tracks
924929
};
925-
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt+NBinsPt+1};
926-
927-
// default values for the cuts m ptXi ptPi chi2PCA dL dLXY cosp cospXY impParXY Xi Pi
928-
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 0 < pT < 1 */
929-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 1 < pT < 2 */
930-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 2 < pT < 3 */
931-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 3 < pT < 4 */
932-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 4 < pT < 5 */
933-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 5 < pT < 6 */
934-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 6 < pT < 8 */
935-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 8 < pT < 12 */
936-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 12 < pT < 24 */
937-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}}; /* 24 < pT < 36 */
938-
939-
// row tables
940-
static const std::vector<std::string> labelsPt = {
941-
"pT bin 0",
942-
"pT bin 1",
943-
"pT bin 2",
944-
"pT bin 3",
945-
"pT bin 4",
946-
"pT bin 5",
947-
"pT bin 6",
948-
"pT bin 7",
949-
"pT bin 8",
950-
"pT bin 9"};
930+
} // namespace hf_cuts_to_xi_pi
951931

952-
// column labels
953-
static const std::vector<std::string> labelsCutVar = {
954-
"m",
955-
"pT Xi",
956-
"pT Pi",
957-
"chi2PCA",
958-
"max decay length",
959-
"max decay length XY",
960-
"cosine pointing angle",
961-
"cosine pointing angle XY",
962-
"max impParXY Xi",
963-
"max impParXY Pi"};
964-
}// namspace hf_cuts_xic0_xicp_to_hadronic::xic0
965-
966-
namespace xicp
932+
namespace hf_cuts_xic_to_xi_pi_pi
967933
{
968-
static constexpr int NBinsPt= 10;
969-
static constexpr int NCutVars= 12;
934+
static constexpr int NBinsPt = 13;
935+
static constexpr int NCutVars = 13;
970936
// default values for the pT bin edges (can be used to configure histogram axis)
971937
// offset by 1 from the bin numbers in cuts array
972-
constexpr double BinsPt[NBinsPt+ 1] = {
938+
constexpr double BinsPt[NBinsPt + 1] = {
973939
0.,
974940
1.,
975941
2.,
976942
3.,
977943
4.,
978944
5.,
979945
6.,
946+
7.,
980947
8.,
948+
9.,
949+
10.,
950+
11.,
981951
12.,
982-
24.,
983-
36.};
984-
auto vecBinsPt= std::vector<double>{BinsPt, BinsPt+ NBinsPt+ 1};
985-
986-
// default values for the cuts m ptXi ptPi0 ptPi1 chi2PCA dL dLXY cosp cospXY impParXY Xi Pi0 Pi1
987-
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 */
988-
{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 */
989-
{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 */
990-
{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 */
991-
{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 */
992-
{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 */
993-
{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 */
994-
{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 */
995-
{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 */
996-
{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 */
952+
20.};
953+
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
954+
955+
// default values for the cuts m Y Eta EtaPi EtaXi pT Pi0 Pi1 Sum chi2SV dL dLXY invMass Xi-Pi pairs
956+
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
957+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
958+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
959+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
960+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
961+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
962+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
963+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
964+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
965+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
966+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
967+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
968+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
969+
997970
// row labels
998-
static const std::vector<std::string> labelsPt= {
971+
static const std::vector<std::string> labelsPt = {
999972
"pT bin 0",
1000973
"pT bin 1",
1001974
"pT bin 2",
@@ -1005,27 +978,14 @@ static const std::vector<std::string> labelsPt= {
1005978
"pT bin 6",
1006979
"pT bin 7",
1007980
"pT bin 8",
1008-
"pT bin 9"};
981+
"pT bin 9",
982+
"pT bin 10",
983+
"pT bin 11",
984+
"pT bin 12"};
1009985

1010986
// column labels
1011-
static const std::vector<std::string> labelsCutVar= {
1012-
"m",
1013-
"pT Xi",
1014-
"pT Pi0",
1015-
"pT Pi1",
1016-
"chi2PCA",
1017-
"max decay length",
1018-
"max decay length XY",
1019-
"cos pointing angle",
1020-
"cos pointing angle XY",
1021-
"max impParXY Xi",
1022-
"max impParXY Pi0",
1023-
"max impParXY Pi1"}; // end Xicp settings
1024-
1025-
}// namspace hf_cuts_xic0_xicp_to_hadronic::xicp
1026-
1027-
}// namespace hf_cuts_xic0_xicp_to_hadronic
1028-
987+
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Pi from XicPlus", "eta Xi Daughters", "pT Pi0", "pT Pi1", "pT Pi0 + Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
988+
} // namespace hf_cuts_xic_to_xi_pi_pi
1029989

1030990
namespace hf_cuts_xicc_to_p_k_pi_pi
1031991
{

‎PWGHF/D2H/Tasks/CMakeLists.txt‎

Lines changed: 0 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -139,11 +139,6 @@ o2physics_add_dpl_workflow(task-xic-to-xi-pi-pi
139139
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
140140
COMPONENT_NAME Analysis)
141141

142-
o2physics_add_dpl_workflow(task-xic0-xicp-to-hadronic
143-
SOURCES taskXic0XicpToHadronic.cxx
144-
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
145-
COMPONENT_NAME Analysis)
146-
147142
o2physics_add_dpl_workflow(task-xicc
148143
SOURCES taskXicc.cxx
149144
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore

0 commit comments

Comments
 (0)