Inheritance diagram for DTGeometryValidate:

Public Member Functions
	DTGeometryValidate (const ParameterSet &)

	~DTGeometryValidate () override

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

void	beginJob () override

void	clearData ()

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

void	validateDTLayerGeometry ()

Private Attributes
vector< float >	bottomWidths_

edm::ESHandle< DTGeometry >	dtGeometry_

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeometryToken_

FWGeometry	fwGeometry_

vector< float >	globalDistances_

string	infileName_

vector< float >	lengths_

TFile *	outFile_

string	outfileName_

vector< float >	thicknesses_

int	tolerance_

vector< float >	topWidths_

Detailed Description

Definition at line 41 of file DTGeometryValidate.cc.

Constructor & Destructor Documentation

DTGeometryValidate::DTGeometryValidate ( const ParameterSet & iConfig )

explicit

Definition at line 87 of file DTGeometryValidate.cc.

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

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

DTGeometryValidate::~DTGeometryValidate ( )

inlineoverride

Definition at line 45 of file DTGeometryValidate.cc.

References analyze(), and bk::beginJob().

45 {}

Member Function Documentation

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

overrideprivate

Definition at line 98 of file DTGeometryValidate.cc.

References dtGeometry_, dtGeometryToken_, edm::EventSetup::getHandle(), edm::ESHandleBase::isValid(), validateDTChamberGeometry(), and validateDTLayerGeometry().

 {
   dtGeometry_ = eventSetup.getHandle(dtGeometryToken_);
 
   if(dtGeometry_.isValid()) {
     LogVerbatim("DTGeometry") << "Validating DT chamber geometry";
     validateDTChamberGeometry();
 
     LogVerbatim("DTGeometry") <<"Validating DT layer geometry";
     validateDTLayerGeometry();
   }
   else
     LogVerbatim("DTGeometry") << "Invalid DT geometry"; 
 }

void DTGeometryValidate::beginJob ( void )

overrideprivate

Definition at line 339 of file DTGeometryValidate.cc.

References outFile_.

                              {
   outFile_->cd();
 }

void DTGeometryValidate::clearData ( )

inlineprivate

Definition at line 64 of file DTGeometryValidate.cc.

Referenced by validateDTChamberGeometry(), and validateDTLayerGeometry().

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

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

private

Definition at line 226 of file DTGeometryValidate.cc.

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

Referenced by validateDTChamberGeometry(), and validateDTLayerGeometry().

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

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

private

Definition at line 213 of file DTGeometryValidate.cc.

References SoftLeptonByDistance_cfi::distance, getDistance(), globalDistances_, and DTRecHitClients_cfi::local.

Referenced by validateDTChamberGeometry(), and validateDTLayerGeometry().

 {
   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);
   globalDistances_.push_back(distance);
 }

void DTGeometryValidate::endJob ( void )

overrideprivate

Definition at line 344 of file DTGeometryValidate.cc.

References DEFINE_FWK_MODULE, outFile_, and outfileName_.

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

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

private

Definition at line 291 of file DTGeometryValidate.cc.

References almost_equal(), and tolerance_.

Referenced by validateDTLayerGeometry().

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

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

private

Definition at line 283 of file DTGeometryValidate.cc.

References mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

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

private

Definition at line 322 of file DTGeometryValidate.cc.

References begin, end, create_public_lumi_plots::hist, and particleFlow_cfi::maxE.

Referenced by makeHistograms(), and validateDTLayerGeometry().

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

void DTGeometryValidate::makeHistograms ( const char * detector )

private

Definition at line 299 of file DTGeometryValidate.cc.

References bottomWidths_, edmIntegrityCheck::d, globalDistances_, lengths_, makeHistogram(), outFile_, thicknesses_, and topWidths_.

Referenced by validateDTChamberGeometry(), and validateDTLayerGeometry().

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

void DTGeometryValidate::validateDTChamberGeometry ( )

private

Definition at line 114 of file DTGeometryValidate.cc.

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

Referenced by analyze().

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

void DTGeometryValidate::validateDTLayerGeometry ( )

private

Definition at line 147 of file DTGeometryValidate.cc.

References clearData(), compareShape(), compareTransform(), dtGeometry_, f, fwGeometry_, getDiff(), FWGeometry::getMatrix(), FWGeometry::getParameters(), FWGeometry::getShapePars(), runTauDisplay::gp, DTGeometry::layers(), makeHistogram(), makeHistograms(), makeMuonMisalignmentScenario::matrix, ecaldqm::nChannels, DetId::rawId(), and ApeEstimator_cff::width.

Referenced by analyze().

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