TrigmuRoI.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigmuRoI.h"
#include "TrigT1Result/RoIBResult.h"
#include "TrigT1Interfaces/RecMuonRoI.h"
#include "CxxUtils/phihelper.h"
 
// ================================================================================
// ================================================================================
 
TrigmuRoI::TrigmuRoI(const std::string& name, ISvcLocator* pSvcLocator)
  : AthReentrantAlgorithm(name, pSvcLocator)
{   
}
 
 
// ================================================================================
// ================================================================================
 
StatusCode TrigmuRoI::initialize()
{
 
  ATH_MSG_DEBUG("MinValueForOutOfTimeBC: "<<m_minValueForOutOfTimeBC);
  ATH_MSG_DEBUG("MaxValueForOutOfTimeBC: "<<m_maxValueForOutOfTimeBC);
 
  // Retrieve the tools/services
  ATH_CHECK(m_trigMuonRoITool.retrieve());
  ATH_CHECK(m_recRPCRoiTool.retrieve());
  ATH_CHECK(m_recTGCRoiTool.retrieve());
  ATH_CHECK(m_roisWriteHandleKey.initialize());
 
  if (!m_monTool.empty()) {
    ATH_MSG_DEBUG("Retrieving monTool");
    CHECK(m_monTool.retrieve());
  } else {
    ATH_MSG_INFO("No monTool configured => NO MONITORING");
  }
 
  return StatusCode::SUCCESS;
}
 
// ================================================================================
// ================================================================================
 
unsigned int TrigmuRoI::getBitMaskValue( const unsigned int uintValue, const unsigned int mask ) const {
  unsigned int result;
  unsigned int maskcopy;
  maskcopy = mask;
  result = uintValue & mask;
  if ( mask != 0 ) {
    while ( ( maskcopy & 0x00000001 ) == 0 ) {
      maskcopy = maskcopy >> 1;
      result = result >> 1;
    }
  }
  return result;
}
 
// ================================================================================
// ================================================================================
 
StatusCode TrigmuRoI::execute(const EventContext& ctx) const
{
 
  std::vector<int> outOfTimeRpc, outOfTimeTgc;
  std::vector<float> etaRoI, phiRoI;
  auto rpcOutOfTime = Monitored::Collection("RpcOutOfTime", outOfTimeRpc);
  auto tgcOutOfTime = Monitored::Collection("TgcOutOfTime", outOfTimeTgc);
  auto etaOutOfTimeRoI = Monitored::Collection("EtaOutOfTime", etaRoI);
  auto phiOutOfTimeRoI = Monitored::Collection("PhiOutOfTime", phiRoI);
 
  auto mon = Monitored::Group(m_monTool, rpcOutOfTime, tgcOutOfTime, etaOutOfTimeRoI, phiOutOfTimeRoI);
 
  int roi_id = 0;
 
  //--------------------------------------------------------------------------
  // Gather the Muon RoIs out of time by the
  //--------------------------------------------------------------------------
      
  SG::WriteHandle<TrigRoiDescriptorCollection> wh_roiCollection(m_roisWriteHandleKey, ctx);
  ATH_CHECK(wh_roiCollection.record(std::make_unique<TrigRoiDescriptorCollection>()));
  auto roiColl = wh_roiCollection.ptr();
 
  //get rois and loop over out of time rois
  auto roiVectors = m_trigMuonRoITool->decodeMuCTPi(ctx);
  if(!roiVectors){
    ATH_MSG_VERBOSE("No RoIs found");
    return StatusCode::SUCCESS;
  }
  for(auto it : roiVectors->outOfTimeRois){
 
    if (msgLvl(MSG::DEBUG)) {
      ATH_MSG_DEBUG(" Difference(RoI(BCID) - Event(BCID)) = " << (it).second);
      ATH_MSG_DEBUG(" ------------------------------------- ");
      ATH_MSG_DEBUG("RoIB word               : 0x" << MSG::hex << ((it).first).roIWord() << MSG::dec);
      ATH_MSG_DEBUG("Threshold               :  pt" << ((it).first).pt());
      ATH_MSG_DEBUG("Sector ID               :  " << ((it).first).getSectorID());
      ATH_MSG_DEBUG("Sector addr             :  0x" << MSG::hex << ((it).first).getSectorAddress() << MSG::dec);
      ATH_MSG_DEBUG("Sector overflow         :  " << ((it).first).getSectorOverflow());
      ATH_MSG_DEBUG("RoI overflow            :  " << ((it).first).getRoiOverflow());
      ATH_MSG_DEBUG("RoI number              :  " << ((it).first).getRoiNumber());
      ATH_MSG_DEBUG("IsHighestPt             :  " << ((it).first).getCandidateIsHighestPt());
      ATH_MSG_DEBUG("=================================================");
    }
            
    unsigned int temp_sysID = getBitMaskValue(((it).first).getSectorAddress(), LVL1::SysIDMask );
    unsigned int sysID = 0;                // Barrel
    if( temp_sysID & 0x2 ) sysID = 1;      // Endcap
    else if( temp_sysID & 0x1 ) sysID = 2; // Forward
 
 
    if ( sysID == 0 ) outOfTimeRpc.push_back((it).second);
    else              outOfTimeTgc.push_back((it).second);
    
    LVL1::TrigT1MuonRecRoiData roiData;
    //std::string region = "";
    if( sysID == 0 ) {
      ATH_CHECK( m_recRPCRoiTool->roiData(((it).first).roIWord(),roiData) );
      //region = "Barrel region";
    } else if ( sysID == 1 ){
      ATH_CHECK( m_recTGCRoiTool->roiData(((it).first).roIWord(),roiData) );
      //region = "Endcap region";
    } else {
      ATH_CHECK( m_recTGCRoiTool->roiData(((it).first).roIWord(),roiData) );
      //region = "Forward region";
    }
       
    // create new trigger element for this out of time RoI
    double eta = roiData.eta();
    double phi = roiData.phi();
    etaRoI.push_back(eta);
    phiRoI.push_back(phi);
     
    double etamin = eta - 0.2;
    double etamax = eta + 0.2;
    double phimin = CxxUtils::wrapToPi(phi - 0.2);
    double phimax = CxxUtils::wrapToPi(phi + 0.2);
 
    if ((it).second >= m_minValueForOutOfTimeBC &&
    (it).second <= m_maxValueForOutOfTimeBC    ) {
 
      // generic TrigRoiDescriptor
      auto roiDescriptor = new TrigRoiDescriptor( ((it).first).roIWord(), 0, roi_id, eta, etamin, etamax, phi, phimin, phimax,0,-255,255);
      roiColl->push_back(roiDescriptor);
 
      ATH_MSG_DEBUG("Created RoI descriptor with id, eta, phi: "<<roi_id<<" "<<eta<<" "<<phi);
    }
    roi_id += 1;
  }
  
  return StatusCode::SUCCESS;
}