Skip to content
Snippets Groups Projects

Draft: calo E, ET and ADC 1D distributions back to being done by default

Closed Aniol Lobo Salvia requested to merge calo_updates into master
1 unresolved thread
2 files
+ 5
4
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -57,6 +57,7 @@ private:
mutable std::optional<Gaudi::Accumulators::Histogram<1>> m_histoDigitX;
mutable std::optional<Gaudi::Accumulators::Histogram<1>> m_histoDigitY;
mutable std::optional<Gaudi::Accumulators::Histogram<2>> m_histoDigitXY;
mutable std::optional<Gaudi::Accumulators::Histogram<1>> m_histoE;
mutable std::optional<Gaudi::Accumulators::Histogram<1>> m_histoET;
mutable std::optional<Gaudi::Accumulators::Histogram<1>> m_histoADC;
mutable Gaudi::Accumulators::StatCounter<> m_digits{this, "# digits over threshold"};
@@ -92,6 +93,7 @@ StatusCode CaloFutureDigitMonitor::initialize() {
m_histoDigitY.emplace( this, "DigitY", "Digit Y", axis1D{m_yBin, m_yMin, m_yMax} );
m_histoDigitXY.emplace( this, "DigitXY", "Digit position x vs y", axis1D{m_xBin, m_xMin, m_xMax},
axis1D{m_yBin, m_yMin, m_yMax} );
m_histoE.emplace( this, "E", "E", axis1D{m_etBin, m_etMin, m_etMax} );
m_histoET.emplace( this, "ET", "ET", axis1D{m_etBin, m_etMin, m_etMax} );
m_histoADC.emplace( this, "ADC", "ADC", axis1D{m_adcBin, m_adcMin, m_adcMax} );
} );
@@ -142,6 +144,9 @@ void CaloFutureDigitMonitor::operator()( const Input& digits, const LHCb::ODIN&
const auto e = digit.energy();
const auto adc = digit.adc();
const auto et = e * calo.cellSine( id );
++m_histoE.value()[e];
++m_histoET.value()[et];
++m_histoADC.value()[adc];
if ( e < m_eFilter ) continue;
if ( et < m_etFilter ) continue;
if ( adc < m_adcFilter ) continue;
@@ -179,9 +184,6 @@ void CaloFutureDigitMonitor::operator()( const Input& digits, const LHCb::ODIN&
if ( msgLevel( MSG::VERBOSE ) ) {
verbose() << " et " << et << " cell " << unit_et << endmsg;
verbose() << " adc " << adc << " cell " << unit_adc << endmsg;
} else {
++m_histoET.value()[et];
++m_histoADC.value()[adc];
}
}
}
Loading