LArDigitalTriggMonAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
 
// ********************************************************************
//
// NAME:     LArDigitalTriggMonAlg.cxx
// PACKAGE:  LArMonitoring
//
// AUTHOR:   Pavol Strizenec (pavol@mail.cern.ch)
//           Ljiljana Morvaj (ljiljana.morvaj@cern.ch)
//           Yesenia Hernadez (yesenia@cern.ch)
//           Based on LAtDigitMon tool by L. Hellary and LArOddCellsMonTool.cxx  by Benjamin Trocme
//
// Monitor a few things in the LArDigit...
//
//  1) Check that the highest value of the LArDigit is contained in an interval. 
//         If it is not the case increment 3 histograms for each subdetector:
//        a) Out of range histograms
//        b) The average histograms: give the average value of the highest digit sample 
//        c) Channel VS FEB histograms: gives wich slot on wich FEB has his highest digit sample ou of the range
//  2) Check if a digits samples are in saturation. If it's the case increment the saturation histograms. 
//
// Available cuts in the jo file:
//
//   a) SampleRangeLow-SampleRangeUp: range to check the digit sample.
//   b) ADCcut : To select Digits Samples with signal.
//   c) ADCsature: lowest value to check if a Digit sample is in saturation.
// ********************************************************************
 
#include "LArDigitalTriggMonAlg.h"
 
//Histograms
//LAr infos:
#include "Identifier/HWIdentifier.h"
#include "LArIdentifier/LArOnlineID.h"
#include "LArRawEvent/LArDigit.h"
#include "CaloGeoHelpers/CaloSampling.h"
#include "LArRawEvent/LArSCDigit.h"
#include "LArRawEvent/LArLATOMEHeaderContainer.h"
#include "CaloIdentifier/CaloCell_SuperCell_ID.h"
#include "LArTrigStreamMatching.h"
 
#include "LArCOOLConditions/LArPedestalSC.h"
 
//STL:
#include <sstream>
#include <iomanip>
#include <cmath>
#include <vector>
#include <algorithm>
#include <set>
 
template<typename C>
bool isEmptyCont(C& c) {
  return (c.key().empty() || !c.isValid() || c->size()==0); 
}
 
 
/*---------------------------------------------------------*/
LArDigitalTriggMonAlg::~LArDigitalTriggMonAlg()
{
}
 
/*---------------------------------------------------------*/
StatusCode LArDigitalTriggMonAlg::initialize()
{
  ATH_MSG_INFO( "Initialize LArDigitalTriggMonAlg" );
 
  ATH_CHECK(detStore()->retrieve( m_LArOnlineIDHelper, "LArOnline_SuperCellID" ));
 
  ATH_CHECK(detStore()->retrieve( m_SCID_helper, "CaloCell_SuperCell_ID" ).isSuccess());
 
  ATH_MSG_INFO("Building tool map");
  m_toolmapLayerNames_digi = Monitored::buildToolMap<int>( m_tools, "LArDigitalTriggerMon_digi", m_layerNames);
  m_toolmapLayerNames_sc = Monitored::buildToolMap<int>( m_tools, "LArDigitalTriggerMon_sc", m_layerNames);
 
  ATH_MSG_INFO("Done building tool map");
 
  ATH_MSG_INFO("Input containers=" << m_rawSCContainerKey << " / " << m_rawSCEtRecoContainerKey);
 
  
  ATH_CHECK(m_bcContKey.initialize());
  ATH_CHECK(m_bcMask.buildBitMask(m_problemsToMask,msg()));
 
  ATH_CHECK(m_digitContainerKey.initialize());
  ATH_CHECK(m_keyPedestalSC.initialize());
  ATH_CHECK(m_caloSuperCellMgrKey.initialize());
  ATH_CHECK(m_rawSCContainerKey.initialize(SG::AllowEmpty));
  ATH_CHECK(m_rawSCEtRecoContainerKey.initialize(SG::AllowEmpty));
  ATH_CHECK(m_cablingKey.initialize());
  ATH_CHECK(m_actualMuKey.initialize());
  ATH_CHECK(m_LATOMEHeaderContainerKey.initialize());
 
  // Property check:
  constexpr unsigned expSize=MAXLYRNS*2+1;
  if (m_layerNames.size() != expSize) {
    ATH_MSG_ERROR("Unexpected size of LayerNames property. Expect "<< expSize << " entries, found "<< m_layerNames.size() << " entries");
    return StatusCode::FAILURE;
  }
 
  if (m_isADCBaseline) {
    ATH_MSG_INFO("IsADCBas set to true");
  }
  return AthMonitorAlgorithm::initialize();
}
 
 
struct Digi_MonValues {
  float digi_eta;
  float digi_phi;
  int digi_ieta;
  int digi_iphi;
  int digi_sampos;
  int digi_adc;
  int digi_latomesourceidbin;
  float digi_pedestal;
  int digi_maxpos;
  int digi_partition;
  float digi_diff_adc_ped_norm;
  float digi_diff_adc_ped;
  float digi_diff_adc0_ped;
  float digi_adc_rms;
  int digi_bcid;
  unsigned int digi_lb;
  bool digi_passDigiNom;
  bool digi_badNotMasked;
 
};
 
 
struct SC_MonValues {
  float sc_eta;
  float sc_phi;
  int sc_ieta;
  int sc_iphi;
  int sc_latomesourceidbin;
  float sc_et_ofl;
  int sc_et_diff;
  float sc_et_onl;
  float sc_et_onl_muscaled;
  float sc_time;
  int sc_bcid;
  unsigned int sc_lb;
  bool sc_zeroET; 
  bool sc_passSCNomInvalid;
  bool sc_passSCNom0_0p325;
  bool sc_passSCNom0p325_1;
  bool sc_passSCNom1;
  bool sc_passSCNom10;
  bool sc_passSCNom10tauGt3;
  bool sc_saturNotMasked;
  bool sc_OFCbOFNotMasked;
  bool sc_notMaskedEoflNe0;
  bool sc_notMaskedEoflGt1;
};
 
 
StatusCode LArDigitalTriggMonAlg::fillHistograms(const EventContext& ctx) const
{
 
  ATH_MSG_DEBUG("in fillHists()" );
 
  // General Monitored variables
  auto lumi_block = Monitored::Scalar<unsigned int>("lumi_block", 0);
  auto time_stamp = Monitored::Scalar<unsigned int>("time_stamp", 0);
  auto BCID = Monitored::Scalar<int>("BCID",0);
  auto Pedestal = Monitored::Scalar<float>("Pedestal",0.0);
  auto PedestalRMS = Monitored::Scalar<float>("PedestalRMS",0.0);
 
  // From digi loop
  auto Digi_Nsamples = Monitored::Scalar<int>("Digi_Nsamples",-1);  // MNsamples
  auto Digi_SCChannel = Monitored::Scalar<int>("Digi_SCChannel",-1); // MSCChannel
  auto Digi_latomeSourceId = Monitored::Scalar<int>("Digi_latomeSourceId",-1); // MlatomeSourceId
  auto Digi_latomeSourceIdBIN = Monitored::Scalar<int>("Digi_latomeSourceIdBIN",1); // MlatomeSourceIdBIN
  auto Digi_phi = Monitored::Scalar<float>("Digi_phi",0.0); // MSCphi
  auto Digi_eta = Monitored::Scalar<float>("Digi_eta",0.0); // MSCeta
  auto Digi_iphi = Monitored::Scalar<int>("Digi_iphi",0.0);
  auto Digi_ieta = Monitored::Scalar<int>("Digi_ieta",0.0);
  auto Digi_maxpos = Monitored::Scalar<int>("Digi_maxpos",-1); // Mmaxpos
  auto Digi_partition = Monitored::Scalar<int>("Digi_partition",-1); // Mpartition
  auto Digi_sampos = Monitored::Scalar<int>("Digi_sampos",-1); // Msampos
  auto Digi_ADC = Monitored::Scalar<int>("Digi_ADC",-1); // MADC
  auto Digi_Diff_ADC_Ped = Monitored::Scalar<float>("Digi_Diff_ADC_Ped", -999); // Diff_ADC_Pedestal
  auto Digi_Diff_ADC0_Ped = Monitored::Scalar<float>("Digi_Diff_ADC0_Ped", -999); // Pedestal diff
  auto Digi_Diff_ADC_Ped_Norm = Monitored::Scalar<float>("Digi_Diff_ADC_Ped_Norm",-999); // Diff_ADC_Pedestal_Norm
  auto Digi_ADC_RMS = Monitored::Scalar<float>("Digi_ADC_RMS",-1); // Digi_ADC_RMS
 
  // cuts
  auto notBadQual = Monitored::Scalar<bool>("notBadQual",false);
  auto ADCped10RMS = Monitored::Scalar<bool>("ADCped10RMS",false);
  auto passDigiNom = Monitored::Scalar<bool>("passDigiNom",false);
  auto badNotMasked = Monitored::Scalar<bool>("badNotMasked",false);
 
  // cuts which are used in both loops
  auto notMasked = Monitored::Scalar<bool>("notMasked",false);
 
  // From SC loop
  auto SC_SCChannel = Monitored::Scalar<int>("SC_SCChannel",-1); // MSCChannel
  auto SC_latomeSourceId = Monitored::Scalar<int>("SC_latomeSourceId",-1); // MlatomeSourceId
  auto SC_partition = Monitored::Scalar<int>("SC_partition",-1); // Mpartition
  auto SC_phi = Monitored::Scalar<float>("SC_phi",0.0); // MSCphi
  auto SC_eta = Monitored::Scalar<float>("SC_eta",0.0); // MSCeta
  auto SC_iphi = Monitored::Scalar<int>("SC_iphi",0.0);
  auto SC_ieta = Monitored::Scalar<int>("SC_ieta",0.0);
  auto SC_energy_onl = Monitored::Scalar<int>("SC_energy_onl",0.0); // Menergy_onl
  auto SC_ET_onl = Monitored::Scalar<float>("SC_ET_onl",0.0); // Menergy_onl
  auto SC_ET_onl_muscaled = Monitored::Scalar<float>("SC_ET_onl_muscaled",0.0); // Menergy_onl
  auto SC_energy_ofl = Monitored::Scalar<int>("SC_energy_ofl",0.0); // Menergy_ofl
  auto SC_ET_ofl = Monitored::Scalar<float>("SC_ET_ofl",0.0); // Menergy_onl
  auto SC_ET_diff = Monitored::Scalar<int>("SC_ET_diff",0.0); // MSCEt_diff
  auto SC_time = Monitored::Scalar<float>("SC_time",0.0); // MSCtime
  auto SC_latomeSourceIdBIN = Monitored::Scalar<int>("SC_latomeSourceIdBIN",1); // MlatomeSourceIdBIN
  auto SC_AvEnergyOverMu = Monitored::Scalar<float>("SC_AvEnergyOverMu",0); // LMAvEnergyOverMu
  // cuts
  auto passTauSel = Monitored::Scalar<bool>("passTauSel",false);
  auto nonZeroET = Monitored::Scalar<bool>("nonZeroET",false); // eTgt0GeV
  auto zeroET = Monitored::Scalar<bool>("zeroET",false); // eTgt0GeV
  auto nonZeroETofl = Monitored::Scalar<bool>("nonZeroETofl",false); // eTgt0GeV
  auto onlofflEmismatch = Monitored::Scalar<bool>("onlofflEmismatch",false);
  auto notSatur = Monitored::Scalar<bool>("notSatur",false);
  auto notOFCbOF = Monitored::Scalar<bool>("notOFCbOF",false);
  auto tauGt3 = Monitored::Scalar<bool>("tauGt3",false);
  auto nonZeroEtau = Monitored::Scalar<bool>("nonZeroEtau",false);
  auto eTgt1GeV = Monitored::Scalar<bool>("eTgt1GeV",false);
  auto eTlt1GeV = Monitored::Scalar<bool>("eTlt1GeV",false);
  auto eTgt0p325GeV = Monitored::Scalar<bool>("eTgt0p325GeV",false);
  auto eTgt0lt0p325GeV = Monitored::Scalar<bool>("eTgt0lt0p325GeV",false);
  auto eTgt10GeV = Monitored::Scalar<bool>("eTgt10GeV",false);
  auto eTlt10GeV = Monitored::Scalar<bool>("eTlt10GeV",false);
  auto eToflGt1GeV = Monitored::Scalar<bool>("eToflGt1GeV",false);
 
 
  //auto passSCNom = Monitored::Scalar<bool>("passSCNom",false);  // pass tau, not satur, not OFCb OF, not masked  nonZeroET < 10 GeV
  auto passSCNomInvalid = Monitored::Scalar<bool>("passSCNomInvalid",false);  // pass tau, not satur, not OFCb OF, not masked and raw E = -99999
  auto passSCNom0_0p325 = Monitored::Scalar<bool>("passSCNom0_0p325",false);  // pass tau, not satur, not OFCb OF, not masked  nonZeroET < 0.2 GeV
  auto passSCNom0p325_1 = Monitored::Scalar<bool>("passSCNom0p325_1",false);  // pass tau, not satur, not OFCb OF, not masked  0.2 < ET < 1 GeV
  auto passSCNom1 = Monitored::Scalar<bool>("passSCNom1",false);  // pass tau, not satur, not OFCb OF, not masked eTgt1GeV
  auto passSCNom10 = Monitored::Scalar<bool>("passSCNom10",false);  // pass tau, not satur, not OFCb OF, not masked eTgt10GeV
  auto passSCNom10tauGt3 = Monitored::Scalar<bool>("passSCNom10tauGt3",false);  // pass tau, not satur, not OFCb OF, not masked eTgt10GeV  tauGt3
  auto saturNotMasked = Monitored::Scalar<bool>("saturNotMasked",false);  // notSatur is false, notMasked is false
  auto OFCbOFNotMasked = Monitored::Scalar<bool>("OFCbOFNotMasked",false);  // notOFCbOF is false, notMasked is false
  auto notMaskedEoflNe0 = Monitored::Scalar<bool>("notMaskedEoflNe0",false);  // not masked OSUM, not satur, not OFCb OF, ET ofl != 0
  auto notMaskedEoflGt1 = Monitored::Scalar<bool>("notMaskedEoflGt1",false);  // not masked OSUM, not satur, not OFCb OF, ET ofl > 1
 
  // From LATOME header loop
  auto thisEvent=this->GetEventInfo(ctx);
 
  const std::vector<unsigned> streamsThisEvent=LArMon::trigStreamMatching(m_streams,thisEvent->streamTags());
 
  SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey,ctx};
  const LArOnOffIdMapping* cabling=*cablingHdl;
 
  SG::ReadHandle<LArDigitContainer> hLArDigitContainer;
  if (!m_digitContainerKey.empty()) {
    hLArDigitContainer= SG::ReadHandle<LArDigitContainer>{m_digitContainerKey,ctx}; //"SC"  
    if (!hLArDigitContainer.isValid()) {
      ATH_MSG_WARNING("The requested digit container key could not be retrieved. Was there a problem retrieving information from the run logger?");
    }
  }
  else {
    ATH_MSG_DEBUG("hLArDigitContainer.size() " << hLArDigitContainer->size());
  }
 
  
  SG::ReadHandle<LArRawSCContainer> hSCetContainer;
  if (!m_rawSCContainerKey.empty()) {
    hSCetContainer = SG::ReadHandle<LArRawSCContainer>{m_rawSCContainerKey, ctx};  //"SC_ET"
    if (!hSCetContainer.isValid()) {
      ATH_MSG_WARNING("The requested SC ET container key could not be retrieved. Was there a problem retrieving information from the run logger?");
    }
  } 
  else {
    ATH_MSG_DEBUG("hSCetContainer.size() " << hSCetContainer->size());
  }
 
  SG::ReadHandle<LArRawSCContainer> hSCetRecoContainer;
  if (!m_rawSCEtRecoContainerKey.empty()) {
    hSCetRecoContainer = SG::ReadHandle<LArRawSCContainer>{m_rawSCEtRecoContainerKey, ctx};  //"SC_ET_RECO"
    if (!hSCetRecoContainer.isValid()) {
      ATH_MSG_WARNING("The requested SC ET reco container key could not be retrieved. Was there a problem retrieving information from the run logger?");
    }
  } 
  else {
    ATH_MSG_DEBUG("hSCetRecoContainer.size() " << hSCetRecoContainer->size());
  }
 
  SG::ReadHandle<LArLATOMEHeaderContainer> hLArLATOMEHeaderContainer{m_LATOMEHeaderContainerKey,ctx}; //"SC_LATOME_HEADER"
  if (!hLArLATOMEHeaderContainer.isValid()) {
    ATH_MSG_WARNING("The requested LATOME header container key could not be retrieved. Was there a problem retrieving information from the run logger?");
  }else{
    ATH_MSG_DEBUG("hLArLATOMEHeaderContainer.size() " << hLArLATOMEHeaderContainer->size());
  }
 
  if (isEmptyCont(hLArDigitContainer) && isEmptyCont(hSCetContainer) && isEmptyCont(hSCetRecoContainer) && isEmptyCont(hLArLATOMEHeaderContainer)) {
    //Make this only warning, come CI tests use the runs without DT info
    ATH_MSG_WARNING("All of the requested containers are empty. Was there a problem retrieving information from the run logger?");
    return StatusCode::SUCCESS;
  }
 
    BCID = thisEvent->bcid();  // - 88)%((36+7)*4 + 36 + 31);
    lumi_block = thisEvent->lumiBlock();
    time_stamp = thisEvent->timeStamp();
 
    float mu = lbInteractionsPerCrossing(ctx);
    float event_mu = lbLuminosityPerBCID(ctx);
 
    ATH_MSG_DEBUG("mu (LB): " << mu);
    ATH_MSG_DEBUG("mu (BCID): " << event_mu);
    ATH_MSG_DEBUG("Event number: " << thisEvent->eventNumber());
    ATH_MSG_DEBUG("LB number: " << thisEvent->lumiBlock());
    ATH_MSG_DEBUG("BCID: " << thisEvent->bcid());
    SG::ReadCondHandle<ILArPedestal> pedestalHdl{m_keyPedestalSC, ctx};
    const ILArPedestal* pedestals = *pedestalHdl;
 
    SG::ReadCondHandle<CaloSuperCellDetDescrManager> caloSuperCellMgrHandle{m_caloSuperCellMgrKey, ctx};
    const CaloSuperCellDetDescrManager* ddman = *caloSuperCellMgrHandle;
 
    // retrieve BadChannel info:
    const LArBadChannelCont* bcCont = nullptr;
    SG::ReadCondHandle<LArBadChannelCont> bcContHdl{m_bcContKey, ctx};
    bcCont = (*bcContHdl);
 
    if ((hLArDigitContainer.isValid())) {
      std::vector<std::vector<Digi_MonValues>> digiMonValueVec(m_layerNames.size());
      for (auto& innerVec : digiMonValueVec) {
        innerVec.reserve(1600);  // (m_layerNcells[ilayer]) * nsamples;
      }
 
      // Loop over digits
      for (const LArDigit* pLArDigit : *hLArDigitContainer) {
        HWIdentifier id = pLArDigit->hardwareID();  // gives online ID
        // skip disconnected channels:
        if (!cabling->isOnlineConnected(id))
          continue;
 
        const unsigned trueNSamples = pLArDigit->nsamples();
        if (trueNSamples == 0)
          continue;
        Digi_Nsamples = trueNSamples;  // Fill the monitored variable
        const int cgain = pLArDigit->gain();
 
        const Identifier offlineID = cabling->cnvToIdentifier(id);
        const CaloDetDescrElement* caloDetElement = ddman->get_element(offlineID);
        if (caloDetElement == 0) {
          ATH_MSG_ERROR("Cannot retrieve caloDetElement");
          continue;
        }
        Digi_eta = caloDetElement->eta_raw();
        Digi_phi = caloDetElement->phi_raw();
 
        Digi_ieta = m_SCID_helper->eta(offlineID);
        Digi_iphi = m_SCID_helper->phi(offlineID);
 
        const int calosample = caloDetElement->getSampling();
 
        const unsigned iLyrNS = m_caloSamplingToLyrNS[calosample];
        const int side = m_LArOnlineIDHelper->pos_neg(id);
        const unsigned iLyr = iLyrNS * 2 + side;
        auto& lvaluemap_digi = digiMonValueVec[iLyr];
        auto& lvaluemap_digi_ALL = digiMonValueVec.back();
 
        // Determine to which partition this channel belongs to
        const int ThisPartition = whatPartition(id, side);
        Digi_partition = ThisPartition;  // Fill the monitored variable
 
        fill(m_scMonGroupName, Digi_Nsamples);
 
        // Check if this is a maskedOSUM SC
        notMasked = true;
        if (m_bcMask.cellShouldBeMasked(bcCont, id)) {
          notMasked = false;
        }
 
        if (pedestals) {
          Pedestal = pedestals->pedestal(id, cgain);
          PedestalRMS = pedestals->pedestalRMS(id, cgain);
        } else
          ATH_MSG_INFO("Pedestal values not received");
 
        const LArSCDigit* scdigi = dynamic_cast<const LArSCDigit*>(pLArDigit);
        if (!scdigi) {
          ATH_MSG_DEBUG(" CAN'T CAST ");
        } else {
          Digi_SCChannel = scdigi->Channel();
          Digi_latomeSourceId = scdigi->SourceId();
          Digi_latomeSourceIdBIN = getXbinFromSourceID(Digi_latomeSourceId);
        }
        // Retrieve samples
        const std::vector<short>* digito = &pLArDigit->samples();
 
        // retrieve the max sample digit ie digitot.back().
        std::vector<short>::const_iterator maxSam = std::max_element(digito->begin(), digito->end());
        int thismaxPos = std::distance(digito->begin(), maxSam);
        Digi_maxpos = thismaxPos + 1;  // count samples [1,5]
        float ADC_max = pLArDigit->samples().at(Digi_maxpos - 1);
 
        float ADC_0 = pLArDigit->samples().at(0);
        if (m_isADCBaseline) {  // SC_ADC_BAS, have to divide by 8
          ADC_0 = ADC_0 / 8;
        }
 
        float samp_sum = std::accumulate(digito->begin(), digito->end(), 0.0);
        float samp_mean = samp_sum / static_cast<float>(trueNSamples);
        float sq_sum = std::inner_product(digito->begin(), digito->end(), digito->begin(), 0.0);
        float rms_arg = sq_sum / static_cast<float>(trueNSamples) - samp_mean * samp_mean;
 
        if (rms_arg < 0)
          Digi_ADC_RMS = -1;
        else
          Digi_ADC_RMS = std::sqrt(rms_arg);
        
 
        // Start Loop over samples
        Digi_Diff_ADC0_Ped = ADC_0 - Pedestal;
        for (unsigned i = 0; i < trueNSamples; ++i) {
          badNotMasked = false;
          notBadQual = false;
          ADCped10RMS = false;
          passDigiNom = false;
          Digi_sampos = i + 1;
          Digi_ADC = pLArDigit->samples().at(i);
          if (m_isADCBaseline) {  // SC_ADC_BAS, have to divide by 8
            Digi_ADC = Digi_ADC / 8;
          }
 
          Digi_Diff_ADC_Ped = Digi_ADC - Pedestal;
          if (ADC_max != Pedestal) {
            Digi_Diff_ADC_Ped_Norm = (Digi_ADC - Pedestal) / std::abs(ADC_max - Pedestal);
          }
 
          // Some selections
          if (Digi_ADC != -1) {
            notBadQual = true;
          } else {
            if (notMasked) {
              badNotMasked = true;
            }
          }
          if (ADC_max - Pedestal > 10 * PedestalRMS) {
            ADCped10RMS = true;
          }
          if (notMasked && notBadQual && ADCped10RMS) {
            passDigiNom = true;
          }
 
          lvaluemap_digi.emplace_back(Digi_eta, Digi_phi, Digi_ieta, Digi_iphi, Digi_sampos, Digi_ADC, Digi_latomeSourceIdBIN, Pedestal, Digi_maxpos,
                                      Digi_partition, Digi_Diff_ADC_Ped_Norm, Digi_Diff_ADC_Ped, Digi_Diff_ADC0_Ped, Digi_ADC_RMS, BCID, lumi_block, passDigiNom,
                                      badNotMasked);
          lvaluemap_digi_ALL.emplace_back(Digi_eta, Digi_phi, Digi_ieta, Digi_iphi, Digi_sampos, Digi_ADC, Digi_latomeSourceIdBIN, Pedestal, Digi_maxpos,
                                          Digi_partition, Digi_Diff_ADC_Ped_Norm, Digi_Diff_ADC_Ped, Digi_Diff_ADC0_Ped, Digi_ADC_RMS, BCID, lumi_block, passDigiNom,
                                          badNotMasked);
 
        }  // End loop over samples
 
      }  // End of loop on LArDigit
 
      // fill, for every layer/threshold
      for (size_t ilayer = 0; ilayer < digiMonValueVec.size(); ++ilayer) {
        const auto& tool = digiMonValueVec[ilayer];
        auto digi_part_eta = Monitored::Collection("Digi_part_eta", tool, [](const auto& v) { return v.digi_eta; });
        auto digi_part_phi = Monitored::Collection("Digi_part_phi", tool, [](const auto& v) { return v.digi_phi; });
        auto digi_part_ieta = Monitored::Collection("Digi_part_ieta", tool, [](const auto& v) { return v.digi_ieta; });
        auto digi_part_iphi = Monitored::Collection("Digi_part_iphi", tool, [](const auto& v) { return v.digi_iphi; });
        auto digi_part_sampos = Monitored::Collection("Digi_part_sampos", tool, [](const auto& v) { return v.digi_sampos; });
        auto digi_part_adc = Monitored::Collection("Digi_part_adc", tool, [](const auto& v) { return v.digi_adc; });
        auto digi_part_latomesourceidbin = Monitored::Collection("Digi_part_latomesourceidbin", tool, [](const auto& v) { return v.digi_latomesourceidbin; });
        auto digi_part_pedestal = Monitored::Collection("Digi_part_pedestal", tool, [](const auto& v) { return v.digi_pedestal; });
        auto digi_part_maxpos = Monitored::Collection("Digi_part_maxpos", tool, [](const auto& v) { return v.digi_maxpos; });
        auto digi_part_partition = Monitored::Collection("Digi_part_partition", tool, [](const auto& v) { return v.digi_partition; });
        auto digi_part_diff_adc_ped_norm = Monitored::Collection("Digi_part_diff_adc_ped_norm", tool, [](const auto& v) { return v.digi_diff_adc_ped_norm; });
        auto digi_part_diff_adc_ped = Monitored::Collection("Digi_part_diff_adc_ped", tool, [](const auto& v) { return v.digi_diff_adc_ped; });
        auto digi_part_diff_adc0_ped = Monitored::Collection("Digi_part_diff_adc0_ped", tool, [](const auto& v) { return v.digi_diff_adc0_ped; });
        auto digi_part_bcid = Monitored::Collection("Digi_part_BCID", tool, [](const auto& v) { return v.digi_bcid; });
        auto digi_part_lb = Monitored::Collection("Digi_part_LB", tool, [](const auto& v) { return v.digi_lb; });
        auto digi_part_passDigiNom = Monitored::Collection("Digi_part_passDigiNom", tool, [](const auto& v) { return v.digi_passDigiNom; });
        auto digi_part_badNotMasked = Monitored::Collection("Digi_part_badNotMasked", tool, [](const auto& v) { return v.digi_badNotMasked; });
        auto digi_part_adc_rms = Monitored::Collection("Digi_part_adc_rms", tool, [](const auto& v) { return v.digi_adc_rms; });
 
        fill(m_tools[m_toolmapLayerNames_digi.at(m_layerNames[ilayer])], digi_part_eta, digi_part_phi, digi_part_ieta, digi_part_iphi, digi_part_sampos,
             digi_part_adc, digi_part_latomesourceidbin, digi_part_pedestal, digi_part_maxpos, digi_part_diff_adc_ped_norm, digi_part_diff_adc_ped,
             digi_part_diff_adc0_ped, digi_part_adc_rms, digi_part_bcid, digi_part_lb, digi_part_passDigiNom, digi_part_badNotMasked);
      }
 
    }  // End if(LArDigitContainer is valid)
 
 
    if (!isEmptyCont(hSCetContainer) && !isEmptyCont(hSCetRecoContainer)) {
      LArRawSCContainer::const_iterator itSC = hSCetContainer->begin();
      LArRawSCContainer::const_iterator itSC_e= hSCetContainer->end();
      LArRawSCContainer::const_iterator itSCReco = hSCetRecoContainer->begin();
      const LArRawSC* rawSC = 0;
      const LArRawSC* rawSCReco = 0;
 
      std::vector<std::vector<SC_MonValues>> scMonValueVec(m_layerNames.size());
      for (auto& innerVec : scMonValueVec) {
        innerVec.reserve(1600);  // (m_layerNcells[ilayer]) * nsamples;
      }
 
      // Loop over SCs
      for (; itSC != itSC_e; ++itSC, ++itSCReco) {
        rawSC = *itSC;
        if (itSCReco < hSCetRecoContainer->end()) {
          rawSCReco = *itSCReco;
        } else {
        //temporarily removed
        //ATH_MSG_WARNING("Looping SC ET container, but we have reached the end of the SC ET Reco iterator. Check the sizes of these containers. Is SC ET Reco size zero? Is there a problem with the digit container name sent by the run logger?");
          rawSCReco = 0;
        }
        SC_SCChannel = rawSC->chan();
        HWIdentifier id = rawSC->hardwareID();  // gives online ID
        // skip disconnected channels:
        if (!cabling->isOnlineConnected(id))
          continue;
 
        const Identifier offlineID = cabling->cnvToIdentifier(id);  // converts online to offline ID
        // Get Physical Coordinates
        const CaloDetDescrElement* caloDetElement = ddman->get_element(offlineID);
        if (caloDetElement == 0) {
          ATH_MSG_ERROR("Cannot retrieve (eta,phi) coordinates for raw channels");
          ATH_MSG_ERROR("  ==============> " << std::hex << ";  offlineID = " << offlineID << "online ID =" << m_LArOnlineIDHelper->channel_name(id)
                                             << "; rawSC->SourceId() = " << rawSC->SourceId());
          continue;
        }
        SC_eta = caloDetElement->eta_raw();
        SC_phi = caloDetElement->phi_raw();
 
        SC_ieta = m_SCID_helper->eta(offlineID);
        SC_iphi = m_SCID_helper->phi(offlineID);
        int calosample = caloDetElement->getSampling();
 
        const unsigned iLyrNS = m_caloSamplingToLyrNS[calosample];
        const unsigned side = m_LArOnlineIDHelper->pos_neg(id);
        const unsigned iLyr = iLyrNS * 2 + side;
 
        auto& lvaluemap_sc = scMonValueVec[iLyr];
        auto& lvaluemap_sc_ALL = scMonValueVec.back();
 
        SC_latomeSourceIdBIN = getXbinFromSourceID(rawSC->SourceId());
 
        // initialise cuts
        notMasked = false;
        passTauSel = false;
        nonZeroET = false;
        zeroET = false;
        notSatur = false;
        nonZeroEtau = false;
        eTgt0p325GeV = false;
        eTgt0lt0p325GeV = false;
        eTgt1GeV = false;
        eTlt1GeV = false;
        eTgt10GeV = false;
        eTlt10GeV = false;
        notOFCbOF = false;
        tauGt3 = false;
        onlofflEmismatch = false;
        passSCNom0_0p325 = false;
        passSCNom0p325_1 = false;
        passSCNom1 = false;
        passSCNom10 = false;
        passSCNom10tauGt3 = false;
        saturNotMasked = false;
        OFCbOFNotMasked = false;
        notMaskedEoflNe0 = false;
        notMaskedEoflGt1 = false;
        nonZeroETofl = false;
        eToflGt1GeV = false;
        passSCNomInvalid = false;
        // Check if this is a maskedOSUM SC
        if (!m_bcMask.cellShouldBeMasked(bcCont, id)) {
          notMasked = true;
        }
        // popopopo
        if (rawSCReco != 0 && rawSCReco->passTauSelection().size() > 0) {  // only compare Et if tau selection is passed
          if (rawSCReco->passTauSelection().at(0) == true)
            passTauSel = true;
        }
        unsigned int bcid_ind = 0;
        if (rawSC->energies().size() > 0) {
          for (auto& SCe : rawSC->bcids()) {
            if (SCe == BCID)
              break;
            bcid_ind++;
          }
          if (bcid_ind >= rawSC->bcids().size()) {
            ATH_MSG_WARNING("BCID not found in SC bcids list!! " << BCID << " BCIDs size: " << rawSC->bcids().size() << ", proposed index: " << bcid_ind);
          } else if (rawSC->bcids().at(bcid_ind) != BCID) {
            ATH_MSG_WARNING("BCID not found in SC bcids list!! " << BCID << " " << rawSC->bcids().at(bcid_ind));
          }
 
          if (rawSC->energies().size() > bcid_ind) {
            SC_energy_onl = rawSC->energies().at(bcid_ind);
          } else {
            ATH_MSG_WARNING("rawSC energies vector is too small for the requested BCID index " << bcid_ind << " (size is " << rawSC->energies().size()
                                                                                               << ", bcid vec size is " << rawSC->bcids().size() << ")");
            SC_energy_onl = 0;
          }
        } else {
          ATH_MSG_WARNING("rawSC energies vector is empty!");
          SC_energy_onl = 0;
        }
 
        if (rawSCReco != 0) {
          if (rawSCReco->energies().size() > 0) {
            SC_energy_ofl = rawSCReco->energies().at(0);  // algorithm already selects the correct energy
          } else {
            ATH_MSG_WARNING("rawSCReco energies vector is empty!");
            SC_energy_ofl = 0;
          }
        }
        SC_ET_diff = SC_energy_onl - SC_energy_ofl;
        SC_ET_onl = (SC_energy_onl * 12.5) / 1000;  // Converted to GeV
        SC_ET_ofl = (SC_energy_ofl * 12.5) / 1000;  // Converted to GeV
        SC_ET_onl_muscaled = event_mu > 0. ? SC_ET_onl / event_mu : SC_ET_onl;
        int Etau = 0;
        if (rawSCReco != 0) {
          if (rawSCReco->tauEnergies().size() > 0) {
            Etau = rawSCReco->tauEnergies().at(0);
          }
        }
        SC_time = (SC_energy_ofl != 0) ? (float)Etau / (float)SC_energy_ofl : Etau;
 
        ATH_MSG_DEBUG("Energy onl - Energy ofl: " << SC_energy_onl << ",  " << SC_energy_ofl << std::endl);
        if (SC_ET_onl != 0) {
          nonZeroET = true;
        } else {
          zeroET = true;
        }
        if (SC_ET_onl > 0.325) {
          eTgt0p325GeV = true;
        } else if (SC_ET_onl < 0.325 &&  SC_ET_onl > 0.) {
          eTgt0lt0p325GeV = true;
        }
        if (SC_ET_onl > 1) {
          eTgt1GeV = true;
        }
        if (SC_ET_onl < 1) {
          eTlt1GeV = true;
        }
        if (SC_ET_onl > 10) {
          eTgt10GeV = true;
        }
        if (SC_ET_onl < 10) {
          eTlt10GeV = true;
        }
        if (SC_ET_ofl != 0) {
          nonZeroETofl = true;
        }
        if (SC_ET_ofl > 1) {
          eToflGt1GeV = true;
        }
        if (rawSC->satur().size() > bcid_ind) {
          if (rawSC->satur().at(bcid_ind)) {
            if (notMasked) {
              saturNotMasked = true;
            }
          } else {
            notSatur = true;
          }
        }
        if (Etau != 0) {
          nonZeroEtau = true;
        }
        if (rawSCReco != 0 && rawSCReco->ofcbOverflow() == false) {
          notOFCbOF = true;
        } else {
          if (notMasked) {
            OFCbOFNotMasked = true;
          }
        }
        if (std::abs(SC_time) > 3) {
          tauGt3 = true;
        }
 
        if (notMasked && notSatur && notOFCbOF) {
          if (nonZeroETofl) {
            notMaskedEoflNe0 = true;
          }
          if (eToflGt1GeV) {
            notMaskedEoflGt1 = true;
          }
          if  (SC_energy_onl == -99999) {
            passSCNomInvalid = true;
          }
 
          //if ( passTauSel ){
          //if ( nonZeroET && eTgt0lt0p325GeV ){
          if (eTgt0lt0p325GeV) {
              passSCNom0_0p325 = true;
            }
            if (eTgt0p325GeV && eTlt1GeV) {
              passSCNom0p325_1 = true;
            }
            if (eTgt1GeV) {
              passSCNom1 = true;
            }
            if (eTgt10GeV) {
              passSCNom10 = true;
              if (tauGt3) {
                passSCNom10tauGt3 = true;
              }
            }
            if (SC_energy_onl != SC_energy_ofl) {
              onlofflEmismatch = true;
            }
        //}
        }  // end nominal selections
 
        lvaluemap_sc.emplace_back(SC_eta, SC_phi, SC_ieta, SC_iphi, SC_latomeSourceIdBIN, SC_ET_ofl, SC_ET_diff, SC_ET_onl, SC_ET_onl_muscaled, SC_time, BCID,
                                  lumi_block, zeroET, passSCNomInvalid, passSCNom0_0p325, passSCNom0p325_1, passSCNom1, passSCNom10, passSCNom10tauGt3, saturNotMasked, OFCbOFNotMasked, notMaskedEoflNe0,
                                  notMaskedEoflGt1);
        lvaluemap_sc_ALL.emplace_back(SC_eta, SC_phi, SC_ieta, SC_iphi, SC_latomeSourceIdBIN, SC_ET_ofl, SC_ET_diff, SC_ET_onl, SC_ET_onl_muscaled, SC_time,
                                      BCID, lumi_block, zeroET, passSCNomInvalid, passSCNom0_0p325, passSCNom0p325_1, passSCNom1, passSCNom10, passSCNom10tauGt3, saturNotMasked, OFCbOFNotMasked,
                                      notMaskedEoflNe0, notMaskedEoflGt1);
 
      }  // end loop over SCs
 
      // fill, for every layer/threshold
      for (size_t ilayer = 0; ilayer < scMonValueVec.size(); ++ilayer) {
        const auto& tool = scMonValueVec[ilayer];
        auto sc_part_eta = Monitored::Collection("SC_part_eta", tool, [](const auto& v) { return v.sc_eta; });
        auto sc_part_phi = Monitored::Collection("SC_part_phi", tool, [](const auto& v) { return v.sc_phi; });
        auto sc_part_ieta = Monitored::Collection("SC_part_ieta", tool, [](const auto& v) { return v.sc_ieta; });
        auto sc_part_iphi = Monitored::Collection("SC_part_iphi", tool, [](const auto& v) { return v.sc_iphi; });
        auto sc_part_latomesourceidbin = Monitored::Collection("SC_part_latomesourceidbin", tool, [](const auto& v) { return v.sc_latomesourceidbin; });
        auto sc_part_et_ofl = Monitored::Collection("SC_part_et_ofl", tool, [](const auto& v) { return v.sc_et_ofl; });
        auto sc_part_et_diff = Monitored::Collection("SC_part_et_diff", tool, [](const auto& v) { return v.sc_et_diff; });
        auto sc_part_et_onl = Monitored::Collection("SC_part_et_onl", tool, [](const auto& v) { return v.sc_et_onl; });
        auto sc_part_et_onl_muscaled = Monitored::Collection("SC_part_et_onl_muscaled", tool, [](const auto& v) { return v.sc_et_onl_muscaled; });
        auto sc_part_time = Monitored::Collection("SC_part_time", tool, [](const auto& v) { return v.sc_time; });
        auto sc_part_bcid = Monitored::Collection("SC_part_BCID", tool, [](const auto& v) { return v.sc_bcid; });
        auto sc_part_lb = Monitored::Collection("SC_part_LB", tool, [](const auto& v) { return v.sc_lb; });
        // auto sc_part_passSCNom = Monitored::Collection("SC_part_passSCNom", tool, [](const auto& v) { return v.sc_passSCNom; });
        auto sc_zeroET = Monitored::Collection("SC_part_zeroET", tool, [](const auto& v) { return v.sc_zeroET; });
        auto sc_part_passSCNomInvalid = Monitored::Collection("SC_part_passSCNomInvalid", tool, [](const auto& v) { return v.sc_passSCNomInvalid; });
        auto sc_part_passSCNom0_0p325 = Monitored::Collection("SC_part_passSCNom0_0p325", tool, [](const auto& v) { return v.sc_passSCNom0_0p325; });
        auto sc_part_passSCNom0p325_1 = Monitored::Collection("SC_part_passSCNom0p325_1", tool, [](const auto& v) { return v.sc_passSCNom0p325_1; });
        auto sc_part_passSCNom1 = Monitored::Collection("SC_part_passSCNom1", tool, [](const auto& v) { return v.sc_passSCNom1; });
        auto sc_part_passSCNom10 = Monitored::Collection("SC_part_passSCNom10", tool, [](const auto& v) { return v.sc_passSCNom10; });
        auto sc_part_passSCNom10tauGt3 = Monitored::Collection("SC_part_passSCNom10tauGt3", tool, [](const auto& v) { return v.sc_passSCNom10tauGt3; });
        auto sc_part_saturNotMasked = Monitored::Collection("SC_part_saturNotMasked", tool, [](const auto& v) { return v.sc_saturNotMasked; });
        auto sc_part_OFCbOFNotMasked = Monitored::Collection("SC_part_OFCbOFNotMasked", tool, [](const auto& v) { return v.sc_OFCbOFNotMasked; });
        auto sc_part_notMaskedEoflNe0 = Monitored::Collection("SC_part_notMaskedEoflNe0", tool, [](const auto& v) { return v.sc_notMaskedEoflNe0; });
        auto sc_part_notMaskedEoflGt1 = Monitored::Collection("SC_part_notMaskedEoflGt1", tool, [](const auto& v) { return v.sc_notMaskedEoflGt1; });
 
 
        fill(m_tools[m_toolmapLayerNames_sc.at(m_layerNames[ilayer])], sc_part_eta, sc_part_phi, sc_part_ieta, sc_part_iphi, sc_part_latomesourceidbin,
             sc_part_et_ofl, sc_part_et_diff, sc_part_et_onl, sc_part_et_onl_muscaled, sc_part_time, sc_part_bcid, sc_part_lb, sc_zeroET, sc_part_passSCNomInvalid,
             sc_part_passSCNom0_0p325, sc_part_passSCNom0p325_1, sc_part_passSCNom1, sc_part_passSCNom10, sc_part_passSCNom10tauGt3, sc_part_saturNotMasked, sc_part_OFCbOFNotMasked, sc_part_notMaskedEoflNe0,
             sc_part_notMaskedEoflGt1);
      }
 
 
    }  // End if(LArSCContainer is valid)
 
 
 
    // LATOME event size
    if ((hLArLATOMEHeaderContainer.isValid())) {
      auto event_size = Monitored::Scalar<float>("event_size", 0);
      for (const LArLATOMEHeader* pLArLATOMEHeader : *hLArLATOMEHeaderContainer) {
        event_size += pLArLATOMEHeader->ROBFragSize() + 48;  // 48 is the offset between rod_ndata and ROB fragment size
      }
      event_size /= (1024 * 1024 / 4);
      fill(m_scMonGroupName, lumi_block, event_size);
    }
 
    // end LATOME event size
 
    return StatusCode::SUCCESS;
  }
 
/*---------------------------------------------------------*/
int LArDigitalTriggMonAlg::whatPartition(HWIdentifier id, int side) const
{
  if (m_LArOnlineIDHelper->isEMBchannel(id)) {
    if(side==0) return 0;
    else return 1;
  } else if (m_LArOnlineIDHelper->isEMECchannel(id)) {
    if(side==0) return 2;
    else return 3;
  } else if (m_LArOnlineIDHelper->isHECchannel(id)) {
    if(side==0) return 4;
    else return 5;
  } else {
    if(side==0) return 6;
    else return 7;
  }
}
 
 
 
unsigned LArDigitalTriggMonAlg::getXbinFromSourceID(const unsigned sourceID) const
{
  //  int NLatomeBins = 117;
  int detStartingBin=m_NLatomeBins;
  const unsigned detID = sourceID >> 16;
  const unsigned value = sourceID & 0xF;
  auto mapit=m_LatomeDetBinMappingQ.find(detID);
  if (mapit!=m_LatomeDetBinMappingQ.end()) {
    detStartingBin=mapit->second;
  }
 
  unsigned binx = detStartingBin+value;
  if (binx>m_NLatomeBins){
    ATH_MSG_WARNING("something wrong with binning, filling overflowbin");
    binx=m_NLatomeBins;
  }
 
  return binx;
}