#include <GEMCheckGeometry.h>

Inheritance diagram for GEMCheckGeometry:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	GEMCheckGeometry (const edm::ParameterSet &gc)

	~GEMCheckGeometry () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Attributes
bool	detailPlot_

double	GE11PhiBegin_

double	GE11PhiStep_

double	maxPhi_

double	minPhi_

std::map< UInt_t, MonitorElement * >	the_st_dphi

std::map< UInt_t, MonitorElement * >	theStdPlots

Detailed Description

Definition at line 13 of file GEMCheckGeometry.h.

Constructor & Destructor Documentation

GEMCheckGeometry::GEMCheckGeometry ( const edm::ParameterSet & gc )

explicit

Definition at line 12 of file GEMCheckGeometry.cc.

References detailPlot_, GE11PhiBegin_, GE11PhiStep_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), maxPhi_, and minPhi_.

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

GEMCheckGeometry::~GEMCheckGeometry ( )

override

Definition at line 89 of file GEMCheckGeometry.cc.

89 {}

Member Function Documentation

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

override

Definition at line 91 of file GEMCheckGeometry.cc.

91 {}

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

override

Definition at line 20 of file GEMCheckGeometry.cc.

References DQMStore::IBooker::book2D(), detailPlot_, MillePedeFileConverter_cfg::e, edm::EventSetup::get(), cond::hash, triggerObjects_cff::id, StripTopology::localPosition(), LogDebug, dataset::name, phi, GEMGeometry::regions(), relativeConstraints::ring, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), relativeConstraints::station, digitizers_cfi::strip, theStdPlots, runGCPTkAlMap::title, and ecaldqm::topology().

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