T2CaloEgammaReFastAlgo.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
/*
 NAME:     T2CaloEgammaReFastAlgo.cxx
 PACKAGE:  Trigger/TrigAlgorithms/TrigT2CaloEgamma
 
 AUTHOR:   Denis Oliveira Damazio
           Carlos Chavez Barajas
 
*/
 
#include "GaudiKernel/StatusCode.h"
 
#include "AthLinks/ElementLink.h"
#include "T2CaloEgammaReFastAlgo.h"
#include "TrigT2CaloCommon/IReAlgToolCalo.h"
#include "CxxUtils/phihelper.h"
#include "xAODTrigCalo/TrigEMClusterAuxContainer.h"
#include "xAODTrigCalo/TrigEMClusterContainer.h"
#include "AthenaMonitoringKernel/Monitored.h"
 
class ISvcLocator;
 
T2CaloEgammaReFastAlgo::T2CaloEgammaReFastAlgo(const std::string& name, ISvcLocator* pSvcLocator) :
    AthReentrantAlgorithm(name, pSvcLocator)
{}
 
StatusCode T2CaloEgammaReFastAlgo::initialize()
{
  m_emAlgTools.retrieve().ignore();
  ATH_CHECK(m_clusterContainerKey.initialize());
  ATH_CHECK(m_roiCollectionKey.initialize());
  ATH_CHECK( m_bcidAvgKey.initialize() );
  if (! m_monTool.empty() ) ATH_CHECK( m_monTool.retrieve() );
  if(m_doCalibWithRings) ATH_CHECK( m_calibWRingsTool.retrieve() );
 
  return StatusCode::SUCCESS;
}
 
StatusCode T2CaloEgammaReFastAlgo::execute(const EventContext& context) const
{
 
  auto timer = Monitored::Timer("TIME_exec");
  auto clET = Monitored::Scalar("TrigEMCluster_eT",-999.0);
  auto clHET = Monitored::Scalar("TrigEMCluster_had1",-999.0);
  auto clEta = Monitored::Scalar("TrigEMCluster_eta",-999.0);
  auto clPhi = Monitored::Scalar("TrigEMCluster_phi",-999.0);
  auto clReta = Monitored::Scalar("TrigEMCluster_rEta",-999.0);
  auto clETrings = Monitored::Scalar("TrigEMCluster_et_rings",-999.0); 
  auto res_et = Monitored::Scalar("Resolution_et",-999.0);
  auto monitoring = Monitored::Group( m_monTool, timer, clET, clHET, clEta, clPhi, clReta,clETrings,res_et);
 
 
  SG::WriteHandle<xAOD::TrigEMClusterContainer> trigEmClusterCollection(m_clusterContainerKey, context);
  ATH_CHECK( trigEmClusterCollection.record(std::make_unique<xAOD::TrigEMClusterContainer>(),
                                            std::make_unique<xAOD::TrigEMClusterAuxContainer>()) );
 
 
  auto roisHandle = SG::makeHandle(m_roiCollectionKey, context);
  if (!roisHandle.isValid()) {
    ATH_MSG_DEBUG("no RoI");
    return StatusCode::SUCCESS;
  }
 
 
  trigEmClusterCollection->reserve(roisHandle->size());
  
  ATH_MSG_DEBUG("RoI descriptor size is " << roisHandle->size() );
 
  for (const TrigRoiDescriptor* roiDescriptor : *roisHandle) {
    float etaL1, phiL1;
    double etamin, etamax, phimin, phimax;
 
    ATH_MSG_DEBUG( "RoI eta = " << roiDescriptor->eta() << " RoI phi = " << roiDescriptor->phi());
 
    if ((m_l1eta < -9.9) && (m_l1phi < -9.9)) {
      
      etamin = std::max(-2.5, roiDescriptor->eta() - m_etaWidth);
      etamax = std::min( 2.5, roiDescriptor->eta() + m_etaWidth);
 
      phimin = CxxUtils::wrapToPi(roiDescriptor->phi() - m_phiWidth);
      phimax = CxxUtils::wrapToPi(roiDescriptor->phi() + m_phiWidth);
 
      etaL1 = roiDescriptor->eta();
      phiL1 = roiDescriptor->phi();
    }
    else {
      etamin = std::max(-2.5, m_l1eta - m_etaWidth);
      etamax = std::min( 2.5, static_cast<double>(m_l1eta) + m_etaWidth);
 
      phimin = CxxUtils::wrapToPi(m_l1phi - m_phiWidth);
      phimax = CxxUtils::wrapToPi(static_cast<double>(m_l1phi) + m_phiWidth);
 
      etaL1 = m_l1eta;
      phiL1 = m_l1phi;
    }
 
    TrigRoiDescriptor newroi(roiDescriptor->eta(), etamin, etamax,
                             roiDescriptor->phi(), phimin, phimax);
 
 
    ATH_MSG_DEBUG(" etamin = " << etamin << " etamax = " << etamax <<
                  " phimin = " << phimin << " phimax = " << phimax);
 
    xAOD::TrigEMCluster* ptrigEmCluster = new xAOD::TrigEMCluster();
    trigEmClusterCollection->push_back(ptrigEmCluster);
    ptrigEmCluster->setEnergy(0.0);
    ptrigEmCluster->setEt(0.0);
    ptrigEmCluster->setRawEnergy(0.0);
    ptrigEmCluster->setRawEt(0.0);
    ptrigEmCluster->setE277(0);
    ptrigEmCluster->setEmaxs1(0);
    ptrigEmCluster->setE2tsts1(0);
    ptrigEmCluster->setEhad1(-999);
    ptrigEmCluster->setWeta2(-999);
    ptrigEmCluster->setFracs1(-999);
    ptrigEmCluster->setE233(-999);
    ptrigEmCluster->setE237(-999);
    ptrigEmCluster->setWstot(-999);
    ptrigEmCluster->setEta1(-999);
    ptrigEmCluster->setNCells(0);
    ptrigEmCluster->setRawEta(-999);
    ptrigEmCluster->setRawPhi(-999);
 
    // It is a good idea to clear the energies
    for (int i = 0; i < CaloSampling::CaloSample::MINIFCAL0; i++) {
      ptrigEmCluster->setEnergy((CaloSampling::CaloSample)i, 0.);
      ptrigEmCluster->setRawEnergy((CaloSampling::CaloSample)i, 0.);
    }
    // Initial cluster position is the LVL1 position
    ptrigEmCluster->setEta(etaL1);
    ptrigEmCluster->setPhi(phiL1);
 
    // Add RoI word to TrigEMCluster
    // Dangerous !!!! we need to define a *new* roiDescriptor if we want to
    // change the size, so we should be careful about *which* roi descriptor
    // we save, and *which* "roiWord" (if any) we store if we need to use it
    // again
    ptrigEmCluster->setRoIword(roiDescriptor->roiWord());
    const CaloDetDescrElement* caloDDE = nullptr;
 
    uint32_t error = 0;
    for (const auto& tool : m_emAlgTools) {
      ATH_CHECK( tool->execute(*ptrigEmCluster, newroi, caloDDE, context) );
    }
    //  // support to new monitoring
    // Cluster quality is a collection of possible errors
    // No error quality=0
    ptrigEmCluster->setClusterQuality(error);
    if ((error & 0xC0000000) || ptrigEmCluster->phi() < -M_PI || ptrigEmCluster->phi() > +M_PI ||
        std::abs(ptrigEmCluster->eta()) > 10.0) {
      // Clustering failed. Transmit ahead L1
      ptrigEmCluster->setEta(etaL1);
      ptrigEmCluster->setPhi(phiL1);
      ptrigEmCluster->setEnergy(0.0);
      ptrigEmCluster->setEt(0.0);
    }
 
 
      if ( caloDDE != nullptr ){
        if ( caloDDE->is_lar_em_barrel() ){
          for( ToolHandleArray<IEgammaCalibration>::const_iterator
                      ical=m_calibsBarrel.begin();
                      ical != m_calibsBarrel.end(); ++ical )
              (*ical)->makeCorrection(ptrigEmCluster,caloDDE);
        }else{
        for( ToolHandleArray<IEgammaCalibration>::const_iterator
                    ical=m_calibsEndcap.begin();
                    ical != m_calibsEndcap.end(); ++ical )
            (*ical)->makeCorrection(ptrigEmCluster,caloDDE);
        }
      }
     
     float et_calib = -999.0;
     float et_uncalib = ptrigEmCluster->et();
     if ( m_doCalibWithRings ){
        ATH_CHECK(m_calibWRingsTool->checkRings(context));
        et_calib = m_calibWRingsTool->makeCalibWRings(context);
        ptrigEmCluster->setEt(et_calib);
 
     }
 
 
    float calZ0 = 0;
 
    // Print out Cluster produced
    if (msgLvl(MSG::DEBUG)) {
      ATH_MSG_DEBUG(" Values of Cluster produced: ");
      ATH_MSG_DEBUG(" REGTEST: emEnergy = " << ptrigEmCluster->energy());
      ATH_MSG_DEBUG(" REGTEST: hadEnergy = " << ptrigEmCluster->ehad1());
      ATH_MSG_DEBUG(" REGTEST: e237= " << ptrigEmCluster->e237());
      ATH_MSG_DEBUG(" REGTEST: e277= " << ptrigEmCluster->e277());
      ATH_MSG_DEBUG(" REGTEST: clusterWidth = " << ptrigEmCluster->weta2());
      ATH_MSG_DEBUG(" REGTEST: frac73 = " << ptrigEmCluster->fracs1());
      ATH_MSG_DEBUG(" REGTEST: e233 = " << ptrigEmCluster->e233());
      ATH_MSG_DEBUG(" REGTEST: wstot = " << ptrigEmCluster->wstot());
      ATH_MSG_DEBUG(" REGTEST: eta = " << ptrigEmCluster->eta());
      ATH_MSG_DEBUG(" REGTEST: phi = " << ptrigEmCluster->phi());
      ATH_MSG_DEBUG(" REGTEST: Eta1 = " << ptrigEmCluster->eta1());
      ATH_MSG_DEBUG(" REGTEST: calZ0 = " << calZ0);
      ATH_MSG_DEBUG(" REGTEST: quality = " << ptrigEmCluster->clusterQuality());
      ATH_MSG_DEBUG(std::hex << " REGTEST: roiWord = 0x" << ptrigEmCluster->RoIword()
                             << std::dec);
    }
    // my monitoring
    clET = et_uncalib*1e-3;
    clHET = ptrigEmCluster->ehad1()*1e-3;
    clEta = ptrigEmCluster->eta();
    clPhi = ptrigEmCluster->phi();
    clETrings = et_calib*1e-3;
   
    if (et_uncalib != 0.) {
      res_et = (et_calib - et_uncalib)/et_uncalib;
    }
 
    if ( ptrigEmCluster->e277() > 0.01 ) clReta = ptrigEmCluster->e237()/ptrigEmCluster->e277();
 
  } // end of roiCollection iterator
 
  return StatusCode::SUCCESS;
}