Project CMSSW displayed by LXR CMSSW_12_5_X_2022-07-01-2300/​HLTrigger/​Timer/​plugins/​FastTimerService.cc	Source navigation Diff markup Identifier search General search
	Version: CMSSW_12_5_X_2022-07-01-2300 CMSSW_12_5_X_2022-06-30-2300 CMSSW_12_5_X_2022-06-29-2300 CMSSW_12_5_X_2022-06-28-2300 CMSSW_12_5_X_2022-06-27-2300 CMSSW_12_5_X_2022-06-26-2300 CMSSW_11_2_1 CMSSW_11_1_7 CMSSW_11_0_3 CMSSW_10_6_20 Revert   Change

File indexing completed on 2022-05-07 00:41:25

 // C++ headers
 #include <cmath>
 #include <fstream>
 #include <iomanip>
 #include <iostream>
 #include <limits>
 #include <mutex>
 #include <sstream>
 #include <string>
 #include <unordered_map>
 #include <unordered_set>
 
 // boost headers
 #include <boost/range/irange.hpp>
 
 // {fmt} headers
 #include <fmt/printf.h>
 
 // JSON headers
 #include <nlohmann/json.hpp>
 using json = nlohmann::json;
 
 // CMSSW headers
 #include "DQMServices/Core/interface/DQMStore.h"
 #include "DataFormats/Common/interface/HLTPathStatus.h"
 #include "DataFormats/Provenance/interface/EventID.h"
 #include "DataFormats/Provenance/interface/ModuleDescription.h"
 #include "DataFormats/Provenance/interface/Timestamp.h"
 #include "DataFormats/Scalers/interface/LumiScalers.h"
 #include "FWCore/Framework/interface/Event.h"
 #include "FWCore/Framework/interface/LuminosityBlock.h"
 #include "FWCore/Framework/interface/TriggerNamesService.h"
 #include "FWCore/ServiceRegistry/interface/PlaceInPathContext.h"
 #include "FWCore/MessageLogger/interface/MessageLogger.h"
 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
 #include "FWCore/ParameterSet/interface/ParameterSet.h"
 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
 #include "FWCore/Utilities/interface/Exception.h"
 #include "FWCore/Utilities/interface/StreamID.h"
 #include "HLTrigger/Timer/interface/memory_usage.h"
 #include "HLTrigger/Timer/interface/processor_model.h"
 #include "FastTimerService.h"
 
 using namespace std::literals;
 
 ///////////////////////////////////////////////////////////////////////////////
 
 namespace {
 
   // convert any std::chrono::duration to milliseconds
   template <class Rep, class Period>
   double ms(std::chrono::duration<Rep, Period> duration) {
     return std::chrono::duration_cast<std::chrono::duration<double, std::milli>>(duration).count();
   }
 
   // convert any boost::chrono::duration to milliseconds
   template <class Rep, class Period>
   double ms(boost::chrono::duration<Rep, Period> duration) {
     return boost::chrono::duration_cast<boost::chrono::duration<double, boost::milli>>(duration).count();
   }
 
   // convert a std::atomic<boost::chrono::nanoseconds::rep> to milliseconds
   double ms(std::atomic<boost::chrono::nanoseconds::rep> const& duration) {
     return boost::chrono::duration_cast<boost::chrono::duration<double, boost::milli>>(
                boost::chrono::nanoseconds(duration.load()))
         .count();
   }
 
   // convert from bytes to kilobytes, rounding down
   uint64_t kB(uint64_t bytes) { return bytes / 1024; }
 
   // convert from bytes to kilobytes, rounding down
   uint64_t kB(std::atomic<uint64_t> const& bytes) { return bytes.load() / 1024; }
 
 }  // namespace
 
 ///////////////////////////////////////////////////////////////////////////////
 
 // resources being monitored by the service
 
 // Resources
 
 FastTimerService::Resources::Resources()
     : time_thread(boost::chrono::nanoseconds::zero()),
       time_real(boost::chrono::nanoseconds::zero()),
       allocated(0ul),
       deallocated(0ul) {}
 
 void FastTimerService::Resources::reset() {
   time_thread = boost::chrono::nanoseconds::zero();
   time_real = boost::chrono::nanoseconds::zero();
   allocated = 0ul;
   deallocated = 0ul;
 }
 
 FastTimerService::Resources& FastTimerService::Resources::operator+=(Resources const& other) {
   time_thread += other.time_thread;
   time_real += other.time_real;
   allocated += other.allocated;
   deallocated += other.deallocated;
   return *this;
 }
 
 FastTimerService::Resources& FastTimerService::Resources::operator+=(AtomicResources const& other) {
   time_thread += boost::chrono::nanoseconds(other.time_thread.load());
   time_real += boost::chrono::nanoseconds(other.time_real.load());
   allocated += other.allocated.load();
   deallocated += other.deallocated.load();
   return *this;
 }
 
 FastTimerService::Resources FastTimerService::Resources::operator+(Resources const& other) const {
   Resources result(*this);
   result += other;
   return result;
 }
 
 FastTimerService::Resources FastTimerService::Resources::operator+(AtomicResources const& other) const {
   Resources result(*this);
   result += other;
   return result;
 }
 
 // AtomicResources
 // operation on the whole object are not atomic, as the operations
 // on the individual fields could be interleaved; however, accumulation
 // of results should yield the correct result.
 
 FastTimerService::AtomicResources::AtomicResources()
     : time_thread(0ul), time_real(0ul), allocated(0ul), deallocated(0ul) {}
 
 FastTimerService::AtomicResources::AtomicResources(AtomicResources const& other)
     : time_thread(other.time_thread.load()),
       time_real(other.time_real.load()),
       allocated(other.allocated.load()),
       deallocated(other.deallocated.load()) {}
 
 void FastTimerService::AtomicResources::reset() {
   time_thread = 0ul;
   time_real = 0ul;
   allocated = 0ul;
   deallocated = 0ul;
 }
 
 FastTimerService::AtomicResources& FastTimerService::AtomicResources::operator=(AtomicResources const& other) {
   time_thread = other.time_thread.load();
   time_real = other.time_real.load();
   allocated = other.allocated.load();
   deallocated = other.deallocated.load();
   return *this;
 }
 
 FastTimerService::AtomicResources& FastTimerService::AtomicResources::operator+=(AtomicResources const& other) {
   time_thread += other.time_thread.load();
   time_real += other.time_real.load();
   allocated += other.allocated.load();
   deallocated += other.deallocated.load();
   return *this;
 }
 
 FastTimerService::AtomicResources& FastTimerService::AtomicResources::operator+=(Resources const& other) {
   time_thread += other.time_thread.count();
   time_real += other.time_real.count();
   allocated += other.allocated;
   deallocated += other.deallocated;
   return *this;
 }
 
 FastTimerService::AtomicResources FastTimerService::AtomicResources::operator+(AtomicResources const& other) const {
   AtomicResources result(*this);
   result += other;
   return result;
 }
 
 FastTimerService::Resources FastTimerService::AtomicResources::operator+(Resources const& other) const {
   return other + *this;
 }
 
 // ResourcesPerModule
 
 FastTimerService::ResourcesPerModule::ResourcesPerModule() noexcept : total(), events(0), has_acquire(false) {}
 
 void FastTimerService::ResourcesPerModule::reset() noexcept {
   total.reset();
   events = 0;
   has_acquire = false;
 }
 
 FastTimerService::ResourcesPerModule& FastTimerService::ResourcesPerModule::operator+=(ResourcesPerModule const& other) {
   total += other.total;
   events += other.events;
   has_acquire = has_acquire or other.has_acquire;
   return *this;
 }
 
 FastTimerService::ResourcesPerModule FastTimerService::ResourcesPerModule::operator+(
     ResourcesPerModule const& other) const {
   ResourcesPerModule result(*this);
   result += other;
   return result;
 }
 
 // ResourcesPerPath
 
 void FastTimerService::ResourcesPerPath::reset() {
   active.reset();
   total.reset();
   last = 0;
   status = false;
 }
 
 FastTimerService::ResourcesPerPath& FastTimerService::ResourcesPerPath::operator+=(ResourcesPerPath const& other) {
   active += other.active;
   total += other.total;
   last = 0;  // summing these makes no sense, reset them instead
   status = false;
   return *this;
 }
 
 FastTimerService::ResourcesPerPath FastTimerService::ResourcesPerPath::operator+(ResourcesPerPath const& other) const {
   ResourcesPerPath result(*this);
   result += other;
   return result;
 }
 
 // ResourcesPerProcess
 
 FastTimerService::ResourcesPerProcess::ResourcesPerProcess(ProcessCallGraph::ProcessType const& process)
     : total(), paths(process.paths_.size()), endpaths(process.endPaths_.size()) {}
 
 void FastTimerService::ResourcesPerProcess::reset() {
   total.reset();
   for (auto& path : paths)
     path.reset();
   for (auto& path : endpaths)
     path.reset();
 }
 
 FastTimerService::ResourcesPerProcess& FastTimerService::ResourcesPerProcess::operator+=(
     ResourcesPerProcess const& other) {
   total += other.total;
   assert(paths.size() == other.paths.size());
   for (unsigned int i : boost::irange(0ul, paths.size()))
     paths[i] += other.paths[i];
   assert(endpaths.size() == other.endpaths.size());
   for (unsigned int i : boost::irange(0ul, endpaths.size()))
     endpaths[i] += other.endpaths[i];
   return *this;
 }
 
 FastTimerService::ResourcesPerProcess FastTimerService::ResourcesPerProcess::operator+(
     ResourcesPerProcess const& other) const {
   ResourcesPerProcess result(*this);
   result += other;
   return result;
 }
 
 // ResourcesPerJob
 
 FastTimerService::ResourcesPerJob::ResourcesPerJob(ProcessCallGraph const& job,
                                                    std::vector<GroupOfModules> const& groups)
     : total(), overhead(), eventsetup(), event(), highlight(groups.size()), modules(job.size()), processes(), events(0) {
   processes.reserve(job.processes().size());
   for (auto const& process : job.processes())
     processes.emplace_back(process);
 }
 
 void FastTimerService::ResourcesPerJob::reset() {
   total.reset();
   overhead.reset();
   eventsetup.reset();
   event.reset();
   for (auto& module : highlight)
     module.reset();
   for (auto& module : modules)
     module.reset();
   for (auto& process : processes)
     process.reset();
   events = 0;
 }
 
 FastTimerService::ResourcesPerJob& FastTimerService::ResourcesPerJob::operator+=(ResourcesPerJob const& other) {
   total += other.total;
   overhead += other.overhead;
   eventsetup += other.eventsetup;
   event += other.event;
   assert(highlight.size() == other.highlight.size());
   for (unsigned int i : boost::irange(0ul, highlight.size()))
     highlight[i] += other.highlight[i];
   assert(modules.size() == other.modules.size());
   for (unsigned int i : boost::irange(0ul, modules.size()))
     modules[i] += other.modules[i];
   assert(processes.size() == other.processes.size());
   for (unsigned int i : boost::irange(0ul, processes.size()))
     processes[i] += other.processes[i];
   events += other.events;
   return *this;
 }
 
 FastTimerService::ResourcesPerJob FastTimerService::ResourcesPerJob::operator+(ResourcesPerJob const& other) const {
   ResourcesPerJob result(*this);
   result += other;
   return result;
 }
 
 // per-thread measurements
 
 // Measurement
 
 FastTimerService::Measurement::Measurement() noexcept { measure(); }
 
 void FastTimerService::Measurement::measure() noexcept {
 #ifdef DEBUG_THREAD_CONCURRENCY
   id = std::this_thread::get_id();
 #endif  // DEBUG_THREAD_CONCURRENCY
   time_thread = boost::chrono::thread_clock::now();
   time_real = boost::chrono::high_resolution_clock::now();
   allocated = memory_usage::allocated();
   deallocated = memory_usage::deallocated();
 }
 
 void FastTimerService::Measurement::measure_and_store(Resources& store) noexcept {
 #ifdef DEBUG_THREAD_CONCURRENCY
   assert(std::this_thread::get_id() == id);
 #endif  // DEBUG_THREAD_CONCURRENCY
   auto new_time_thread = boost::chrono::thread_clock::now();
   auto new_time_real = boost::chrono::high_resolution_clock::now();
   auto new_allocated = memory_usage::allocated();
   auto new_deallocated = memory_usage::deallocated();
   store.time_thread = new_time_thread - time_thread;
   store.time_real = new_time_real - time_real;
   store.allocated = new_allocated - allocated;
   store.deallocated = new_deallocated - deallocated;
   time_thread = new_time_thread;
   time_real = new_time_real;
   allocated = new_allocated;
   deallocated = new_deallocated;
 }
 
 void FastTimerService::Measurement::measure_and_accumulate(Resources& store) noexcept {
 #ifdef DEBUG_THREAD_CONCURRENCY
   assert(std::this_thread::get_id() == id);
 #endif  // DEBUG_THREAD_CONCURRENCY
   auto new_time_thread = boost::chrono::thread_clock::now();
   auto new_time_real = boost::chrono::high_resolution_clock::now();
   auto new_allocated = memory_usage::allocated();
   auto new_deallocated = memory_usage::deallocated();
   store.time_thread += new_time_thread - time_thread;
   store.time_real += new_time_real - time_real;
   store.allocated += new_allocated - allocated;
   store.deallocated += new_deallocated - deallocated;
   time_thread = new_time_thread;
   time_real = new_time_real;
   allocated = new_allocated;
   deallocated = new_deallocated;
 }
 
 void FastTimerService::Measurement::measure_and_accumulate(AtomicResources& store) noexcept {
 #ifdef DEBUG_THREAD_CONCURRENCY
   assert(std::this_thread::get_id() == id);
 #endif  // DEBUG_THREAD_CONCURRENCY
   auto new_time_thread = boost::chrono::thread_clock::now();
   auto new_time_real = boost::chrono::high_resolution_clock::now();
   auto new_allocated = memory_usage::allocated();
   auto new_deallocated = memory_usage::deallocated();
   store.time_thread += boost::chrono::duration_cast<boost::chrono::nanoseconds>(new_time_thread - time_thread).count();
   store.time_real += boost::chrono::duration_cast<boost::chrono::nanoseconds>(new_time_real - time_real).count();
   store.allocated += new_allocated - allocated;
   store.deallocated += new_deallocated - deallocated;
   time_thread = new_time_thread;
   time_real = new_time_real;
   allocated = new_allocated;
   deallocated = new_deallocated;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void FastTimerService::PlotsPerElement::book(dqm::reco::DQMStore::IBooker& booker,
                                              std::string const& name,
                                              std::string const& title,
                                              PlotRanges const& ranges,
                                              unsigned int lumisections,
                                              bool byls) {
   int time_bins = (int)std::ceil(ranges.time_range / ranges.time_resolution);
   int mem_bins = (int)std::ceil(ranges.memory_range / ranges.memory_resolution);
   std::string y_title_ms = fmt::sprintf("events / %.1f ms", ranges.time_resolution);
   std::string y_title_kB = fmt::sprintf("events / %.1f kB", ranges.memory_resolution);
 
   // MonitorElement::setStatOverflows(kTRUE) includes underflows and overflows in the computation of mean and RMS
   time_thread_ =
       booker.book1D(name + " time_thread", title + " processing time (cpu)", time_bins, 0., ranges.time_range);
   time_thread_->setXTitle("processing time [ms]");
   time_thread_->setYTitle(y_title_ms);
   time_thread_->setStatOverflows(kTRUE);
 
   time_real_ = booker.book1D(name + " time_real", title + " processing time (real)", time_bins, 0., ranges.time_range);
   time_real_->setXTitle("processing time [ms]");
   time_real_->setYTitle(y_title_ms);
   time_real_->setStatOverflows(kTRUE);
 
   if (memory_usage::is_available()) {
     allocated_ = booker.book1D(name + " allocated", title + " allocated memory", mem_bins, 0., ranges.memory_range);
     allocated_->setXTitle("memory [kB]");
     allocated_->setYTitle(y_title_kB);
     allocated_->setStatOverflows(kTRUE);
 
     deallocated_ =
         booker.book1D(name + " deallocated", title + " deallocated memory", mem_bins, 0., ranges.memory_range);
     deallocated_->setXTitle("memory [kB]");
     deallocated_->setYTitle(y_title_kB);
     deallocated_->setStatOverflows(kTRUE);
   }
 
   if (not byls)
     return;
 
   time_thread_byls_ = booker.bookProfile(name + " time_thread_byls",
                                          title + " processing time (cpu) vs. lumisection",
                                          lumisections,
                                          0.5,
                                          lumisections + 0.5,
                                          time_bins,
                                          0.,
                                          std::numeric_limits<double>::infinity(),
                                          " ");
   time_thread_byls_->setXTitle("lumisection");
   time_thread_byls_->setYTitle("processing time [ms]");
   time_thread_byls_->setStatOverflows(kTRUE);
 
   time_real_byls_ = booker.bookProfile(name + " time_real_byls",
                                        title + " processing time (real) vs. lumisection",
                                        lumisections,
                                        0.5,
                                        lumisections + 0.5,
                                        time_bins,
                                        0.,
                                        std::numeric_limits<double>::infinity(),
                                        " ");
   time_real_byls_->setXTitle("lumisection");
   time_real_byls_->setYTitle("processing time [ms]");
   time_real_byls_->setStatOverflows(kTRUE);
 
   if (memory_usage::is_available()) {
     allocated_byls_ = booker.bookProfile(name + " allocated_byls",
                                          title + " allocated memory vs. lumisection",
                                          lumisections,
                                          0.5,
                                          lumisections + 0.5,
                                          mem_bins,
                                          0.,
                                          std::numeric_limits<double>::infinity(),
                                          " ");
     allocated_byls_->setXTitle("lumisection");
     allocated_byls_->setYTitle("memory [kB]");
     allocated_byls_->setStatOverflows(kTRUE);
 
     deallocated_byls_ = booker.bookProfile(name + " deallocated_byls",
                                            title + " deallocated memory vs. lumisection",
                                            lumisections,
                                            0.5,
                                            lumisections + 0.5,
                                            mem_bins,
                                            0.,
                                            std::numeric_limits<double>::infinity(),
                                            " ");
     deallocated_byls_->setXTitle("lumisection");
     deallocated_byls_->setYTitle("memory [kB]");
     deallocated_byls_->setStatOverflows(kTRUE);
   }
 }
 
 void FastTimerService::PlotsPerElement::fill(Resources const& data, unsigned int lumisection) {
   if (time_thread_)
     time_thread_->Fill(ms(data.time_thread));
 
   if (time_thread_byls_)
     time_thread_byls_->Fill(lumisection, ms(data.time_thread));
 
   if (time_real_)
     time_real_->Fill(ms(data.time_real));
 
   if (time_real_byls_)
     time_real_byls_->Fill(lumisection, ms(data.time_real));
 
   if (allocated_)
     allocated_->Fill(kB(data.allocated));
 
   if (allocated_byls_)
     allocated_byls_->Fill(lumisection, kB(data.allocated));
 
   if (deallocated_)
     deallocated_->Fill(kB(data.deallocated));
 
   if (deallocated_byls_)
     deallocated_byls_->Fill(lumisection, kB(data.deallocated));
 }
 
 void FastTimerService::PlotsPerElement::fill(AtomicResources const& data, unsigned int lumisection) {
   if (time_thread_)
     time_thread_->Fill(ms(boost::chrono::nanoseconds(data.time_thread.load())));
 
   if (time_thread_byls_)
     time_thread_byls_->Fill(lumisection, ms(boost::chrono::nanoseconds(data.time_thread.load())));
 
   if (time_real_)
     time_real_->Fill(ms(boost::chrono::nanoseconds(data.time_real.load())));
 
   if (time_real_byls_)
     time_real_byls_->Fill(lumisection, ms(boost::chrono::nanoseconds(data.time_real.load())));
 
   if (allocated_)
     allocated_->Fill(kB(data.allocated));
 
   if (allocated_byls_)
     allocated_byls_->Fill(lumisection, kB(data.allocated));
 
   if (deallocated_)
     deallocated_->Fill(kB(data.deallocated));
 
   if (deallocated_byls_)
     deallocated_byls_->Fill(lumisection, kB(data.deallocated));
 }
 
 void FastTimerService::PlotsPerElement::fill_fraction(Resources const& data,
                                                       Resources const& part,
                                                       unsigned int lumisection) {
   float total;
   float fraction;
 
   total = ms(data.time_thread);
   fraction = (total > 0.) ? (ms(part.time_thread) / total) : 0.;
   if (time_thread_)
     time_thread_->Fill(total, fraction);
 
   if (time_thread_byls_)
     time_thread_byls_->Fill(lumisection, total, fraction);
 
   total = ms(data.time_real);
   fraction = (total > 0.) ? (ms(part.time_real) / total) : 0.;
   if (time_real_)
     time_real_->Fill(total, fraction);
 
   if (time_real_byls_)
     time_real_byls_->Fill(lumisection, total, fraction);
 
   total = kB(data.allocated);
   fraction = (total > 0.) ? (kB(part.allocated) / total) : 0.;
   if (allocated_)
     allocated_->Fill(total, fraction);
 
   if (allocated_byls_)
     allocated_byls_->Fill(lumisection, total, fraction);
 
   total = kB(data.deallocated);
   fraction = (total > 0.) ? (kB(part.deallocated) / total) : 0.;
   if (deallocated_)
     deallocated_->Fill(total, fraction);
 
   if (deallocated_byls_)
     deallocated_byls_->Fill(lumisection, total, fraction);
 }
 
 void FastTimerService::PlotsPerPath::book(dqm::reco::DQMStore::IBooker& booker,
                                           std::string const& prefixDir,
                                           ProcessCallGraph const& job,
                                           ProcessCallGraph::PathType const& path,
                                           PlotRanges const& ranges,
                                           unsigned int lumisections,
                                           bool byls) {
   const std::string basedir = booker.pwd();
   booker.setCurrentFolder(basedir + "/" + prefixDir + path.name_);
 
   total_.book(booker, "path", path.name_, ranges, lumisections, byls);
 
   // MonitorElement::setStatOverflows(kTRUE) includes underflows and overflows in the computation of mean and RMS
   unsigned int bins = path.modules_and_dependencies_.size();
   module_counter_ = booker.book1DD("module_counter", "module counter", bins + 1, -0.5, bins + 0.5);
   module_counter_->setYTitle("events");
   module_counter_->setStatOverflows(kTRUE);
   module_time_thread_total_ =
       booker.book1DD("module_time_thread_total", "total module time (cpu)", bins, -0.5, bins - 0.5);
   module_time_thread_total_->setYTitle("processing time [ms]");
   module_time_thread_total_->setStatOverflows(kTRUE);
   module_time_real_total_ =
       booker.book1DD("module_time_real_total", "total module time (real)", bins, -0.5, bins - 0.5);
   module_time_real_total_->setYTitle("processing time [ms]");
   module_time_real_total_->setStatOverflows(kTRUE);
   if (memory_usage::is_available()) {
     module_allocated_total_ =
         booker.book1DD("module_allocated_total", "total allocated memory", bins, -0.5, bins - 0.5);
     module_allocated_total_->setYTitle("memory [kB]");
     module_allocated_total_->setStatOverflows(kTRUE);
     module_deallocated_total_ =
         booker.book1DD("module_deallocated_total", "total deallocated memory", bins, -0.5, bins - 0.5);
     module_deallocated_total_->setYTitle("memory [kB]");
     module_deallocated_total_->setStatOverflows(kTRUE);
   }
   for (unsigned int bin : boost::irange(0u, bins)) {
     auto const& module = job[path.modules_and_dependencies_[bin]];
     std::string const& label =
         module.scheduled_ ? module.module_.moduleLabel() : module.module_.moduleLabel() + " (unscheduled)";
     module_counter_->setBinLabel(bin + 1, label);
     module_time_thread_total_->setBinLabel(bin + 1, label);
     module_time_real_total_->setBinLabel(bin + 1, label);
     if (memory_usage::is_available()) {
       module_allocated_total_->setBinLabel(bin + 1, label);
       module_deallocated_total_->setBinLabel(bin + 1, label);
     }
   }
   module_counter_->setBinLabel(bins + 1, "");
 
   booker.setCurrentFolder(basedir);
 }
 
 void FastTimerService::PlotsPerPath::fill(ProcessCallGraph::PathType const& description,
                                           ResourcesPerJob const& data,
                                           ResourcesPerPath const& path,
                                           unsigned int ls) {
   // fill the total path time
   total_.fill(path.total, ls);
 
   // fill the modules that actually ran and the total time spent in each od them
   for (unsigned int i = 0; i < path.last; ++i) {
     auto const& module = data.modules[description.modules_and_dependencies_[i]];
     if (module_counter_)
       module_counter_->Fill(i);
 
     if (module_time_thread_total_)
       module_time_thread_total_->Fill(i, ms(module.total.time_thread));
 
     if (module_time_real_total_)
       module_time_real_total_->Fill(i, ms(module.total.time_real));
 
     if (module_allocated_total_)
       module_allocated_total_->Fill(i, kB(module.total.allocated));
 
     if (module_deallocated_total_)
       module_deallocated_total_->Fill(i, kB(module.total.deallocated));
   }
   if (module_counter_ and path.status)
     module_counter_->Fill(path.last);
 }
 
 FastTimerService::PlotsPerProcess::PlotsPerProcess(ProcessCallGraph::ProcessType const& process)
     : event_(), paths_(process.paths_.size()), endpaths_(process.endPaths_.size()) {}
 
 void FastTimerService::PlotsPerProcess::book(dqm::reco::DQMStore::IBooker& booker,
                                              ProcessCallGraph const& job,
                                              ProcessCallGraph::ProcessType const& process,
                                              PlotRanges const& event_ranges,
                                              PlotRanges const& path_ranges,
                                              unsigned int lumisections,
                                              bool bypath,
                                              bool byls) {
   const std::string basedir = booker.pwd();
   event_.book(booker, "process " + process.name_, "process " + process.name_, event_ranges, lumisections, byls);
   if (bypath) {
     booker.setCurrentFolder(basedir + "/process " + process.name_ + " paths");
     for (unsigned int id : boost::irange(0ul, paths_.size())) {
       paths_[id].book(booker, "path ", job, process.paths_[id], path_ranges, lumisections, byls);
     }
     for (unsigned int id : boost::irange(0ul, endpaths_.size())) {
       endpaths_[id].book(booker, "endpath ", job, process.endPaths_[id], path_ranges, lumisections, byls);
     }
     booker.setCurrentFolder(basedir);
   }
 }
 
 void FastTimerService::PlotsPerProcess::fill(ProcessCallGraph::ProcessType const& description,
                                              ResourcesPerJob const& data,
                                              ResourcesPerProcess const& process,
                                              unsigned int ls) {
   // fill process event plots
   event_.fill(process.total, ls);
 
   // fill all paths plots
   for (unsigned int id : boost::irange(0ul, paths_.size()))
     paths_[id].fill(description.paths_[id], data, process.paths[id], ls);
 
   // fill all endpaths plots
   for (unsigned int id : boost::irange(0ul, endpaths_.size()))
     endpaths_[id].fill(description.endPaths_[id], data, process.endpaths[id], ls);
 }
 
 FastTimerService::PlotsPerJob::PlotsPerJob(ProcessCallGraph const& job, std::vector<GroupOfModules> const& groups)
     : event_(), event_ex_(), overhead_(), highlight_(groups.size()), modules_(job.size()), processes_() {
   processes_.reserve(job.processes().size());
   for (auto const& process : job.processes())
     processes_.emplace_back(process);
 }
 
 void FastTimerService::PlotsPerJob::book(dqm::reco::DQMStore::IBooker& booker,
                                          ProcessCallGraph const& job,
                                          std::vector<GroupOfModules> const& groups,
                                          PlotRanges const& event_ranges,
                                          PlotRanges const& path_ranges,
                                          PlotRanges const& module_ranges,
                                          unsigned int lumisections,
                                          bool bymodule,
                                          bool bypath,
                                          bool byls,
                                          bool transitions) {
   const std::string basedir = booker.pwd();
 
   // event summary plots
   event_.book(booker, "event", "Event", event_ranges, lumisections, byls);
 
   event_ex_.book(booker, "explicit", "Event (explicit)", event_ranges, lumisections, byls);
 
   overhead_.book(booker, "overhead", "Overhead", event_ranges, lumisections, byls);
 
   modules_[job.source().id()].book(booker, "source", "Source", module_ranges, lumisections, byls);
 
   if (transitions) {
     lumi_.book(booker, "lumi", "LumiSection transitions", event_ranges, lumisections, byls);
 
     run_.book(booker, "run", "Run transtions", event_ranges, lumisections, false);
   }
 
   // plot the time spent in few given groups of modules
   for (unsigned int group : boost::irange(0ul, groups.size())) {
     auto const& label = groups[group].label;
     highlight_[group].book(booker, "highlight " + label, "Highlight " + label, event_ranges, lumisections, byls);
   }
 
   // plots per subprocess (event, modules, paths and endpaths)
   for (unsigned int pid : boost::irange(0ul, job.processes().size())) {
     auto const& process = job.processDescription(pid);
     processes_[pid].book(booker, job, process, event_ranges, path_ranges, lumisections, bypath, byls);
 
     if (bymodule) {
       booker.setCurrentFolder(basedir + "/process " + process.name_ + " modules");
       for (unsigned int id : process.modules_) {
         auto const& module_name = job.module(id).moduleLabel();
         modules_[id].book(booker, module_name, module_name, module_ranges, lumisections, byls);
       }
       booker.setCurrentFolder(basedir);
     }
   }
 }
 
 void FastTimerService::PlotsPerJob::fill(ProcessCallGraph const& job, ResourcesPerJob const& data, unsigned int ls) {
   // fill total event plots
   event_.fill(data.total, ls);
   event_ex_.fill(data.event, ls);
   overhead_.fill(data.overhead, ls);
 
   // fill highltight plots
   for (unsigned int group : boost::irange(0ul, highlight_.size()))
     highlight_[group].fill_fraction(data.total, data.highlight[group], ls);
 
   // fill modules plots
   for (unsigned int id : boost::irange(0ul, modules_.size()))
     modules_[id].fill(data.modules[id].total, ls);
 
   for (unsigned int pid : boost::irange(0ul, processes_.size()))
     processes_[pid].fill(job.processDescription(pid), data, data.processes[pid], ls);
 }
 
 void FastTimerService::PlotsPerJob::fill_run(AtomicResources const& data) {
   // fill run transition plots
   run_.fill(data, 0);
 }
 
 void FastTimerService::PlotsPerJob::fill_lumi(AtomicResources const& data, unsigned int ls) {
   // fill lumisection transition plots
   lumi_.fill(data, ls);
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 FastTimerService::FastTimerService(const edm::ParameterSet& config, edm::ActivityRegistry& registry)
     :  // configuration
       callgraph_(),
       // job configuration
       concurrent_lumis_(0),
       concurrent_runs_(0),
       concurrent_streams_(0),
       concurrent_threads_(0),
       print_event_summary_(config.getUntrackedParameter<bool>("printEventSummary")),
       print_run_summary_(config.getUntrackedParameter<bool>("printRunSummary")),
       print_job_summary_(config.getUntrackedParameter<bool>("printJobSummary")),
       // JSON configuration
       //write_json_per_event_(config.getUntrackedParameter<bool>("writeJSONByEvent")),
       //write_json_per_ls_(config.getUntrackedParameter<bool>("writeJSONByLumiSection")),
       //write_json_per_run_(config.getUntrackedParameter<bool>("writeJSONByRun")),
       write_json_summary_(config.getUntrackedParameter<bool>("writeJSONSummary")),
       json_filename_(config.getUntrackedParameter<std::string>("jsonFileName")),
       // DQM configuration
       enable_dqm_(config.getUntrackedParameter<bool>("enableDQM")),
       enable_dqm_bymodule_(config.getUntrackedParameter<bool>("enableDQMbyModule")),
       enable_dqm_bypath_(config.getUntrackedParameter<bool>("enableDQMbyPath")),
       enable_dqm_byls_(config.getUntrackedParameter<bool>("enableDQMbyLumiSection")),
       enable_dqm_bynproc_(config.getUntrackedParameter<bool>("enableDQMbyProcesses")),
       enable_dqm_transitions_(config.getUntrackedParameter<bool>("enableDQMTransitions")),
       dqm_event_ranges_({config.getUntrackedParameter<double>("dqmTimeRange"),                  // ms
                          config.getUntrackedParameter<double>("dqmTimeResolution"),             // ms
                          config.getUntrackedParameter<double>("dqmMemoryRange"),                // kB
                          config.getUntrackedParameter<double>("dqmMemoryResolution")}),         // kB
       dqm_path_ranges_({config.getUntrackedParameter<double>("dqmPathTimeRange"),               // ms
                         config.getUntrackedParameter<double>("dqmPathTimeResolution"),          // ms
                         config.getUntrackedParameter<double>("dqmPathMemoryRange"),             // kB
                         config.getUntrackedParameter<double>("dqmPathMemoryResolution")}),      // kB
       dqm_module_ranges_({config.getUntrackedParameter<double>("dqmModuleTimeRange"),           // ms
                           config.getUntrackedParameter<double>("dqmModuleTimeResolution"),      // ms
                           config.getUntrackedParameter<double>("dqmModuleMemoryRange"),         // kB
                           config.getUntrackedParameter<double>("dqmModuleMemoryResolution")}),  // kB
       dqm_lumisections_range_(config.getUntrackedParameter<unsigned int>("dqmLumiSectionsRange")),
       dqm_path_(config.getUntrackedParameter<std::string>("dqmPath")),
       // highlight configuration
       highlight_module_psets_(config.getUntrackedParameter<std::vector<edm::ParameterSet>>("highlightModules")),
       highlight_modules_(highlight_module_psets_.size())  // filled in postBeginJob()
 {
   // start observing when a thread enters or leaves the TBB global thread arena
   tbb::task_scheduler_observer::observe();
 
   // register EDM call backs
   registry.watchPreallocate(this, &FastTimerService::preallocate);
   registry.watchPreBeginJob(this, &FastTimerService::preBeginJob);
   registry.watchPostBeginJob(this, &FastTimerService::postBeginJob);
   registry.watchPostEndJob(this, &FastTimerService::postEndJob);
   registry.watchPreGlobalBeginRun(this, &FastTimerService::preGlobalBeginRun);
   //registry.watchPostGlobalBeginRun(         this, & FastTimerService::postGlobalBeginRun );
   //registry.watchPreGlobalEndRun(            this, & FastTimerService::preGlobalEndRun );
   registry.watchPostGlobalEndRun(this, &FastTimerService::postGlobalEndRun);
   registry.watchPreStreamBeginRun(this, &FastTimerService::preStreamBeginRun);
   //registry.watchPostStreamBeginRun(         this, & FastTimerService::postStreamBeginRun );
   //registry.watchPreStreamEndRun(            this, & FastTimerService::preStreamEndRun );
   registry.watchPostStreamEndRun(this, &FastTimerService::postStreamEndRun);
   registry.watchPreGlobalBeginLumi(this, &FastTimerService::preGlobalBeginLumi);
   //registry.watchPostGlobalBeginLumi(        this, & FastTimerService::postGlobalBeginLumi );
   //registry.watchPreGlobalEndLumi(           this, & FastTimerService::preGlobalEndLumi );
   registry.watchPostGlobalEndLumi(this, &FastTimerService::postGlobalEndLumi);
   registry.watchPreStreamBeginLumi(this, &FastTimerService::preStreamBeginLumi);
   //registry.watchPostStreamBeginLumi(        this, & FastTimerService::postStreamBeginLumi );
   //registry.watchPreStreamEndLumi(           this, & FastTimerService::preStreamEndLumi );
   registry.watchPostStreamEndLumi(this, &FastTimerService::postStreamEndLumi);
   registry.watchPreEvent(this, &FastTimerService::preEvent);
   registry.watchPostEvent(this, &FastTimerService::postEvent);
   registry.watchPrePathEvent(this, &FastTimerService::prePathEvent);
   registry.watchPostPathEvent(this, &FastTimerService::postPathEvent);
   registry.watchPreSourceConstruction(this, &FastTimerService::preSourceConstruction);
   //registry.watchPostSourceConstruction(     this, & FastTimerService::postSourceConstruction);
   registry.watchPreSourceRun(this, &FastTimerService::preSourceRun);
   registry.watchPostSourceRun(this, &FastTimerService::postSourceRun);
   registry.watchPreSourceLumi(this, &FastTimerService::preSourceLumi);
   registry.watchPostSourceLumi(this, &FastTimerService::postSourceLumi);
   registry.watchPreSourceEvent(this, &FastTimerService::preSourceEvent);
   registry.watchPostSourceEvent(this, &FastTimerService::postSourceEvent);
   //registry.watchPreModuleConstruction(      this, & FastTimerService::preModuleConstruction);
   //registry.watchPostModuleConstruction(     this, & FastTimerService::postModuleConstruction);
   //registry.watchPreModuleBeginJob(          this, & FastTimerService::preModuleBeginJob );
   //registry.watchPostModuleBeginJob(         this, & FastTimerService::postModuleBeginJob );
   //registry.watchPreModuleEndJob(            this, & FastTimerService::preModuleEndJob );
   //registry.watchPostModuleEndJob(           this, & FastTimerService::postModuleEndJob );
   //registry.watchPreModuleBeginStream(       this, & FastTimerService::preModuleBeginStream );
   //registry.watchPostModuleBeginStream(      this, & FastTimerService::postModuleBeginStream );
   //registry.watchPreModuleEndStream(         this, & FastTimerService::preModuleEndStream );
   //registry.watchPostModuleEndStream(        this, & FastTimerService::postModuleEndStream );
   registry.watchPreModuleGlobalBeginRun(this, &FastTimerService::preModuleGlobalBeginRun);
   registry.watchPostModuleGlobalBeginRun(this, &FastTimerService::postModuleGlobalBeginRun);
   registry.watchPreModuleGlobalEndRun(this, &FastTimerService::preModuleGlobalEndRun);
   registry.watchPostModuleGlobalEndRun(this, &FastTimerService::postModuleGlobalEndRun);
   registry.watchPreModuleGlobalBeginLumi(this, &FastTimerService::preModuleGlobalBeginLumi);
   registry.watchPostModuleGlobalBeginLumi(this, &FastTimerService::postModuleGlobalBeginLumi);
   registry.watchPreModuleGlobalEndLumi(this, &FastTimerService::preModuleGlobalEndLumi);
   registry.watchPostModuleGlobalEndLumi(this, &FastTimerService::postModuleGlobalEndLumi);
   registry.watchPreModuleStreamBeginRun(this, &FastTimerService::preModuleStreamBeginRun);
   registry.watchPostModuleStreamBeginRun(this, &FastTimerService::postModuleStreamBeginRun);
   registry.watchPreModuleStreamEndRun(this, &FastTimerService::preModuleStreamEndRun);
   registry.watchPostModuleStreamEndRun(this, &FastTimerService::postModuleStreamEndRun);
   registry.watchPreModuleStreamBeginLumi(this, &FastTimerService::preModuleStreamBeginLumi);
   registry.watchPostModuleStreamBeginLumi(this, &FastTimerService::postModuleStreamBeginLumi);
   registry.watchPreModuleStreamEndLumi(this, &FastTimerService::preModuleStreamEndLumi);
   registry.watchPostModuleStreamEndLumi(this, &FastTimerService::postModuleStreamEndLumi);
   //registry.watchPreModuleEventPrefetching(  this, & FastTimerService::preModuleEventPrefetching );
   //registry.watchPostModuleEventPrefetching( this, & FastTimerService::postModuleEventPrefetching );
   registry.watchPreModuleEventAcquire(this, &FastTimerService::preModuleEventAcquire);
   registry.watchPostModuleEventAcquire(this, &FastTimerService::postModuleEventAcquire);
   registry.watchPreModuleEvent(this, &FastTimerService::preModuleEvent);
   registry.watchPostModuleEvent(this, &FastTimerService::postModuleEvent);
   registry.watchPreModuleEventDelayedGet(this, &FastTimerService::preModuleEventDelayedGet);
   registry.watchPostModuleEventDelayedGet(this, &FastTimerService::postModuleEventDelayedGet);
   registry.watchPreEventReadFromSource(this, &FastTimerService::preEventReadFromSource);
   registry.watchPostEventReadFromSource(this, &FastTimerService::postEventReadFromSource);
   registry.watchPreESModule(this, &FastTimerService::preESModule);
   registry.watchPostESModule(this, &FastTimerService::postESModule);
 }
 
 void FastTimerService::ignoredSignal(const std::string& signal) const {
   LogDebug("FastTimerService") << "The FastTimerService received is currently not monitoring the signal \"" << signal
                                << "\".\n";
 }
 
 void FastTimerService::unsupportedSignal(const std::string& signal) const {
   // warn about each signal only once per job
   if (unsupported_signals_.insert(signal).second)
     edm::LogWarning("FastTimerService") << "The FastTimerService received the unsupported signal \"" << signal
                                         << "\".\n"
                                         << "Please report how to reproduce the issue to cms-hlt@cern.ch .";
 }
 
 void FastTimerService::preGlobalBeginRun(edm::GlobalContext const& gc) {
   ignoredSignal(__func__);
 
   // reset the run counters only during the main process being run
   if (isFirstSubprocess(gc)) {
     auto index = gc.runIndex();
     subprocess_global_run_check_[index] = 0;
     run_transition_[index].reset();
     run_summary_[index].reset();
 
     // book the DQM plots
     if (enable_dqm_) {
       // define a callback to book the MonitorElements
       auto bookTransactionCallback = [&, this](dqm::reco::DQMStore::IBooker& booker, dqm::reco::DQMStore::IGetter&) {
         auto scope = dqm::reco::DQMStore::IBooker::UseRunScope(booker);
         // we should really do this, but only DQMStore is allowed to touch it
         // We could move to postGlobalBeginRun, then the DQMStore has sure set it up.
         //booker.setRunLumi(gc.luminosityBlockID());
         booker.setCurrentFolder(dqm_path_);
         plots_->book(booker,
                      callgraph_,
                      highlight_modules_,
                      dqm_event_ranges_,
                      dqm_path_ranges_,
                      dqm_module_ranges_,
                      dqm_lumisections_range_,
                      enable_dqm_bymodule_,
                      enable_dqm_bypath_,
                      enable_dqm_byls_,
                      enable_dqm_transitions_);
       };
 
       // book MonitorElements for this stream
       edm::Service<dqm::legacy::DQMStore>()->meBookerGetter(bookTransactionCallback);
     }
   }
 }
 
 void FastTimerService::postGlobalBeginRun(edm::GlobalContext const& gc) { ignoredSignal(__func__); }
 
 void FastTimerService::preStreamBeginRun(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::preallocate(edm::service::SystemBounds const& bounds) {
   concurrent_lumis_ = bounds.maxNumberOfConcurrentLuminosityBlocks();
   concurrent_runs_ = bounds.maxNumberOfConcurrentRuns();
   concurrent_streams_ = bounds.maxNumberOfStreams();
   concurrent_threads_ = bounds.maxNumberOfThreads();
 
   if (enable_dqm_bynproc_)
     dqm_path_ += fmt::sprintf(
         "/Running on %s with %d streams on %d threads", processor_model, concurrent_streams_, concurrent_threads_);
 
   // clean characters that are deemed unsafe for DQM
   // see the definition of `s_safe` in DQMServices/Core/src/DQMStore.cc
   auto safe_for_dqm = "/ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789-+=_()# "s;
   for (auto& c : dqm_path_)
     if (safe_for_dqm.find(c) == std::string::npos)
       c = '_';
 
   // allocate atomic variables to keep track of the completion of each step, process by process
   subprocess_event_check_ = std::make_unique<std::atomic<unsigned int>[]>(concurrent_streams_);
   for (unsigned int i = 0; i < concurrent_streams_; ++i)
     subprocess_event_check_[i] = 0;
   subprocess_global_run_check_ = std::make_unique<std::atomic<unsigned int>[]>(concurrent_runs_);
   for (unsigned int i = 0; i < concurrent_runs_; ++i)
     subprocess_global_run_check_[i] = 0;
   subprocess_global_lumi_check_ = std::make_unique<std::atomic<unsigned int>[]>(concurrent_lumis_);
   for (unsigned int i = 0; i < concurrent_lumis_; ++i)
     subprocess_global_lumi_check_[i] = 0;
 
   // allocate buffers to keep track of the resources spent in the lumi and run transitions
   lumi_transition_.resize(concurrent_lumis_);
   run_transition_.resize(concurrent_runs_);
 }
 
 void FastTimerService::preSourceConstruction(edm::ModuleDescription const& module) {
   callgraph_.preSourceConstruction(module);
 }
 
 void FastTimerService::preBeginJob(edm::PathsAndConsumesOfModulesBase const& pathsAndConsumes,
                                    edm::ProcessContext const& context) {
   callgraph_.preBeginJob(pathsAndConsumes, context);
 }
 
 void FastTimerService::postBeginJob() {
   unsigned int modules = callgraph_.size();
 
   // module highlights
   for (unsigned int group : boost::irange(0ul, highlight_module_psets_.size())) {
     // copy and sort for faster search via std::binary_search
     auto labels = highlight_module_psets_[group].getUntrackedParameter<std::vector<std::string>>("modules");
     std::sort(labels.begin(), labels.end());
 
     highlight_modules_[group].label = highlight_module_psets_[group].getUntrackedParameter<std::string>("label");
     highlight_modules_[group].modules.reserve(labels.size());
     // convert the module labels in module ids
     for (unsigned int i = 0; i < modules; ++i) {
       auto const& label = callgraph_.module(i).moduleLabel();
       if (std::binary_search(labels.begin(), labels.end(), label))
         highlight_modules_[group].modules.push_back(i);
     }
   }
   highlight_module_psets_.clear();
 
   // allocate the resource counters for each stream, process, path and module
   ResourcesPerJob temp(callgraph_, highlight_modules_);
   streams_.resize(concurrent_streams_, temp);
   run_summary_.resize(concurrent_runs_, temp);
   job_summary_ = temp;
 
   // check that the DQMStore service is available
   if (enable_dqm_ and not edm::Service<dqm::legacy::DQMStore>().isAvailable()) {
     // the DQMStore is not available, disable all DQM plots
     enable_dqm_ = false;
     edm::LogWarning("FastTimerService") << "The DQMStore is not avalable, the DQM plots will not be generated";
   }
 
   // allocate the structures to hold pointers to the DQM plots
   if (enable_dqm_) {
     plots_ = std::make_unique<PlotsPerJob>(callgraph_, highlight_modules_);
   }
 }
 
 void FastTimerService::postStreamBeginRun(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::preStreamEndRun(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::postStreamEndRun(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::preGlobalBeginLumi(edm::GlobalContext const& gc) {
   ignoredSignal(__func__);
 
   // reset the lumi counters only during the main process being run
   if (isFirstSubprocess(gc)) {
     auto index = gc.luminosityBlockIndex();
     subprocess_global_lumi_check_[index] = 0;
     lumi_transition_[index].reset();
   }
 }
 
 void FastTimerService::postGlobalBeginLumi(edm::GlobalContext const& gc) { ignoredSignal(__func__); }
 
 void FastTimerService::preGlobalEndLumi(edm::GlobalContext const& gc) { ignoredSignal(__func__); }
 
 void FastTimerService::postGlobalEndLumi(edm::GlobalContext const& gc) {
   ignoredSignal(__func__);
 
   // handle the summaries only after the last subprocess has run
   auto index = gc.luminosityBlockIndex();
   bool last = isLastSubprocess(subprocess_global_lumi_check_[index]);
   if (not last)
     return;
 
   edm::LogVerbatim out("FastReport");
   auto const& label =
       fmt::sprintf("run %d, lumisection %d", gc.luminosityBlockID().run(), gc.luminosityBlockID().luminosityBlock());
   printTransition(out, lumi_transition_[index], label);
 
   if (enable_dqm_transitions_) {
     plots_->fill_lumi(lumi_transition_[index], gc.luminosityBlockID().luminosityBlock());
   }
 }
 
 void FastTimerService::preStreamBeginLumi(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::postStreamBeginLumi(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::preStreamEndLumi(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::postStreamEndLumi(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::preGlobalEndRun(edm::GlobalContext const& gc) { ignoredSignal(__func__); }
 
 void FastTimerService::postGlobalEndRun(edm::GlobalContext const& gc) {
   ignoredSignal(__func__);
 
   // handle the summaries only after the last subprocess has run
   auto index = gc.runIndex();
   bool last = isLastSubprocess(subprocess_global_run_check_[index]);
   if (not last)
     return;
 
   edm::LogVerbatim out("FastReport");
   auto const& label = fmt::sprintf("Run %d", gc.luminosityBlockID().run());
   if (print_run_summary_) {
     printSummary(out, run_summary_[index], label);
   }
   printTransition(out, run_transition_[index], label);
 
   if (enable_dqm_transitions_) {
     plots_->fill_run(run_transition_[index]);
   }
 }
 
 void FastTimerService::preSourceRun(edm::RunIndex index) { thread().measure_and_accumulate(overhead_); }
 
 void FastTimerService::postSourceRun(edm::RunIndex index) { thread().measure_and_accumulate(run_transition_[index]); }
 
 void FastTimerService::preSourceLumi(edm::LuminosityBlockIndex index) { thread().measure_and_accumulate(overhead_); }
 
 void FastTimerService::postSourceLumi(edm::LuminosityBlockIndex index) {
   thread().measure_and_accumulate(lumi_transition_[index]);
 }
 
 void FastTimerService::postEndJob() {
   // stop observing to avoid potential race conditions at exit
   tbb::task_scheduler_observer::observe(false);
   guard_.finalize();
   if (print_job_summary_) {
     edm::LogVerbatim out("FastReport");
     printSummary(out, job_summary_, "Job");
   }
   if (write_json_summary_) {
     writeSummaryJSON(job_summary_, json_filename_);
   }
 }
 
 template <typename T>
 void FastTimerService::printHeader(T& out, std::string const& label) const {
   out << "FastReport ";
   if (label.size() < 60)
     for (unsigned int i = (60 - label.size()) / 2; i > 0; --i)
       out << '-';
   out << ' ' << label << " Summary ";
   if (label.size() < 60)
     for (unsigned int i = (59 - label.size()) / 2; i > 0; --i)
       out << '-';
   out << '\n';
 }
 
 template <typename T>
 void FastTimerService::printEventHeader(T& out, std::string const& label) const {
   out << "FastReport       CPU time      Real time      Allocated    Deallocated  " << label << "\n";
   //      FastReport  ########.# ms  ########.# ms  +######### kB  -######### kB  ...
 }
 
 template <typename T>
 void FastTimerService::printEventLine(T& out, Resources const& data, std::string const& label) const {
   out << fmt::sprintf("FastReport  %10.1f ms  %10.1f ms  %+10d kB  %+10d kB  %s\n",
                       ms(data.time_thread),
                       ms(data.time_real),
                       +static_cast<int64_t>(kB(data.allocated)),
                       -static_cast<int64_t>(kB(data.deallocated)),
                       label);
 }
 
 template <typename T>
 void FastTimerService::printEventLine(T& out, AtomicResources const& data, std::string const& label) const {
   out << fmt::sprintf("FastReport  %10.1f ms  %10.1f ms  %+10d kB  %+10d kB  %s\n",
                       ms(boost::chrono::nanoseconds(data.time_thread.load())),
                       ms(boost::chrono::nanoseconds(data.time_real.load())),
                       +static_cast<int64_t>(kB(data.allocated)),
                       -static_cast<int64_t>(kB(data.deallocated)),
                       label);
 }
 
 template <typename T>
 void FastTimerService::printEvent(T& out, ResourcesPerJob const& data) const {
   printHeader(out, "Event");
   printEventHeader(out, "Modules");
   auto const& source_d = callgraph_.source();
   auto const& source = data.modules[source_d.id()];
   printEventLine(out, source.total, source_d.moduleLabel());
   for (unsigned int i = 0; i < callgraph_.processes().size(); ++i) {
     auto const& proc_d = callgraph_.processDescription(i);
     auto const& proc = data.processes[i];
     printEventLine(out, proc.total, "process " + proc_d.name_);
     for (unsigned int m : proc_d.modules_) {
       auto const& module_d = callgraph_.module(m);
       auto const& module = data.modules[m];
       printEventLine(out, module.total, "  " + module_d.moduleLabel());
     }
   }
   printEventLine(out, data.total, "total");
   out << '\n';
   printEventHeader(out, "Processes and Paths");
   printEventLine(out, source.total, source_d.moduleLabel());
   for (unsigned int i = 0; i < callgraph_.processes().size(); ++i) {
     auto const& proc_d = callgraph_.processDescription(i);
     auto const& proc = data.processes[i];
     printEventLine(out, proc.total, "process " + proc_d.name_);
     for (unsigned int p = 0; p < proc.paths.size(); ++p) {
       auto const& name = proc_d.paths_[p].name_;
       auto const& path = proc.paths[p];
       printEventLine(out, path.active, name + " (only scheduled modules)");
       printEventLine(out, path.total, name + " (including dependencies)");
     }
     for (unsigned int p = 0; p < proc.endpaths.size(); ++p) {
       auto const& name = proc_d.endPaths_[p].name_;
       auto const& path = proc.endpaths[p];
       printEventLine(out, path.active, name + " (only scheduled modules)");
       printEventLine(out, path.total, name + " (including dependencies)");
     }
   }
   printEventLine(out, data.total, "total");
   out << '\n';
   for (unsigned int group : boost::irange(0ul, highlight_modules_.size())) {
     printEventHeader(out, "Highlighted modules");
     for (unsigned int m : highlight_modules_[group].modules) {
       auto const& module_d = callgraph_.module(m);
       auto const& module = data.modules[m];
       printEventLine(out, module.total, "  " + module_d.moduleLabel());
     }
     printEventLine(out, data.highlight[group], highlight_modules_[group].label);
     out << '\n';
   }
 }
 
 template <typename T>
 void FastTimerService::printSummaryHeader(T& out, std::string const& label, bool detailed) const {
   // clang-format off
   if (detailed)
     out << "FastReport   CPU time avg.      when run  Real time avg.      when run     Alloc. avg.      when run   Dealloc. avg.      when run  ";
   //        FastReport  ########.# ms  ########.# ms  ########.# ms  ########.# ms  +######### kB  +######### kB  -######### kB  -######### kB  ...
   else
     out << "FastReport   CPU time avg.                Real time avg.                   Alloc. avg.                 Dealloc. avg.                ";
   //        FastReport  ########.# ms                 ########.# ms                 +######### kB                 -######### kB                 ...
   out << label << '\n';
   // clang-format on
 }
 
 template <typename T>
 void FastTimerService::printPathSummaryHeader(T& out, std::string const& label) const {
   // clang-format off
   out << "FastReport     CPU time sched. / depend.    Real time sched. / depend.       Alloc. sched. / depend.     Dealloc. sched. / depend.  ";
   //      FastReport  ########.# ms  ########.# ms  ########.# ms  ########.# ms  +######### kB  +######### kB  -######### kB  -######### kB  ...
   out << label << '\n';
   // clang-format on
 }
 
 template <typename T>
 void FastTimerService::printSummaryLine(T& out, Resources const& data, uint64_t events, std::string const& label) const {
   out << fmt::sprintf(
       // clang-format off
       "FastReport  %10.1f ms                 %10.1f ms                 %+10d kB                 %+10d kB                 %s\n",
       // clang-format on
       (events ? ms(data.time_thread) / events : 0),
       (events ? ms(data.time_real) / events : 0),
       (events ? +static_cast<int64_t>(kB(data.allocated) / events) : 0),
       (events ? -static_cast<int64_t>(kB(data.deallocated) / events) : 0),
       label);
 }
 
 template <typename T>
 void FastTimerService::printSummaryLine(
     T& out, AtomicResources const& data, uint64_t events, uint64_t active, std::string const& label) const {
   out << fmt::sprintf(
       // clang-format off
       "FastReport  %10.1f ms  %10.1f ms  %10.1f ms  %10.1f ms  %+10d kB  %+10d kB  %+10d kB  %+10d kB  %s\n",
       // clang-format on
       (events ? ms(data.time_thread) / events : 0),
       (active ? ms(data.time_thread) / active : 0),
       (events ? ms(data.time_real) / events : 0),
       (active ? ms(data.time_real) / active : 0),
       (events ? +static_cast<int64_t>(kB(data.allocated) / events) : 0),
       (active ? +static_cast<int64_t>(kB(data.allocated) / active) : 0),
       (events ? -static_cast<int64_t>(kB(data.deallocated) / events) : 0),
       (active ? -static_cast<int64_t>(kB(data.deallocated) / active) : 0),
       label);
 }
 
 template <typename T>
 void FastTimerService::printSummaryLine(T& out,
                                         AtomicResources const& data,
                                         uint64_t events,
                                         std::string const& label) const {
   out << fmt::sprintf(
       // clang-format off
       "FastReport  %10.1f ms                 %10.1f ms                 %+10d kB                 %+10d kB                 %s\n",
       // clang-format on
       (events ? ms(data.time_thread) / events : 0),
       (events ? ms(data.time_real) / events : 0),
       (events ? +static_cast<int64_t>(kB(data.allocated) / events) : 0),
       (events ? -static_cast<int64_t>(kB(data.deallocated) / events) : 0),
       label);
 }
 
 template <typename T>
 void FastTimerService::printSummaryLine(
     T& out, Resources const& data, uint64_t events, uint64_t active, std::string const& label) const {
   out << fmt::sprintf(
       "FastReport  %10.1f ms  %10.1f ms  %10.1f ms  %10.1f ms  %+10d kB  %+10d kB  %+10d kB  %+10d kB  %s\n",
       (events ? ms(data.time_thread) / events : 0),
       (active ? ms(data.time_thread) / active : 0),
       (events ? ms(data.time_real) / events : 0),
       (active ? ms(data.time_real) / active : 0),
       (events ? +static_cast<int64_t>(kB(data.allocated) / events) : 0),
       (active ? +static_cast<int64_t>(kB(data.allocated) / active) : 0),
       (events ? -static_cast<int64_t>(kB(data.deallocated) / events) : 0),
       (active ? -static_cast<int64_t>(kB(data.deallocated) / active) : 0),
       label);
 }
 
 template <typename T>
 void FastTimerService::printPathSummaryLine(
     T& out, Resources const& data, Resources const& total, uint64_t events, std::string const& label) const {
   out << fmt::sprintf(
       "FastReport  %10.1f ms  %10.1f ms  %10.1f ms  %10.1f ms  %+10d kB  %+10d kB  %+10d kB  %+10d kB  %s\n",
       (events ? ms(data.time_thread) / events : 0),
       (events ? ms(total.time_thread) / events : 0),
       (events ? ms(data.time_real) / events : 0),
       (events ? ms(total.time_real) / events : 0),
       (events ? +static_cast<int64_t>(kB(data.allocated) / events) : 0),
       (events ? +static_cast<int64_t>(kB(total.allocated) / events) : 0),
       (events ? -static_cast<int64_t>(kB(data.deallocated) / events) : 0),
       (events ? -static_cast<int64_t>(kB(total.deallocated) / events) : 0),
       label);
 }
 
 template <typename T>
 void FastTimerService::printSummary(T& out, ResourcesPerJob const& data, std::string const& label) const {
   printHeader(out, label);
   printSummaryHeader(out, "Modules", true);
   auto const& source_d = callgraph_.source();
   auto const& source = data.modules[source_d.id()];
   printSummaryLine(out, source.total, data.events, source.events, source_d.moduleLabel());
   for (unsigned int i = 0; i < callgraph_.processes().size(); ++i) {
     auto const& proc_d = callgraph_.processDescription(i);
     auto const& proc = data.processes[i];
     printSummaryLine(out, proc.total, data.events, "process " + proc_d.name_);
     for (unsigned int m : proc_d.modules_) {
       auto const& module_d = callgraph_.module(m);
       auto const& module = data.modules[m];
       printSummaryLine(out, module.total, data.events, module.events, "  " + module_d.moduleLabel());
     }
   }
   printSummaryLine(out, data.total, data.events, "total");
   printSummaryLine(out, data.overhead, data.events, "other");
   printSummaryLine(out, data.eventsetup, data.events, "eventsetup");
   out << '\n';
   printPathSummaryHeader(out, "Processes and Paths");
   printSummaryLine(out, source.total, data.events, source_d.moduleLabel());
   for (unsigned int i = 0; i < callgraph_.processes().size(); ++i) {
     auto const& proc_d = callgraph_.processDescription(i);
     auto const& proc = data.processes[i];
     printSummaryLine(out, proc.total, data.events, "process " + proc_d.name_);
     for (unsigned int p = 0; p < proc.paths.size(); ++p) {
       auto const& name = proc_d.paths_[p].name_;
       auto const& path = proc.paths[p];
       printPathSummaryLine(out, path.active, path.total, data.events, "  " + name);
     }
     for (unsigned int p = 0; p < proc.endpaths.size(); ++p) {
       auto const& name = proc_d.endPaths_[p].name_;
       auto const& path = proc.endpaths[p];
       printPathSummaryLine(out, path.active, path.total, data.events, "  " + name);
     }
   }
   printSummaryLine(out, data.total, data.events, "total");
   printSummaryLine(out, data.overhead, data.events, "other");
   printSummaryLine(out, data.eventsetup, data.events, "eventsetup");
   out << '\n';
   for (unsigned int group : boost::irange(0ul, highlight_modules_.size())) {
     printSummaryHeader(out, "Highlighted modules", true);
     for (unsigned int m : highlight_modules_[group].modules) {
       auto const& module_d = callgraph_.module(m);
       auto const& module = data.modules[m];
       printSummaryLine(out, module.total, data.events, module.events, module_d.moduleLabel());
     }
     printSummaryLine(out, data.highlight[group], data.events, highlight_modules_[group].label);
     out << '\n';
   }
 }
 
 template <typename T>
 void FastTimerService::printTransition(T& out, AtomicResources const& data, std::string const& label) const {
   printEventHeader(out, "Transition");
   printEventLine(out, data, label);
 }
 
 template <typename T>
 json FastTimerService::encodeToJSON(std::string const& type,
                                     std::string const& label,
                                     unsigned int events,
                                     T const& data) const {
   return json{{"type", type},
               {"label", label},
               {"events", events},
               {"time_thread", ms(data.time_thread)},
               {"time_real", ms(data.time_real)},
               {"mem_alloc", kB(data.allocated)},
               {"mem_free", kB(data.deallocated)}};
 }
 
 json FastTimerService::encodeToJSON(edm::ModuleDescription const& module, ResourcesPerModule const& data) const {
   return encodeToJSON(module.moduleName(), module.moduleLabel(), data.events, data.total);
 }
 
 void FastTimerService::writeSummaryJSON(ResourcesPerJob const& data, std::string const& filename) const {
   json j;
 
   // write a description of the resources
   j["resources"] = json::array({json{{"time_real", "real time"}},
                                 json{{"time_thread", "cpu time"}},
                                 json{{"mem_alloc", "allocated memory"}},
                                 json{{"mem_free", "deallocated memory"}}});
 
   // write the resources used by the job
   j["total"] = encodeToJSON(
       "Job", callgraph_.processDescription(0).name_, data.events, data.total + data.overhead + data.eventsetup);
 
   // write the resources used by every module
   j["modules"] = json::array();
   for (unsigned int i = 0; i < callgraph_.size(); ++i) {
     auto const& module = callgraph_.module(i);
     auto const& data_m = data.modules[i];
     j["modules"].push_back(encodeToJSON(module, data_m));
   }
 
   // add an entry for the "overhead"
   j["modules"].push_back(encodeToJSON("other", "other", data.events, data.overhead));
   j["modules"].push_back(encodeToJSON("eventsetup", "eventsetup", data.events, data.eventsetup));
 
   std::ofstream out(filename);
   out << std::setw(2) << j << std::flush;
 }
 
 // check if this is the first process being signalled
 bool FastTimerService::isFirstSubprocess(edm::StreamContext const& sc) {
   return (not sc.processContext()->isSubProcess());
 }
 
 bool FastTimerService::isFirstSubprocess(edm::GlobalContext const& gc) {
   return (not gc.processContext()->isSubProcess());
 }
 
 // check if this is the last process being signalled
 bool FastTimerService::isLastSubprocess(std::atomic<unsigned int>& check) {
   // release-acquire semantic guarantees that all writes in this and other threads are visible
   // after this operation; full sequentially-consistent ordering is (probably) not needed.
   unsigned int old_value = check.fetch_add(1, std::memory_order_acq_rel);
   return (old_value == callgraph_.processes().size() - 1);
 }
 
 void FastTimerService::preEvent(edm::StreamContext const& sc) { ignoredSignal(__func__); }
 
 void FastTimerService::postEvent(edm::StreamContext const& sc) {
   ignoredSignal(__func__);
 
   unsigned int pid = callgraph_.processId(*sc.processContext());
   unsigned int sid = sc.streamID();
   auto& stream = streams_[sid];
   auto& process = callgraph_.processDescription(pid);
 
   // measure the event resources as the sum of all modules' resources
   auto& data = stream.processes[pid].total;
   for (unsigned int id : process.modules_)
     data += stream.modules[id].total;
   stream.total += data;
 
   // handle the summaries and fill the plots only after the last subprocess has run
   bool last = isLastSubprocess(subprocess_event_check_[sid]);
   if (not last)
     return;
 
   // measure the event resources explicitly
   stream.event_measurement.measure_and_store(stream.event);
 
   // add to the event resources those used by source (which is not part of any process)
   unsigned int id = 0;
   stream.total += stream.modules[id].total;
 
   // highlighted modules
   for (unsigned int group : boost::irange(0ul, highlight_modules_.size()))
     for (unsigned int i : highlight_modules_[group].modules)
       stream.highlight[group] += stream.modules[i].total;
 
   // avoid concurrent access to the summary objects
   {
     std::lock_guard<std::mutex> guard(summary_mutex_);
     job_summary_ += stream;
     run_summary_[sc.runIndex()] += stream;
   }
 
   if (print_event_summary_) {
     edm::LogVerbatim out("FastReport");
     printEvent(out, stream);
   }
 
   if (enable_dqm_) {
     plots_->fill(callgraph_, stream, sc.eventID().luminosityBlock());
   }
 }
 
 void FastTimerService::preSourceEvent(edm::StreamID sid) {
   // clear the event counters
   auto& stream = streams_[sid];
   stream.reset();
   ++stream.events;
 
   subprocess_event_check_[sid] = 0;
 
   // reuse the same measurement for the Source module and for the explicit begin of the Event
   auto& measurement = thread();
   measurement.measure_and_accumulate(stream.overhead);
   stream.event_measurement = measurement;
 }
 
 void FastTimerService::postSourceEvent(edm::StreamID sid) {
   edm::ModuleDescription const& md = callgraph_.source();
   unsigned int id = md.id();
   auto& stream = streams_[sid];
   auto& module = stream.modules[id];
 
   thread().measure_and_store(module.total);
   ++stream.modules[id].events;
 }
 
 void FastTimerService::prePathEvent(edm::StreamContext const& sc, edm::PathContext const& pc) {
   unsigned int sid = sc.streamID().value();
   unsigned int pid = callgraph_.processId(*sc.processContext());
   unsigned int id = pc.pathID();
   auto& stream = streams_[sid];
   auto& data = pc.isEndPath() ? stream.processes[pid].endpaths[id] : stream.processes[pid].paths[id];
   data.status = false;
   data.last = 0;
 }
 
 void FastTimerService::postPathEvent(edm::StreamContext const& sc,
                                      edm::PathContext const& pc,
                                      edm::HLTPathStatus const& status) {
   unsigned int sid = sc.streamID().value();
   unsigned int pid = callgraph_.processId(*sc.processContext());
   unsigned int id = pc.pathID();
   auto& stream = streams_[sid];
   auto& data = pc.isEndPath() ? stream.processes[pid].endpaths[id] : stream.processes[pid].paths[id];
 
   auto const& path =
       pc.isEndPath() ? callgraph_.processDescription(pid).endPaths_[id] : callgraph_.processDescription(pid).paths_[id];
   unsigned int index = path.modules_on_path_.empty() ? 0 : status.index() + 1;
   data.last = path.modules_on_path_.empty() ? 0 : path.last_dependency_of_module_[status.index()];
 
   for (unsigned int i = 0; i < index; ++i) {
     auto const& module = stream.modules[path.modules_on_path_[i]];
     data.active += module.total;
   }
   for (unsigned int i = 0; i < data.last; ++i) {
     auto const& module = stream.modules[path.modules_and_dependencies_[i]];
     data.total += module.total;
   }
 }
 
 void FastTimerService::preModuleEventAcquire(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   unsigned int sid = sc.streamID().value();
   auto& stream = streams_[sid];
   thread().measure_and_accumulate(stream.overhead);
 }
 
 void FastTimerService::postModuleEventAcquire(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   edm::ModuleDescription const& md = *mcc.moduleDescription();
   unsigned int id = md.id();
   unsigned int sid = sc.streamID().value();
   auto& stream = streams_[sid];
   auto& module = stream.modules[id];
 
   module.has_acquire = true;
   thread().measure_and_store(module.total);
 }
 
 void FastTimerService::preModuleEvent(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   unsigned int sid = sc.streamID().value();
   auto& stream = streams_[sid];
   thread().measure_and_accumulate(stream.overhead);
 }
 
 void FastTimerService::postModuleEvent(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   edm::ModuleDescription const& md = *mcc.moduleDescription();
   unsigned int id = md.id();
   unsigned int sid = sc.streamID().value();
   auto& stream = streams_[sid];
   auto& module = stream.modules[id];
 
   if (module.has_acquire) {
     thread().measure_and_accumulate(module.total);
   } else {
     thread().measure_and_store(module.total);
   }
   ++module.events;
 }
 
 void FastTimerService::preModuleEventDelayedGet(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   unsupportedSignal(__func__);
 }
 
 void FastTimerService::postModuleEventDelayedGet(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   unsupportedSignal(__func__);
 }
 
 void FastTimerService::preModuleEventPrefetching(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   ignoredSignal(__func__);
 }
 
 void FastTimerService::postModuleEventPrefetching(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   ignoredSignal(__func__);
 }
 
 void FastTimerService::preEventReadFromSource(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   if (mcc.state() == edm::ModuleCallingContext::State::kPrefetching) {
     auto& stream = streams_[sc.streamID()];
     thread().measure_and_accumulate(stream.overhead);
   }
 }
 
 void FastTimerService::postEventReadFromSource(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   if (mcc.state() == edm::ModuleCallingContext::State::kPrefetching) {
     edm::ModuleDescription const& md = callgraph_.source();
     unsigned int id = md.id();
     auto& stream = streams_[sc.streamID()];
     auto& module = stream.modules[id];
 
     thread().measure_and_accumulate(module.total);
   }
 }
 
 void FastTimerService::preModuleGlobalBeginRun(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleGlobalBeginRun(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   auto index = gc.runIndex();
   thread().measure_and_accumulate(run_transition_[index]);
 }
 
 void FastTimerService::preModuleGlobalEndRun(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleGlobalEndRun(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   auto index = gc.runIndex();
   thread().measure_and_accumulate(run_transition_[index]);
 }
 
 void FastTimerService::preModuleGlobalBeginLumi(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleGlobalBeginLumi(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   auto index = gc.luminosityBlockIndex();
   thread().measure_and_accumulate(lumi_transition_[index]);
 }
 
 void FastTimerService::preModuleGlobalEndLumi(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleGlobalEndLumi(edm::GlobalContext const& gc, edm::ModuleCallingContext const& mcc) {
   auto index = gc.luminosityBlockIndex();
   thread().measure_and_accumulate(lumi_transition_[index]);
 }
 
 void FastTimerService::preModuleStreamBeginRun(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleStreamBeginRun(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   auto index = sc.runIndex();
   thread().measure_and_accumulate(run_transition_[index]);
 }
 
 void FastTimerService::preModuleStreamEndRun(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleStreamEndRun(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   auto index = sc.runIndex();
   thread().measure_and_accumulate(run_transition_[index]);
 }
 
 void FastTimerService::preModuleStreamBeginLumi(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleStreamBeginLumi(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   auto index = sc.luminosityBlockIndex();
   thread().measure_and_accumulate(lumi_transition_[index]);
 }
 
 void FastTimerService::preModuleStreamEndLumi(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   thread().measure_and_accumulate(overhead_);
 }
 
 void FastTimerService::postModuleStreamEndLumi(edm::StreamContext const& sc, edm::ModuleCallingContext const& mcc) {
   auto index = sc.luminosityBlockIndex();
   thread().measure_and_accumulate(lumi_transition_[index]);
 }
 
 void FastTimerService::preESModule(edm::eventsetup::EventSetupRecordKey const&, edm::ESModuleCallingContext const& cc) {
   auto top = cc.getTopModuleCallingContext();
   if (top->type() == edm::ParentContext::Type::kPlaceInPath) {
     //Paths are only used when processing an Event
     unsigned int sid = top->parent().placeInPathContext()->pathContext()->streamContext()->streamID().value();
     auto& stream = streams_[sid];
     thread().measure_and_accumulate(stream.overhead);
   }
 }
 void FastTimerService::postESModule(edm::eventsetup::EventSetupRecordKey const&,
                                     edm::ESModuleCallingContext const& cc) {
   auto top = cc.getTopModuleCallingContext();
   if (top->type() == edm::ParentContext::Type::kPlaceInPath) {
     unsigned int sid = top->parent().placeInPathContext()->pathContext()->streamContext()->streamID().value();
     auto& stream = streams_[sid];
     thread().measure_and_accumulate(stream.eventsetup);
   }
 }
 
 FastTimerService::ThreadGuard::ThreadGuard() {
   auto err = ::pthread_key_create(&key_, retire_thread);
   if (err) {
     throw cms::Exception("FastTimerService") << "ThreadGuard key creation failed: " << ::strerror(err);
   }
 }
 
 // If this is a new thread, register it and return true
 bool FastTimerService::ThreadGuard::register_thread(FastTimerService::AtomicResources& r) {
   auto ptr = ::pthread_getspecific(key_);
 
   if (not ptr) {
     auto p = thread_resources_.emplace_back(std::make_shared<specific_t>(r));
     auto pp = new std::shared_ptr<specific_t>(*p);
     auto err = ::pthread_setspecific(key_, pp);
     if (err) {
       throw cms::Exception("FastTimerService") << "ThreadGuard pthread_setspecific failed: " << ::strerror(err);
     }
     return true;
   }
   return false;
 }
 
 std::shared_ptr<FastTimerService::ThreadGuard::specific_t>* FastTimerService::ThreadGuard::ptr(void* p) {
   return static_cast<std::shared_ptr<specific_t>*>(p);
 }
 
 // called when a thread exits
 void FastTimerService::ThreadGuard::retire_thread(void* p) {
   auto ps = ptr(p);
   auto expected = true;
   if ((*ps)->live_.compare_exchange_strong(expected, false)) {
     // account any resources used or freed by the thread before leaving the TBB pool
     (*ps)->measurement_.measure_and_accumulate((*ps)->resource_);
   }
   delete ps;
 }
 
 // finalize all threads that have not retired
 void FastTimerService::ThreadGuard::finalize() {
   for (auto& p : thread_resources_) {
     auto expected = true;
     if (p->live_.compare_exchange_strong(expected, false)) {
       p->measurement_.measure_and_accumulate(p->resource_);
     }
   }
 }
 
 FastTimerService::Measurement& FastTimerService::ThreadGuard::thread() {
   return (*ptr(::pthread_getspecific(key_)))->measurement_;
 }
 
 void FastTimerService::on_scheduler_entry(bool worker) {
   if (guard_.register_thread(overhead_)) {
     // initialise the measurement point for a thread that has newly joined the TBB pool
     thread().measure();
   }
 }
 
 void FastTimerService::on_scheduler_exit(bool worker) {}
 
 FastTimerService::Measurement& FastTimerService::thread() { return guard_.thread(); }
 
 // describe the module's configuration
 void FastTimerService::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
   edm::ParameterSetDescription desc;
   desc.addUntracked<bool>("printEventSummary", false);
   desc.addUntracked<bool>("printRunSummary", true);
   desc.addUntracked<bool>("printJobSummary", true);
   // JSON configuration
   //desc.addUntracked<bool>("writeJSONByEvent", false);
   //desc.addUntracked<bool>("writeJSONByLumiSection", false);
   //desc.addUntracked<bool>("writeJSONByRun", false);
   desc.addUntracked<bool>("writeJSONSummary", false);
   desc.addUntracked<std::string>("jsonFileName", "resources.json");
   // DQM configuration
   desc.addUntracked<bool>("enableDQM", true);
   desc.addUntracked<bool>("enableDQMbyModule", false);
   desc.addUntracked<bool>("enableDQMbyPath", false);
   desc.addUntracked<bool>("enableDQMbyLumiSection", false);
   desc.addUntracked<bool>("enableDQMbyProcesses", false);
   desc.addUntracked<bool>("enableDQMTransitions", false);
   desc.addUntracked<double>("dqmTimeRange", 1000.);              // ms
   desc.addUntracked<double>("dqmTimeResolution", 5.);            // ms
   desc.addUntracked<double>("dqmMemoryRange", 1000000.);         // kB
   desc.addUntracked<double>("dqmMemoryResolution", 5000.);       // kB
   desc.addUntracked<double>("dqmPathTimeRange", 100.);           // ms
   desc.addUntracked<double>("dqmPathTimeResolution", 0.5);       // ms
   desc.addUntracked<double>("dqmPathMemoryRange", 1000000.);     // kB
   desc.addUntracked<double>("dqmPathMemoryResolution", 5000.);   // kB
   desc.addUntracked<double>("dqmModuleTimeRange", 40.);          // ms
   desc.addUntracked<double>("dqmModuleTimeResolution", 0.2);     // ms
   desc.addUntracked<double>("dqmModuleMemoryRange", 100000.);    // kB
   desc.addUntracked<double>("dqmModuleMemoryResolution", 500.);  // kB
   desc.addUntracked<unsigned>("dqmLumiSectionsRange", 2500);     // ~ 16 hours
   desc.addUntracked<std::string>("dqmPath", "HLT/TimerService");
 
   edm::ParameterSetDescription highlightModulesDescription;
   highlightModulesDescription.addUntracked<std::vector<std::string>>("modules", {});
   highlightModulesDescription.addUntracked<std::string>("label", "producers");
   desc.addVPSetUntracked("highlightModules", highlightModulesDescription, {});
 
   // # OBSOLETE - these parameters are ignored, they are left only not to break old configurations
   // they will not be printed in the generated cfi.py file
   desc.addOptionalNode(edm::ParameterDescription<bool>("useRealTimeClock", true, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableTimingPaths", true, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableTimingModules", true, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableTimingExclusive", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableTimingSummary", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("skipFirstPath", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathActive", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathTotal", true, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathOverhead", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathDetails", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathCounters", true, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyPathExclusive", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMbyModuleType", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
   desc.addOptionalNode(edm::ParameterDescription<bool>("enableDQMSummary", false, false), false)
       ->setComment("This parameter is obsolete and will be ignored.");
 
   descriptions.add("FastTimerService", desc);
 }
 
 // declare FastTimerService as a framework Service
 #include "FWCore/ServiceRegistry/interface/ServiceMaker.h"
 DEFINE_FWK_SERVICE(FastTimerService);