VecHistogramFiller1D.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/  
  
#ifndef AthenaMonitoringKernel_HistogramFiller_VecHistogramFiller1D_h
#define AthenaMonitoringKernel_HistogramFiller_VecHistogramFiller1D_h
 
#include "HistogramFiller1D.h"
 
namespace Monitored {
  class VecHistogramFiller1D : public HistogramFiller1D {
  public:
    VecHistogramFiller1D(const HistogramDef& definition, std::shared_ptr<IHistogramProvider> provider)
      : HistogramFiller1D(definition, std::move(provider)) {}
 
 
    virtual unsigned fill(const HistogramFiller::VariablesPack& vars) const override {
      if ( ATH_UNLIKELY(vars.size() == 0 or vars[0] == nullptr) ) { return 0; }
 
      if (vars.cut) {
        const size_t maskSize = vars.cut->size();
        // Abort if no cut entries or first (and only) entry is false
        if (maskSize == 0 || (maskSize == 1 && !vars.cut->get(0))) { return 0; }
 
        if (ATH_UNLIKELY(maskSize > 1 && maskSize != vars[0]->size())) {
          MsgStream log(Athena::getMessageSvc(), "VecHistogramFiller1D");
          log << MSG::ERROR << "CutMask does not match the size of plotted variable: "
              << maskSize << " " << vars[0]->size() << endmsg;
        }
      }
 
      auto histogram = this->histogram<TH1>();
       
      if ( ATH_UNLIKELY( static_cast<size_t>(histogram->GetNbinsX()) + (m_histDef->kVecUO ? 2 : 0) != vars[0]->size() ) ) {
          MsgStream log(Athena::getMessageSvc(), "VecHistogramFiller1D");
          log << MSG::WARNING << "Histogram " << histogram->GetName()
              << " filled with kVec(UO) option with variable " << vars[0]->name()
              << " have incompatible sizes (histogram) " << histogram->GetNbinsX()
              << " (variable) "  << vars[0]->size()
              << " They ought to match exactly for kVec option or n. hist. bins +2 == var. size fro kVecOU" << endmsg;
              return 0;
      }
 
      const unsigned offset = m_histDef->kVecUO ? 0 : 1;
      for (unsigned i = 0; i < vars[0]->size(); ++i) {
        if (vars.cut == nullptr or vars.cut->get(i)) {
          const double value = vars[0]->get(i);
          histogram->AddBinContent(i+offset, value);
          histogram->SetEntries(histogram->GetEntries() + value);
        }
      }
 
      return vars[0]->size();
    }
  };
}
 
#endif /* AthenaMonitoringKernel_HistogramFiller_VecHistogramFiller1D_h */