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 ( const edm::ParameterSet & cfg )

Definition at line 86 of file GEMDQMSource.cc.

References fGlobXMax_, fGlobXMin_, fGlobYMax_, fGlobYMin_, edm::ParameterSet::getParameter(), nClusterSizeBinNum_, nIdxFirstStrip_, and tagRecHit_.

                                                      {
   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");
 }

GEMDQMSource::~GEMDQMSource ( )

inlineoverride

Definition at line 28 of file GEMDQMSource.cc.

References fillDescriptions().

28 {};

Member Function Documentation

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

overrideprotected

Definition at line 196 of file GEMDQMSource.cc.

References DEFINE_FWK_MODULE, gemChambers_, GEMGeometry_, gemRecHits_cfi::gemRecHits, mps_fire::i, GEMGeometry::idToDet(), initGeometry(), edm::HandleBase::isValid(), GEMDetId::layer(), nClusterSizeBinNum_, nIdxFirstStrip_, recGlobalPos, recHitME_, GEMDetId::region(), rh_vs_eta_, GEMDetId::ring(), GEMDetId::roll(), GEMDetId::station(), StripsFired_vs_eta_, tagRecHit_, GeomDet::toGlobal(), VFAT_vs_ClusterSize_, PV3DBase< T, PVType, FrameType >::x(), hit::x, and PV3DBase< T, PVType, FrameType >::y().

Referenced by dqmBeginRun().

                                                                                  {
   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 (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());
       }
     }
   }
 }

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

overrideprotected

Definition at line 113 of file GEMDQMSource.cc.

References GEMDetId::chamber(), fGlobXMax_, fGlobXMin_, fGlobYMax_, fGlobYMin_, gemChambers_, GEMGeometry_, initGeometry(), cmsLHEtoEOSManager::l, GEMDetId::layer(), nClusterSizeBinNum_, recGlobalPos, recHitME_, GEMDetId::region(), rh_vs_eta_, GEMDetId::ring(), GEMDetId::station(), AlCaHLTBitMon_QueryRunRegistry::string, StripsFired_vs_eta_, GEMGeometry::superChambers(), and VFAT_vs_ClusterSize_.

Referenced by dqmBeginRun().

                                                                                                       {
   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 (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 (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_);
   }
 }

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

inlineoverrideprotected

Definition at line 32 of file GEMDQMSource.cc.

References analyze(), bookHistograms(), and MillePedeFileConverter_cfg::e.

32 {};

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

static

Definition at line 98 of file GEMDQMSource.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_2018_cff::InputTag.

Referenced by ~GEMDQMSource().

                                                                               {
   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);
 }

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

private

Definition at line 62 of file GEMDQMSource.cc.

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

                                                                       {
   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_;
   }
 }

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

private

Definition at line 73 of file GEMDQMSource.cc.

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

Referenced by analyze(), and bookHistograms().

                                                                          {
   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_;
 }