From 2784671732ca5044cabd2324911d8bce6e8fe76d Mon Sep 17 00:00:00 2001 From: vsitaraman Date: Tue, 11 Aug 2026 14:19:29 -0400 Subject: [PATCH] modified: TrackRecon.C PID gate on the suspected proton locus and prudcing the Ex function corresponding to the state. Also implemented a1c0 dithering for the reaction data --- TrackRecon.C | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/TrackRecon.C b/TrackRecon.C index df6940d..8af3d04 100644 --- a/TrackRecon.C +++ b/TrackRecon.C @@ -2077,6 +2077,8 @@ void pcCalibratedHistograms(HistPlotter *plotter, const std::vector &QQQ_ { for (const auto &pcevent : PC_Events_calibrated) { + if (pcevent.multi1 > 2 || pc.multi2 > 4) + continue; const std::string topo = "_a" + std::to_string(pcevent.multi1) + "c" + std::to_string(pcevent.multi2); const bool hasCathode = (pcevent.Cathodech >= 0); if (hasCathode) @@ -4071,7 +4073,7 @@ static void reaction_ax_core(HistPlotter *plotter, const std::vector &Si_ double level_residual = 0.0; // double snapped_level = snapToNearestLevel(Ex, levels, level_residual); double ebeam_kin_MeV = (Ex < ex_gate_MeV) - ? invertBeamEnergyMeV(m_beam, mass_4He, m3, m4, Efix, theta * 180 / M_PI, Ex) + ? invertBeamEnergyMeV(m_beam, mass_4He, m3, m4, Efix, theta * 180 / M_PI, 0) : -1.0; // double ebeam_kin_MeV = invertBeamEnergyMeV(m_beam, mass_4He, m3, m4, Efix, theta * 180 / M_PI, snapped_level);