#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 101 of file GEMCheckGeometry.cc.

                                     {
  
 }

Member Function Documentation

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

override

Definition at line 106 of file GEMCheckGeometry.cc.

                                                                       {
 
 }

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

override

Definition at line 23 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(), mathSSE::return(), relativeConstraints::ring, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), relativeConstraints::station, digitizers_cfi::strip, theStdPlots, fftjetcommon_cfi::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);
                 }
                 
               }
             }
 
 
 }