GEMCheckGeometry.cc

Go to the documentation of this file.

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/GEMGeometry/interface/GEMEtaPartitionSpecs.h"
#include "Geometry/GEMGeometry/interface/GEMGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "Validation/MuonGEMDigis/interface/GEMCheckGeometry.h"

#include <iomanip>

GEMCheckGeometry::GEMCheckGeometry(const edm::ParameterSet &gc) {
  GE11PhiBegin_ = gc.getUntrackedParameter<double>("GE11PhiBegin", -5.);
  GE11PhiStep_ = gc.getUntrackedParameter<double>("GE11PhiStep", 10);
  minPhi_ = gc.getUntrackedParameter<double>("minPhi", -180.);
  maxPhi_ = gc.getUntrackedParameter<double>("maxPhi", +180.);
  detailPlot_ = gc.getParameter<bool>("detailPlot");
}

void GEMCheckGeometry::bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &Run, edm::EventSetup const &iSetup) {
  if (!detailPlot_)
    return;
  const GEMGeometry *GEMGeometry_;

  try {
    edm::ESHandle<GEMGeometry> hGeom;
    iSetup.get<MuonGeometryRecord>().get(hGeom);
    GEMGeometry_ = &*hGeom;
  } catch (edm::eventsetup::NoProxyException<GEMGeometry> &e) {
    edm::LogError("MuonGEMGeometry") << "+++ Error : GEM geometry  is unavailable on event loop. +++\n";
    return;
  }

  ibooker.setCurrentFolder("MuonGEMDigisV/GEMDigisTask");
  LogDebug("GEMCheckGeometry") << "ibooker set current folder\n";

  for (auto region : GEMGeometry_->regions()) {
    TString title = TString::Format("Geometry's phi distribution on Region %d ; #phi(degree); ;", region->region());
    TString name = TString::Format("geo_phi_r%d", region->region());
    auto temp_me = ibooker.book2D(name.Data(), title.Data(), 360000, -180., 180, 12, 1, 13);
    temp_me->setBinLabel(1, "St1,La1_even", 2);
    temp_me->setBinLabel(2, "St1,La1_odd", 2);
    temp_me->setBinLabel(3, "St1,La2_even", 2);
    temp_me->setBinLabel(4, "St1,La2_odd", 2);
    temp_me->setBinLabel(5, "St2,La1_even", 2);
    temp_me->setBinLabel(6, "St2,La1_odd", 2);
    temp_me->setBinLabel(7, "St2,La2_even", 2);
    temp_me->setBinLabel(8, "St2,La2_odd", 2);
    theStdPlots.insert(std::map<UInt_t, MonitorElement *>::value_type(name.Hash(), temp_me));
  }

  for (auto region : GEMGeometry_->regions())
    for (auto station : region->stations())
      for (auto ring : station->rings())
        for (auto sch : ring->superChambers())
          for (auto ch : sch->chambers())
            for (auto roll : ch->etaPartitions()) {
              const StripTopology *topology(&(roll->specificTopology()));
              auto parameters(roll->specs()->parameters());
              float nStrips(parameters[3]);
              for (int strip = 0; strip <= nStrips; strip++) {
                LocalPoint lEdge(topology->localPosition(strip));
                // double x = roll->toGlobal(lEdge).x();

                // double y = roll->toGlobal(lEdge).y();
                // double z = roll->toGlobal(lEdge).z();
                // double eta = roll->toGlobal(lEdge).eta();
                double phi = roll->toGlobal(lEdge).phi().degrees();

                GEMDetId id(roll->id());
                int region_idx = id.region();
                int station_idx = id.station();
                int chamber_idx = id.chamber();
                int layer_idx = id.layer();
                // int roll_idx = id.roll();
                int value = (station_idx - 1) * 4 + (layer_idx - 1) * 2 + (chamber_idx % 2) + 1;

                if (region_idx == 1) {
                  UInt_t hash = TString("geo_phi_r1").Hash();
                  theStdPlots[hash]->Fill(phi, value);
                } else {
                  UInt_t hash = TString("geo_phi_r-1").Hash();
                  theStdPlots[hash]->Fill(phi, value);
                }
              }
            }
}

GEMCheckGeometry::~GEMCheckGeometry() {}

void GEMCheckGeometry::analyze(const edm::Event &e, const edm::EventSetup &) {}