DQM/DTMonitorModule/src/DTTimeEvolutionHisto.cc - master - githubmirror/cmssw

DQM/DTMonitorModule/src/DTTimeEvolutionHisto.cc
223 строки8 KB
iarspider
Update DTTimeEvolutionHisto.cc
18 окт 2024, 14:08
18 окт 2024, 14:084f3dd79

/*
 *  See header file for a description of this class.
 *
 *  \author G. Cerminara - INFN Torino
 */

#include "DQM/DTMonitorModule/interface/DTTimeEvolutionHisto.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;
using namespace edm;

DTTimeEvolutionHisto::DTTimeEvolutionHisto(DQMStore::IBooker& ibooker,
                                           const string& name,
                                           const string& title,
                                           int nbins,
                                           int lsPrescale,
                                           bool sliding,
                                           int mode)
    : DTTimeEvolutionHisto(ibooker, name, title, nbins, 1, lsPrescale, sliding, mode) {
  nBookedBins = histo->getNbinsX();
}

DTTimeEvolutionHisto::DTTimeEvolutionHisto(DQMStore::IBooker& ibooker,
                                           const string& name,
                                           const string& title,
                                           int nbins,
                                           int firstLS,
                                           int lsPrescale,
                                           bool sliding,
                                           int mode)
    : valueLastTimeSlot(0), theFirstLS(firstLS), theLSPrescale(lsPrescale), doSlide(sliding), theMode(mode) {
  // set the number of bins to be booked
  nBookedBins = nbins;
  binLabelCounter = -1;
  if (sliding)
    nBookedBins++;
  if (!sliding && theMode == 0)
    LogWarning("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
        << "[DTTimeEvolutionHisto]***Error: wrong configuration" << endl;

  stringstream realTitle;
  realTitle << title << "/" << theLSPrescale << " LS";

  // book the ME

  histo = ibooker.book1D(name, realTitle.str(), nBookedBins, (float)theFirstLS, nBookedBins + 1.);

  // set the axis label
  if (sliding) {
    histo->setBinLabel(1, "avg. previous", 1);
  } else {
    // loop over bins and

    for (int bin = 1; bin != nBookedBins + 1; ++bin) {
      stringstream label;
      if (theLSPrescale > 1) {
        label << "LS " << ((bin - 1) * theLSPrescale) + theFirstLS << "-" << bin * theLSPrescale + theFirstLS;
      } else {
        label << "LS " << ((bin - 1) * theLSPrescale) + theFirstLS;
      }
      if (bin % (2 * (int)theLSPrescale) == 0)
        histo->setBinLabel(bin, label.str(), 1);  //JF to allow easy reading of x-axis
    }
  }
}

//FR changed previous constructor with 2 arguments:
//no igetter here!! so I get the histo from the client and just instanciate here a DTTimeEvolutionHisto

DTTimeEvolutionHisto::DTTimeEvolutionHisto(MonitorElement* histoGot)
    : valueLastTimeSlot(0),
      theFirstLS(1),
      theLSPrescale(-1),
      doSlide(false),
      theMode(0) {  // FIXME: set other memebers to sensible values
  LogVerbatim("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
      << "[DTTimeEvolutionHisto] Retrieve ME with name: "
      << " " << endl;
  histo = histoGot;
}

DTTimeEvolutionHisto::~DTTimeEvolutionHisto() {}

void DTTimeEvolutionHisto::setTimeSlotValue(float value, int timeSlot) {
  if (!doSlide) {
    histo->Fill(timeSlot, value);
  } else {
    for (int bin = 1; bin != nBookedBins; ++bin) {
      float value = histo->getBinContent(bin);

      if (bin == 1) {  // average of previous time slots (fixme)
        histo->setBinContent(bin, (value + histo->getBinContent(bin + 1)) / 2.);
      } else if (bin != nBookedBins) {
        histo->setBinContent(bin, histo->getBinContent(bin + 1));
        histo->setBinError(bin, histo->getBinError(bin + 1));
        histo->setBinLabel(bin, histo->getTH1F()->GetXaxis()->GetBinLabel(bin + 1), 1);  //slide to left
        histo->setBinLabel(bin + 1, "", 1);  //delete old label to avoid duplication
      }
    }
    histo->setBinContent(nBookedBins, value);
  }
}

void DTTimeEvolutionHisto::accumulateValueTimeSlot(float value) { valueLastTimeSlot += value; }

void DTTimeEvolutionHisto::updateTimeSlot(int ls, int nEventsInLS) {
  if (doSlide) {  // sliding bins
    // count LS in this time-slot
    if (nEventsInLastTimeSlot.find(ls) != nEventsInLastTimeSlot.end()) {
      nEventsInLastTimeSlot[ls] += nEventsInLS;
      nLumiTrInLastTimeSlot[ls]++;
    } else {
      nEventsInLastTimeSlot[ls] = nEventsInLS;
      nLumiTrInLastTimeSlot[ls] = 1;
    }

    if (!nEventsInLastTimeSlot.empty() &&
        nEventsInLastTimeSlot.size() % theLSPrescale == 0) {  // update the value of the slot and reset the counters
      int firstLSinTimeSlot = nEventsInLastTimeSlot.begin()->first;
      int lastLSinTimeSlot = nEventsInLastTimeSlot.rbegin()->first;

      map<int, int>::const_iterator nEventsIt = nEventsInLastTimeSlot.begin();
      map<int, int>::const_iterator nEventsEnd = nEventsInLastTimeSlot.end();

      int nEvents = 0;
      for (; nEventsIt != nEventsEnd; ++nEventsIt)
        nEvents += nEventsIt->second;

      LogVerbatim("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
          << "[DTTimeEvolutionHisto] Update time-slot, # entries: " << valueLastTimeSlot << " # events: " << nEvents
          << endl;
      // set the bin content

      float value = 0;
      if (theMode == 0) {
        if (nEvents != 0)
          value = valueLastTimeSlot / (float)nEvents;
      } else if (theMode == 1) {
        value = valueLastTimeSlot;
      } else if (theMode == 2) {
        value = nEvents;
      } else if (theMode == 3) {
        map<int, int>::const_iterator nLumiTrIt = nLumiTrInLastTimeSlot.begin();
        map<int, int>::const_iterator nLumiTrEnd = nLumiTrInLastTimeSlot.end();

        float nLumiTr = 0.;
        for (; nLumiTrIt != nLumiTrEnd; ++nLumiTrIt)
          nLumiTr += nLumiTrIt->second;

        value = valueLastTimeSlot / nLumiTr;
      }
      setTimeSlotValue(value, nBookedBins);
      LogVerbatim("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
          << "       updated value: " << histo->getBinContent(nBookedBins) << endl;

      // set the bin label
      stringstream binLabel;
      binLabel << "LS " << firstLSinTimeSlot;
      if (nEventsInLastTimeSlot.size() > 1) {
        binLabel << "-" << lastLSinTimeSlot;
        binLabelCounter++;
      }

      // Set only labels which can be seen in the plot without zooming
      if (nBookedBins >= 25 && binLabelCounter % (nBookedBins / 25) == 0)  //around 25 labels fit in a full size hist
        histo->setBinLabel(nBookedBins, binLabel.str(), 1);
      //first label, LS=1 ideally, is guaranteed by setting binLabelCounter=-1 in constructor

      // reset the counters for the time slot
      nEventsInLastTimeSlot.clear();
      nLumiTrInLastTimeSlot.clear();
      valueLastTimeSlot = 0;
    }

  } else {
    int binN = (int)ls - (theFirstLS - 1) / (int)theLSPrescale;
    // set the bin content
    float value = 0;
    if (theMode == 1) {
      value = valueLastTimeSlot;
    } else if (theMode == 2) {
      value = nEventsInLS;
    } else if (theMode == 3) {
      value = valueLastTimeSlot / theLSPrescale;
    }
    LogVerbatim("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
        << "[DTTimeEvolutionHisto] Update time-slot: " << binN << " with value: " << value << endl;
    setTimeSlotValue(value, binN);
  }
}

void DTTimeEvolutionHisto::normalizeTo(const MonitorElement* histForNorm) {
  if (histo == nullptr) {
    LogWarning("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
        << "[DTTimeEvolutionHisto]***Error: pointer to ME is NULL" << endl;
    return;
  }
  int nBins = histo->getNbinsX();
  if (histForNorm->getNbinsX() != nBins) {
    LogWarning("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
        << "[DTTimeEvolutionHisto]***Error: normalizing histos with != # of bins" << endl;
    return;
  }
  for (int bin = 1; bin <= nBins; ++bin) {  // loop over bins
    if (histForNorm->getBinContent(bin) != 0) {
      double normValue = histo->getBinContent(bin) / histForNorm->getBinContent(bin);
      LogVerbatim("DTDQM|DTMonitorModule|DTMonitorClient|DTTimeEvolutionHisto")
          << "[DTTimeEvolutionHisto] Normalizing bin: " << bin << " to: " << histo->getBinContent(bin) << " / "
          << histForNorm->getBinContent(bin) << " = " << normValue << endl;
      histo->setBinContent(bin, normValue);
    } else {
      histo->setBinContent(bin, 0);
    }
  }
}

// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End: