Skip to content

Commit 3957158

Browse files
jinhyunnijinhyunni
committed
[Xic0Omegac0 Update] First draft of updated Xic0Omegac0 DCAFitte & KFPartice integrated workflow
1 parent fda2ca5 commit 3957158

8 files changed

Lines changed: 2837 additions & 2211 deletions

‎PWGHF/Core/SelectorCuts.h‎

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

788-
namespace hf_cuts_xic_to_xi_pi_pi
788+
namespace hf_cuts_to_xi_pi
789789
{
790790
static constexpr int NBinsPt = 13;
791-
static constexpr int NCutVars = 11;
791+
static constexpr int NCutVars = 26;
792792
// default values for the pT bin edges (can be used to configure histogram axis)
793793
// offset by 1 from the bin numbers in cuts array
794794
constexpr double BinsPt[NBinsPt + 1] = {
@@ -808,20 +808,22 @@ constexpr double BinsPt[NBinsPt + 1] = {
808808
20.};
809809
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
810810

811-
// default values for the cuts m Y Eta EtaDau ptPi0 ptPi1 chi2SV dL dLXY invMass Xi-Pi pairs
812-
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
813-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
814-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
815-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
816-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
817-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
818-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
819-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
820-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
821-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
822-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
823-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
824-
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
811+
// default values for the cuts
812+
constexpr double Cuts[NBinsPt][NCutVars] = {
813+
{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},
814+
{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},
815+
{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},
816+
{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},
817+
{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},
818+
{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},
819+
{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},
820+
{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},
821+
{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},
822+
{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},
823+
{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},
824+
{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},
825+
{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},
826+
};
825827

826828
// row labels
827829
static const std::vector<std::string> labelsPt = {
@@ -840,105 +842,76 @@ static const std::vector<std::string> labelsPt = {
840842
"pT bin 12"};
841843

842844
// column labels
843-
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Daughters", "pT Pi0", "pT Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
844-
} // namespace hf_cuts_xic_to_xi_pi_pi
845-
846-
namespace hf_cuts_xic0_xicp_to_hadronic
847-
{
848-
// cuts for xic0
849-
namespace xic0
850-
{
851-
static constexpr int NBinsPt = 10;
852-
static constexpr int NCutVars = 10;
853-
// default values for the pT bin edges
854-
// offset by 1 from the bin numbers in cuts array
855-
constexpr double BinsPt[NBinsPt+1] = {
856-
0.,
857-
1.,
858-
2.,
859-
3.,
860-
4.,
861-
5.,
862-
6.,
863-
8.,
864-
12.,
865-
24.,
866-
36.
845+
static const std::vector<std::string> labelsCutVar = {
846+
"mCharmBaryonMin",
847+
"mCharmBaryonMax",
848+
"massWindowCascade",
849+
"massWindowV0", // Invariant mass
850+
"ptCharmBachelorMin",
851+
"ptBachelorMin", // pt of daughter
852+
"dcaCascDauMax",
853+
"dcaV0DauMax",
854+
"dcaCharmBaryonDauMax", // dca between daughters
855+
"dcaxyV0PosDauToPvMin",
856+
"dcaxyV0NegDauToPvMin",
857+
"dcaxyBachToPvMin", // DCAXY(impact parameter)
858+
"impactParXYCharmBachelorMin",
859+
"impactParXYCharmBachelorMax",
860+
"impactParZCharmBachelorMin",
861+
"impactParZCharmBachelorMax", // DCAXY(impact parameter)
862+
"impactParXYCascMin",
863+
"impactParXYCascMax",
864+
"impactParZCascMin",
865+
"impactParZCascMax", // DCAXY(impact parameter)
866+
"cosPaCascMin",
867+
"cosPaV0Min", // CosPa
868+
"radiusCascMin",
869+
"radiusV0Min", // radius
870+
"etaTrackLFDauMax",
871+
"etaTrackCharmBachMax" // eta cut to daughter tracks
867872
};
868-
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt+NBinsPt+1};
869-
870-
// default values for the cuts m ptXi ptPi chi2PCA dL dLXY cosp cospXY impParXY Xi Pi
871-
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 */
872-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 1 < pT < 2 */
873-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 2 < pT < 3 */
874-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 3 < pT < 4 */
875-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 4 < pT < 5 */
876-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 5 < pT < 6 */
877-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 6 < pT < 8 */
878-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 8 < pT < 12 */
879-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}, /* 12 < pT < 24 */
880-
{0.4, 0.4, 0.4, 1e-5, 0.5, 0.5, 0.9, 0.9, 0.1, 0.1}}; /* 24 < pT < 36 */
881-
882-
// row tables
883-
static const std::vector<std::string> labelsPt = {
884-
"pT bin 0",
885-
"pT bin 1",
886-
"pT bin 2",
887-
"pT bin 3",
888-
"pT bin 4",
889-
"pT bin 5",
890-
"pT bin 6",
891-
"pT bin 7",
892-
"pT bin 8",
893-
"pT bin 9"};
873+
} // namespace hf_cuts_to_xi_pi
894874

895-
// column labels
896-
static const std::vector<std::string> labelsCutVar = {
897-
"m",
898-
"pT Xi",
899-
"pT Pi",
900-
"chi2PCA",
901-
"max decay length",
902-
"max decay length XY",
903-
"cosine pointing angle",
904-
"cosine pointing angle XY",
905-
"max impParXY Xi",
906-
"max impParXY Pi"};
907-
}// namspace hf_cuts_xic0_xicp_to_hadronic::xic0
908-
909-
namespace xicp
875+
namespace hf_cuts_xic_to_xi_pi_pi
910876
{
911-
static constexpr int NBinsPt= 10;
912-
static constexpr int NCutVars= 12;
877+
static constexpr int NBinsPt = 13;
878+
static constexpr int NCutVars = 11;
913879
// default values for the pT bin edges (can be used to configure histogram axis)
914880
// offset by 1 from the bin numbers in cuts array
915-
constexpr double BinsPt[NBinsPt+ 1] = {
881+
constexpr double BinsPt[NBinsPt + 1] = {
916882
0.,
917883
1.,
918884
2.,
919885
3.,
920886
4.,
921887
5.,
922888
6.,
889+
7.,
923890
8.,
891+
9.,
892+
10.,
893+
11.,
924894
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 */
895+
20.};
896+
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
897+
898+
// default values for the cuts m Y Eta EtaDau ptPi0 ptPi1 chi2SV dL dLXY invMass Xi-Pi pairs
899+
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
900+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
901+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
902+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
903+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
904+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
905+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
906+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
907+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
908+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
909+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
910+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
911+
{0.4, 0.8, 0.8, 0.8, 0.4, 0.4, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
912+
940913
// row labels
941-
static const std::vector<std::string> labelsPt= {
914+
static const std::vector<std::string> labelsPt = {
942915
"pT bin 0",
943916
"pT bin 1",
944917
"pT bin 2",
@@ -948,27 +921,14 @@ static const std::vector<std::string> labelsPt= {
948921
"pT bin 6",
949922
"pT bin 7",
950923
"pT bin 8",
951-
"pT bin 9"};
924+
"pT bin 9",
925+
"pT bin 10",
926+
"pT bin 11",
927+
"pT bin 12"};
952928

953929
// 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
969-
970-
}// namespace hf_cuts_xic0_xicp_to_hadronic
971-
930+
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Daughters", "pT Pi0", "pT Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
931+
} // namespace hf_cuts_xic_to_xi_pi_pi
972932

973933
namespace hf_cuts_xicc_to_p_k_pi_pi
974934
{

‎PWGHF/D2H/Tasks/CMakeLists.txt‎

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

137-
o2physics_add_dpl_workflow(task-xic0-xicp-to-hadronic
138-
SOURCES taskXic0XicpToHadronic.cxx
139-
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
140-
COMPONENT_NAME Analysis)
141-
142137
o2physics_add_dpl_workflow(task-xicc
143138
SOURCES taskXicc.cxx
144139
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore

0 commit comments

Comments
 (0)