L1TPUM.cc

Go to the documentation of this file.

/*
 * \file L1TPUM.cc
 *
 * N. Smith <nick.smith@cern.ch>
 */

#include "DQM/L1TMonitor/interface/L1TPUM.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Provenance/interface/EventAuxiliary.h"

#include "CondFormats/RunInfo/interface/RunInfo.h"
#include "CondFormats/DataRecord/interface/RunSummaryRcd.h"

// TODO: move to configuration?
namespace {
  const unsigned int R10BINS = 1024;
  const float R10MIN = -0.5;
  const float R10MAX = 1023.5;

  const unsigned int PUMETABINS = 22;
  const unsigned int PUMNORMALIZE = 22;

  const unsigned int PUMBINS = 18;
  const float PUMMIN = -0.5;
  const float PUMMAX = 17.5;
}  // namespace

L1TPUM::L1TPUM(const edm::ParameterSet& ps)
    : regionSource_(consumes<L1CaloRegionCollection>(ps.getParameter<edm::InputTag>("regionSource"))),
      histFolder_(ps.getParameter<std::string>("histFolder")) {}

L1TPUM::~L1TPUM() {}

void L1TPUM::analyze(const edm::Event& event, const edm::EventSetup& es) {
  edm::Handle<L1CaloRegionCollection> regionCollection;
  event.getByToken(regionSource_, regionCollection);

  int nonzeroRegionsBxP2{0};
  int nonzeroRegionsBx0{0};
  int nonzeroRegionsBxM2{0};

  float regionsTotalEtBxP2{0.};
  float regionsTotalEtBx0{0.};
  float regionsTotalEtBxM2{0.};

  for (const auto& region : *regionCollection) {
    if (region.et() > 0) {
      if (region.bx() == 0) {
        nonzeroRegionsBx0++;
        regionsTotalEtBx0 += region.et();
      } else if (region.bx() == 2) {
        nonzeroRegionsBxP2++;
        regionsTotalEtBxP2 += region.et();
      } else if (region.bx() == -2) {
        nonzeroRegionsBxM2++;
        regionsTotalEtBxM2 += region.et();
      }
    }
  }

  regionsTotalEtBxP2_->Fill(regionsTotalEtBxP2);
  regionsTotalEtBx0_->Fill(regionsTotalEtBx0);
  regionsTotalEtBxM2_->Fill(regionsTotalEtBxM2);

  regionsAvgEtBxP2_->Fill(regionsTotalEtBxP2 / 396.);
  regionsAvgEtBx0_->Fill(regionsTotalEtBx0 / 396.);
  regionsAvgEtBxM2_->Fill(regionsTotalEtBxM2 / 396.);

  regionsAvgNonZeroEtBxP2_->Fill(regionsTotalEtBxP2 / nonzeroRegionsBxP2);
  regionsAvgNonZeroEtBx0_->Fill(regionsTotalEtBx0 / nonzeroRegionsBx0);
  regionsAvgNonZeroEtBxM2_->Fill(regionsTotalEtBxM2 / nonzeroRegionsBxM2);

  nonZeroRegionsBxP2_->Fill(nonzeroRegionsBxP2);
  nonZeroRegionsBx0_->Fill(nonzeroRegionsBx0);
  nonZeroRegionsBxM2_->Fill(nonzeroRegionsBxM2);

  for (const auto& region : *regionCollection) {
    size_t etaBin = region.gctEta();
    regionBxPopulation_->Fill(etaBin * 18 + region.gctPhi(), region.bx());
    regionBxEtSum_->Fill(etaBin * 18 + region.gctPhi(), region.bx(), region.et());
    if (region.bx() == 0)
      regionsPUMEtaBx0_[etaBin]->Fill(nonzeroRegionsBx0 / PUMNORMALIZE, region.et());
    else if (region.bx() == 2)
      regionsPUMEtaBxP2_[etaBin]->Fill(nonzeroRegionsBxP2 / PUMNORMALIZE, region.et());
    else if (region.bx() == -2)
      regionsPUMEtaBxM2_[etaBin]->Fill(nonzeroRegionsBxM2 / PUMNORMALIZE, region.et());
  }
}

void L1TPUM::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run& run, const edm::EventSetup& es) {
  ibooker.setCurrentFolder(histFolder_ + "/BX0");
  regionsPUMEtaBx0_.resize(PUMETABINS);
  for (size_t ieta = 0; ieta < PUMETABINS; ++ieta) {
    regionsTotalEtBx0_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
    regionsAvgEtBx0_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
    regionsAvgNonZeroEtBx0_ =
        ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
    nonZeroRegionsBx0_ = ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
    regionsPUMEtaBx0_[ieta] = ibooker.book2D("regionsPUMEta" + std::to_string(ieta),
                                             "PUM Bin rank distribution;PU bin;Rank",
                                             PUMBINS,
                                             PUMMIN,
                                             PUMMAX,
                                             R10BINS,
                                             R10MIN,
                                             R10MAX);
  }

  ibooker.setCurrentFolder(histFolder_ + "/BXP2");
  regionsPUMEtaBxP2_.resize(PUMETABINS);
  for (size_t ieta = 0; ieta < PUMETABINS; ++ieta) {
    regionsTotalEtBxP2_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
    regionsAvgEtBxP2_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
    regionsAvgNonZeroEtBxP2_ =
        ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
    nonZeroRegionsBxP2_ =
        ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
    regionsPUMEtaBxP2_[ieta] = ibooker.book2D("regionsPUMEta" + std::to_string(ieta),
                                              "PUM Bin rank distribution;PU bin;Rank",
                                              PUMBINS,
                                              PUMMIN,
                                              PUMMAX,
                                              R10BINS,
                                              R10MIN,
                                              R10MAX);
  }

  ibooker.setCurrentFolder(histFolder_ + "/BXM2");
  regionsPUMEtaBxM2_.resize(PUMETABINS);
  for (size_t ieta = 0; ieta < PUMETABINS; ++ieta) {
    regionsTotalEtBxM2_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
    regionsAvgEtBxM2_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
    regionsAvgNonZeroEtBxM2_ =
        ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
    nonZeroRegionsBxM2_ =
        ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
    regionsPUMEtaBxM2_[ieta] = ibooker.book2D("regionsPUMEta" + std::to_string(ieta),
                                              "PUM Bin rank distribution;PU bin;Rank",
                                              PUMBINS,
                                              PUMMIN,
                                              PUMMAX,
                                              R10BINS,
                                              R10MIN,
                                              R10MAX);
  }

  ibooker.setCurrentFolder(histFolder_ + "/RegionBxInfo");
  regionBxPopulation_ =
      ibooker.book2D("regionBxPopulation",
                     "Event counts per region per bunch crossing;Region index (18*eta+phi);BX index;Counts",
                     396,
                     -0.5,
                     395.5,
                     5,
                     -2.5,
                     2.5);
  regionBxEtSum_ = ibooker.book2D("regionBxEtSum",
                                  "Et per region per bunch crossing;Region index (18*eta+phi);BX index;Counts*et",
                                  396,
                                  -0.5,
                                  395.5,
                                  5,
                                  -2.5,
                                  2.5);
}