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

modified: QQQ_Calcheck.C made changes to bring it in line with the new consistent Wedge-Ring mapping

Browse Source 
modified:   TrackRecon.C
This commit is contained in:
Vignesh Sitaraman 2026-01-25 14:15:45 -05:00
parent c4b543bdeb
commit ac035370b4
2 changed files with 14 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
QQQ_Calcheck.C
View File

@ -133,10 +133,10 @@ Bool_t QQQ_Calcheck::Process(Long64_t entry)
else else
continue; continue;
// plug in calibrations // plug in calibrations
if (qqqCalibValid[qqq.id[i]][chRing][chWedge]) if (qqqCalibValid[qqq.id[i]][chWedge][chRing])
{ {
eWedgeMeV = eWedge * qqqCalib[qqq.id[i]][chRing][chWedge] / 1000; eWedgeMeV = eWedge * qqqCalib[qqq.id[i]][chWedge][chRing] / 1000;
eRingMeV = eRing * qqqCalib[qqq.id[i]][chRing][chWedge] / 1000; eRingMeV = eRing * qqqCalib[qqq.id[i]][chWedge][chRing] / 1000;
} }
else else
continue; continue;

19
TrackRecon.C
View File

@ -146,8 +146,8 @@ void TrackRecon::Begin(TTree * /*tree*/)
double gainw, gainr; double gainw, gainr;
while (infile >> det >> wedge >> ring >> gainw >> gainr) while (infile >> det >> wedge >> ring >> gainw >> gainr)
{ {
qqqGain[det][ring][wedge] = gainw; qqqGain[det][wedge][ring] = gainw;
qqqGainValid[det][ring][wedge] = (gainw > 0); qqqGainValid[det][wedge][ring] = (gainw > 0);
// std::cout << "QQQ Gain Loaded: Det " << det << " Ring " << ring << " Wedge " << wedge << " GainW " << gainw << " GainR " << gainr << std::endl; // std::cout << "QQQ Gain Loaded: Det " << det << " Ring " << ring << " Wedge " << wedge << " GainW " << gainw << " GainR " << gainr << std::endl;
} }
infile.close(); infile.close();
@ -162,8 +162,8 @@ void TrackRecon::Begin(TTree * /*tree*/)
double slope; double slope;
while (infile >> det >> wedge >> ring >> slope) while (infile >> det >> wedge >> ring >> slope)
{ {
qqqCalib[det][ring][wedge] = slope; qqqCalib[det][wedge][ring] = slope;
qqqCalibValid[det][ring][wedge] = (slope > 0); qqqCalibValid[det][wedge][ring] = (slope > 0);
// std::cout << "QQQ Calib Loaded: Det " << det << " Ring " << ring << " Wedge " << wedge << " Slope " << slope << std::endl; // std::cout << "QQQ Calib Loaded: Det " << det << " Ring " << ring << " Wedge " << wedge << " Slope " << slope << std::endl;
} }
infile.close(); infile.close();
@ -278,14 +278,14 @@ Bool_t TrackRecon::Process(Long64_t entry)
else else
continue; continue;
plotter->Fill1D("Wedgetime_Vs_Ringtime", 200, -1000, 1000, tWedge - tRing, "hCalQQQ"); plotter->Fill1D("Wedgetime_Vs_Ringtime", 100, -1000, 1000, tWedge - tRing, "hCalQQQ");
plotter->Fill2D("RingE_vs_Index", 16 * 4, 0, 16 * 4, 1000, 0, 16000, chRing + qqq.id[i] * 16, eRing, "hRawQQQ"); plotter->Fill2D("RingE_vs_Index", 16 * 4, 0, 16 * 4, 1000, 0, 16000, chRing + qqq.id[i] * 16, eRing, "hRawQQQ");
plotter->Fill2D("WedgeE_vs_Index", 16 * 4, 0, 16 * 4, 1000, 0, 16000, chWedge + qqq.id[i] * 16, eWedge, "hRawQQQ"); plotter->Fill2D("WedgeE_vs_Index", 16 * 4, 0, 16 * 4, 1000, 0, 16000, chWedge + qqq.id[i] * 16, eWedge, "hRawQQQ");
if (qqqCalibValid[qqq.id[i]][chRing][chWedge]) if (qqqCalibValid[qqq.id[i]][chWedge][chRing])
{ {
eWedgeMeV = eWedge * qqqCalib[qqq.id[i]][chRing][chWedge] / 1000; eWedgeMeV = eWedge * qqqCalib[qqq.id[i]][chWedge][chRing] / 1000;
eRingMeV = eRing * qqqCalib[qqq.id[i]][chRing][chWedge] / 1000; eRingMeV = eRing * qqqCalib[qqq.id[i]][chWedge][chRing] / 1000;
} }
else else
continue; continue;
@ -317,6 +317,8 @@ Bool_t TrackRecon::Process(Long64_t entry)
// // } // // }
plotter->Fill2D("Timing_Difference_QQQ_PC", 500, -1000, 1000, 16, 0, 16, tRing - static_cast<double>(pc.t[k]), chRing, "hCalQQQ"); plotter->Fill2D("Timing_Difference_QQQ_PC", 500, -1000, 1000, 16, 0, 16, tRing - static_cast<double>(pc.t[k]), chRing, "hCalQQQ");
plotter->Fill2D("DelT_Vs_QQQRingECal", 500, -1000, 1000, 1000, 0, 10, tRing - static_cast<double>(pc.t[k]), eRingMeV, "hCalQQQ"); plotter->Fill2D("DelT_Vs_QQQRingECal", 500, -1000, 1000, 1000, 0, 10, tRing - static_cast<double>(pc.t[k]), eRingMeV, "hCalQQQ");
plotter->Fill2D("CalibratedQQQEvsPCE", 1000, 0, 10, 2000, 0, 30000, eRingMeV, pc.e[k], "hGMPC");
plotter->Fill2D("CalibratedQQQEvsPCE", 1000, 0, 10, 2000, 0, 30000, eWedgeMeV, pc.e[k], "hGMPC");
if (tRing - static_cast<double>(pc.t[k]) < -150 && tRing - static_cast<double>(pc.t[k]) > -450) // 27Al if (tRing - static_cast<double>(pc.t[k]) < -150 && tRing - static_cast<double>(pc.t[k]) > -450) // 27Al
// if (tRing - static_cast<double>(pc.t[k]) < -70 && tRing - static_cast<double>(pc.t[k]) > -150) // 17F // if (tRing - static_cast<double>(pc.t[k]) < -70 && tRing - static_cast<double>(pc.t[k]) > -150) // 17F
{ {
@ -326,6 +328,7 @@ Bool_t TrackRecon::Process(Long64_t entry)
if (pc.index[k] >= 24 && pc.e[k] > 50) if (pc.index[k] >= 24 && pc.e[k] > 50)
{ {
plotter->Fill2D("Timing_Difference_QQQ_PC_Cathode", 500, -1000, 1000, 16, 0, 16, tRing - static_cast<double>(pc.t[k]), chRing, "hCalQQQ"); plotter->Fill2D("Timing_Difference_QQQ_PC_Cathode", 500, -1000, 1000, 16, 0, 16, tRing - static_cast<double>(pc.t[k]), chRing, "hCalQQQ");
} }
} }