GEMBaseValidation.cc

Go to the documentation of this file.

#include "Validation/MuonGEMHits/interface/GEMBaseValidation.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/GEMGeometry/interface/GEMEtaPartitionSpecs.h"

#include <memory>

using namespace dqm::impl;

GEMBaseValidation::GEMBaseValidation(const edm::ParameterSet& ps, std::string log_category)
    : kLogCategory_(log_category) {
  pid_list_ = ps.getUntrackedParameter<std::vector<Int_t> >("pidList");
  zr_occ_num_bins_ = ps.getUntrackedParameter<std::vector<Int_t> >("ZROccNumBins");
  zr_occ_range_ = ps.getUntrackedParameter<std::vector<Double_t> >("ZROccRange");
  xy_occ_num_bins_ = ps.getUntrackedParameter<Int_t>("XYOccNumBins", 360);
  // TODO depends on the station.. for detail plots..
  eta_range_ = ps.getUntrackedParameter<std::vector<Double_t> >("EtaOccRange");

  detail_plot_ = ps.getParameter<Bool_t>("detailPlot");
}

GEMBaseValidation::~GEMBaseValidation() {}

Int_t GEMBaseValidation::getDetOccBinX(Int_t num_layers, Int_t chamber_id, Int_t layer_id) {
  return num_layers * chamber_id + layer_id - num_layers;
}

Bool_t GEMBaseValidation::isMuonSimHit(const PSimHit& simhit) { return std::abs(simhit.particleType()) == kMuonPDGId_; }

Float_t GEMBaseValidation::toDegree(Float_t radian) {
  Float_t degree = radian / M_PI * 180;
  if (degree < -5)
    return degree + 360;
  else
    return degree;
}

Int_t GEMBaseValidation::getPidIdx(Int_t pid) {
  return std::find(pid_list_.begin(), pid_list_.end(), pid) - pid_list_.begin();
}

MonitorElement* GEMBaseValidation::bookZROccupancy(DQMStore::IBooker& booker,
                                                   Int_t region_id,
                                                   const char* name_prefix,
                                                   const char* title_prefix) {
  auto name_suffix = GEMUtils::getSuffixName(region_id);
  auto title_suffix = GEMUtils::getSuffixTitle(region_id);

  TString name = TString::Format("%s_occ_zr%s", name_prefix, name_suffix.Data());
  TString title = TString::Format("%s ZR Occupancy :%s;|Z| [cm];R [cm]", title_prefix, title_suffix.Data());

  Double_t station0_xmin = zr_occ_range_[0];
  Double_t station0_xmax = zr_occ_range_[1];
  Double_t station1_xmin = zr_occ_range_[4];
  Double_t station1_xmax = zr_occ_range_[5];
  Double_t station2_xmin = zr_occ_range_[8];
  Double_t station2_xmax = zr_occ_range_[9];

  std::vector<Double_t> xbins_vector;
  for (Double_t i = station0_xmin - 1; i < station2_xmax + 1; i += 0.25) {
    if (i > station0_xmax + 1 and i < station1_xmin - 1)
      continue;
    if (i > station1_xmax + 1 and i < station2_xmin - 1)
      continue;
    xbins_vector.push_back(i);
  }

  Int_t nbinsx = xbins_vector.size() - 1;

  Int_t nbinsy = zr_occ_num_bins_[2];
  Double_t ylow = std::min(zr_occ_range_[2], std::min(zr_occ_range_[6], zr_occ_range_[10]));
  Double_t yup = std::max(zr_occ_range_[3], std::max(zr_occ_range_[7], zr_occ_range_[11]));

  auto hist = new TH2F(name, title, nbinsx, &xbins_vector[0], nbinsy, ylow, yup);
  return booker.book2D(name, hist);
}