RoiUpdaterTool.cxx

Go to the documentation of this file.

 
/*
Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigSteeringEvent/TrigRoiDescriptorCollection.h"
#include "RoiUpdaterTool.h"
 
 
 
RoiUpdaterTool::RoiUpdaterTool(const std::string& type, const std::string& name, const IInterface* parent)
  : base_class(type, name, parent), m_update(false) 
{}
 
 
 
StatusCode RoiUpdaterTool::initialize()  {
 
  ATH_MSG_DEBUG("initialize()");
 
  if ( m_etaWidth!=-999 || m_phiWidth!=-999 || m_zedWidth!=-999 ) m_update = true;
 
  if ( m_useBeamspot ) { 
    m_update = true;
    ATH_CHECK( m_beamspotKey.initialize( m_useBeamspot ) );
  }
 
  return StatusCode::SUCCESS;
}
 
 
 
void zrange( double nsigma, double fence, double& zed, double& zedMinus, double& zedPlus, const SG::ReadCondHandleKey<InDet::BeamSpotData>& bs, const EventContext& ctx ) { 
 
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle( bs, ctx );
 
  uint32_t bitmap = beamSpotHandle->beamStatus();
  
  if ( ( (bitmap & 0x7) == 0x7 ) || !(bitmap & 0x4) ) { 
    double zsigma = beamSpotHandle->beamSigma(2);
    double zpos   = beamSpotHandle->beamPos()[2];
    
    zed      = zpos; 
    zedMinus = zpos - nsigma*zsigma - fence;
    zedPlus  = zpos + nsigma*zsigma + fence;
  }
  else { 
    zedMinus = -RoiDescriptor::zedWidthDefault();     
    zedPlus  =  RoiDescriptor::zedWidthDefault();   
  }
  
}
 
 
 
 
std::unique_ptr<TrigRoiDescriptor> RoiUpdaterTool::execute( const EventContext& ctx ) const {  
  RoiDescriptor troi( RoiDescriptor::FULLSCAN );
  return execute( &troi, ctx );
}
 
 
std::unique_ptr<TrigRoiDescriptor> RoiUpdaterTool::execute( const IRoiDescriptor* iroi, const EventContext& ctx ) const {  
 
  ATH_MSG_DEBUG( "execute()" );
  
  if ( !m_update ) return std::make_unique<TrigRoiDescriptor>( *iroi );
 
  if ( iroi->composite() )  { 
    ATH_MSG_WARNING( "Will not update composite Roi parameters: " << *iroi );
    return std::make_unique<TrigRoiDescriptor>( *iroi );
  }  
 
  double eta      = iroi->eta();
  double etaMinus = iroi->etaMinus();
  double etaPlus  = iroi->etaPlus();
  
  double phi      = iroi->phi();
  double phiMinus = iroi->phiMinus();
  double phiPlus  = iroi->phiPlus();
  
  double zed      = iroi->zed();
  double zedMinus = iroi->zedMinus();
  double zedPlus  = iroi->zedPlus();
  
  if ( m_etaWidth!=-999 ) { 
    etaMinus = eta - m_etaWidth;  
    etaPlus  = eta + m_etaWidth;  
  }
  
  if ( m_phiWidth!=-999 ) { 
    phiMinus = phi - m_phiWidth;  
    phiPlus  = phi + m_phiWidth;  
  }
  
  if ( m_zedWidth!=-999 ) { 
    zedMinus = zed - m_zedWidth;
    zedPlus  = zed + m_zedWidth;
  }
  
  if ( m_useBeamspot ) zrange( m_nsigma, m_fence, zed, zedMinus, zedPlus, m_beamspotKey, ctx );
  
  if ( zedMinus<-RoiDescriptor::zedWidthDefault() ) zedMinus = -RoiDescriptor::zedWidthDefault();     
  if (  zedPlus> RoiDescriptor::zedWidthDefault() ) zedPlus  =  RoiDescriptor::zedWidthDefault();  
  
  return std::make_unique<TrigRoiDescriptor>( eta, etaMinus, etaPlus, phi, phiMinus, phiPlus, zed, zedMinus, zedPlus );
 
}