EMTauInputProvider.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "./EMTauInputProvider.h"
 
#include <math.h>
 
#include "TrigT1CaloEvent/EmTauROI_ClassDEF.h"
 
#include "TrigT1Interfaces/CPRoIDecoder.h"
 
#include "L1TopoEvent/ClusterTOB.h"
#include "L1TopoEvent/TopoInputEvent.h"
 
using namespace std;
using namespace LVL1;
 
EMTauInputProvider::EMTauInputProvider(const std::string& type, const std::string& name, 
                                       const IInterface* parent) :
   base_class(type, name, parent),
   m_emTauLocation( TrigT1CaloDefs::EmTauTopoTobLocation )
{
   declareInterface<LVL1::IInputTOBConverter>( this );
   declareProperty( "EmTauROILocation", m_emTauLocation, "Storegate key for the EMTAU info from CMX" );
}
 
EMTauInputProvider::~EMTauInputProvider()
{}
 
StatusCode
EMTauInputProvider::initialize() {
 
   CHECK(m_emTauLocation.initialize());
 
   if (!m_monTool.empty()) ATH_CHECK(m_monTool.retrieve());
 
   return StatusCode::SUCCESS;
}
 
StatusCode
EMTauInputProvider::fillTopoInputEvent(TCS::TopoInputEvent& inputEvent) const {
 
   // Retrieve EMTAU RoIs (they are built by EMTAUTrigger)
 
   SG::ReadHandle<DataVector < CPCMXTopoData> > emtau(m_emTauLocation);
   if( !emtau.isValid() ) {
      ATH_MSG_WARNING("No CPCMXTopoDataCollection with SG key '" << m_emTauLocation.key() << "' found in the event. No EM or TAU input for the L1Topo simulation.");
      return StatusCode::RECOVERABLE;
   }
 
   ATH_MSG_DEBUG("Filling the input event. Number of emtau topo data objects: " << emtau->size());
   // topoData is read in in reverse in order to obtain the same crate order as in the hardware (3->0, not 0->3)
   for(auto iTopoData = emtau->rbegin(); iTopoData != emtau->rend(); ++iTopoData) {
      const CPCMXTopoData *topoData = *iTopoData;
 
      // fill the vector of TOBs
      std::vector< CPTopoTOB > tobs;
      topoData->tobs(tobs);
      ATH_MSG_DEBUG("Emtau topo data object has # TOBs: " << tobs.size());
      for(const CPTopoTOB & tob : tobs ) {
         ATH_MSG_DEBUG( "EMTAU TOB with cmx = " << tob.cmx() << "[" << (tob.cmx()==0?"EM":"TAU") << "]"
                        << " : e = " << setw(3) << tob.et() << ", isolation " << tob.isolation()
                        << ", eta = " << setw(2) << tob.eta() << ", phi = " << tob.phi()
                        << ", ieta = " << setw(2) << tob.ieta() << ", iphi = " << tob.iphi()
                        << ", word = " << hex << tob.roiWord() << dec
                        );
 
         TCS::ClusterTOB cl(tob.et(), tob.isolation(), tob.ieta(), tob.iphi(), tob.cmx()==0 ? TCS::CLUSTER : TCS::TAU, tob.roiWord() );
         cl.setEtaDouble( tob.eta() );
         cl.setPhiDouble( tob.phi() );
    
         if(tob.cmx()==0) {
            inputEvent.addCluster( cl );
            auto mon_hEMEt = Monitored::Scalar("EMTOBEt", cl.Et());
            auto mon_hEMPhi = Monitored::Scalar("EMTOBPhi", cl.phi());
            auto mon_hEMEta = Monitored::Scalar("EMTOBEta", cl.eta());  
            Monitored::Group(m_monTool, mon_hEMEt, mon_hEMPhi, mon_hEMEta);
         } else {
            inputEvent.addTau( cl );            
            auto mon_hTauEt = Monitored::Scalar("TauTOBEt", cl.Et());
            auto mon_hTauPhi = Monitored::Scalar("TauTOBPhi", cl.phi());
            auto mon_hTauEta = Monitored::Scalar("TauTOBEta", cl.eta());  
            Monitored::Group(m_monTool, mon_hTauEt, mon_hTauPhi, mon_hTauEta);
         }
      }
      if(topoData->overflow()){
          inputEvent.setOverflowFromEmtauInput(true);
          ATH_MSG_DEBUG("setOverflowFromEmtauInput : true");
      }
   }
   return StatusCode::SUCCESS;
}
 
void 
EMTauInputProvider::CalculateCoordinates(int32_t roiWord, double & eta, double & phi) const {
   CPRoIDecoder get;
   constexpr double TwoPI = 2 * M_PI;
   CoordinateRange coordRange = get.coordinate( roiWord );
   
   eta = coordRange.eta();
   phi = coordRange.phi();
   if( phi > M_PI ) phi -= TwoPI;
}