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

modified: TrackRecon.C

Browse Source 
modified:   run_tr.sh
This commit is contained in:
Vignesh Sitaraman 2026-09-06 10:43:55 -04:00
parent fde75b8a77
commit d3c4c47e41
2 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
TrackRecon.C
View File

@ -3611,7 +3611,7 @@ void protonAlphaElastic_core(HistPlotter *plotter, const std::vector<Event> &Si_
std::string ejtag;
if (alphaHyp && sievent.Energy1 < 5)
ejtag = "_a";
else if (sievent.Energy1 >= 6.4 && sievent.Energy1 < 7.0)
else if (sievent.Energy1 >= 5.0 && sievent.Energy1 < 7.0)
ejtag = "_p";
else
ejtag = "_maybep";
@ -3632,6 +3632,16 @@ void protonAlphaElastic_core(HistPlotter *plotter, const std::vector<Event> &Si_
plotter->Fill2D(rx + "_Ex_vs_theta" + ejtag + sfx, 360, 0, 180, 800, -10, 10, theta * 180 / M_PI, Ex, pmlabel);
plotter->Fill2D(rx + "_Ex_vs_phi" + ejtag + sfx, 180, -180, 180, 800, -10, 10, sievent.pos.Phi() * 180 / M_PI, Ex, pmlabel);
for (const auto &pcevent : PC_Events)
{
if (!(pcevent.multi1 == 1 && (pcevent.multi2 == 1 || pcevent.multi2 == 2)))
continue;
if (TMath::Abs(sievent.pos.DeltaPhi(pcevent.pos)) > phi_win)
continue;
plotter->Fill2D(rx + "_Ex_vs_dT" + ejtag + sfx, 500, -2000, 2000, 600, -10, 20, (sievent.Time1 - pcevent.Time1), Ex, pmlabel);
plotter->Fill2D(rx + "_dEgasCalib_vs_dT" + ejtag + sfx, 500, -2000, 2000, 800, 0, 0.6, (sievent.Time1 - pcevent.Time1), anodeE_MeV, pmlabel);
}
// Ground-state beam-energy consistency check -- elastic scattering has
// no excited levels, so there's only ever a "ground state" hypothesis
// here, unlike the (a,p) reaction branch's snapped levels.

2
run_tr.sh
View File

@ -55,7 +55,7 @@ if [[ 1 -eq 0 ]]; then
fi
# --- Block 3: 27Al Alpha+Gas Runs (9, 12) ---
if [[ 1 -eq 1 ]]; then
if [[ 1 -eq 0 ]]; then
export DATASET="27Al"
export PREFIX="Run_"
export OUT_DIR="Output_a"