L1CaloTriggerTowerSelector.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigT1CaloCalibUtils/L1CaloTriggerTowerSelector.h"
 
#include "TrigT1CaloEvent/TriggerTowerCollection.h"
 
L1CaloTriggerTowerSelector::L1CaloTriggerTowerSelector(const std::string& name, ISvcLocator* pSvcLocator): AthAlgorithm(name, pSvcLocator), m_towerKey(0)
{
    //setting selection range to the whole TTs by default
    m_vEtaSelectionRange.push_back(-5.);
    m_vEtaSelectionRange.push_back(5.);
 
    m_vPhiSelectionRange.push_back(0.);
    m_vPhiSelectionRange.push_back(7.);
 
    declareProperty("OutputTriggerTowerLocation", m_inputTriggerTowerLocation = "LVL1TriggerTowers");
    declareProperty("InputTriggerTowerLocation", m_outputTriggerTowerLocation = "SelectedLVL1TriggerTowers");
 
    declareProperty("EtaSelectionRange", m_vEtaSelectionRange);
    declareProperty("PhiSelectionRange", m_vPhiSelectionRange);
}
 
L1CaloTriggerTowerSelector::~L1CaloTriggerTowerSelector()
{
     delete m_towerKey;
}
 
//-------------------------------------------
// Initialize
//-------------------------------------------
StatusCode L1CaloTriggerTowerSelector::initialize()
{
    ATH_MSG_INFO("From Initialize...");
 
    m_towerKey = new LVL1::TriggerTowerKey();
 
    return StatusCode::SUCCESS;
}
 
//----------------------------------------
// Execute
//----------------------------------------
StatusCode L1CaloTriggerTowerSelector::execute(const EventContext& /*ctx*/)
{
    StatusCode sc;
 
    // retrieve triggertowers container from storegate
    const TriggerTowerCollection* ttCollection = 0;
    sc = evtStore()->retrieve(ttCollection, m_inputTriggerTowerLocation);
    if (sc.isFailure() || !ttCollection) {
        ATH_MSG_ERROR( "No Trigger Towers found" );
        return StatusCode::SUCCESS;
    } else ATH_MSG_DEBUG("TriggerTowerCollection retrieved from StoreGate, size: "<< ttCollection->size() );
 
    // create TriggerTowerCollection that will hold the selected TT
    TriggerTowerCollection* outputTTCollection = new TriggerTowerCollection();
 
    //---
    //loop on trigger towers
    TriggerTowerCollection::const_iterator pos = ttCollection->begin();
    TriggerTowerCollection::const_iterator pos_end = ttCollection->end();
 
        for(;pos!=pos_end;++pos) {
        const LVL1::TriggerTower* ptt = (*pos);
 
            // Compute TowerTrigger Key out of its coordinates
            double eta = ptt->eta();
        double phi = ptt->phi();
 
        if (this->isSelected(eta, phi)) {
            LVL1::TriggerTower* selectedTriggerTower = new LVL1::TriggerTower(*ptt);
            outputTTCollection->push_back(selectedTriggerTower);
        }
 
    } // -- End of trigger tower loop --
 
    sc = evtStore()->record(outputTTCollection, m_outputTriggerTowerLocation);
    if (sc.isFailure()) {
        ATH_MSG_ERROR( "Failed to record output TriggerTowerCollection to TDS with key: "<< m_outputTriggerTowerLocation );
        return StatusCode::SUCCESS;
    } else ATH_MSG_DEBUG("TriggerTowerCollection recorded to TDS with key: "<< m_outputTriggerTowerLocation<<", size: "<< outputTTCollection->size() );
 
    return StatusCode::SUCCESS;
}
 
//-----------------------------------------
// Finalize
//-----------------------------------------
StatusCode L1CaloTriggerTowerSelector::finalize()
{
    return StatusCode::SUCCESS;
}
 
bool L1CaloTriggerTowerSelector::isSelected(double eta, double phi) {
    return (eta<m_vEtaSelectionRange.at(1)&&eta>m_vEtaSelectionRange.at(0))&&(phi<m_vPhiSelectionRange.at(1)&&phi>m_vPhiSelectionRange.at(0));
}