Skip to content

Commit 2a5ffc7

Browse files
committed
Add missing THnSparse for delta Eta, Phi and qinv, kT
1 parent 06e3aa1 commit 2a5ffc7

File tree

1 file changed

+22
-10
lines changed

1 file changed

+22
-10
lines changed

PWGEM/PhotonMeson/Tasks/photonhbt.cxx

Lines changed: 22 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -633,12 +633,13 @@ struct photonhbt {
633633
}
634634

635635
fRegistry.add("Pair/same/hDeltaRCosOA", "distance between 2 conversion points / cos(#theta_{op}/2);#Delta r / cos(#theta_{op}/2) (cm);counts", kTH1D, {{100, 0, 100}}, true);
636-
fRegistry.add("Pair/same/hSparse_DEtaDPhi_kT", "same-event (#Delta#eta,#Delta#phi,q_{inv},k_{T}) for efficiency reweighting;" "#Delta#eta;#Delta#phi (rad);q_{inv} (GeV/c);k_{T} (GeV/c)",
636+
fRegistry.add("Pair/same/hSparse_DEtaDPhi_kT",
637+
"same-event (#Delta#eta,#Delta#phi,q_{inv},k_{T}) for efficiency reweighting;"
638+
"#Delta#eta;#Delta#phi (rad);q_{inv} (GeV/c);k_{T} (GeV/c)",
637639
kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
638640

639-
fRegistry.add("Pair/same/hPhi_lowerPtV0",
640-
"azimuthal angle of lower-p_{T} V0 in pair;#phi (rad);counts",
641-
kTH1D, {axisPhi}, true);
641+
fRegistry.add("Pair/same/hPhi_lowerPtV0", "azimuthal angle of lower-p_{T} V0 in pair;#phi (rad);counts", kTH1D, {axisPhi}, true);
642+
fRegistry.add("Pair/same/hSparse_DEtaDPhi_qinv_kT", "azimuthal angle of lower-p_{T} V0 in pair;#phi (rad);counts", kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisQinv, axisKt}, true);
642643

643644
addQAHistogramsForStep("Pair/same/QA/Before/");
644645
addQAHistogramsForStep("Pair/same/QA/AfterDRCosOA/");
@@ -651,10 +652,22 @@ struct photonhbt {
651652
addFullRangeHistograms("Pair/same/FullRange/");
652653
fRegistry.addClone("Pair/same/", "Pair/mix/");
653654

654-
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_sameSideV0_sameSideLegs", "both V0 same #eta-side, all legs same side;" "#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)", kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
655-
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_sameSideV0_mixedLegs", "both V0 same #eta-side, legs mixed sides;" "#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)", kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
656-
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_diffSideV0_matchingLegs", "V0 on opposite #eta-sides, legs match their V0 side;" "#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)", kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
657-
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_diffSideV0_crossedLegs", "V0 on opposite #eta-sides, legs cross their V0 side;" "#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)", kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
655+
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_sameSideV0_sameSideLegs",
656+
"both V0 same #eta-side, all legs same side;"
657+
"#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)",
658+
kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
659+
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_sameSideV0_mixedLegs",
660+
"both V0 same #eta-side, legs mixed sides;"
661+
"#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)",
662+
kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
663+
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_diffSideV0_matchingLegs",
664+
"V0 on opposite #eta-sides, legs match their V0 side;"
665+
"#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)",
666+
kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
667+
fRegistry.add("Pair/same/EtaTopology/hSparse_DEtaDPhi_kT_diffSideV0_crossedLegs",
668+
"V0 on opposite #eta-sides, legs cross their V0 side;"
669+
"#Delta#eta;#Delta#phi (rad);k_{T} (GeV/c)",
670+
kTHnSparseD, {axisDeltaEta, axisDeltaPhi, axisKt}, true);
658671
}
659672

660673
void DefineEMEventCut()
@@ -1073,7 +1086,7 @@ struct photonhbt {
10731086

10741087
// fRegistryPairMC — reco-level pair histograms, per MC truth type
10751088

1076-
// Per-type CF + observables
1089+
// Per-type CF + observables
10771090
static constexpr std::array<std::string_view, 6> kTypes = {
10781091
"TrueTrueDistinct/", "TrueTrueSamePhoton/", "SharedMcLeg/",
10791092
"TrueFake/", "FakeFake/", "Pi0Daughters/"};
@@ -1094,7 +1107,6 @@ struct photonhbt {
10941107
kTH2D, {hbtanalysis.cfgUseLCMS ? axisQabsLcms : axisQinv, axisKt}, true);
10951108
}
10961109

1097-
10981110
// 1D observables
10991111
fRegistryPairMC.add((base + "hQinv").c_str(), "q_{inv};q_{inv} (GeV/c);counts", kTH1D, {axisQinv}, true);
11001112
fRegistryPairMC.add((base + "hDeltaR").c_str(), "|R_{1}-R_{2}|;|R_{1}-R_{2}| (cm);counts", kTH1D, {axisDeltaR}, true);

0 commit comments

Comments
 (0)