CaloCellEnergyRescaler.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "CaloCellEnergyRescaler.h" 
#include "CaloEvent/CaloCellContainer.h"
#include "CaloCondBlobObjs/CaloCondBlobFlt.h"
#include <memory>
 
CaloCellEnergyRescaler::CaloCellEnergyRescaler (const std::string& type, 
                        const std::string& name, 
                        const IInterface* parent) :
  base_class (type, name, parent) {}
 
CaloCellEnergyRescaler::~CaloCellEnergyRescaler() = default;
 
 
StatusCode CaloCellEnergyRescaler::initialize() {
  ATH_CHECK(m_eneShiftFldr.initialize());
  return StatusCode::SUCCESS;
}
 
 
StatusCode
CaloCellEnergyRescaler::process (CaloCellContainer* theCaloCellContainer,
                                 const EventContext& ctx) const {
 
  SG::ReadCondHandle<AthenaAttributeList> eneShiftHdl(m_eneShiftFldr,ctx);
 
  //Possible optimization: If the following lines turn out to be slow to be executed on every event, 
  //put them into a conditions algo
  const coral::Blob& blob = (**eneShiftHdl)["CaloCondBlob16M"].data<coral::Blob>();
  if (blob.size()<3) {
    ATH_MSG_DEBUG("Found empty blob, no corretion needed");
    return StatusCode::SUCCESS;
  }
 
  std::unique_ptr<const CaloCondBlobFlt> corrValues(CaloCondBlobFlt::getInstance(blob));
  ATH_MSG_DEBUG("Database folder has values for " << corrValues->getNChans() << " channels and " << corrValues->getNGains() << " gains.");
  
 
  CaloCellContainer::iterator it=theCaloCellContainer->begin();
  CaloCellContainer::iterator it_e=theCaloCellContainer->end();
  for(;it!=it_e;++it) {
    CaloCell* theCell=(*it);
    const IdentifierHash& hash_id=theCell->caloDDE()->calo_hash();
    if (hash_id<corrValues->getNChans()) {
      const float& scale= corrValues->getData(hash_id);
      theCell->setEnergy(theCell->energy()*scale);
    }//end if hash_id<NChans
  }//end loop over cells
  return StatusCode::SUCCESS;
}