From 75125f22561a18aab1a21cb40f199885a431c4ce Mon Sep 17 00:00:00 2001
From: Mark Jones <jones@jlab.org>
Date: Tue, 20 Nov 2018 09:57:19 -0500
Subject: [PATCH] Modify THcAerogel.cxx
Modify THcAerogel::CoarseProcess
For each PMT, fGoodPosAdcMult.at(npmt) was being set to the total number of positive PMTs with an
adc hit and fGoodNegAdcMult.at(npmt) was being set to the total number of negative PMTs with an
adc hit.
To be like the other dectectors make fGoodPosAdcMult.at(npmt)
and GoodNegAdcMult.at(npmt) be the number of adc hits for that particular
PMT.
---
src/THcAerogel.cxx | 21 ++++++++++++++-------
1 file changed, 14 insertions(+), 7 deletions(-)
diff --git a/src/THcAerogel.cxx b/src/THcAerogel.cxx
index 1dc6285..3b53d46 100644
--- a/src/THcAerogel.cxx
+++ b/src/THcAerogel.cxx
@@ -716,10 +716,14 @@ Int_t THcAerogel::CoarseProcess( TClonesArray& ) //tracks
Bool_t pulseTimeCut = adctdcdiffTime > fAdcPosTimeWindowMin[npmt] && adctdcdiffTime < fAdcPosTimeWindowMax[npmt];
// By default, the last hit within the timing cut will be considered "good"
- if (!errorFlag && pulseTimeCut) {
+ if (!errorFlag)
+ {
+ fGoodPosAdcMult.at(npmt) += 1;
+ }
+
+ if (!errorFlag && pulseTimeCut) {
fGoodPosAdcPed.at(npmt) = pulsePed;
- fGoodPosAdcMult.at(npmt) = frPosAdcPulseInt->GetEntries();
- // cout << " out = " << npmt << " " << frPosAdcPulseInt->GetEntries() << " " <<fGoodPosAdcMult.at(npmt);
+ // cout << " out = " << npmt << " " << frPosAdcPulseInt->GetEntries() << " " <<fGoodPosAdcMult.at(npmt);
fGoodPosAdcPulseInt.at(npmt) = pulseInt;
fGoodPosAdcPulseIntRaw.at(npmt) = pulseIntRaw;
fGoodPosAdcPulseAmp.at(npmt) = pulseAmp;
@@ -748,15 +752,18 @@ Int_t THcAerogel::CoarseProcess( TClonesArray& ) //tracks
Bool_t errorFlag = ((THcSignalHit*) fNegAdcErrorFlag->ConstructedAt(ielem))->GetData();
//// Bool_t pulseTimeCut = adctdcdiffTime > fAdcTimeWindowMin && adctdcdiffTime < fAdcTimeWindowMax;
Bool_t pulseTimeCut = adctdcdiffTime > fAdcNegTimeWindowMin[npmt] && adctdcdiffTime < fAdcNegTimeWindowMax[npmt];
-
+ if (!errorFlag)
+ {
+ fGoodNegAdcMult.at(npmt) += 1;
+ }
+
// By default, the last hit within the timing cut will be considered "good"
if (!errorFlag && pulseTimeCut) {
fGoodNegAdcPed.at(npmt) = pulsePed;
- fGoodNegAdcMult.at(npmt) = frNegAdcPulseInt->GetEntries();
- fGoodNegAdcPulseInt.at(npmt) = pulseInt;
fGoodNegAdcPulseIntRaw.at(npmt) = pulseIntRaw;
fGoodNegAdcPulseAmp.at(npmt) = pulseAmp;
- fGoodNegAdcPulseTime.at(npmt) = pulseTime;
+ fGoodNegAdcPulseInt.at(npmt) = pulseInt;
+ fGoodNegAdcPulseTime.at(npmt) = pulseTime;
fGoodNegAdcTdcDiffTime.at(npmt) = adctdcdiffTime;
fNegNpe.at(npmt) = fNegGain[npmt]*fGoodNegAdcPulseInt.at(npmt);
--
GitLab