rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

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

Browse Source
This commit is contained in:
Vignesh Sitaraman 2026-08-11 14:19:29 -04:00
parent 774050dfdf
commit 2784671732
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
TrackRecon.C
View File

@ -2077,6 +2077,8 @@ void pcCalibratedHistograms(HistPlotter *plotter, const std::vector<Event> &QQQ_
{ {
for (const auto &pcevent : PC_Events_calibrated) 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 std::string topo = "_a" + std::to_string(pcevent.multi1) + "c" + std::to_string(pcevent.multi2);
const bool hasCathode = (pcevent.Cathodech >= 0); const bool hasCathode = (pcevent.Cathodech >= 0);
if (hasCathode) if (hasCathode)
@ -4071,7 +4073,7 @@ static void reaction_ax_core(HistPlotter *plotter, const std::vector<Event> &Si_
double level_residual = 0.0; double level_residual = 0.0;
// double snapped_level = snapToNearestLevel(Ex, levels, level_residual); // double snapped_level = snapToNearestLevel(Ex, levels, level_residual);
double ebeam_kin_MeV = (Ex < ex_gate_MeV) 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; : -1.0;
// double ebeam_kin_MeV = invertBeamEnergyMeV(m_beam, mass_4He, m3, m4, Efix, theta * 180 / M_PI, snapped_level); // double ebeam_kin_MeV = invertBeamEnergyMeV(m_beam, mass_4He, m3, m4, Efix, theta * 180 / M_PI, snapped_level);