GEMDQMSource Class Reference

Inheritance diagram for GEMDQMSource:

Public Member Functions
	GEMDQMSource (const edm::ParameterSet &cfg)
	~GEMDQMSource () override

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

Private Member Functions
int	findVFAT (float min_, float max_, float x_, int roll_)
const GEMGeometry *	initGeometry (edm::EventSetup const &iSetup)

Private Attributes
float	fGlobXMax_
float	fGlobXMin_
float	fGlobYMax_
float	fGlobYMin_
std::vector< GEMChamber >	gemChambers_
const GEMGeometry *	GEMGeometry_
int	nClusterSizeBinNum_
int	nIdxFirstStrip_
std::unordered_map< UInt_t, MonitorElement * >	recGlobalPos
std::unordered_map< UInt_t, MonitorElement * >	recHitME_
std::unordered_map< UInt_t, MonitorElement * >	rh_vs_eta_
std::unordered_map< UInt_t, MonitorElement * >	StripsFired_vs_eta_
edm::EDGetToken	tagRecHit_
std::unordered_map< UInt_t, MonitorElement * >	VFAT_vs_ClusterSize_

Detailed Description

Definition at line 25 of file GEMDQMSource.cc.

Constructor & Destructor Documentation

GEMDQMSource()

GEMDQMSource::GEMDQMSource

(

const edm::ParameterSet &

cfg

)

Definition at line 86 of file GEMDQMSource.cc.

                                                     {
  tagRecHit_ = consumes<GEMRecHitCollection>(cfg.getParameter<edm::InputTag>("recHitsInputLabel"));
 
  nIdxFirstStrip_ = cfg.getParameter<int>("idxFirstStrip");
  nClusterSizeBinNum_ = cfg.getParameter<int>("ClusterSizeBinNum");
 
  fGlobXMin_ = cfg.getParameter<double>("global_x_bound_min");
  fGlobXMax_ = cfg.getParameter<double>("global_x_bound_max");
  fGlobYMin_ = cfg.getParameter<double>("global_y_bound_min");
  fGlobYMax_ = cfg.getParameter<double>("global_y_bound_max");
}

References looper::cfg.

~GEMDQMSource()

GEMDQMSource::~GEMDQMSource ( )

inlineoverride

Definition at line 28 of file GEMDQMSource.cc.

{};

Member Function Documentation

analyze()

void GEMDQMSource::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotected

Definition at line 196 of file GEMDQMSource.cc.

                                                                                 {
  const GEMGeometry* GEMGeometry_ = initGeometry(eventSetup);
  if (GEMGeometry_ == nullptr)
    return;
 
  edm::Handle<GEMRecHitCollection> gemRecHits;
  event.getByToken(this->tagRecHit_, gemRecHits);
  if (!gemRecHits.isValid()) {
    edm::LogError("GEMDQMSource") << "GEM recHit is not valid.\n";
    return;
  }
  for (const auto& ch : gemChambers_) {
    GEMDetId cId = ch.id();
    for (auto roll : ch.etaPartitions()) {
      GEMDetId rId = roll->id();
      const auto& recHitsRange = gemRecHits->get(rId);
      auto gemRecHit = recHitsRange.first;
      for (auto hit = gemRecHit; hit != recHitsRange.second; ++hit) {
        Int_t nIdxStrip = hit->firstClusterStrip() + 0.5 * hit->clusterSize() - nIdxFirstStrip_;
        Int_t nVfat = 8 * ((Int_t)(nIdxStrip / (roll->nstrips() / 3)) + 1) - rId.roll();  // Strip:Start at 0
        recHitME_[cId]->Fill(nVfat);
        rh_vs_eta_[cId]->Fill(hit->localPosition().x(), rId.roll());
        Int_t nCLSOverflow = (hit->clusterSize() <= nClusterSizeBinNum_ ? hit->clusterSize() : nClusterSizeBinNum_);
        VFAT_vs_ClusterSize_[cId]->Fill(nCLSOverflow, nVfat);
        for (int i = hit->firstClusterStrip(); i < (hit->firstClusterStrip() + hit->clusterSize()); i++) {
          StripsFired_vs_eta_[cId]->Fill(i, rId.roll());
        }
 
        GlobalPoint recHitGP = GEMGeometry_->idToDet(hit->gemId())->surface().toGlobal(hit->localPosition());
        GEMDetId idLayer(rId.region(), rId.ring(), rId.station(), rId.layer(), 0, 0);
        recGlobalPos[idLayer]->Fill(recHitGP.x(), recHitGP.y());
      }
    }
  }
}

References muonGEMRecHitPSet::gemRecHit, gemRecHits_cfi::gemRecHits, mps_fire::i, GEMGeometry::idToDet(), GEMDetId::layer(), GEMDetId::region(), GEMDetId::ring(), GEMDetId::roll(), GEMDetId::station(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), hit::x, and PV3DBase< T, PVType, FrameType >::y().

bookHistograms()

void GEMDQMSource::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &  iSetup  )

overrideprotected

Definition at line 113 of file GEMDQMSource.cc.

                                                                                                      {
  std::vector<GEMDetId> listLayerOcc;
 
  GEMGeometry_ = initGeometry(iSetup);
  if (GEMGeometry_ == nullptr)
    return;
 
  const std::vector<const GEMSuperChamber*>& superChambers_ = GEMGeometry_->superChambers();
  for (auto sch : superChambers_) {
    int n_lay = sch->nChambers();
    for (int l = 0; l < n_lay; l++) {
      Bool_t bExist = false;
      for (const auto& ch : gemChambers_)
        if (ch.id() == sch->chamber(l + 1)->id())
          bExist = true;
      if (bExist)
        continue;
 
      gemChambers_.push_back(*sch->chamber(l + 1));
    }
  }
 
  ibooker.cd();
  ibooker.setCurrentFolder("GEM/recHit");
 
  for (const auto& ch : gemChambers_) {
    GEMDetId gid = ch.id();
 
    std::string strIdxName = "Gemini_" + to_string(gid.chamber()) + "_GE" + (gid.region() > 0 ? "p" : "m") +
                             to_string(gid.station()) + "_" + to_string(gid.layer());
    std::string strIdxTitle = "GEMINIm" + to_string(gid.chamber()) + " in GE" + (gid.region() > 0 ? "+" : "-") +
                              to_string(gid.station()) + "/" + to_string(gid.layer());
 
    std::string strStId = (gid.region() > 0 ? "p" : "m") + std::to_string(gid.station());
    std::string strStT = (gid.region() > 0 ? "+" : "-") + std::to_string(gid.station());
 
    std::string strLa = std::to_string(gid.layer());
    std::string strCh = std::to_string(gid.chamber());
 
    string hName = "recHit_" + strIdxName;
    string hTitle = "recHit " + strIdxTitle;
    hTitle += ";VFAT;";
    recHitME_[gid] = ibooker.book1D(hName, hTitle, 24, 0, 24);
 
    string hName_2 = "VFAT_vs_ClusterSize_" + strIdxName;
    string hTitle_2 = "VFAT vs ClusterSize " + strIdxTitle;
    hTitle_2 += ";Cluster size;VFAT";
    VFAT_vs_ClusterSize_[gid] =
        ibooker.book2D(hName_2, hTitle_2, nClusterSizeBinNum_, 1, nClusterSizeBinNum_ + 1, 24, 0, 24);
 
    string hName_fired = "StripFired_" + strIdxName;
    string hTitle_fired = "StripsFired " + strIdxTitle;
    hTitle_fired += ";Strip;iEta";
    StripsFired_vs_eta_[gid] = ibooker.book2D(hName_fired, hTitle_fired, 384, 1, 385, 8, 1, 9);
 
    string hName_rh = "recHit_x_" + strIdxName;
    string hTitle_rh = "recHit local x " + strIdxTitle;
    hTitle_rh += ";Local x (cm);iEta";
    rh_vs_eta_[gid] = ibooker.book2D(hName_rh, hTitle_rh, 50, -25, 25, 8, 1, 9);
 
    GEMDetId lid(gid.region(), gid.ring(), gid.station(), gid.layer(), 0, 0);
    Int_t nIdxOcc = 0;
 
    for (; nIdxOcc < (Int_t)listLayerOcc.size(); nIdxOcc++)
      if (listLayerOcc[nIdxOcc] == lid)
        break;
    if (nIdxOcc >= (Int_t)listLayerOcc.size())
      listLayerOcc.push_back(lid);
  }
 
  for (auto gid : listLayerOcc) {
    std::string strIdxName =
        std::string("GE") + (gid.region() > 0 ? "p" : "m") + to_string(gid.station()) + "_" + to_string(gid.layer());
    std::string strIdxTitle =
        std::string("GE") + (gid.region() > 0 ? "+" : "-") + to_string(gid.station()) + "/" + to_string(gid.layer());
 
    string hName_rh = "recHit_globalPos_Gemini_" + strIdxName;
    string hTitle_rh = "recHit global position Gemini chamber : " + strIdxTitle;
    hTitle_rh += ";Global X (cm);Global Y (cm)";
    recGlobalPos[gid] = ibooker.book2D(hName_rh, hTitle_rh, 100, fGlobXMin_, fGlobXMax_, 100, fGlobYMin_, fGlobYMax_);
  }
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), GEMDetId::chamber(), cmsLHEtoEOSManager::l, GEMDetId::layer(), GEMDetId::region(), GEMDetId::ring(), dqm::implementation::NavigatorBase::setCurrentFolder(), GEMDetId::station(), and AlCaHLTBitMon_QueryRunRegistry::string.

dqmBeginRun()

void GEMDQMSource::dqmBeginRun ( edm::Run const &  , edm::EventSetup const &    )

inlineoverrideprotected

Definition at line 32 of file GEMDQMSource.cc.

{};

fillDescriptions()

void GEMDQMSource::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 98 of file GEMDQMSource.cc.

                                                                              {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("recHitsInputLabel", edm::InputTag("gemRecHits", ""));
 
  desc.add<int>("idxFirstStrip", 0);
  desc.add<int>("ClusterSizeBinNum", 9);
 
  desc.add<double>("global_x_bound_min", -350);
  desc.add<double>("global_x_bound_max", 350);
  desc.add<double>("global_y_bound_min", -260);
  desc.add<double>("global_y_bound_max", 260);
 
  descriptions.add("GEMDQMSource", desc);
}

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

findVFAT()

int GEMDQMSource::findVFAT ( float  min_, float  max_, float  x_, int  roll_  )

private

Definition at line 62 of file GEMDQMSource.cc.

                                                                      {
  float step = abs(max_ - min_) / 3.0;
  if (x_ < (min(min_, max_) + step)) {
    return 8 - roll_;
  } else if (x_ < (min(min_, max_) + 2.0 * step)) {
    return 16 - roll_;
  } else {
    return 24 - roll_;
  }
}

References funct::abs(), and min().

initGeometry()

const GEMGeometry * GEMDQMSource::initGeometry ( edm::EventSetup const &  iSetup )

private

Definition at line 73 of file GEMDQMSource.cc.

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

References MillePedeFileConverter_cfg::e, edm::EventSetup::get(), and get.

Member Data Documentation

fGlobXMax_

float GEMDQMSource::fGlobXMax_

private

Definition at line 39 of file GEMDQMSource.cc.

fGlobXMin_

float GEMDQMSource::fGlobXMin_

private

Definition at line 39 of file GEMDQMSource.cc.

fGlobYMax_

float GEMDQMSource::fGlobYMax_

private

Definition at line 40 of file GEMDQMSource.cc.

fGlobYMin_

float GEMDQMSource::fGlobYMin_

private

Definition at line 40 of file GEMDQMSource.cc.

gemChambers_

std::vector<GEMChamber> GEMDQMSource::gemChambers_

private

Definition at line 50 of file GEMDQMSource.cc.

GEMGeometry_

const GEMGeometry* GEMDQMSource::GEMGeometry_

private

Definition at line 48 of file GEMDQMSource.cc.

nClusterSizeBinNum_

int GEMDQMSource::nClusterSizeBinNum_

private

Definition at line 43 of file GEMDQMSource.cc.

nIdxFirstStrip_

int GEMDQMSource::nIdxFirstStrip_

private

Definition at line 42 of file GEMDQMSource.cc.

recGlobalPos

std::unordered_map<UInt_t, MonitorElement*> GEMDQMSource::recGlobalPos

private

Definition at line 56 of file GEMDQMSource.cc.

recHitME_

std::unordered_map<UInt_t, MonitorElement*> GEMDQMSource::recHitME_

private

Definition at line 52 of file GEMDQMSource.cc.

rh_vs_eta_

std::unordered_map<UInt_t, MonitorElement*> GEMDQMSource::rh_vs_eta_

private

Definition at line 55 of file GEMDQMSource.cc.

StripsFired_vs_eta_

std::unordered_map<UInt_t, MonitorElement*> GEMDQMSource::StripsFired_vs_eta_

private

Definition at line 54 of file GEMDQMSource.cc.

tagRecHit_

edm::EDGetToken GEMDQMSource::tagRecHit_

private

Definition at line 37 of file GEMDQMSource.cc.

VFAT_vs_ClusterSize_

std::unordered_map<UInt_t, MonitorElement*> GEMDQMSource::VFAT_vs_ClusterSize_

private

Definition at line 53 of file GEMDQMSource.cc.