Inheritance diagram for GEMGeometryValidate:

Public Member Functions
	GEMGeometryValidate (const ParameterSet &)

	~GEMGeometryValidate () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	clearData ()

void	clearData2 ()

void	compareShape (const GeomDet , const float )

void	compareTransform (const GlobalPoint &, const TGeoMatrix *)

void	endJob () override

float	getDiff (const float, const float)

float	getDistance (const GlobalPoint &, const GlobalPoint &)

void	makeHistogram (const string &, vector< float > &)

void	makeHistograms (const char *)

void	makeHistograms2 (const char *)

void	validateGEMChamberGeometry ()

void	validateGEMEtaPartitionGeometry ()

Private Attributes
vector< float >	bottomWidths_

FWGeometry	fwGeometry_

edm::ESHandle< GEMGeometry >	gemGeometry_

vector< float >	globalDistances_

string	infileName_

vector< float >	lengths_

vector< float >	nstrips_

TFile *	outFile_

string	outfileName_

vector< float >	pitch_

vector< float >	stripslen_

vector< float >	thicknesses_

int	tolerance_

vector< float >	topWidths_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 50 of file GEMGeometryValidate.cc.

Constructor & Destructor Documentation

◆ GEMGeometryValidate()

GEMGeometryValidate::GEMGeometryValidate ( const ParameterSet & iConfig )

explicit

Definition at line 103 of file GEMGeometryValidate.cc.

     : infileName_(iConfig.getUntrackedParameter<string>("infileName", "cmsRecoGeom-2021.root")),
       outfileName_(iConfig.getUntrackedParameter<string>("outfileName", "validateGEMGeometry.root")),
       tolerance_(iConfig.getUntrackedParameter<int>("tolerance", 6)) {
   fwGeometry_.loadMap(infileName_.c_str());
   outFile_ = TFile::Open(outfileName_.c_str(), "RECREATE");
 }

References fwGeometry_, infileName_, FWGeometry::loadMap(), outFile_, and outfileName_.

◆ ~GEMGeometryValidate()

GEMGeometryValidate::~GEMGeometryValidate ( )

inlineoverride

Definition at line 53 of file GEMGeometryValidate.cc.

53 {}

Member Function Documentation

◆ analyze()

void GEMGeometryValidate::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 111 of file GEMGeometryValidate.cc.

                                                                                         {
   eventSetup.get<MuonGeometryRecord>().get(gemGeometry_);
  
   if (gemGeometry_.isValid()) {
     LogVerbatim("GEMGeometry") << "Validating GEM chamber geometry";
  
     validateGEMChamberGeometry();
  
     validateGEMEtaPartitionGeometry();
  
   } else
     LogVerbatim("GEMGeometry") << "Invalid GEM geometry";
 }

References gemGeometry_, edm::EventSetup::get(), get, edm::ESHandleBase::isValid(), validateGEMChamberGeometry(), and validateGEMEtaPartitionGeometry().

◆ beginJob()

void GEMGeometryValidate::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 302 of file GEMGeometryValidate.cc.

302 { outFile_->cd(); }

References outFile_.

◆ clearData()

void GEMGeometryValidate::clearData ( )

inlineprivate

Definition at line 73 of file GEMGeometryValidate.cc.

                    {
     globalDistances_.clear();
     topWidths_.clear();
     bottomWidths_.clear();
     lengths_.clear();
     thicknesses_.clear();
   }

Referenced by validateGEMChamberGeometry().

◆ clearData2()

void GEMGeometryValidate::clearData2 ( )

inlineprivate

Definition at line 81 of file GEMGeometryValidate.cc.

                     {
     nstrips_.clear();
     pitch_.clear();
     stripslen_.clear();
   }

Referenced by validateGEMEtaPartitionGeometry().

◆ compareShape()

void GEMGeometryValidate::compareShape	(	const GeomDet *	det,
		const float *	shape
	)

private

Definition at line 188 of file GEMGeometryValidate.cc.

                                                                              {
   float shapeTopWidth;
   float shapeBottomWidth;
   float shapeLength;
   float shapeThickness;
  
   if (shape[0] == 1) {
     shapeTopWidth = shape[2];
     shapeBottomWidth = shape[1];
     shapeLength = shape[4];
     shapeThickness = shape[3];
   } else if (shape[0] == 2) {
     shapeTopWidth = shape[1];
     shapeBottomWidth = shape[1];
     shapeLength = shape[2];
     shapeThickness = shape[3];
   } else {
     LogVerbatim("GEMGeometry") << "Failed to get box or trapezoid from shape";
  
     return;
   }
  
   float topWidth, bottomWidth;
   float length, thickness;
  
   const Bounds* bounds = &(det->surface().bounds());
   if (const TrapezoidalPlaneBounds* tpbs = dynamic_cast<const TrapezoidalPlaneBounds*>(bounds)) {
     array<const float, 4> const& ps = tpbs->parameters();
  
     assert(ps.size() == 4);
  
     bottomWidth = ps[0];
     topWidth = ps[1];
     thickness = ps[2];
     length = ps[3];
   } else if ((dynamic_cast<const RectangularPlaneBounds*>(bounds))) {
     length = det->surface().bounds().length() * 0.5;
     topWidth = det->surface().bounds().width() * 0.5;
     bottomWidth = topWidth;
     thickness = det->surface().bounds().thickness() * 0.5;
   } else {
     LogVerbatim("GEMGeometry") << "Failed to get bounds";
  
     return;
   }
   topWidths_.push_back(fabs(shapeTopWidth - topWidth));
   bottomWidths_.push_back(fabs(shapeBottomWidth - bottomWidth));
   lengths_.push_back(fabs(shapeLength - length));
   thicknesses_.push_back(fabs(shapeThickness - thickness));
 }

References cms::cuda::assert(), bottomWidths_, Surface::bounds(), Bounds::length(), lengths_, GeomDet::surface(), Bounds::thickness(), Calorimetry_cff::thickness, thicknesses_, topWidths_, and Bounds::width().

Referenced by validateGEMChamberGeometry().

◆ compareTransform()

void GEMGeometryValidate::compareTransform	(	const GlobalPoint &	gp,
		const TGeoMatrix *	matrix
	)

private

Definition at line 176 of file GEMGeometryValidate.cc.

                                                                                           {
   double local[3] = {0.0, 0.0, 0.0};
   double global[3];
  
   matrix->LocalToMaster(local, global);
  
   float distance = getDistance(GlobalPoint(global[0], global[1], global[2]), gp);
   if (abs(distance) < 1.0e-7)
     distance = 0.0;  // set a tollerance for the distance inside Histos
   globalDistances_.push_back(distance);
 }

References funct::abs(), HLT_2018_cff::distance, MillePedeFileConverter_cfg::e, getDistance(), globalDistances_, runTauDisplay::gp, DTRecHitClients_cfi::local, and makeMuonMisalignmentScenario::matrix.

Referenced by validateGEMChamberGeometry().

◆ endJob()

void GEMGeometryValidate::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 304 of file GEMGeometryValidate.cc.

                                  {
   LogVerbatim("GEMGeometry") << "Done.";
   LogVerbatim("GEMGeometry") << "Results written to " << outfileName_;
   outFile_->Close();
 }

References outFile_, and outfileName_.

◆ getDiff()

float GEMGeometryValidate::getDiff	(	const float	val1,
		const float	val2
	)

private

Definition at line 244 of file GEMGeometryValidate.cc.

                                                                      {
   if (almost_equal(val1, val2, tolerance_))
     return 0.0f;
   else
     return (val1 - val2);
 }

References almost_equal(), and tolerance_.

◆ getDistance()

float GEMGeometryValidate::getDistance	(	const GlobalPoint &	p1,
		const GlobalPoint &	p2
	)

private

Definition at line 239 of file GEMGeometryValidate.cc.

                                                                                    {
   return sqrt((p1.x() - p2.x()) * (p1.x() - p2.x()) + (p1.y() - p2.y()) * (p1.y() - p2.y()) +
               (p1.z() - p2.z()) * (p1.z() - p2.z()));
 }

References p1, p2, and mathSSE::sqrt().

Referenced by compareTransform().

◆ makeHistogram()

void GEMGeometryValidate::makeHistogram	(	const string &	name,
		vector< float > &	data
	)

private

Definition at line 287 of file GEMGeometryValidate.cc.

                                                                                {
   if (data.empty())
     return;
  
   const auto [minE, maxE] = minmax_element(begin(data), end(data));
  
   TH1D hist(name.c_str(), name.c_str(), 100, *minE * (1 + 0.10), *maxE * (1 + 0.10));
  
   for (auto const& it : data)
     hist.Fill(it);
  
   hist.GetXaxis()->SetTitle("[cm]");
   hist.Write();
 }

References begin, data, end, compare::hist, HLT_2018_cff::maxE, and Skims_PA_cff::name.

Referenced by makeHistograms(), and makeHistograms2().

◆ makeHistograms()

void GEMGeometryValidate::makeHistograms ( const char * detector )

private

Definition at line 251 of file GEMGeometryValidate.cc.

                                                              {
   outFile_->cd();
  
   string d(detector);
  
   string gdn = d + ": distance between points in global coordinates";
   makeHistogram(gdn, globalDistances_);
  
   string twn = d + ": absolute difference between top widths (along X)";
   makeHistogram(twn, topWidths_);
  
   string bwn = d + ": absolute difference between bottom widths (along X)";
   makeHistogram(bwn, bottomWidths_);
  
   string ln = d + ": absolute difference between lengths (along Y)";
   makeHistogram(ln, lengths_);
  
   string tn = d + ": absolute difference between thicknesses (along Z)";
   makeHistogram(tn, thicknesses_);
 }

References bottomWidths_, ztail::d, hgcalTestNeighbor_cfi::detector, globalDistances_, lengths_, makeHistogram(), outFile_, thicknesses_, and topWidths_.

Referenced by validateGEMChamberGeometry().

◆ makeHistograms2()

void GEMGeometryValidate::makeHistograms2 ( const char * detector )

private

Definition at line 272 of file GEMGeometryValidate.cc.

                                                               {
   outFile_->cd();
  
   string d(detector);
  
   string ns = d + ": absolute difference between nStrips";
   makeHistogram(ns, nstrips_);
  
   string pi = d + ": absolute difference between Strips Pitch";
   makeHistogram(pi, pitch_);
  
   string pl = d + ": absolute difference between Strips Length";
   makeHistogram(pl, stripslen_);
 }

References ztail::d, hgcalTestNeighbor_cfi::detector, makeHistogram(), nstrips_, outFile_, pi, pitch_, and stripslen_.

Referenced by validateGEMEtaPartitionGeometry().

◆ validateGEMChamberGeometry()

void GEMGeometryValidate::validateGEMChamberGeometry ( )

private

Definition at line 125 of file GEMGeometryValidate.cc.

                                                      {
   clearData();
  
   for (auto const& it : gemGeometry_->chambers()) {
     GEMDetId chId = it->id();
     GlobalPoint gp = it->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
  
     const TGeoMatrix* matrix = fwGeometry_.getMatrix(chId.rawId());
  
     if (!matrix) {
       LogVerbatim("GEMGeometry") << "Failed to get matrix of GEM chamber with detid: " << chId.rawId();
       continue;
     }
     compareTransform(gp, matrix);
  
     auto const& shape = fwGeometry_.getShapePars(chId.rawId());
  
     if (!shape) {
       LogVerbatim("GEMGeometry") << "Failed to get shape of GEM chamber with detid: " << chId.rawId();
       continue;
     }
     compareShape(it, shape);
   }
   makeHistograms("GEM Chamber");
 }

References GEMGeometry::chambers(), clearData(), compareShape(), compareTransform(), fwGeometry_, gemGeometry_, FWGeometry::getMatrix(), FWGeometry::getShapePars(), runTauDisplay::gp, makeHistograms(), makeMuonMisalignmentScenario::matrix, and DetId::rawId().

Referenced by analyze().

◆ validateGEMEtaPartitionGeometry()

void GEMGeometryValidate::validateGEMEtaPartitionGeometry ( )

private

Definition at line 151 of file GEMGeometryValidate.cc.

                                                           {
   clearData2();
  
   for (auto const& it : gemGeometry_->etaPartitions()) {
     GEMDetId chId = it->id();
     const int n_strips = it->nstrips();
     const float n_pitch = it->pitch();
     const StripTopology& topo = it->specificTopology();
     const float stripLen = topo.stripLength();
     const float* parameters = fwGeometry_.getParameters(chId.rawId());
     nstrips_.push_back(abs(n_strips - parameters[0]));
  
     if (n_strips) {
       for (int istrips = 1; istrips <= n_strips; istrips++) {
         pitch_.push_back(fabs(n_pitch - parameters[2]));
  
         stripslen_.push_back(fabs(stripLen - parameters[1]));
       }
     } else {
       LogVerbatim("GEMGeometry") << "ATTENTION! nStrips == 0";
     }
   }
   makeHistograms2("GEM Eta Partition");
 }

References funct::abs(), clearData2(), GEMGeometry::etaPartitions(), fwGeometry_, gemGeometry_, FWGeometry::getParameters(), makeHistograms2(), nstrips_, pitch_, DetId::rawId(), StripTopology::stripLength(), and stripslen_.

Referenced by analyze().

Member Data Documentation

◆ bottomWidths_

vector<float> GEMGeometryValidate::bottomWidths_

private

Definition at line 92 of file GEMGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ fwGeometry_

FWGeometry GEMGeometryValidate::fwGeometry_

private

Definition at line 88 of file GEMGeometryValidate.cc.

Referenced by GEMGeometryValidate(), validateGEMChamberGeometry(), and validateGEMEtaPartitionGeometry().

◆ gemGeometry_

edm::ESHandle<GEMGeometry> GEMGeometryValidate::gemGeometry_

private

Definition at line 87 of file GEMGeometryValidate.cc.

Referenced by analyze(), validateGEMChamberGeometry(), and validateGEMEtaPartitionGeometry().

◆ globalDistances_

vector<float> GEMGeometryValidate::globalDistances_

private

Definition at line 90 of file GEMGeometryValidate.cc.

Referenced by compareTransform(), and makeHistograms().

◆ infileName_

string GEMGeometryValidate::infileName_

private

Definition at line 98 of file GEMGeometryValidate.cc.

Referenced by GEMGeometryValidate().

◆ lengths_

vector<float> GEMGeometryValidate::lengths_

private

Definition at line 93 of file GEMGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ nstrips_

vector<float> GEMGeometryValidate::nstrips_

private

Definition at line 95 of file GEMGeometryValidate.cc.

Referenced by makeHistograms2(), and validateGEMEtaPartitionGeometry().

◆ outFile_

TFile* GEMGeometryValidate::outFile_

private

Definition at line 89 of file GEMGeometryValidate.cc.

Referenced by beginJob(), endJob(), GEMGeometryValidate(), makeHistograms(), and makeHistograms2().

◆ outfileName_

string GEMGeometryValidate::outfileName_

private

Definition at line 99 of file GEMGeometryValidate.cc.

Referenced by endJob(), and GEMGeometryValidate().

◆ pitch_

vector<float> GEMGeometryValidate::pitch_

private

Definition at line 96 of file GEMGeometryValidate.cc.

Referenced by makeHistograms2(), and validateGEMEtaPartitionGeometry().

◆ stripslen_

vector<float> GEMGeometryValidate::stripslen_

private

Definition at line 97 of file GEMGeometryValidate.cc.

Referenced by makeHistograms2(), and validateGEMEtaPartitionGeometry().

◆ thicknesses_

vector<float> GEMGeometryValidate::thicknesses_

private

Definition at line 94 of file GEMGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ tolerance_

int GEMGeometryValidate::tolerance_

private

Definition at line 100 of file GEMGeometryValidate.cc.

Referenced by getDiff().

◆ topWidths_

vector<float> GEMGeometryValidate::topWidths_

private

Definition at line 91 of file GEMGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ GEMGeometryValidate()

◆ ~GEMGeometryValidate()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ clearData()

◆ clearData2()

◆ compareShape()

◆ compareTransform()

◆ endJob()

◆ getDiff()

◆ getDistance()

◆ makeHistogram()

◆ makeHistograms()

◆ makeHistograms2()

◆ validateGEMChamberGeometry()

◆ validateGEMEtaPartitionGeometry()

Member Data Documentation

◆ bottomWidths_

◆ fwGeometry_

◆ gemGeometry_

◆ globalDistances_

◆ infileName_

◆ lengths_

◆ nstrips_

◆ outFile_

◆ outfileName_

◆ pitch_

◆ stripslen_

◆ thicknesses_

◆ tolerance_

◆ topWidths_