From d0b07c900b4fa1d3329eea45e88ecd2e212a5099 Mon Sep 17 00:00:00 2001
From: Carlos Yero <nuclear@deuteron.local>
Date: Mon, 21 Dec 2020 11:17:49 -0500
Subject: [PATCH] minor bug fix on the beam current BCM saturation correction
on the 1st standard scaler event
---
src/THcScalerEvtHandler.cxx | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/THcScalerEvtHandler.cxx b/src/THcScalerEvtHandler.cxx
index d4cf849..358cffb 100644
--- a/src/THcScalerEvtHandler.cxx
+++ b/src/THcScalerEvtHandler.cxx
@@ -424,7 +424,8 @@ Int_t THcScalerEvtHandler::AnalyzeBuffer(UInt_t* rdata, Bool_t onlysync)
dvars[ivar]=0.;
if (bcm_ind != -1) {
dvars[ivar]=((scalers[idx]->GetData(ichan))/fDeltaTime-fBCM_Offset[bcm_ind])/fBCM_Gain[bcm_ind];
- dvars[ivar]=dvars[ivar]+fBCM_SatOffset[bcm_ind]*TMath::Max(dvars[ivar]-fBCM_SatOffset[i],0.0);
+ dvars[ivar]=dvars[ivar]+fBCM_SatQuadratic[bcm_ind]*TMath::Power(TMath::Max(dvars[ivar]-fBCM_SatOffset[bcm_ind],0.0),2.0);
+
}
if (bcm_ind == fbcm_Current_Threshold_Index) scal_current= dvars[ivar];
}
--
GitLab