Skip to content

Commit 0efcf10

Browse files
Refactor decay radius calculation and add histogram print
Updated the decay radius calculation to use std::sqrt for better precision and added a histogram print statement for inspection.
1 parent 4974dc9 commit 0efcf10

1 file changed

Lines changed: 3 additions & 1 deletion

File tree

PWGLF/Tasks/QC/v0cascadesqa.cxx

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -445,6 +445,8 @@ struct v0cascadesQA {
445445
histos.add("histos_Casc/InvMassOmegaPlusTrue", "InvMassOmegaPlusTrue", {HistType::kTH3F, {{100, 0.f, 10.f}, {100, 0.f, 50.f}, {80, 1.63f, 1.71f}}});
446446
histos.add("histos_Casc/InvMassOmegaMinusTrue", "InvMassOmegaMinusTrue", {HistType::kTH3F, {{100, 0.f, 10.f}, {100, 0.f, 50.f}, {80, 1.63f, 1.71f}}});
447447
}
448+
// inspect histogram sizes, please
449+
histos.print();
448450

449451
// Initialise the RCTFlagsChecker
450452
rctFlagsChecker.init(rctConfigurations.cfgRCTLabel.value, rctConfigurations.cfgCheckZDC, rctConfigurations.cfgTreatLimitedAcceptanceAsBad);
@@ -1398,7 +1400,7 @@ struct v0cascadesQA {
13981400
if (vx != 0 && vy != 0)
13991401
break;
14001402
}
1401-
double R_Decay = TMath::Sqrt(vx * vx + vy * vy);
1403+
double R_Decay = std::sqrt(vx * vx + vy * vy);
14021404

14031405
if (mcparticle.isPhysicalPrimary() && std::abs(mcparticle.y()) < v0Selections.rapidityCut) {
14041406
if (mcparticle.pdgCode() == PDG_t::kK0Short)

0 commit comments

Comments
 (0)