Inheritance diagram for ValidateGeometry:

Public Member Functions
	ValidateGeometry (const edm::ParameterSet &)

	~ValidateGeometry () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	clearData ()

void	compareShape (const GeomDet det, const float shape)

void	compareTransform (const GlobalPoint &point, const TGeoMatrix *matrix)

bool	dataEmpty ()

void	endJob () override

double	getDistance (const GlobalPoint &point1, const GlobalPoint &point2)

void	makeHistogram (const std::string &name, std::vector< double > &data)

void	makeHistograms (const char *detector)

void	validateCaloGeometry (DetId::Detector detector, int subdetector, const char *detname)

void	validateCSChamberGeometry (const int endcap, const char *detname)

void	validateCSCLayerGeometry (const int endcap, const char *detname)

void	validateDTChamberGeometry ()

void	validateDTLayerGeometry ()

void	validatePixelTopology (const TrackerGeometry::DetContainer &dets, const char *detname)

void	validateRPCGeometry (const int regionNumber, const char *regionName)

void	validateStripTopology (const TrackerGeometry::DetContainer &dets, const char *detname)

void	validateTrackerGeometry (const TrackerGeometry::DetContainer &dets, const char *detname)

Private Attributes
std::vector< double >	bottomWidths_

edm::ESHandle< CaloGeometry >	caloGeometry_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

edm::ESHandle< CSCGeometry >	cscGeometry_

edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	cscGeometryToken_

bool	doCalo_

bool	doMuon_

bool	doTracker_

edm::ESHandle< DTGeometry >	dtGeometry_

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeometryToken_

FWGeometry	fwGeometry_

std::vector< double >	globalDistances_

std::string	infileName_

std::vector< double >	lengths_

TFile *	outFile_

std::string	outfileName_

edm::ESHandle< RPCGeometry >	rpcGeometry_

edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeometryToken_

std::vector< double >	thicknesses_

std::vector< double >	topWidths_

edm::ESHandle< TrackerGeometry >	trackerGeometry_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 62 of file ValidateGeometry.cc.

Constructor & Destructor Documentation

◆ ValidateGeometry()

ValidateGeometry::ValidateGeometry ( const edm::ParameterSet & iConfig )

explicit

Definition at line 141 of file ValidateGeometry.cc.

References caloGeometryToken_, cscGeometryToken_, doCalo_, doMuon_, doTracker_, dtGeometryToken_, edm::EDConsumerBase::esConsumes(), fwGeometry_, edm::ParameterSet::getUntrackedParameter(), infileName_, FWGeometry::loadMap(), outFile_, outfileName_, rpcGeometryToken_, and trackerGeometryToken_.

     : infileName_(iConfig.getUntrackedParameter<std::string>("infileName")),
       outfileName_(iConfig.getUntrackedParameter<std::string>("outfileName")) {
   doTracker_ = iConfig.getUntrackedParameter<bool>("Tracker", true);
   doMuon_ = iConfig.getUntrackedParameter<bool>("Muon", true);
   doCalo_ = iConfig.getUntrackedParameter<bool>("Calo", true);
 
   if (doMuon_) {
     rpcGeometryToken_ = esConsumes();
     dtGeometryToken_ = esConsumes();
     cscGeometryToken_ = esConsumes();
   }
   if (doCalo_) {
     caloGeometryToken_ = esConsumes();
   }
   if (doTracker_) {
     trackerGeometryToken_ = esConsumes();
   }
 
   fwGeometry_.loadMap(infileName_.c_str());
 
   outFile_ = new TFile(outfileName_.c_str(), "RECREATE");
 }

◆ ~ValidateGeometry()

ValidateGeometry::~ValidateGeometry ( )

override

Definition at line 165 of file ValidateGeometry.cc.

165 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 167 of file ValidateGeometry.cc.

References DetId::Calo, caloGeometry_, caloGeometryToken_, gather_cfg::cout, cscGeometry_, cscGeometryToken_, TrackerGeometry::detsPXB(), TrackerGeometry::detsPXF(), TrackerGeometry::detsTEC(), TrackerGeometry::detsTIB(), TrackerGeometry::detsTID(), TrackerGeometry::detsTOB(), doCalo_, doMuon_, doTracker_, dtGeometry_, dtGeometryToken_, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, options_cfi::eventSetup, fwLog, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, edm::ESHandleBase::isValid(), fwlog::kWarning, rpcGeometry_, rpcGeometryToken_, HcalCastorDetId::SubdetectorId, HcalZDCDetId::SubdetectorId, trackerGeometry_, trackerGeometryToken_, validateCaloGeometry(), validateCSChamberGeometry(), validateCSCLayerGeometry(), validateDTChamberGeometry(), validateDTLayerGeometry(), validatePixelTopology(), validateRPCGeometry(), validateStripTopology(), and validateTrackerGeometry().

                                                                                      {
   if (doMuon_) {
     rpcGeometry_ = eventSetup.getHandle(rpcGeometryToken_);
 
     if (rpcGeometry_.isValid()) {
       std::cout << "Validating RPC -z endcap geometry" << std::endl;
       validateRPCGeometry(-1, "RPC -z endcap");
 
       std::cout << "Validating RPC +z endcap geometry" << std::endl;
       validateRPCGeometry(+1, "RPC +z endcap");
 
       std::cout << "Validating RPC barrel geometry" << std::endl;
       validateRPCGeometry(0, "RPC barrel");
     } else
       fwLog(fwlog::kWarning) << "Invalid RPC geometry" << std::endl;
 
     dtGeometry_ = eventSetup.getHandle(dtGeometryToken_);
 
     if (dtGeometry_.isValid()) {
       std::cout << "Validating DT chamber geometry" << std::endl;
       validateDTChamberGeometry();
 
       std::cout << "Validating DT layer geometry" << std::endl;
       validateDTLayerGeometry();
     } else
       fwLog(fwlog::kWarning) << "Invalid DT geometry" << std::endl;
 
     cscGeometry_ = eventSetup.getHandle(cscGeometryToken_);
 
     if (cscGeometry_.isValid()) {
       std::cout << "Validating CSC -z geometry" << std::endl;
       validateCSChamberGeometry(-1, "CSC chamber -z endcap");
 
       std::cout << "Validating CSC +z geometry" << std::endl;
       validateCSChamberGeometry(+1, "CSC chamber +z endcap");
 
       std::cout << "Validating CSC layer -z geometry" << std::endl;
       validateCSCLayerGeometry(-1, "CSC layer -z endcap");
 
       std::cout << "Validating CSC layer +z geometry" << std::endl;
       validateCSCLayerGeometry(+1, "CSC layer +z endcap");
     } else
       fwLog(fwlog::kWarning) << "Invalid CSC geometry" << std::endl;
   }
 
   if (doTracker_) {
     trackerGeometry_ = eventSetup.getHandle(trackerGeometryToken_);
 
     if (trackerGeometry_.isValid()) {
       std::cout << "Validating TIB geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsTIB(), "TIB");
       validateStripTopology(trackerGeometry_->detsTIB(), "TIB");
 
       std::cout << "Validating TOB geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsTOB(), "TOB");
       validateStripTopology(trackerGeometry_->detsTOB(), "TOB");
 
       std::cout << "Validating TEC geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsTEC(), "TEC");
       validateStripTopology(trackerGeometry_->detsTEC(), "TEC");
 
       std::cout << "Validating TID geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsTID(), "TID");
       validateStripTopology(trackerGeometry_->detsTID(), "TID");
 
       std::cout << "Validating PXB geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsPXB(), "PXB");
       validatePixelTopology(trackerGeometry_->detsPXB(), "PXB");
 
       std::cout << "Validating PXF geometry and topology" << std::endl;
       validateTrackerGeometry(trackerGeometry_->detsPXF(), "PXF");
       validatePixelTopology(trackerGeometry_->detsPXF(), "PXF");
     } else
       fwLog(fwlog::kWarning) << "Invalid Tracker geometry" << std::endl;
   }
 
   if (doCalo_) {
     caloGeometry_ = eventSetup.getHandle(caloGeometryToken_);
 
     if (caloGeometry_.isValid()) {
       std::cout << "Validating EB geometry" << std::endl;
       validateCaloGeometry(DetId::Ecal, EcalBarrel, "EB");
 
       std::cout << "Validating EE geometry" << std::endl;
       validateCaloGeometry(DetId::Ecal, EcalEndcap, "EE");
 
       std::cout << "Validating ES geometry" << std::endl;
       validateCaloGeometry(DetId::Ecal, EcalPreshower, "ES");
 
       std::cout << "Validating HB geometry" << std::endl;
       validateCaloGeometry(DetId::Hcal, HcalBarrel, "HB");
 
       std::cout << "Validating HE geometry" << std::endl;
       validateCaloGeometry(DetId::Hcal, HcalEndcap, "HE");
 
       std::cout << "Validating HO geometry" << std::endl;
       validateCaloGeometry(DetId::Hcal, HcalOuter, "HO");
 
       std::cout << "Validating HF geometry" << std::endl;
       validateCaloGeometry(DetId::Hcal, HcalForward, "HF");
 
       std::cout << "Validating Castor geometry" << std::endl;
       validateCaloGeometry(DetId::Calo, HcalCastorDetId::SubdetectorId, "Castor");
 
       std::cout << "Validating ZDC geometry" << std::endl;
       validateCaloGeometry(DetId::Calo, HcalZDCDetId::SubdetectorId, "ZDC");
     } else
       fwLog(fwlog::kWarning) << "Invalid Calo geometry" << std::endl;
   }
 }

◆ beginJob()

void ValidateGeometry::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1050 of file ValidateGeometry.cc.

References outFile_.

1050 { outFile_->cd(); }

ValidateGeometry::outFile_

TFile * outFile_

Definition: ValidateGeometry.cc:115

◆ clearData()

void ValidateGeometry::clearData ( )

inlineprivate

Definition at line 123 of file ValidateGeometry.cc.

References bottomWidths_, globalDistances_, lengths_, thicknesses_, and topWidths_.

Referenced by validateCaloGeometry(), validateCSChamberGeometry(), validateCSCLayerGeometry(), validateDTChamberGeometry(), validateDTLayerGeometry(), validateRPCGeometry(), and validateTrackerGeometry().

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

◆ compareShape()

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

private

Definition at line 923 of file ValidateGeometry.cc.

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

Referenced by validateCSChamberGeometry(), validateCSCLayerGeometry(), validateDTChamberGeometry(), validateDTLayerGeometry(), validateRPCGeometry(), and validateTrackerGeometry().

                                                                           {
   double shape_topWidth;
   double shape_bottomWidth;
   double shape_length;
   double shape_thickness;
 
   if (shape[0] == 1) {
     shape_topWidth = shape[2];
     shape_bottomWidth = shape[1];
     shape_length = shape[4];
     shape_thickness = shape[3];
   }
 
   else if (shape[0] == 2) {
     shape_topWidth = shape[1];
     shape_bottomWidth = shape[1];
     shape_length = shape[2];
     shape_thickness = shape[3];
   }
 
   else {
     std::cout << "Failed to get box or trapezoid from shape" << std::endl;
     return;
   }
 
   double topWidth, bottomWidth;
   double length, thickness;
 
   const Bounds* bounds = &(det->surface().bounds());
 
   if (const TrapezoidalPlaneBounds* tpbs = dynamic_cast<const TrapezoidalPlaneBounds*>(bounds)) {
     std::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 {
     std::cout << "Failed to get bounds" << std::endl;
     return;
   }
 
   //assert((tgeotrap && trapezoid) || (tgeobbox && rectangle));
 
   /*
   std::cout<<"topWidth: "<< shape_topWidth <<" "<< topWidth <<std::endl;
   std::cout<<"bottomWidth: "<< shape_bottomWidth <<" "<< bottomWidth <<std::endl;
   std::cout<<"length: "<< shape_length <<" "<< length <<std::endl;
   std::cout<<"thickness: "<< shape_thickness <<" "<< thickness <<std::endl;
   */
 
   topWidths_.push_back(fabs(shape_topWidth - topWidth));
   bottomWidths_.push_back(fabs(shape_bottomWidth - bottomWidth));
   lengths_.push_back(fabs(shape_length - length));
   thicknesses_.push_back(fabs(shape_thickness - thickness));
 
   return;
 }

◆ compareTransform()

void ValidateGeometry::compareTransform	(	const GlobalPoint &	point,
		const TGeoMatrix *	matrix
	)

private

Definition at line 912 of file ValidateGeometry.cc.

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

Referenced by validateCSChamberGeometry(), validateCSCLayerGeometry(), validateDTChamberGeometry(), validateDTLayerGeometry(), validateRPCGeometry(), and validateTrackerGeometry().

                                                                                        {
   double local[3] = {0.0, 0.0, 0.0};
 
   double global[3];
 
   matrix->LocalToMaster(local, global);
 
   double distance = getDistance(GlobalPoint(global[0], global[1], global[2]), gp);
   globalDistances_.push_back(distance);
 }

◆ dataEmpty()

bool ValidateGeometry::dataEmpty ( )

inlineprivate

Definition at line 131 of file ValidateGeometry.cc.

References bottomWidths_, globalDistances_, lengths_, thicknesses_, and topWidths_.

                    {
     return (globalDistances_.empty() && topWidths_.empty() && bottomWidths_.empty() && lengths_.empty() &&
             thicknesses_.empty());
   }

◆ endJob()

void ValidateGeometry::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1052 of file ValidateGeometry.cc.

References gather_cfg::cout, outFile_, and outfileName_.

                               {
   std::cout << "Done. " << std::endl;
   std::cout << "Results written to " << outfileName_ << std::endl;
   outFile_->Close();
 }

◆ getDistance()

double ValidateGeometry::getDistance	(	const GlobalPoint &	point1,
		const GlobalPoint &	point2
	)

private

Definition at line 993 of file ValidateGeometry.cc.

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

Referenced by compareTransform(), and validateCaloGeometry().

                                                                                  {
   /*
   std::cout<<"X: "<< p1.x() <<" "<< p2.x() <<std::endl;
   std::cout<<"Y: "<< p1.y() <<" "<< p2.y() <<std::endl;
   std::cout<<"Z: "<< p1.z() <<" "<< p2.z() <<std::endl;
   */
 
   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 ValidateGeometry::makeHistogram	(	const std::string &	name,
		std::vector< double > &	data
	)

private

Definition at line 1027 of file ValidateGeometry.cc.

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

Referenced by makeHistograms(), validateCSCLayerGeometry(), validateDTLayerGeometry(), validatePixelTopology(), validateRPCGeometry(), and validateStripTopology().

                                                                                    {
   if (data.empty())
     return;
 
   std::vector<double>::iterator it;
 
   it = std::min_element(data.begin(), data.end());
   double minE = *it;
 
   it = std::max_element(data.begin(), data.end());
   double maxE = *it;
 
   std::vector<double>::iterator itEnd = data.end();
 
   TH1D hist(name.c_str(), name.c_str(), 100, minE * (1 + 0.10), maxE * (1 + 0.10));
 
   for (it = data.begin(); it != itEnd; ++it)
     hist.Fill(*it);
 
   hist.GetXaxis()->SetTitle("[cm]");
   hist.Write();
 }

◆ makeHistograms()

void ValidateGeometry::makeHistograms ( const char * detector )

private

Definition at line 1004 of file ValidateGeometry.cc.

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

Referenced by validateCaloGeometry(), validateCSChamberGeometry(), validateCSCLayerGeometry(), validateDTChamberGeometry(), validateDTLayerGeometry(), validateRPCGeometry(), and validateTrackerGeometry().

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

◆ validateCaloGeometry()

void ValidateGeometry::validateCaloGeometry	(	DetId::Detector	detector,
		int	subdetector,
		const char *	detname
	)

private

Definition at line 722 of file ValidateGeometry.cc.

References cms::cuda::assert(), caloGeometry_, clearData(), gather_cfg::cout, hgcalTestNeighbor_cfi::detector, HLT_2024v14_cff::distance, fwGeometry_, FWGeometry::getCorners(), getDistance(), CaloGeometry::getSubdetectorGeometry(), globalDistances_, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistograms(), hltrates_dqm_sourceclient-live_cfg::offset, nano_mu_digi_cff::rawId, EZArrayFL< T >::size(), subdetector, x, y, and z.

Referenced by analyze().

                                                                                                         {
   clearData();
 
   const CaloSubdetectorGeometry* geometry = caloGeometry_->getSubdetectorGeometry(detector, subdetector);
 
   const std::vector<DetId>& ids = geometry->getValidDetIds(detector, subdetector);
 
   for (auto it = ids.begin(), iEnd = ids.end(); it != iEnd; ++it) {
     unsigned int rawId = (*it).rawId();
 
     const float* points = fwGeometry_.getCorners(rawId);
 
     if (points == nullptr) {
       std::cout << "Failed to get points of " << detname << " element with detid: " << rawId << std::endl;
       continue;
     }
 
     auto cellGeometry = geometry->getGeometry(*it);
     const CaloCellGeometry::CornersVec& corners = cellGeometry->getCorners();
 
     assert(corners.size() == 8);
 
     for (unsigned int i = 0, offset = 0; i < 8; ++i) {
       offset = 2 * i;
 
       double distance = getDistance(GlobalPoint(points[i + offset], points[i + 1 + offset], points[i + 2 + offset]),
                                     GlobalPoint(corners[i].x(), corners[i].y(), corners[i].z()));
 
       globalDistances_.push_back(distance);
     }
   }
 
   makeHistograms(detname);
 }

◆ validateCSChamberGeometry()

void ValidateGeometry::validateCSChamberGeometry	(	const int	endcap,
		const char *	detname
	)

private

Definition at line 456 of file ValidateGeometry.cc.

References relativeConstraints::chamber, CSCGeometry::chambers(), chambers, clearData(), compareShape(), compareTransform(), gather_cfg::cout, cscGeometry_, hcalRecHitTable_cff::detId, makeMuonMisalignmentScenario::endcap, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getShapePars(), runTauDisplay::gp, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistograms(), makeMuonMisalignmentScenario::matrix, and l1trig_cff::shape.

Referenced by analyze().

                                                                                       {
   clearData();
 
   auto const& chambers = cscGeometry_->chambers();
 
   for (auto it = chambers.begin(), itEnd = chambers.end(); it != itEnd; ++it) {
     const CSCChamber* chamber = *it;
 
     if (chamber && chamber->id().endcap() == endcap) {
       DetId detId = chamber->geographicalId();
       GlobalPoint gp = chamber->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
 
       const TGeoMatrix* matrix = fwGeometry_.getMatrix(detId.rawId());
 
       if (!matrix) {
         std::cout << "Failed to get matrix of CSC chamber with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       compareTransform(gp, matrix);
 
       const float* shape = fwGeometry_.getShapePars(detId.rawId());
 
       if (!shape) {
         std::cout << "Failed to get shape of CSC chamber with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       compareShape(chamber, shape);
     }
   }
 
   makeHistograms(detname);
 }

◆ validateCSCLayerGeometry()

void ValidateGeometry::validateCSCLayerGeometry	(	const int	endcap,
		const char *	detname
	)

private

Definition at line 491 of file ValidateGeometry.cc.

References cms::cuda::assert(), clearData(), compareShape(), compareTransform(), gather_cfg::cout, cscGeometry_, hcalRecHitTable_cff::detId, makeMuonMisalignmentScenario::endcap, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getParameters(), FWGeometry::getShapePars(), runTauDisplay::gp, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, CSCGeometry::layers(), makeHistogram(), makeHistograms(), makeMuonMisalignmentScenario::matrix, relativeConstraints::ring, l1trig_cff::shape, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

                                                                                      {
   clearData();
   std::vector<double> strip_positions;
   std::vector<double> wire_positions;
 
   std::vector<double> me11_wiresLocal;
   std::vector<double> me12_wiresLocal;
   std::vector<double> me13_wiresLocal;
   std::vector<double> me14_wiresLocal;
   std::vector<double> me21_wiresLocal;
   std::vector<double> me22_wiresLocal;
   std::vector<double> me31_wiresLocal;
   std::vector<double> me32_wiresLocal;
   std::vector<double> me41_wiresLocal;
   std::vector<double> me42_wiresLocal;
 
   auto const& layers = cscGeometry_->layers();
 
   for (auto it = layers.begin(), itEnd = layers.end(); it != itEnd; ++it) {
     const CSCLayer* layer = *it;
 
     if (layer && layer->id().endcap() == endcap) {
       DetId detId = layer->geographicalId();
       GlobalPoint gp = layer->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
 
       const TGeoMatrix* matrix = fwGeometry_.getMatrix(detId.rawId());
 
       if (!matrix) {
         std::cout << "Failed to get matrix of CSC layer with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       compareTransform(gp, matrix);
 
       const float* shape = fwGeometry_.getShapePars(detId.rawId());
 
       if (!shape) {
         std::cout << "Failed to get shape of CSC layer with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       compareShape(layer, shape);
 
       double length;
 
       if (shape[0] == 1) {
         length = shape[4];
       }
 
       else {
         std::cout << "Failed to get trapezoid from shape for CSC layer with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       const float* parameters = fwGeometry_.getParameters(detId.rawId());
 
       if (parameters == nullptr) {
         std::cout << "Parameters empty for CSC layer with detid: " << detId.rawId() << std::endl;
         continue;
       }
 
       int yAxisOrientation = layer->geometry()->topology()->yAxisOrientation();
       assert(yAxisOrientation == parameters[0]);
 
       float centreToIntersection = layer->geometry()->topology()->centreToIntersection();
       assert(centreToIntersection == parameters[1]);
 
       float yCentre = layer->geometry()->topology()->yCentreOfStripPlane();
       assert(yCentre == parameters[2]);
 
       float phiOfOneEdge = layer->geometry()->topology()->phiOfOneEdge();
       assert(phiOfOneEdge == parameters[3]);
 
       float stripOffset = layer->geometry()->topology()->stripOffset();
       assert(stripOffset == parameters[4]);
 
       float angularWidth = layer->geometry()->topology()->angularWidth();
       assert(angularWidth == parameters[5]);
 
       for (int nStrip = 1; nStrip <= layer->geometry()->numberOfStrips(); ++nStrip) {
         float xOfStrip1 = layer->geometry()->xOfStrip(nStrip);
 
         double stripAngle = phiOfOneEdge + yAxisOrientation * (nStrip - (0.5 - stripOffset)) * angularWidth;
         double xOfStrip2 = yAxisOrientation * (centreToIntersection - yCentre) * tan(stripAngle);
 
         strip_positions.push_back(xOfStrip1 - xOfStrip2);
       }
 
       int station = layer->id().station();
       int ring = layer->id().ring();
 
       double wireSpacingInGroup = layer->geometry()->wireTopology()->wireSpacing();
       assert(wireSpacingInGroup == parameters[6]);
 
       double wireSpacing = 0.0;
       // we calculate an average wire group
       // spacing from radialExtentOfTheWirePlane / numOfWireGroups
 
       double extentOfWirePlane = 0.0;
 
       if (ring == 2) {
         if (station == 1)
           extentOfWirePlane = 174.81;  //wireSpacing = 174.81/64;
         else
           extentOfWirePlane = 323.38;  //wireSpacing = 323.38/64;
       } else if (station == 1 && (ring == 1 || ring == 4))
         extentOfWirePlane = 150.5;  //wireSpacing = 150.5/48;
       else if (station == 1 && ring == 3)
         extentOfWirePlane = 164.47;  //wireSpacing = 164.47/32;
       else if (station == 2 && ring == 1)
         extentOfWirePlane = 189.97;  //wireSpacing = 189.97/112;
       else if (station == 3 && ring == 1)
         extentOfWirePlane = 170.01;  //wireSpacing = 170.01/96;
       else if (station == 4 && ring == 1)
         extentOfWirePlane = 149.73;  //wireSpacing = 149.73/96;
 
       float wireAngle = layer->geometry()->wireTopology()->wireAngle();
       assert(wireAngle == parameters[7]);
 
       //float cosWireAngle = cos(wireAngle);
 
       /* NOTE
          Some parameters don't seem available in a public interface
          so have to perhaps hard-code. This may not be too bad as there
          seems to be a lot of degeneracy. 
       */
 
       double alignmentPinToFirstWire;
       double yAlignmentFrame = 3.49;
 
       if (station == 1) {
         if (ring == 1 || ring == 4) {
           alignmentPinToFirstWire = 1.065;
           yAlignmentFrame = 0.0;
         }
 
         else  // ME12, ME 13
           alignmentPinToFirstWire = 2.85;
       }
 
       else if (station == 4 && ring == 1)
         alignmentPinToFirstWire = 3.04;
 
       else if (station == 3 && ring == 1)
         alignmentPinToFirstWire = 2.84;
 
       else  // ME21, ME22, ME32, ME42
         alignmentPinToFirstWire = 2.87;
 
       double yOfFirstWire = (yAlignmentFrame - length) + alignmentPinToFirstWire;
 
       int nWireGroups = layer->geometry()->numberOfWireGroups();
       double E = extentOfWirePlane / nWireGroups;
 
       for (int nWireGroup = 1; nWireGroup <= nWireGroups; ++nWireGroup) {
         LocalPoint centerOfWireGroup = layer->geometry()->localCenterOfWireGroup(nWireGroup);
         double yOfWire1 = centerOfWireGroup.y();
 
         //double yOfWire2 = (-0.5 - (nWireGroups*0.5 - 1) + (nWireGroup-1))*E;
         //yOfWire2 += 0.5*E;
         double yOfWire2 = yOfFirstWire + ((nWireGroup - 1) * E);
         yOfWire2 += wireSpacing * 0.5;
 
         double ydiff_local = yOfWire1 - yOfWire2;
         wire_positions.push_back(ydiff_local);
 
         //GlobalPoint globalPoint = layer->surface().toGlobal(LocalPoint(0.0,yOfWire1,0.0));
 
         /*
         float fwLocalPoint[3] = 
         {
           0.0, yOfWire2, 0.0
         };
         
         float fwGlobalPoint[3]; 
         fwGeometry_.localToGlobal(detId.rawId(), fwLocalPoint, fwGlobalPoint);
         double ydiff_global = globalPoint.y() - fwGlobalPoint[1]; 
         */
 
         if (station == 1) {
           if (ring == 1) {
             me11_wiresLocal.push_back(ydiff_local);
           } else if (ring == 2) {
             me12_wiresLocal.push_back(ydiff_local);
           } else if (ring == 3) {
             me13_wiresLocal.push_back(ydiff_local);
           } else if (ring == 4) {
             me14_wiresLocal.push_back(ydiff_local);
           }
         } else if (station == 2) {
           if (ring == 1) {
             me21_wiresLocal.push_back(ydiff_local);
           } else if (ring == 2) {
             me22_wiresLocal.push_back(ydiff_local);
           }
         } else if (station == 3) {
           if (ring == 1) {
             me31_wiresLocal.push_back(ydiff_local);
           } else if (ring == 2) {
             me32_wiresLocal.push_back(ydiff_local);
           }
         } else if (station == 4) {
           if (ring == 1) {
             me41_wiresLocal.push_back(ydiff_local);
           } else if (ring == 2) {
             me42_wiresLocal.push_back(ydiff_local);
           }
         }
       }
     }
   }
 
   std::string hn(detname);
   makeHistogram(hn + ": xOfStrip", strip_positions);
 
   makeHistogram(hn + ": local yOfWire", wire_positions);
 
   makeHistogram("ME11: local yOfWire", me11_wiresLocal);
   makeHistogram("ME12: local yOfWire", me12_wiresLocal);
   makeHistogram("ME13: local yOfWire", me13_wiresLocal);
   makeHistogram("ME14: local yOfWire", me14_wiresLocal);
   makeHistogram("ME21: local yOfWire", me21_wiresLocal);
   makeHistogram("ME22: local yOfWire", me22_wiresLocal);
   makeHistogram("ME31: local yOfWire", me31_wiresLocal);
   makeHistogram("ME32: local yOfWire", me32_wiresLocal);
   makeHistogram("ME41: local yOfWire", me41_wiresLocal);
   makeHistogram("ME42: local yOfWire", me42_wiresLocal);
 
   makeHistograms(detname);
 }

◆ validateDTChamberGeometry()

void ValidateGeometry::validateDTChamberGeometry ( )

private

Definition at line 350 of file ValidateGeometry.cc.

References relativeConstraints::chamber, DTGeometry::chambers(), chambers, clearData(), compareShape(), compareTransform(), gather_cfg::cout, dtGeometry_, 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();
 
   auto const& chambers = dtGeometry_->chambers();
 
   for (auto it = chambers.begin(), itEnd = chambers.end(); it != itEnd; ++it) {
     const DTChamber* chamber = *it;
 
     if (chamber) {
       DTChamberId chId = chamber->id();
       GlobalPoint gp = chamber->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
 
       const TGeoMatrix* matrix = fwGeometry_.getMatrix(chId.rawId());
 
       if (!matrix) {
         std::cout << "Failed to get matrix of DT chamber with detid: " << chId.rawId() << std::endl;
         continue;
       }
 
       compareTransform(gp, matrix);
 
       const float* shape = fwGeometry_.getShapePars(chId.rawId());
 
       if (!shape) {
         std::cout << "Failed to get shape of DT chamber with detid: " << chId.rawId() << std::endl;
         continue;
       }
 
       compareShape(chamber, shape);
     }
   }
 
   makeHistograms("DT chamber");
 }

◆ validateDTLayerGeometry()

void ValidateGeometry::validateDTLayerGeometry ( )

private

Definition at line 385 of file ValidateGeometry.cc.

References cms::cuda::assert(), clearData(), compareShape(), compareTransform(), gather_cfg::cout, dtGeometry_, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getParameters(), FWGeometry::getShapePars(), runTauDisplay::gp, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, DTGeometry::layers(), makeHistogram(), makeHistograms(), makeMuonMisalignmentScenario::matrix, ecaldqm::nChannels, DetId::rawId(), l1trig_cff::shape, Calorimetry_cff::thickness, and ApeEstimator_cff::width.

Referenced by analyze().

                                                {
   clearData();
 
   std::vector<double> wire_positions;
 
   auto const& layers = dtGeometry_->layers();
 
   for (auto it = layers.begin(), itEnd = layers.end(); it != itEnd; ++it) {
     const DTLayer* layer = *it;
 
     if (layer) {
       DTLayerId layerId = layer->id();
       GlobalPoint gp = layer->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
 
       const TGeoMatrix* matrix = fwGeometry_.getMatrix(layerId.rawId());
 
       if (!matrix) {
         std::cout << "Failed to get matrix of DT layer with detid: " << layerId.rawId() << std::endl;
         continue;
       }
 
       compareTransform(gp, matrix);
 
       const float* shape = fwGeometry_.getShapePars(layerId.rawId());
 
       if (!shape) {
         std::cout << "Failed to get shape of DT layer with detid: " << layerId.rawId() << std::endl;
         continue;
       }
 
       compareShape(layer, shape);
 
       const float* parameters = fwGeometry_.getParameters(layerId.rawId());
 
       if (parameters == nullptr) {
         std::cout << "Parameters empty for DT layer with detid: " << layerId.rawId() << std::endl;
         continue;
       }
 
       float width = layer->surface().bounds().width();
       assert(width == parameters[6]);
 
       float thickness = layer->surface().bounds().thickness();
       assert(thickness == parameters[7]);
 
       float length = layer->surface().bounds().length();
       assert(length == parameters[8]);
 
       int firstChannel = layer->specificTopology().firstChannel();
       assert(firstChannel == parameters[3]);
 
       int lastChannel = layer->specificTopology().lastChannel();
       int nChannels = parameters[5];
       assert(nChannels == (lastChannel - firstChannel) + 1);
 
       for (int wireN = firstChannel; wireN - lastChannel <= 0; ++wireN) {
         double localX1 = layer->specificTopology().wirePosition(wireN);
         double localX2 = (wireN - (firstChannel - 1) - 0.5) * parameters[0] - nChannels / 2.0 * parameters[0];
 
         wire_positions.push_back(localX1 - localX2);
 
         //std::cout<<"wireN, localXpos: "<< wireN <<" "<< localX1 <<" "<< localX2 <<std::endl;
       }
     }
   }
 
   makeHistogram("DT layer wire localX", wire_positions);
 
   makeHistograms("DT layer");
 }

◆ validatePixelTopology()

void ValidateGeometry::validatePixelTopology	(	const TrackerGeometry::DetContainer &	dets,
		const char *	detname
	)

private

Definition at line 787 of file ValidateGeometry.cc.

References gather_cfg::cout, fwGeometry_, FWGeometry::getParameters(), TrackerGeometry::idToDetUnit(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistogram(), fireworks::pixelLocalX(), fireworks::pixelLocalY(), nano_mu_digi_cff::rawId, AlCaHLTBitMon_QueryRunRegistry::string, trackerGeometry_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

                                                                                                          {
   std::vector<double> pixelLocalXs;
   std::vector<double> pixelLocalYs;
 
   for (TrackerGeometry::DetContainer::const_iterator it = dets.begin(), itEnd = dets.end(); it != itEnd; ++it) {
     unsigned int rawId = (*it)->geographicalId().rawId();
 
     const float* parameters = fwGeometry_.getParameters(rawId);
 
     if (parameters == nullptr) {
       std::cout << "Parameters empty for " << detname << " element with detid: " << rawId << std::endl;
       continue;
     }
 
     if (const PixelGeomDetUnit* det =
             dynamic_cast<const PixelGeomDetUnit*>(trackerGeometry_->idToDetUnit((*it)->geographicalId()))) {
       if (const PixelTopology* rpt = &det->specificTopology()) {
         int nrows = rpt->nrows();
         int ncolumns = rpt->ncolumns();
 
         for (int row = 1; row <= nrows; ++row) {
           for (int column = 1; column <= ncolumns; ++column) {
             LocalPoint localPoint = rpt->localPosition(MeasurementPoint(row, column));
 
             pixelLocalXs.push_back(localPoint.x() - fireworks::pixelLocalX(row, parameters));
             pixelLocalYs.push_back(localPoint.y() - fireworks::pixelLocalY(column, parameters));
           }
         }
       }
 
       else
         std::cout << "No topology for " << detname << " " << rawId << std::endl;
     }
 
     else
       std::cout << "No geomDetUnit for " << detname << " " << rawId << std::endl;
   }
 
   std::string hn(detname);
   makeHistogram(hn + " pixelLocalX", pixelLocalXs);
   makeHistogram(hn + " pixelLocalY", pixelLocalYs);
 }

◆ validateRPCGeometry()

void ValidateGeometry::validateRPCGeometry	(	const int	regionNumber,
		const char *	regionName
	)

private

Definition at line 278 of file ValidateGeometry.cc.

References cms::cuda::assert(), clearData(), compareShape(), compareTransform(), gather_cfg::cout, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getParameters(), FWGeometry::getShapePars(), runTauDisplay::gp, RPCGeometry::idToDetUnit(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistogram(), makeHistograms(), makeMuonMisalignmentScenario::matrix, me0TriggerPseudoDigis_cff::nStrips, hltrates_dqm_sourceclient-live_cfg::offset, DetId::rawId(), RPCDetId::region(), nano_mu_digi_cff::roll, RPCGeometry::rolls(), rpcGeometry_, l1trig_cff::shape, AlCaHLTBitMon_QueryRunRegistry::string, nano_mu_digi_cff::strip, GeomDet::surface(), Surface::toGlobal(), and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                                                          {
   clearData();
 
   std::vector<double> centers;
 
   auto const& rolls = rpcGeometry_->rolls();
 
   for (auto it = rolls.begin(), itEnd = rolls.end(); it != itEnd; ++it) {
     const RPCRoll* roll = *it;
 
     if (roll) {
       RPCDetId rpcDetId = roll->id();
 
       if (rpcDetId.region() == regionNumber) {
         const GeomDetUnit* det = rpcGeometry_->idToDetUnit(rpcDetId);
         GlobalPoint gp = det->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
 
         const TGeoMatrix* matrix = fwGeometry_.getMatrix(rpcDetId.rawId());
 
         if (!matrix) {
           std::cout << "Failed to get matrix of RPC with detid: " << rpcDetId.rawId() << std::endl;
           continue;
         }
 
         compareTransform(gp, matrix);
 
         const float* shape = fwGeometry_.getShapePars(rpcDetId.rawId());
 
         if (!shape) {
           std::cout << "Failed to get shape of RPC with detid: " << rpcDetId.rawId() << std::endl;
           continue;
         }
 
         compareShape(det, shape);
 
         const float* parameters = fwGeometry_.getParameters(rpcDetId.rawId());
 
         if (parameters == nullptr) {
           std::cout << "Parameters empty for RPC with detid: " << rpcDetId.rawId() << std::endl;
           continue;
         }
 
         // Yes, I know that below I'm comparing the equivalence
         // of floating point numbers
 
         int nStrips = roll->nstrips();
         assert(nStrips == parameters[0]);
 
         float stripLength = roll->specificTopology().stripLength();
         assert(stripLength == parameters[1]);
 
         float pitch = roll->specificTopology().pitch();
         assert(pitch == parameters[2]);
 
         float offset = -0.5 * nStrips * pitch;
 
         for (int strip = 1; strip <= roll->nstrips(); ++strip) {
           LocalPoint centreOfStrip1 = roll->centreOfStrip(strip);
           LocalPoint centreOfStrip2 = LocalPoint((strip - 0.5) * pitch + offset, 0.0);
 
           centers.push_back(centreOfStrip1.x() - centreOfStrip2.x());
         }
       }
     }
   }
 
   std::string hn(regionName);
   makeHistogram(hn + ": centreOfStrip", centers);
 
   makeHistograms(regionName);
 }

◆ validateStripTopology()

void ValidateGeometry::validateStripTopology	(	const TrackerGeometry::DetContainer &	dets,
		const char *	detname
	)

private

Definition at line 830 of file ValidateGeometry.cc.

References cms::cuda::assert(), gather_cfg::cout, fwGeometry_, FWGeometry::getParameters(), TrackerGeometry::idToDet(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, StripTopology::localPosition(), makeHistogram(), StripTopology::nstrips(), StripTopology::pitch(), nano_mu_digi_cff::rawId, AlCaHLTBitMon_QueryRunRegistry::string, nano_mu_digi_cff::strip, StripTopology::stripLength(), funct::tan(), trackerGeometry_, and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                                                                          {
   std::vector<double> radialStripLocalXs;
   std::vector<double> rectangularStripLocalXs;
 
   for (TrackerGeometry::DetContainer::const_iterator it = dets.begin(), itEnd = dets.end(); it != itEnd; ++it) {
     unsigned int rawId = (*it)->geographicalId().rawId();
 
     const float* parameters = fwGeometry_.getParameters(rawId);
 
     if (parameters == nullptr) {
       std::cout << "Parameters empty for " << detname << " element with detid: " << rawId << std::endl;
       continue;
     }
 
     if (const StripGeomDetUnit* det =
             dynamic_cast<const StripGeomDetUnit*>(trackerGeometry_->idToDet((*it)->geographicalId()))) {
       // NOTE: why the difference in dets vs. units between these and pixels? The dynamic cast above
       // fails for many of the detids...
 
       const StripTopology* st = dynamic_cast<const StripTopology*>(&det->specificTopology());
 
       if (st) {
         //assert(parameters[0] == 0);
         int nstrips = st->nstrips();
         assert(parameters[1] == nstrips);
         assert(parameters[2] == st->stripLength());
 
         if (const RadialStripTopology* rst =
                 dynamic_cast<const RadialStripTopology*>(&(det->specificType().specificTopology()))) {
           assert(parameters[0] == 1);
           assert(parameters[3] == rst->yAxisOrientation());
           assert(parameters[4] == rst->originToIntersection());
           assert(parameters[5] == rst->phiOfOneEdge());
           assert(parameters[6] == rst->angularWidth());
 
           for (uint16_t strip = 1; strip <= nstrips; ++strip) {
             float stripAngle1 = rst->stripAngle(strip);
             float stripAngle2 = parameters[3] * (parameters[5] + strip * parameters[6]);
 
             assert((stripAngle1 - stripAngle2) == 0);
 
             LocalPoint stripPosition = st->localPosition(strip);
 
             float stripX = parameters[4] * tan(stripAngle2);
             radialStripLocalXs.push_back(stripPosition.x() - stripX);
           }
         }
 
         else if (dynamic_cast<const RectangularStripTopology*>(&(det->specificType().specificTopology()))) {
           assert(parameters[0] == 2);
           assert(parameters[3] == st->pitch());
 
           for (uint16_t strip = 1; strip <= nstrips; ++strip) {
             LocalPoint stripPosition = st->localPosition(strip);
             float stripX = -parameters[1] * 0.5 * parameters[3];
             stripX += strip * parameters[3];
             rectangularStripLocalXs.push_back(stripPosition.x() - stripX);
           }
         }
 
         else if (dynamic_cast<const TrapezoidalStripTopology*>(&(det->specificType().specificTopology()))) {
           assert(parameters[0] == 3);
           assert(parameters[3] == st->pitch());
         }
 
         else
           std::cout << "Failed to get pitch for " << detname << " " << rawId << std::endl;
       }
 
       else
         std::cout << "Failed cast to StripTopology for " << detname << " " << rawId << std::endl;
     }
 
     //else
     //  std::cout<<"Failed cast to StripGeomDetUnit for "<< detname <<" "<< rawId <<std::endl;
   }
 
   std::string hn(detname);
   makeHistogram(hn + " radial strip localX", radialStripLocalXs);
   makeHistogram(hn + " rectangular strip localX", rectangularStripLocalXs);
 }

◆ validateTrackerGeometry()

void ValidateGeometry::validateTrackerGeometry	(	const TrackerGeometry::DetContainer &	dets,
		const char *	detname
	)

private

Definition at line 757 of file ValidateGeometry.cc.

References clearData(), compareShape(), compareTransform(), gather_cfg::cout, fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getShapePars(), runTauDisplay::gp, TrackerGeometry::idToDet(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistograms(), makeMuonMisalignmentScenario::matrix, nano_mu_digi_cff::rawId, l1trig_cff::shape, and trackerGeometry_.

Referenced by analyze().

                                                                                                            {
   clearData();
 
   for (TrackerGeometry::DetContainer::const_iterator it = dets.begin(), itEnd = dets.end(); it != itEnd; ++it) {
     GlobalPoint gp =
         (trackerGeometry_->idToDet((*it)->geographicalId()))->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
     unsigned int rawId = (*it)->geographicalId().rawId();
 
     const TGeoMatrix* matrix = fwGeometry_.getMatrix(rawId);
 
     if (!matrix) {
       std::cout << "Failed to get matrix of " << detname << " element with detid: " << rawId << std::endl;
       continue;
     }
 
     compareTransform(gp, matrix);
 
     const float* shape = fwGeometry_.getShapePars(rawId);
 
     if (!shape) {
       std::cout << "Failed to get shape of " << detname << " element with detid: " << rawId << std::endl;
       continue;
     }
 
     compareShape(*it, shape);
   }
 
   makeHistograms(detname);
 }