TrigTLAMonitorAlgorithm.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRIGTLAMONITORING_TRIGTLAMONITORALGORITHM_H
#define TRIGTLAMONITORING_TRIGTLAMONITORALGORITHM_H
 
#include "AthenaMonitoring/AthMonitorAlgorithm.h"
#include "AthenaMonitoringKernel/Monitored.h"
 
#include "TrigDecisionTool/TrigDecisionTool.h"
 
#include "xAODJet/JetContainer.h"
#include "xAODMuon/MuonContainer.h"
#include "xAODEgamma/PhotonContainer.h"
#include "xAODTracking/TrackParticleContainer.h"
#include "xAODTrigger/TrigCompositeContainer.h"
 
class TrigTLAMonitorAlgorithm : public AthMonitorAlgorithm {
 public:
  TrigTLAMonitorAlgorithm( const std::string& name, ISvcLocator* pSvcLocator );
  virtual ~TrigTLAMonitorAlgorithm();
  virtual StatusCode initialize() override;
  virtual StatusCode fillHistograms( const EventContext& ctx ) const override;
 
 private:
  Gaudi::Property<std::vector<std::string>> m_allChains{this,"AllChains",{}};
 
  SG::ReadHandleKey<xAOD::JetContainer   > m_jetContainerKey   {this,"JetContainerName"   ,"HLT_AntiKt4EMTopoJets_subjesIS_TLA"           , "Jet Container Name"};
  SG::ReadHandleKey<xAOD::JetContainer   > m_pfjetContainerKey {this,"PFJetContainerName" ,"HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf_TLA", "Particle Flow Jet Container Name"};
  SG::ReadHandleKey<xAOD::PhotonContainer> m_photonContainerKey{this,"PhotonContainerName","HLT_egamma_Photons_TLA"                       , "Photon Container Name"};
  SG::ReadHandleKey<xAOD::MuonContainer  > m_muonContainerKey  {this,"MuonContainerName"  ,"HLT_MuonsCB_RoI_TLA"                          , "Muon Container Name"};
  SG::ReadHandleKey<xAOD::TrackParticleContainer> m_trackParticleContainerKey{this,"TrackContainerName"      ,"HLT_IDTrack_FS_FTF"        , "Track Container Name"};
  SG::ReadHandleKey<xAOD::TrigCompositeContainer> m_tcEventInfoContainerKey  {this,"TCEventInfoContainerName","HLT_TCEventInfo_TLA"       , "TCEventInfo Container Name"};
 
  PublicToolHandle<Trig::TrigDecisionTool> m_trigDecisionTool{this,"TriggerDecisionTool","Trigger decision tool"};
 
 
  template <typename T>
  StatusCode readContainer(SG::ReadHandle<T> & container, SG::ReadHandleKey<T> key, const EventContext& ctx) const;
 
 
  template <typename T>
  StatusCode fillEventInfoHistogram(SG::ReadHandle<DataVector<xAOD::TrigComposite_v1> > tcEventInfo, const std::string& varname, const std::string& prefix) const;
 
 
  template <class T>
  StatusCode fillParticleHistograms(SG::ReadHandle<DataVector<T>> particles, const std::string& prefix, const std::string& trigName) const;
 
 
  template<class C, typename T>
  StatusCode fillObjectVariableHistogram(SG::ReadHandle<DataVector<C>> container, const std::string& varname, const std::string& prefix , const std::string& trigName, T default_val=-1) const;
 
  //  The first element (corresponding to the hard scatter vertex) is filled.
  template <typename T>
  StatusCode fillJetTrackVariableHistogram(SG::ReadHandle<xAOD::JetContainer> jets,  const std::string& varname,  const std::string& prefix,  const std::string& trigName, T default_val=-1) const;
  
 
  template <typename U, typename T>
  StatusCode fill2DHistogram(const std::string& varname1, U var1, const std::string& varname2, T var2, const std::string& prefix) const;
 
 
  StatusCode fillJetPtCalibStatesHistograms(SG::ReadHandle<xAOD::JetContainer> jets, const std::string& calibState, const std::string& prefix,  const std::string& trigName) const;
 
 
  StatusCode fillDeltaRHistograms(const xAOD::IParticle* p0, const xAOD::IParticle* p1, const std::string& prefix, const std::string& trigName) const;
};
 
 
/*************************************************
* Implementation of templated methods            *
**************************************************/
 
template <typename T>
StatusCode TrigTLAMonitorAlgorithm::readContainer(SG::ReadHandle<T> & container, SG::ReadHandleKey<T> key, const EventContext& ctx) const{
 
  container = SG::makeHandle( key, ctx );
  if (! container.isValid() ) {
    ATH_MSG_ERROR("evtStore() does not contain Collection with name "<< key);
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}
 
 
template <typename T>
StatusCode TrigTLAMonitorAlgorithm::fillEventInfoHistogram(SG::ReadHandle<DataVector<xAOD::TrigComposite_v1> > tcEventInfo, const std::string& varname, const std::string& prefix) const {
 
  Monitored::Scalar<T> variable (prefix+"_"+varname, -1);
 
  for(auto tcEI : *tcEventInfo) {
    auto status = tcEI->getDetail<T>(varname, variable);
    if (!status) ATH_MSG_WARNING("No "<<varname<<" for this event");
    else {
      fill("TrigTLAMonitor", variable);
      ATH_MSG_DEBUG("Retrieved EventInfo variable: "<<varname<<" = "<<variable);
    }
  }
 
  return StatusCode::SUCCESS;
}
 
 
template <class T>
StatusCode TrigTLAMonitorAlgorithm::fillParticleHistograms(SG::ReadHandle<DataVector<T>> particles, const std::string& prefix, const std::string& trigName) const {
 
  // histograms
  Monitored::Scalar<int   > n   ("n"+prefix+"_"+trigName,0  );
  Monitored::Scalar<double> pt  (prefix+"pt_"  +trigName,0.0);
  Monitored::Scalar<double> eta (prefix+"eta_" +trigName,0.0);
  Monitored::Scalar<double> phi (prefix+"phi_" +trigName,0.0);
 
  Monitored::Scalar<double> pt0 (prefix+"0pt_" +trigName,0.0);
  Monitored::Scalar<double> eta0(prefix+"0eta_"+trigName,0.0);
  Monitored::Scalar<double> phi0(prefix+"0phi_"+trigName,0.0);
 
  // fill
  n = particles->size();
  fill("TrigTLAMonitor", n);
 
  for(const T* particle : *particles) {
    pt  = particle->pt()/1e3;
    eta = particle->eta();
    phi = particle->phi();
    fill("TrigTLAMonitor",pt );
    fill("TrigTLAMonitor",eta);
    fill("TrigTLAMonitor",phi);
  }
 
  // fill leading particle
  if(particles->size()>0) {
    pt0 =particles->at(0)->pt ()/1e3;
    eta0=particles->at(0)->eta();
    phi0=particles->at(0)->phi();
    fill("TrigTLAMonitor",pt0 );
    fill("TrigTLAMonitor",eta0);
    fill("TrigTLAMonitor",phi0);
  }
 
  return StatusCode::SUCCESS;
 
}
 
 
template<class C,typename T>
StatusCode TrigTLAMonitorAlgorithm::fillObjectVariableHistogram(SG::ReadHandle<DataVector<C>> container, const std::string& varname, const std::string& prefix , const std::string& trigName, T default_val) const {
  
  // change var name to avoid confilcts with Monitoring standards ("var_trigname")
  std::string hname = varname;
  std::replace( hname.begin(), hname.end(), '_', '-' );
 
  Monitored::Scalar<T> mon_var (prefix+hname+"_"+trigName,default_val);
  SG::AuxElement::ConstAccessor<T> accessor(varname);
  
  unsigned cnt(0);
  for(const auto element : *container) {
    mon_var = accessor(*element);
    fill("TrigTLAMonitor", mon_var);
    if (cnt < 5) ATH_MSG_DEBUG(prefix<<" "<<varname<<" = "<<mon_var);
    cnt++;
  }
  return StatusCode::SUCCESS;
}
 
 
template <typename T>
StatusCode TrigTLAMonitorAlgorithm::fillJetTrackVariableHistogram(SG::ReadHandle<xAOD::JetContainer> jets,  const std::string& varname,  const std::string& prefix,  const std::string& trigName, T default_val) const {
 
  std::vector<T> variable_vec;
  Monitored::Scalar<T> variable  (prefix+varname+"_"+trigName,default_val);
 
  unsigned cnt(0);
  for(auto jet : *jets) {
    auto status = jet->getAttribute<std::vector<T>>(varname, variable_vec);
    if (!status){
      ATH_MSG_WARNING("Failed retrieving "<<varname<<" for "<<prefix);
    }
    variable = variable_vec.at(0);
    fill("TrigTLAMonitor",variable);
    if (cnt < 3) ATH_MSG_DEBUG(prefix<<" "<<varname<<" = "<<variable);
    cnt++;
  }
 
  return StatusCode::SUCCESS;
}
 
 
template <typename U, typename T>
StatusCode TrigTLAMonitorAlgorithm::fill2DHistogram(const std::string& varname1, U var1, const std::string& varname2, T var2, const std::string& prefix) const{
 
  Monitored::Scalar<U> mon_var1(prefix+"_"+varname1, var1);
  Monitored::Scalar<T> mon_var2(prefix+"_"+varname2, var2);
  ATH_MSG_DEBUG("Filling 2D histogram. "<<varname1<<" = "<<var1<<" ; "<<varname2<<" = "<<var2);
  fill("TrigTLAMonitor", mon_var1, mon_var2);
 
  return StatusCode::SUCCESS;
}
 
 
#endif