TriggerRatesMonitor.cc

Go to the documentation of this file.

// Note to self: the implementation uses TH1F's to store the L1 and HLT rates.
// Assuming a maximum rate of 100 kHz times a period of 23.31 s, one needs to store counts up to ~2.3e6.
// A "float" has 24 bits of precision, so it can store up to 2**24 ~ 16.7e6 without loss of precision.


// C++ headers
#include <string>
#include <cstring>

// boost headers
#include <boost/regex.hpp>

// Root headers
#include <TH1F.h>

// CMSSW headers
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/Registry.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/Provenance/interface/ProcessHistory.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMaskAlgoTrigRcd.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMaskTechTrigRcd.h"
#include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
#include "CondFormats/L1TObjects/interface/L1GtTriggerMask.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/MonitorElement.h"

// length of a lumisections, corresponding to 2**18 LHC orbits, or 23.31 seconds
static const double SECS_PER_LUMI = 23.31040958083832;

// helper functions
template <typename T>
static 
const T * get(const edm::Event & event, const edm::EDGetTokenT<T> & token) {
  edm::Handle<T> handle;
  event.getByToken(token, handle);
  if (not handle.isValid())
    throw * handle.whyFailed();
  return handle.product();
}   
    
template <typename R, typename T>
static
const T * get(const edm::EventSetup & setup) {
  edm::ESHandle<T> handle;
  setup.get<R>().get(handle);
  return handle.product();
}


class TriggerRatesMonitor : public DQMEDAnalyzer {
public:
  explicit TriggerRatesMonitor(edm::ParameterSet const &);
  ~TriggerRatesMonitor();

  static void fillDescriptions(edm::ConfigurationDescriptions & descriptions);

private:
  virtual void dqmBeginRun(edm::Run const &, edm::EventSetup const &) override;
  virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  virtual void analyze(edm::Event const &, edm::EventSetup const &) override;

  // module configuration
  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord>    m_l1t_results;
  edm::EDGetTokenT<edm::TriggerResults>             m_hlt_results;
  std::string                                       m_dqm_path;
  uint32_t                                          m_lumisections_range;

  // L1T and HLT configuration
  struct HLTIndices {
    unsigned int index_l1_seed;
    unsigned int index_prescale;

    HLTIndices() :
      index_l1_seed(  (unsigned int) -1),
      index_prescale( (unsigned int) -1)
    { }
  };

  L1GtTriggerMenu const *       m_l1tMenu;
  L1GtTriggerMask const *       m_l1tAlgoMask;
  L1GtTriggerMask const *       m_l1tTechMask;
  HLTConfigProvider             m_hltConfig;
  std::vector<HLTIndices>       m_hltIndices;

  std::vector<std::vector<unsigned int>>    m_datasets;


  struct HLTRatesPlots {
    TH1F * wasrun;
    TH1F * pass_l1_seed;
    TH1F * pass_prescale;
    TH1F * accept;
    TH1F * reject;
    TH1F * error;
  };

  // overall event count and event types
  TH1F *                        m_events_processed;
  TH1F *                        m_events_physics;
  TH1F *                        m_events_calibration;
  TH1F *                        m_events_random;

  // L1T triggers
  std::vector<TH1F *>           m_l1t_algo_counts;
  std::vector<TH1F *>           m_l1t_tech_counts;

  // HLT triggers
  std::vector<HLTRatesPlots>    m_hlt_counts;

  // datasets
  std::vector<TH1F *>           m_dataset_counts;

};



void TriggerRatesMonitor::fillDescriptions(edm::ConfigurationDescriptions & descriptions)
{
  edm::ParameterSetDescription desc;
  desc.addUntracked<edm::InputTag>( "l1tResults",       edm::InputTag("gtDigis"));
  desc.addUntracked<edm::InputTag>( "hltResults",       edm::InputTag("TriggerResults"));
  desc.addUntracked<std::string>(   "dqmPath",          "HLT/TriggerRates" );
  desc.addUntracked<uint32_t>(      "lumisectionRange", 2500 );             // ~16 hours
  descriptions.add("triggerRatesMonitor", desc);
}


TriggerRatesMonitor::TriggerRatesMonitor(edm::ParameterSet const & config) :
  // module configuration
  m_l1t_results( consumes<L1GlobalTriggerReadoutRecord>( config.getUntrackedParameter<edm::InputTag>( "l1tResults" ) ) ),
  m_hlt_results( consumes<edm::TriggerResults>(          config.getUntrackedParameter<edm::InputTag>( "hltResults" ) ) ),
  m_dqm_path(                                            config.getUntrackedParameter<std::string>(   "dqmPath" ) ),
  m_lumisections_range(                                  config.getUntrackedParameter<uint32_t>(      "lumisectionRange" ) ),
  // L1T and HLT configuration
  m_l1tMenu( nullptr ),
  m_l1tAlgoMask( nullptr ),
  m_l1tTechMask( nullptr),
  m_hltConfig(),
  m_hltIndices(),
  m_datasets(),
  // overall event count and event types
  m_events_processed( nullptr ),
  m_events_physics( nullptr ),
  m_events_calibration( nullptr ),
  m_events_random( nullptr ),
  // L1T triggers
  m_l1t_algo_counts(),
  m_l1t_tech_counts(),
  // HLT triggers
  m_hlt_counts(),
  // datasets
  m_dataset_counts()
{
}

TriggerRatesMonitor::~TriggerRatesMonitor()
{
}

void TriggerRatesMonitor::dqmBeginRun(edm::Run const & run, edm::EventSetup const & setup)
{
  m_events_processed    = nullptr;
  m_events_physics      = nullptr;
  m_events_calibration  = nullptr;
  m_events_random       = nullptr;

  // cache the L1 trigger menu
  m_l1tMenu     = get<L1GtTriggerMenuRcd, L1GtTriggerMenu>(setup);
  // FIXME - do we really need this ?
  //(const_cast<L1GtTriggerMenu *>(m_l1tMenu))->buildGtConditionMap();
  m_l1tAlgoMask = get<L1GtTriggerMaskAlgoTrigRcd, L1GtTriggerMask>(setup);
  m_l1tTechMask = get<L1GtTriggerMaskTechTrigRcd, L1GtTriggerMask>(setup);
  if (m_l1tMenu and m_l1tAlgoMask and m_l1tTechMask) {
    m_l1t_algo_counts.clear();
    m_l1t_algo_counts.resize( m_l1tAlgoMask->gtTriggerMask().size(), nullptr );
    m_l1t_tech_counts.clear();
    m_l1t_tech_counts.resize( m_l1tTechMask->gtTriggerMask().size(), nullptr );
  } else {
    // L1GtUtils not initialised, skip the the L1T monitoring
    edm::LogError("TriggerRatesMonitor") << "failed to read the L1 menu or masks from the EventSetup, the L1 trigger rates will not be monitored";
  }

  // initialise the HLTConfigProvider
  bool changed = true;
  edm::EDConsumerBase::Labels labels;
  labelsForToken(m_hlt_results, labels);
  if (m_hltConfig.init(run, setup, labels.process, changed)) {
    m_hlt_counts.clear();
    m_hlt_counts.resize( m_hltConfig.size(), HLTRatesPlots() );
    m_hltIndices.resize( m_hltConfig.size(), HLTIndices() );

    unsigned int datasets = m_hltConfig.datasetNames().size();
    m_datasets.clear();
    m_datasets.resize( datasets, {} );
    for (unsigned int i = 0; i < datasets; ++i) {
      auto const & paths = m_hltConfig.datasetContent(i);
      m_datasets[i].reserve(paths.size());
      for (auto const & path: paths)
        m_datasets[i].push_back(m_hltConfig.triggerIndex(path));
    }
    m_dataset_counts.clear();
    m_dataset_counts.resize( datasets, nullptr );
  } else {
    // HLTConfigProvider not initialised, skip the the HLT monitoring
    edm::LogError("TriggerRatesMonitor") << "failed to initialise HLTConfigProvider, the HLT trigger and datasets rates will not be monitored";
  }
}

void TriggerRatesMonitor::bookHistograms(DQMStore::IBooker & booker, edm::Run const & run, edm::EventSetup const & setup)
{
  // book the overall event count and event types histograms
  booker.setCurrentFolder( m_dqm_path );
  m_events_processed    = booker.book1D("processed",    "Processed events",     m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
  m_events_physics      = booker.book1D("physics",      "Physics evenst",       m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
  m_events_calibration  = booker.book1D("calibration",  "Calibration events",   m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
  m_events_random       = booker.book1D("random",       "Random events",        m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
  
  if (m_l1tMenu and m_l1tAlgoMask and m_l1tTechMask) {
    // book the L1T rate histograms
    booker.setCurrentFolder( m_dqm_path + "/L1" );

    for (auto const & keyval: m_l1tMenu->gtAlgorithmAliasMap()) {
      int bit = keyval.second.algoBitNumber();
      std::string name = keyval.first.substr( 0, keyval.first.find_first_of(".") );
      // check if the trigger is unmasked in *any* partition
      if ((m_l1tAlgoMask->gtTriggerMask()[bit] & 0xff) != 0xff)
        m_l1t_algo_counts[bit] = booker.book1D(name, keyval.first,                  m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      else
        m_l1t_algo_counts[bit] = booker.book1D(name, keyval.first + " (masked)",    m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
    }

    for (auto const & keyval: m_l1tMenu->gtTechnicalTriggerMap()) {
      int bit = keyval.second.algoBitNumber();
      std::string name = keyval.first.substr( 0, keyval.first.find_first_of(".") );
      // check if the trigger is unmasked in *any* partition
      if ((m_l1tTechMask->gtTriggerMask()[bit] & 0xff) != 0xff)
        m_l1t_tech_counts[bit] = booker.book1D(name, keyval.first,                  m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      else
        m_l1t_tech_counts[bit] = booker.book1D(name, keyval.first + " (masked)",    m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
    }
  }

  if (m_hltConfig.inited()) {
    // book the HLT triggers rate histograms
    booker.setCurrentFolder( m_dqm_path + "/HLT" );
    for (unsigned int i = 0; i < m_hltConfig.size(); ++i) {
      std::string const & name = m_hltConfig.triggerName(i);
      m_hlt_counts[i].wasrun        = booker.book1D(name + "_wasrun",           name + " counts",           m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      m_hlt_counts[i].pass_l1_seed  = booker.book1D(name + "_pass_l1_seed",     name + " pass L1 seed",     m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      m_hlt_counts[i].pass_prescale = booker.book1D(name + "_pass_prescale",    name + " pass prescaler",   m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      m_hlt_counts[i].accept        = booker.book1D(name + "_accept",           name + " trigger",          m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      m_hlt_counts[i].reject        = booker.book1D(name + "_reject",           name + " reject",           m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      m_hlt_counts[i].error         = booker.book1D(name + "_error",            name + " error count",      m_lumisections_range + 1,   -0.5,   m_lumisections_range + 0.5)->getTH1F();
      // look for the index of the (last) L1 seed and prescale module in each path
      m_hltIndices[i].index_l1_seed  = m_hltConfig.size(i);
      m_hltIndices[i].index_prescale = m_hltConfig.size(i);
      for (unsigned int j = 0; j < m_hltConfig.size(i); ++j) {
        std::string const & label = m_hltConfig.moduleLabel(i, j);
        std::string const & type  = m_hltConfig.moduleType(label);
        if (type == "HLTLevel1GTSeed" or type == "HLTLevel1Activity" or type == "HLTLevel1Pattern") {
          // there might be more L1 seed filters in sequence
          // keep looking and store the index of the last one
          m_hltIndices[i].index_l1_seed  = j;
        } else if (type == "HLTPrescaler") {
          // there should be only one prescaler in a path, and it should follow all L1 seed filters
          m_hltIndices[i].index_prescale = j;
          break;
        }

      }
    }

    // book the HLT datasets rate histograms
    booker.setCurrentFolder( m_dqm_path + "/Datasets" );
    auto const & datasets = m_hltConfig.datasetNames();
    for (unsigned int i = 0; i < datasets.size(); ++i)
      m_dataset_counts[i] = booker.book1D(datasets[i], datasets[i], m_lumisections_range + 1, -0.5, m_lumisections_range + 0.5)->getTH1F();
  }
}


void TriggerRatesMonitor::analyze(edm::Event const & event, edm::EventSetup const & setup)
{
  unsigned int lumisection = event.luminosityBlock();

  // book the overall event count and event types rates
  m_events_processed->Fill(lumisection);
  switch (event.experimentType()) {
    case edm::EventAuxiliary::PhysicsTrigger :
      m_events_physics->Fill(lumisection);
      break;
    case edm::EventAuxiliary::CalibrationTrigger :
      m_events_calibration->Fill(lumisection);
      break;
    case edm::EventAuxiliary::RandomTrigger :
      m_events_random->Fill(lumisection);
      break;
    case edm::EventAuxiliary::Undefined :
    case edm::EventAuxiliary::Reserved :
    case edm::EventAuxiliary::TracedEvent :
    case edm::EventAuxiliary::TestTrigger :
    case edm::EventAuxiliary::ErrorTrigger :
      // ignore these event types
      break;
  }

  // monitor the L1 triggers rates
  if (m_l1tMenu and m_l1tAlgoMask and m_l1tTechMask) {
    L1GlobalTriggerReadoutRecord const & l1tResults = * get<L1GlobalTriggerReadoutRecord>(event, m_l1t_results);

    const std::vector<bool> & algoword = l1tResults.decisionWord();
    assert(algoword.size() == m_l1t_algo_counts.size());
    for (unsigned int i = 0; i < m_l1t_algo_counts.size(); ++i)
      if (algoword[i])
        m_l1t_algo_counts[i]->Fill(lumisection);

    const std::vector<bool> & techword = l1tResults.technicalTriggerWord();
    assert(techword.size() == m_l1t_tech_counts.size());
    for (unsigned int i = 0; i < m_l1t_tech_counts.size(); ++i)
      if (techword[i])
        m_l1t_tech_counts[i]->Fill(lumisection);
  }

  // monitor the HLT triggers and datsets rates
  if (m_hltConfig.inited()) {
    edm::TriggerResults const & hltResults = * get<edm::TriggerResults>(event, m_hlt_results);
    assert(hltResults.size() == m_hlt_counts.size());
    for (unsigned int i = 0; i < m_hlt_counts.size(); ++i) {
      edm::HLTPathStatus const & path = hltResults.at(i);
      if (path.wasrun())
        m_hlt_counts[i].wasrun->Fill(lumisection);
      if (path.index() > m_hltIndices[i].index_l1_seed)
        m_hlt_counts[i].pass_l1_seed->Fill(lumisection);
      if  (path.index() > m_hltIndices[i].index_prescale)
        m_hlt_counts[i].pass_prescale->Fill(lumisection);
      if (path.accept())
        m_hlt_counts[i].accept->Fill(lumisection);
      else if (path.error())
        m_hlt_counts[i].error ->Fill(lumisection);
      else
        m_hlt_counts[i].reject->Fill(lumisection);
    }

    for (unsigned int i = 0; i < m_datasets.size(); ++i)
      for (unsigned int j: m_datasets[i])
        if (hltResults.at(j).accept()) {
          m_dataset_counts[i]->Fill(lumisection);
          // ensure each dataset is incremented only once per event 
          break;
        }
  }
}


//define this as a plug-in
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(TriggerRatesMonitor);