Skip to content
Snippets Groups Projects
Commit 0a503a85 authored by John Derek Chapman's avatar John Derek Chapman Committed by Atlas Nightlybuild
Browse files

Merge branch 'fix_SCAutoCorr_21.3' into '21.3'

Fix SuperCell AutoCorrTotal

See merge request atlas/athena!32582

(cherry picked from commit dcb77655)

0a86c42d fixing MinBias key
parent 8942a71d
No related branches found
No related tags found
No related merge requests found
...@@ -32,6 +32,7 @@ def LArAutoCorrTotalSCToolDefault (name="LArAutoCorrTotalSCToolDefault", **kw): ...@@ -32,6 +32,7 @@ def LArAutoCorrTotalSCToolDefault (name="LArAutoCorrTotalSCToolDefault", **kw):
tool.keyNoise = "LArNoiseSC" tool.keyNoise = "LArNoiseSC"
tool.keyPedestal = "LArPedestalSC" tool.keyPedestal = "LArPedestalSC"
tool.keyfSampl = "LARfSamplSC" tool.keyfSampl = "LARfSamplSC"
tool.keyMinBias = "LArMinBiasSC"
from AthenaCommon.BeamFlags import jobproperties from AthenaCommon.BeamFlags import jobproperties
tool.deltaBunch = int(jobproperties.Beam.bunchSpacing()/( 25.*ns)+0.5) tool.deltaBunch = int(jobproperties.Beam.bunchSpacing()/( 25.*ns)+0.5)
......
...@@ -291,9 +291,11 @@ StatusCode LArAutoCorrTotalTool::getTerms() ...@@ -291,9 +291,11 @@ StatusCode LArAutoCorrTotalTool::getTerms()
// m_nsamples_AC_OFC = size of AC matrix & OFC vector (e.g 5 in Atlas) // m_nsamples_AC_OFC = size of AC matrix & OFC vector (e.g 5 in Atlas)
//::::::::::::::::::::::::::::::: //:::::::::::::::::::::::::::::::
float fSigma2=0.; float fSigma2=0.;
float MinBiasRMS=0.;
float SigmaNoise=0.;
float Adc2MeV=0.;
if(!m_NoPile) if(!m_NoPile)
{ {
float SigmaNoise;
if(m_isMC) if(m_isMC)
SigmaNoise = m_dd_noise->noise(id,igain); SigmaNoise = m_dd_noise->noise(id,igain);
else else
...@@ -306,11 +308,10 @@ StatusCode LArAutoCorrTotalTool::getTerms() ...@@ -306,11 +308,10 @@ StatusCode LArAutoCorrTotalTool::getTerms()
SigmaNoise = 0.;//(we will have the ERROR message below) SigmaNoise = 0.;//(we will have the ERROR message below)
} }
float fSampl = m_dd_fSampl->FSAMPL(id); float fSampl = m_dd_fSampl->FSAMPL(id);
float MinBiasRMS = m_dd_minbias->minBiasRMS(id); MinBiasRMS = m_dd_minbias->minBiasRMS(id);
if(fSampl!=0) MinBiasRMS/=fSampl; if(fSampl!=0) MinBiasRMS/=fSampl;
const std::vector<float> * const std::vector<float> *
polynom_adc2mev = &(m_adc2mevTool->ADC2MEV(id,igain)); polynom_adc2mev = &(m_adc2mevTool->ADC2MEV(id,igain));
float Adc2MeV=0.;
if (polynom_adc2mev->size()>0) { if (polynom_adc2mev->size()>0) {
Adc2MeV=(*polynom_adc2mev)[1]; Adc2MeV=(*polynom_adc2mev)[1];
} }
......
...@@ -455,9 +455,9 @@ LArOFCTool::computeOFC(int aORb, const HWIdentifier& CellID, ...@@ -455,9 +455,9 @@ LArOFCTool::computeOFC(int aORb, const HWIdentifier& CellID,
for(c=0;c<nsamples_AC_OFC;++c) for(c=0;c<nsamples_AC_OFC;++c)
std::cout<<Shape[c+iBeginOfNSamples]<<" "; std::cout<<Shape[c+iBeginOfNSamples]<<" ";
std::cout<<" <- "<<iBeginOfNSamples<<std::endl; std::cout<<" <- "<<iBeginOfNSamples<<std::endl;
for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PS[i]<<" "; for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PS[i]<<" ";
std::cout<<std::endl; std::cout<<std::endl;
for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PSD[i]<<" "; for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PSD[i]<<" ";
std::cout<<std::endl; std::cout<<std::endl;
std::cout<<" Q1="<<Q1<<" Q2="<<Q2<<" Q3="<<Q3 std::cout<<" Q1="<<Q1<<" Q2="<<Q2<<" Q3="<<Q3
<<" DELTA="<<DELTA<<std::endl; <<" DELTA="<<DELTA<<std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment