LArSuperCellMonAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
// NAME:     LArSuperCellMonAlg.cxx
//           based on LArCellMonAlg:
//                L. Morvaj, P.Strizenec, D. Oliveira Damazio - develop for Digital Trigger monitoring (2021)
// ********************************************************************
#include "LArSuperCellMonAlg.h"
 
 
#include "CaloDetDescr/CaloDetDescrElement.h"
#include "CaloIdentifier/CaloGain.h"
#include "Identifier/Identifier.h"
#include "Identifier/HWIdentifier.h"
#include "CaloIdentifier/CaloCell_ID.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
 
//#include "AthenaMonitoring/DQBadLBFilterTool.h"
//#include "AthenaMonitoring/DQAtlasReadyFilterTool.h"
 
#include "CaloEvent/CaloCellContainer.h"
#include "AthenaKernel/Units.h"
 
#include "LArElecCalib/LArProvenance.h"
 
#include <cassert>
#include <algorithm>
 
LArSuperCellMonAlg::LArSuperCellMonAlg(const std::string& name, ISvcLocator* pSvcLocator) 
  :AthMonitorAlgorithm(name, pSvcLocator),
   //   m_badChannelMask("BadLArRawChannelMask",this),
   m_calo_id(nullptr)
{    
}
 
 
 
LArSuperCellMonAlg::~LArSuperCellMonAlg() {
}
 
StatusCode LArSuperCellMonAlg::initialize() {
 
  ATH_MSG_DEBUG("LArSuperCellMonAlg::initialize() start");
 
 
  //Identfier-helpers 
  ATH_CHECK( detStore()->retrieve(m_calo_id) );
 
  ATH_CHECK( m_superCellContainerKey.initialize() );
  ATH_CHECK( m_superCellContainerRefKey.initialize() );
  ATH_CHECK( m_noiseCDOKey.initialize() );
  ATH_CHECK( m_bcDataKey.initialize(SG::AllowEmpty) );
 
  if(m_superCellContainerRecoKey.empty()) m_doSCReco=false;
  ATH_CHECK(m_superCellContainerRecoKey.initialize(m_doSCReco));  
 
  ATH_MSG_DEBUG("LArSuperCellMonAlg::initialize() is done!");
 
  return AthMonitorAlgorithm::initialize();
}
 
 
 
StatusCode LArSuperCellMonAlg::initThresh() {
 
  return StatusCode::SUCCESS;
}
 
 
StatusCode LArSuperCellMonAlg::fillHistograms(const EventContext& ctx) const{  
 
  ATH_MSG_DEBUG("LArSuperCellMonAlg::fillHistograms() starts");
 
  SG::ReadHandle<CaloCellContainer> superCellHdl{m_superCellContainerKey, ctx};
  const CaloCellContainer* superCellCont = superCellHdl.cptr();
  if(!superCellCont){
     ATH_MSG_ERROR("The requested SC container key " << m_superCellContainerKey.key() << " could not be retrieved. !!!");
     return StatusCode::SUCCESS;
  }
  
  SG::ReadHandle<CaloCellContainer> superCellRefHdl{m_superCellContainerRefKey, ctx};
  const CaloCellContainer* superCellRefCont = superCellRefHdl.cptr();
  if(!superCellRefCont){
     ATH_MSG_ERROR("The requested SC ref container key " << m_superCellContainerRefKey.key() << " could not be retrieved. !!!");
     return StatusCode::SUCCESS;
  }
  
  const CaloCellContainer *superCellRecoCont = nullptr;
  if(m_doSCReco){
     SG::ReadHandle<CaloCellContainer > hSCetRecoContainer{m_superCellContainerRecoKey,ctx}; 
     if (!hSCetRecoContainer.isValid()) {
        ATH_MSG_ERROR("The requested SC ET reco container key could not be retrieved. !!!");
     }else{
        superCellRecoCont = hSCetRecoContainer.cptr(); 
        ATH_MSG_DEBUG("SCetRecoContainer.size() " << hSCetRecoContainer->size());
     }
  }
     
  if (ctx.evt()==0) {
    SG::ReadCondHandle<CaloNoise> noiseHdl{m_noiseCDOKey, ctx};
    const CaloNoise *noisep = *noiseHdl;
    ATH_CHECK(createPerJobHistograms(superCellCont, noisep));
  }
 
  //get LB
  auto lumiBlock = Monitored::Scalar<unsigned int>("lumiBlock",0);
  lumiBlock = ctx.eventID().lumi_block();
  auto bcid = Monitored::Scalar<unsigned int>("bcid",0);
  bcid = ctx.eventID().bunch_crossing_id();
  int bcidFFB = bcid;
  if (!m_bcDataKey.empty()){ 
    SG::ReadCondHandle<BunchCrossingCondData> bccd (m_bcDataKey,ctx);
    bcid=bccd->distanceFromFront(bcid,BunchCrossingCondData::BunchCrossings);
  }
 
  // create local variables to speed up things
  // per layer
  std::vector<std::vector<std::string> > nameHistos;
  for ( auto& layerName : m_layerNames){
        std::vector<std::string> vec;
    vec.push_back("superCellEt_"+layerName);
    vec.push_back("superCelltime_"+layerName);
    vec.push_back("superCellprovenance_"+layerName);
    vec.push_back("superCellEta_"+layerName);
    vec.push_back("superCellPhi_"+layerName);
    vec.push_back("resolution_"+layerName);
    vec.push_back("resolutionPass_"+layerName);
    vec.push_back("resolutionHET_"+layerName);
    vec.push_back("superCellEtRef_"+layerName);
    vec.push_back("superCelltimeRef_"+layerName);
    vec.push_back("superCellprovenanceRef_"+layerName);
    vec.push_back("superCellEtDiff_"+layerName);
        nameHistos.push_back(vec);
  }
 
 
 
 
  CaloCellContainer::const_iterator SCit = superCellCont->begin(); 
  CaloCellContainer::const_iterator SCit_e = superCellCont->end(); 
 
  std::vector<std::reference_wrapper<Monitored::IMonitoredVariable>> variables;
 
  for ( ; SCit!=SCit_e;++SCit) {
 
    const CaloCell* superCell = *SCit; 
    variables.clear();
    // Discard masked cells from monitoring
    int SCprov = superCell->provenance()&0xFFF;
    if (m_removeMasked && ((SCprov&0x80)==0x80)) continue;
    float SCet = superCell->et();
    const CaloDetDescrElement* SCcaloDDE = superCell->caloDDE();
    float SCeta,SCphi;
    unsigned iLyr, iLyrNS;
    float SCt = superCell->time();
    
    getHistoCoordinates(SCcaloDDE, SCeta, SCphi, iLyr, iLyrNS);
    
 
    bool SCpassTime = LArProv::test(SCprov,LArProv::SCTIMEPASS);//SCprov & 0x200;
    bool SCpassPF =   LArProv::test(SCprov,LArProv::SCPASSBCIDMAX);// SCprov & 0x40;
 
    const CaloCell* superCellRef = superCellRefCont->findCell( SCcaloDDE->identifyHash() );
    float SCetRef = superCellRef->et();
    float SCetDiff = SCet - SCetRef;
    float resolution = -1000;
    float resolutionPass = -1000;
    float resolutionHET = -1000;
    if ( SCetRef > m_thresholdsForResolution) resolution = 100.0*(SCet-SCetRef)/SCetRef;
    if ( SCpassTime || SCpassPF ) resolutionPass = resolution;
    if ( SCet > 4e3 ) resolutionHET=resolution;
 
    // real monitoring business
    auto MSCet = Monitored::Scalar<float>("superCellEt",SCet);
    auto MSCt = Monitored::Scalar<float>("superCelltime",SCt);
    auto MSCprov = Monitored::Scalar<int>("superCellprovenance",SCprov);
    auto MSCeta = Monitored::Scalar<float>("superCellEta",SCeta);
    auto MSCphi = Monitored::Scalar<float>("superCellPhi",SCphi);
    auto MSCres = Monitored::Scalar<float>("resolution",resolution);
    auto MSCresPass = Monitored::Scalar<float>("resolutionPass",resolutionPass);
    auto MSCresHET = Monitored::Scalar<float>("resolutionHET",resolutionHET);
    auto MSCetRef = Monitored::Scalar<float>("superCellEtRef",SCetRef);
    auto MSCtRef = Monitored::Scalar<float>("superCelltimeRef",superCellRef->time());
    auto MSCprovRef = Monitored::Scalar<int>("superCellprovenanceRef",(superCellRef->provenance()&0xFFF));
    auto MSCetDiff = Monitored::Scalar<float>("superCellEtDiff",SCetDiff);
    variables.push_back(MSCet);
    variables.push_back(MSCt);
    variables.push_back(MSCprov);
    variables.push_back(MSCeta);
    if (  SCetRef > m_thresholdsForResolution ) variables.push_back(MSCres);
    if ( (SCetRef > m_thresholdsForResolution ) && (SCpassTime || SCpassPF ) ) variables.push_back(MSCresPass);
    if ( (SCetRef > m_thresholdsForResolution ) && (SCet > 4e3 ) ) variables.push_back(MSCresHET);
    variables.push_back(MSCphi);
    variables.push_back(MSCetRef);
    // let us put conditional to force building the linearity plot
    // only when the new signal passes BCID
    variables.push_back(MSCtRef);
    variables.push_back(MSCprovRef);
    variables.push_back(MSCetDiff);
 
    // per layer
    auto layerName=m_layerNames[iLyr];
    auto LMSCet = Monitored::Scalar<float>(nameHistos[iLyr][0],SCet);
    auto LMSCt = Monitored::Scalar<float>(nameHistos[iLyr][1],SCt);
    auto LMSCprov = Monitored::Scalar<int>(nameHistos[iLyr][2],SCprov);
    auto LMSCeta = Monitored::Scalar<float>(nameHistos[iLyr][3],SCeta);
    auto LMSCphi = Monitored::Scalar<float>(nameHistos[iLyr][4],SCphi);
    auto LMSCres = Monitored::Scalar<float>(nameHistos[iLyr][5],resolution);
    auto LMSCresPass = Monitored::Scalar<float>(nameHistos[iLyr][6],resolutionPass);
    auto LMSCresHET = Monitored::Scalar<float>(nameHistos[iLyr][7],resolutionHET);
    auto LMSCetRef = Monitored::Scalar<float>(nameHistos[iLyr][8],SCetRef);
    auto LMSCtRef = Monitored::Scalar<float>(nameHistos[iLyr][9],superCellRef->time());
    auto LMSCprovRef = Monitored::Scalar<int>(nameHistos[iLyr][10],(superCellRef->provenance()&0xFFF));
 
    auto MBCIDFFB = Monitored::Scalar<int>("BCID",bcidFFB);
    auto LMSCetDiff = Monitored::Scalar<float>(nameHistos[iLyr][11],SCetDiff);
    variables.push_back(LMSCet);
    variables.push_back(LMSCt);
    variables.push_back(LMSCprov);
    variables.push_back(LMSCeta);
    if (  SCetRef > m_thresholdsForResolution ) variables.push_back(LMSCres);
    if ( (SCetRef > m_thresholdsForResolution ) && (SCpassTime || SCpassPF ) ) variables.push_back(LMSCresPass);
    if ( (SCetRef > m_thresholdsForResolution ) && (SCet > 4e3 ) ) variables.push_back(LMSCresHET);
    variables.push_back(LMSCphi);
    variables.push_back(LMSCetRef);
    if ( SCpassTime || SCpassPF ) variables.push_back(LMSCtRef);
    variables.push_back(LMSCprovRef);
    variables.push_back(MBCIDFFB);
    variables.push_back(LMSCetDiff);
 
    if(m_doSCReco){
       auto MSCtReco = Monitored::Scalar<float>("superCelltimeReco",0.);
       auto MSCetReco = Monitored::Scalar<float>("superCellEtReco",0.);
       auto LMSCtReco = Monitored::Scalar<float>("superCelltimeReco_"+layerName,0.);
       const CaloCell* superCellReco = superCellRecoCont->findCell( SCcaloDDE->identifyHash() );
       if(superCellReco) {
          float SCetReco = superCellReco->et();
          float SCtimeReco = superCellReco->time();
          MSCtReco = SCtimeReco;
          MSCetReco = SCetReco;
          LMSCtReco = SCtimeReco;
          variables.push_back(MSCtReco);
          variables.push_back(LMSCtReco);
          variables.push_back(MSCetReco);
       }
    } 
    fill(m_MonGroupName,variables);
 
  } // end loop over SC
  
 
 
 
  ATH_MSG_DEBUG("LArSuperCellMonAlg::fillLarHists() is done");
  return StatusCode::SUCCESS;
}
 
 
StatusCode LArSuperCellMonAlg::createPerJobHistograms(const CaloCellContainer* cellCont, const CaloNoise *noisep ) const {
 
  ATH_MSG_INFO("Creating the once-per-job histograms");
 
  if(!noisep){
    ATH_MSG_ERROR("Do not have DB noise, doing nothing !!!");
    return StatusCode::SUCCESS;
  }
 
  //The following histograms can be considered constants for one job
  //(in fact, they are constant for an entire run or even run-periode)
  //ActiveCells in eta/phi (to normalize 1D occupancy plots)
  //BadChannel word
  //Database noise
   
  auto doDatabaseNoisePlot = Monitored::Scalar<bool>("doDatabaseNoisePlot",m_doDatabaseNoiseVsEtaPhi);
 
  //  if(!doDatabaseNoisePlot && !doCellsActiveEtaPlot && !doCellsActivePhiPlot) {
  if(!doDatabaseNoisePlot) {
    ATH_MSG_INFO("No once-per-job histogram requested");
    return StatusCode::SUCCESS;
  }
 
  
  //filling:
 
  CaloCellContainer::const_iterator it = cellCont->begin(); 
  CaloCellContainer::const_iterator it_e = cellCont->end(); 
  for ( ; it!=it_e;++it) { 
    const CaloCell* cell = *it; 
    Identifier id = cell->ID();
    bool is_lar=m_calo_id->is_lar(id);
    if(!is_lar) continue;
    const CaloDetDescrElement* caloDDEl=cell->caloDDE();
    float celleta, cellphi;
    unsigned iLyr, iLyrNS;
 
 
    getHistoCoordinates(caloDDEl, celleta, cellphi, iLyr, iLyrNS);
    
    auto mon_eta = Monitored::Scalar<float>("celleta_"+m_layerNames[iLyr],celleta);
    auto mon_phi = Monitored::Scalar<float>("cellphi_"+m_layerNames[iLyr],cellphi);   
    auto cellnoisedb = Monitored::Scalar<float>("cellnoisedb_"+m_layerNames[iLyr],noisep->getNoise(id,cell->gain()));
 
    fill(m_MonGroupName,cellnoisedb,mon_eta,mon_phi);
 
    
  }//end loop over cells
 
  return StatusCode::SUCCESS;
}
 
 
 
std::string  LArSuperCellMonAlg::strToLower(const std::string& input) const {
  std::string output;
  for (const auto& c : input) {
    output.push_back(std::tolower(c));
  }
  return output;
}
 
 
 
void LArSuperCellMonAlg::getHistoCoordinates(const CaloDetDescrElement* dde, float& celleta, float& cellphi, unsigned& iLyr, unsigned& iLyrNS) const {
  
  celleta=dde->eta_raw();
  cellphi=dde->phi_raw();
  
  int calosample=dde->getSampling();
  if (dde->is_lar_em_endcap_inner()) calosample-=1; //Here, we consider the two layers of the EMEC IW as EME1 and EME2 instad of layer 2 and 3
  iLyrNS=m_caloSamplingToLyrNS.at(calosample); //will throw if out of bounds
  if ((iLyrNS==EMB1NS || iLyrNS==EMB2NS) && m_calo_id->region(dde->identify())==1) {
    //We are in the awkward region 1 of the EM Barrel
    //Looking at the image at http://atlas.web.cern.ch/Atlas/GROUPS/LIQARGEXT/TDR/figures6/figure6-17.eps
    //may be useful to understand what's going on here. 
 
    //In brief: Region 1, layer 1 is closer related ot the middle compartment (aka layer 2)
    //          and region 1, layer 2 closer related to the back compartment (aka layer 3)
    iLyrNS+=1;
 
    //Layer 2: 0<eta<1.4 + 1.4<eta<1.475, deta = 0.025. 3 rows of cells from region 1
    //Layer 3: 0<eta<1.35 (deta=0.050) +  1.4<eta<1.475 (deta = 0.075).  1 row of cell from region 1 with different dEta
  }
  
  const unsigned side=(celleta>0) ? 0 : 1; //Value >0 means A-side
  iLyr=iLyrNS*2+side;  //Getting LayerEnum value. This logic works because of the way the enums LayerEnum and LayerEnumNoSides are set up. 
  return;
}