From 4f234aa8bddf7d1cd2817c55647f24a0367136d9 Mon Sep 17 00:00:00 2001
From: "Stephen A. Wood" <saw@jlab.org>
Date: Thu, 19 Feb 2015 13:58:54 -0500
Subject: [PATCH] Correct off by one errors in time histogram to replicate
ENGINE
---
src/THcHodoscope.cxx | 2 +-
src/THcScintillatorPlane.cxx | 4 ++--
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/THcHodoscope.cxx b/src/THcHodoscope.cxx
index c49a318..8ab6845 100644
--- a/src/THcHodoscope.cxx
+++ b/src/THcHodoscope.cxx
@@ -927,7 +927,7 @@ Int_t THcHodoscope::FineProcess( TClonesArray& tracks )
}
}
- if ( fJMax >= 0 ){ // Line 248. Here I followed the code of THcSCintilaltorPlane::PulseHeightCorrection
+ if ( fJMax > 0 ){ // Line 248. Here I followed the code of THcSCintilaltorPlane::PulseHeightCorrection
fTMin = 0.5 * fJMax;
for(Int_t iphit = 0; iphit < fNScinHits[ip]; iphit++) { // Loop over sinc. hits. in plane
if ( ( fTOFPInfo[iphit].time_pos > fTMin ) && ( fTOFPInfo[iphit].time_pos < ( fTMin + fTofTolerance ) ) ) {
diff --git a/src/THcScintillatorPlane.cxx b/src/THcScintillatorPlane.cxx
index 69be89f..afca7a3 100644
--- a/src/THcScintillatorPlane.cxx
+++ b/src/THcScintillatorPlane.cxx
@@ -444,7 +444,7 @@ Int_t THcScintillatorPlane::PulseHeightCorrection()
time_neg[i]=negtime[i]-(fZpos+(j%2)*fDzpos)/(29.979*hbeta_pcent);
nfound++;
for (int k=0;k<200;k++) {
- tmin=0.5*k;
+ tmin=0.5*(k+1);
if ((time_pos[i]> tmin) && (time_pos[i] < tmin+toftolerance)) {
timehist[k]++;
}
@@ -471,7 +471,7 @@ Int_t THcScintillatorPlane::PulseHeightCorrection()
(((THcSignalHit*) fPosTDCHits->At(i))->GetData()<=maxtdc) &&
(((THcSignalHit*) fNegTDCHits->At(i))->GetData()>=mintdc) &&
(((THcSignalHit*) fNegTDCHits->At(i))->GetData()<=maxtdc)) {
- if(jmax>=0) {
+ if(jmax>0) {
tmin = 0.5*jmax;
if ((time_pos[i]>tmin) && (time_pos[i]<tmin+toftolerance) &&
(time_neg[i]>tmin) && (time_neg[i]<tmin+toftolerance))
--
GitLab