Inheritance diagram for CSCGeometryValidate:

Public Member Functions
	CSCGeometryValidate (const ParameterSet &)

	~CSCGeometryValidate () override

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	validateCSCChamberGeometry ()

void	validateCSCLayerGeometry ()

Private Attributes
vector< float >	angularWidth_

vector< float >	bottomWidths_

vector< float >	centreToIntersection_

const CSCGeometry *	cscGeometry_

FWGeometry	fwGeometry_

vector< float >	globalDistances_

string	infileName_

vector< float >	lengths_

TFile *	outFile_

string	outfileName_

vector< float >	phiOfOneEdge_

vector< float >	sOffset_

vector< float >	thicknesses_

const edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	tokCSC_

int	tolerance_

vector< float >	topWidths_

vector< float >	wireAngle_

vector< float >	wireSpacing_

vector< float >	yAxisOrientation_

vector< float >	yCentreOfStripPlane_

Detailed Description

Definition at line 52 of file CSCGeometryValidate.cc.

Constructor & Destructor Documentation

◆ CSCGeometryValidate()

CSCGeometryValidate::CSCGeometryValidate ( const ParameterSet & iConfig )

explicit

Definition at line 120 of file CSCGeometryValidate.cc.

     : tokCSC_{esConsumes<CSCGeometry, MuonGeometryRecord>(edm::ESInputTag{})},
       infileName_(iConfig.getUntrackedParameter<string>("infileName", "cmsGeom10.root")),
       outfileName_(iConfig.getUntrackedParameter<string>("outfileName", "validateCSCGeometry.root")),
       tolerance_(iConfig.getUntrackedParameter<int>("tolerance", 6)) {
   fwGeometry_.loadMap(infileName_.c_str());
   outFile_ = TFile::Open(outfileName_.c_str(), "RECREATE");
 }

◆ ~CSCGeometryValidate()

CSCGeometryValidate::~CSCGeometryValidate ( )

inlineoverride

Definition at line 55 of file CSCGeometryValidate.cc.

55 {}

Member Function Documentation

◆ analyze()

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

overrideprivate

Definition at line 129 of file CSCGeometryValidate.cc.

References cscGeometry_, options_cfi::eventSetup, tokCSC_, validateCSCChamberGeometry(), and validateCSCLayerGeometry().

                                                                                         {
   cscGeometry_ = &eventSetup.getData(tokCSC_);
   LogVerbatim("CSCGeometry") << "Validating CSC chamber geometry";
   validateCSCChamberGeometry();
   validateCSCLayerGeometry();
 }

◆ beginJob()

void CSCGeometryValidate::beginJob ( void )

overrideprivate

Definition at line 348 of file CSCGeometryValidate.cc.

References outFile_.

348 { outFile_->cd(); }

CSCGeometryValidate::outFile_

TFile * outFile_

Definition: CSCGeometryValidate.cc:97

◆ clearData()

void CSCGeometryValidate::clearData ( )

inlineprivate

Definition at line 75 of file CSCGeometryValidate.cc.

Referenced by validateCSCChamberGeometry().

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

◆ clearData2()

void CSCGeometryValidate::clearData2 ( )

inlineprivate

Definition at line 83 of file CSCGeometryValidate.cc.

Referenced by validateCSCLayerGeometry().

                     {
     yAxisOrientation_.clear();
     sOffset_.clear();
     yCentreOfStripPlane_.clear();
     angularWidth_.clear();
     centreToIntersection_.clear();
     phiOfOneEdge_.clear();
     wireSpacing_.clear();
     wireAngle_.clear();
   }

◆ compareShape()

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

private

Definition at line 219 of file CSCGeometryValidate.cc.

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

Referenced by validateCSCChamberGeometry().

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

◆ compareTransform()

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

private

Definition at line 207 of file CSCGeometryValidate.cc.

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

Referenced by validateCSCChamberGeometry().

                                                                                           {
   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 ((distance >= 0.0) && (distance < 1.0e-7))
     distance = 0.0;  // set a tollerance for the distance inside Histos
   globalDistances_.push_back(distance);
 }

◆ endJob()

void CSCGeometryValidate::endJob ( void )

overrideprivate

Definition at line 350 of file CSCGeometryValidate.cc.

References outFile_, and outfileName_.

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

◆ getDiff()

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

private

Definition at line 275 of file CSCGeometryValidate.cc.

References almost_equal(), and tolerance_.

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

◆ getDistance()

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

private

Definition at line 270 of file CSCGeometryValidate.cc.

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

Referenced by compareTransform().

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

◆ makeHistogram()

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

private

Definition at line 333 of file CSCGeometryValidate.cc.

References data, mps_fire::end, compareTotals::hist, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HLT_2024v14_cff::maxE, and Skims_PA_cff::name.

Referenced by makeHistograms(), and makeHistograms2().

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

◆ makeHistograms()

void CSCGeometryValidate::makeHistograms ( const char * detector )

private

Definition at line 282 of file CSCGeometryValidate.cc.

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

Referenced by validateCSCChamberGeometry().

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

◆ makeHistograms2()

void CSCGeometryValidate::makeHistograms2 ( const char * detector )

private

Definition at line 303 of file CSCGeometryValidate.cc.

References angularWidth_, centreToIntersection_, ztail::d, hgcalTestNeighbor_cfi::detector, makeHistogram(), outFile_, phiOfOneEdge_, pi, sOffset_, wireAngle_, wireSpacing_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::ws, yAxisOrientation_, and yCentreOfStripPlane_.

Referenced by validateCSCLayerGeometry().

                                                               {
   outFile_->cd();
 
   string d(detector);
 
   string ns = d + ": absolute difference between Y Axis Orientation of the Strips";
   makeHistogram(ns, yAxisOrientation_);
 
   string pi = d + ": absolute difference between Strips Offset";
   makeHistogram(pi, sOffset_);
 
   string pl = d + ": absolute difference between 'Y centre' of the Strips Planes";
   makeHistogram(pl, yCentreOfStripPlane_);
 
   string aw = d + ": absolute difference between 'angular width' of the Strips ";
   makeHistogram(aw, angularWidth_);
 
   string ci = d + ": absolute difference between 'centre to intersection' of the Strips ";
   makeHistogram(ci, centreToIntersection_);
 
   string po = d + ": absolute difference between 'phi of one edge' of the Strips ";
   makeHistogram(po, phiOfOneEdge_);
 
   string ws = d + ": absolute difference between 'wire spacing' of the Wires ";
   makeHistogram(ws, wireSpacing_);
 
   string wa = d + ": absolute difference between 'wire angle' of the Wires ";
   makeHistogram(wa, wireAngle_);
 }

◆ validateCSCChamberGeometry()

void CSCGeometryValidate::validateCSCChamberGeometry ( )

private

Definition at line 136 of file CSCGeometryValidate.cc.

References CSCGeometry::chambers(), clearData(), compareShape(), compareTransform(), cscGeometry_, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getShapePars(), runTauDisplay::gp, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistograms(), makeMuonMisalignmentScenario::matrix, DetId::rawId(), and l1trig_cff::shape.

Referenced by analyze().

                                                      {
   clearData();
 
   for (auto const& it : cscGeometry_->chambers()) {
     CSCDetId 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("CSCGeometry") << "Failed to get matrix of CSC chamber with detid: " << chId.rawId();
       continue;
     }
     compareTransform(gp, matrix);
 
     auto const& shape = fwGeometry_.getShapePars(chId.rawId());
 
     if (!shape) {
       LogVerbatim("CSCGeometry") << "Failed to get shape of CSC chamber with detid: " << chId.rawId();
       continue;
     }
     compareShape(it, shape);
   }
   makeHistograms("CSC Chamber");
 }

◆ validateCSCLayerGeometry()

void CSCGeometryValidate::validateCSCLayerGeometry ( )

private

Definition at line 162 of file CSCGeometryValidate.cc.

References CSCRadialStripTopology::angularWidth(), angularWidth_, CSCRadialStripTopology::centreToIntersection(), centreToIntersection_, clearData2(), cscGeometry_, fwGeometry_, FWGeometry::getParameters(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, CSCGeometry::layers(), makeHistograms2(), CSCLayerGeometry::numberOfStrips(), CSCLayerGeometry::numberOfWires(), CSCRadialStripTopology::phiOfOneEdge(), phiOfOneEdge_, DetId::rawId(), sOffset_, CSCLayerGeometry::stripOffset(), CSCLayerGeometry::topology(), CSCWireTopology::wireAngle(), wireAngle_, CSCWireTopology::wireSpacing(), wireSpacing_, CSCLayerGeometry::wireTopology(), CSCRadialStripTopology::yAxisOrientation(), yAxisOrientation_, CSCRadialStripTopology::yCentreOfStripPlane(), and yCentreOfStripPlane_.

Referenced by analyze().

                                                    {
   clearData2();
 
   for (auto const& it : cscGeometry_->layers()) {
     CSCDetId chId = it->id();
     const CSCLayerGeometry* laygeo = it->geometry();
     const int n_strips = laygeo->numberOfStrips();
     const int n_wire = laygeo->numberOfWires();
     const float strips_offset = laygeo->stripOffset();
     const CSCStripTopology* stopo = laygeo->topology();
     const float ycentre_of_strip_plane = stopo->yCentreOfStripPlane();
     const float angular_width = stopo->angularWidth();
     const float y_axis_orientation = stopo->yAxisOrientation();
     const float centre_to_intersection = stopo->centreToIntersection();
     const float phi_of_one_edge = stopo->phiOfOneEdge();
     const float* parameters = fwGeometry_.getParameters(chId.rawId());
     const CSCWireTopology* wiretopo = laygeo->wireTopology();
     const double wire_spacing = wiretopo->wireSpacing();
     const float wire_angle = wiretopo->wireAngle();
 
     if (n_strips) {
       for (int istrips = 1; istrips <= n_strips; istrips++) {
         yAxisOrientation_.push_back(fabs(y_axis_orientation - parameters[0]));
         sOffset_.push_back(fabs(strips_offset - parameters[4]));
         yCentreOfStripPlane_.push_back(fabs(ycentre_of_strip_plane - parameters[2]));
         angularWidth_.push_back(fabs(angular_width - parameters[5]));
         centreToIntersection_.push_back(fabs(centre_to_intersection - parameters[1]));
         phiOfOneEdge_.push_back(fabs(phi_of_one_edge - parameters[3]));
       }
     } else {
       LogVerbatim("CSCGeometry") << "ATTENTION! nStrips == 0";
     }
 
     if (n_wire) {
       for (int iwires = 1; iwires <= n_wire; iwires++) {
         wireSpacing_.push_back(fabs(wire_spacing - parameters[6]));
         wireAngle_.push_back(fabs(wire_angle - parameters[7]));
       }
     } else {
       LogVerbatim("CSCGeometry") << "ATTENTION! nWires == 0";
     }
   }
   makeHistograms2("CSC Layer");
 }

Member Data Documentation

◆ angularWidth_

vector<float> CSCGeometryValidate::angularWidth_

private

Definition at line 108 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ bottomWidths_

vector<float> CSCGeometryValidate::bottomWidths_

private

Definition at line 101 of file CSCGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ centreToIntersection_

vector<float> CSCGeometryValidate::centreToIntersection_

private

Definition at line 109 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ cscGeometry_

const CSCGeometry* CSCGeometryValidate::cscGeometry_

private

Definition at line 95 of file CSCGeometryValidate.cc.

Referenced by analyze(), validateCSCChamberGeometry(), and validateCSCLayerGeometry().

◆ fwGeometry_

FWGeometry CSCGeometryValidate::fwGeometry_

private

Definition at line 96 of file CSCGeometryValidate.cc.

Referenced by validateCSCChamberGeometry(), and validateCSCLayerGeometry().

◆ globalDistances_

vector<float> CSCGeometryValidate::globalDistances_

private

Definition at line 99 of file CSCGeometryValidate.cc.

Referenced by compareTransform(), and makeHistograms().

◆ infileName_

string CSCGeometryValidate::infileName_

private

Definition at line 115 of file CSCGeometryValidate.cc.

◆ lengths_

vector<float> CSCGeometryValidate::lengths_

private

Definition at line 102 of file CSCGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ outFile_

TFile* CSCGeometryValidate::outFile_

private

Definition at line 97 of file CSCGeometryValidate.cc.

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

◆ outfileName_

string CSCGeometryValidate::outfileName_

private

Definition at line 116 of file CSCGeometryValidate.cc.

Referenced by endJob().

◆ phiOfOneEdge_

vector<float> CSCGeometryValidate::phiOfOneEdge_

private

Definition at line 110 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ sOffset_

vector<float> CSCGeometryValidate::sOffset_

private

Definition at line 106 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ thicknesses_

vector<float> CSCGeometryValidate::thicknesses_

private

Definition at line 103 of file CSCGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ tokCSC_

const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> CSCGeometryValidate::tokCSC_

private

Definition at line 94 of file CSCGeometryValidate.cc.

Referenced by analyze().

◆ tolerance_

int CSCGeometryValidate::tolerance_

private

Definition at line 117 of file CSCGeometryValidate.cc.

Referenced by getDiff().

◆ topWidths_

vector<float> CSCGeometryValidate::topWidths_

private

Definition at line 100 of file CSCGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

◆ wireAngle_

vector<float> CSCGeometryValidate::wireAngle_

private

Definition at line 113 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ wireSpacing_

vector<float> CSCGeometryValidate::wireSpacing_

private

Definition at line 112 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ yAxisOrientation_

vector<float> CSCGeometryValidate::yAxisOrientation_

private

Definition at line 105 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

◆ yCentreOfStripPlane_

vector<float> CSCGeometryValidate::yCentreOfStripPlane_

private

Definition at line 107 of file CSCGeometryValidate.cc.

Referenced by makeHistograms2(), and validateCSCLayerGeometry().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CSCGeometryValidate()

◆ ~CSCGeometryValidate()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ clearData()

◆ clearData2()

◆ compareShape()

◆ compareTransform()

◆ endJob()

◆ getDiff()

◆ getDistance()

◆ makeHistogram()

◆ makeHistograms()

◆ makeHistograms2()

◆ validateCSCChamberGeometry()

◆ validateCSCLayerGeometry()

Member Data Documentation

◆ angularWidth_

◆ bottomWidths_

◆ centreToIntersection_

◆ cscGeometry_

◆ fwGeometry_

◆ globalDistances_

◆ infileName_

◆ lengths_

◆ outFile_

◆ outfileName_

◆ phiOfOneEdge_

◆ sOffset_

◆ thicknesses_

◆ tokCSC_

◆ tolerance_

◆ topWidths_

◆ wireAngle_

◆ wireSpacing_

◆ yAxisOrientation_

◆ yCentreOfStripPlane_