LArDSPThresholds2Ntuple.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArCalibTools/LArDSPThresholds2Ntuple.h"
#include "LArIdentifier/LArOnlineID.h"
 
#include "LArRawConditions/LArDSPThresholdsComplete.h"
#include "LArCOOLConditions/LArDSPThresholdsFlat.h"
 
#include "CoralBase/Blob.h"
#include "AthenaPoolUtilities/AthenaAttributeList.h"
 
LArDSPThresholds2Ntuple::LArDSPThresholds2Ntuple(const std::string& name, ISvcLocator* pSvcLocator): 
  LArCond2NtupleBase(name, pSvcLocator){
 
  declareProperty("NtupleTitle",m_ntTitle="Gain Thresholds");
  declareProperty("NtupleName",m_ntname="DSPTHR");
  declareProperty("DumpFlat",m_doFlat=false);
  declareProperty("FlatFolder",m_folder="");
  m_ntpath="/NTUPLES/FILE1/"+m_ntname;
  m_addFEBTemp=false;
}
 
 
StatusCode LArDSPThresholds2Ntuple::initialize() {
 
  StatusCode sc=LArCond2NtupleBase::initialize();
   if (sc!=StatusCode::SUCCESS) {
     ATH_MSG_ERROR( "Base init failed" );
     return StatusCode::FAILURE;
   }
 
   sc=m_nt->addItem("tQThr",m_tQThr);
   if (sc!=StatusCode::SUCCESS) {
     ATH_MSG_ERROR( "addItem tQThr failed" );
     return StatusCode::FAILURE;
   }
   
   sc=m_nt->addItem("samplesThr",m_samplesThr);
   if (sc!=StatusCode::SUCCESS) {
     ATH_MSG_ERROR( "addItem samplesThr failed" );
     return StatusCode::FAILURE;
   }
   
   sc=m_nt->addItem("trigThr",m_trigThr);
   if (sc!=StatusCode::SUCCESS) {
     ATH_MSG_ERROR( "addItem trigThr failed" );
     return StatusCode::FAILURE;
   }
 
   return StatusCode::SUCCESS;
}
 
LArDSPThresholds2Ntuple::~LArDSPThresholds2Ntuple() 
= default;
 
StatusCode LArDSPThresholds2Ntuple::stop() {
  StatusCode sc;
   
   const LArDSPThresholdsComplete *dc = nullptr;
   const LArDSPThresholdsFlat *df = nullptr;
   const AthenaAttributeList* attrList=nullptr;
   if(m_doFlat) {
      sc=detStore()->retrieve(attrList,m_folder);
      if (sc.isFailure()) {
         ATH_MSG_ERROR( "Failed to retrieve AthenaAttributeList with key " << m_folder );
         return sc;
      }
      
      const coral::Blob& blob = (attrList->coralList())["tQThr"].data<coral::Blob>();
      if (blob.size()<3) {
        ATH_MSG_INFO( "Found empty blob, nothing to do");
        return StatusCode::SUCCESS;
      }
      
      df=new LArDSPThresholdsFlat(attrList);
 
   } else {
      sc=m_detStore->retrieve(dc);
      if(sc!=StatusCode::SUCCESS) {
           ATH_MSG_ERROR("Could not retrieve LArDSPThresholdsComplete....");
           return StatusCode::FAILURE;
      }
   }
   
   for (const HWIdentifier hwid: m_onlineId->channel_range()) {
     if(m_doFlat) {
        m_tQThr=df->tQThr(hwid);
        m_samplesThr=df->samplesThr(hwid);
        m_trigThr=df->trigSumThr(hwid);
     } else {
         m_tQThr=dc->tQThr(hwid);
         m_samplesThr=dc->samplesThr(hwid);
         m_trigThr=dc->trigSumThr(hwid);
     }
     
     fillFromIdentifier(hwid);
     //ATH_MSG_INFO("hwid: "<<hwid.getString()<<" "<<tQThr<<" : "<<samplesThr<<" : "<<trigThr);
     
     sc=ntupleSvc()->writeRecord(m_nt);      
     if (sc!=StatusCode::SUCCESS) {
       ATH_MSG_ERROR( "writeRecord failed" );
       return StatusCode::FAILURE;
     }
   }
   
 
   ATH_MSG_INFO( "LArDSPThresholds2Ntuple has finished." );
   return StatusCode::SUCCESS;
   
}// end finalize-method.