CSCGeometryValidate Class Reference

Inheritance diagram for CSCGeometryValidate:

Public Member Functions
	CSCGeometryValidate (const ParameterSet &)
	~CSCGeometryValidate () 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	validateCSCChamberGeometry ()
void	validateCSCLayerGeometry ()

Private Attributes
vector< float >	angularWidth_
vector< float >	bottomWidths_
vector< float >	centreToIntersection_
edm::ESHandle< 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_
int	tolerance_
vector< float >	topWidths_
vector< float >	wireAngle_
vector< float >	wireSpacing_
vector< float >	yAxisOrientation_
vector< float >	yCentreOfStripPlane_

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 53 of file CSCGeometryValidate.cc.

Constructor & Destructor Documentation

CSCGeometryValidate()

CSCGeometryValidate::CSCGeometryValidate ( const ParameterSet &  iConfig )

explicit

Definition at line 120 of file CSCGeometryValidate.cc.

    : 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");
}

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

~CSCGeometryValidate()

CSCGeometryValidate::~CSCGeometryValidate ( )

inlineoverride

Definition at line 56 of file CSCGeometryValidate.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 128 of file CSCGeometryValidate.cc.

                                                                                        {
  eventSetup.get<MuonGeometryRecord>().get(cscGeometry_);
  if (cscGeometry_.isValid()) {
    LogVerbatim("CSCGeometry") << "Validating CSC chamber geometry";
    validateCSCChamberGeometry();
    validateCSCLayerGeometry();
  } else
    LogVerbatim("CSCGeometry") << "Invalid CSC geometry";
}

References cscGeometry_, edm::EventSetup::get(), get, edm::ESHandleBase::isValid(), validateCSCChamberGeometry(), and validateCSCLayerGeometry().

beginJob()

void CSCGeometryValidate::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 350 of file CSCGeometryValidate.cc.

{ outFile_->cd(); }

References outFile_.

clearData()

void CSCGeometryValidate::clearData ( )

inlineprivate

Definition at line 76 of file CSCGeometryValidate.cc.

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

Referenced by validateCSCChamberGeometry().

clearData2()

void CSCGeometryValidate::clearData2 ( )

inlineprivate

Definition at line 84 of file CSCGeometryValidate.cc.

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

Referenced by validateCSCLayerGeometry().

compareShape()

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

private

Definition at line 221 of file CSCGeometryValidate.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("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));
}

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

Referenced by validateCSCChamberGeometry().

compareTransform()

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

private

Definition at line 209 of file CSCGeometryValidate.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 ((distance >= 0.0) && (distance < 1.0e-7))
    distance = 0.0;  // set a tollerance for the distance inside Histos
  globalDistances_.push_back(distance);
}

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

Referenced by validateCSCChamberGeometry().

endJob()

void CSCGeometryValidate::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 352 of file CSCGeometryValidate.cc.

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

References outFile_, and outfileName_.

getDiff()

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

private

Definition at line 277 of file CSCGeometryValidate.cc.

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

References almost_equal(), and tolerance_.

getDistance()

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

private

Definition at line 272 of file CSCGeometryValidate.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 CSCGeometryValidate::makeHistogram ( const string &  name, vector< float > &  data  )

private

Definition at line 335 of file CSCGeometryValidate.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 CSCGeometryValidate::makeHistograms ( const char *  detector )

private

Definition at line 284 of file CSCGeometryValidate.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 validateCSCChamberGeometry().

makeHistograms2()

void CSCGeometryValidate::makeHistograms2 ( const char *  detector )

private

Definition at line 305 of file CSCGeometryValidate.cc.

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

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

Referenced by validateCSCLayerGeometry().

validateCSCChamberGeometry()

void CSCGeometryValidate::validateCSCChamberGeometry ( )

private

Definition at line 138 of file CSCGeometryValidate.cc.

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

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

Referenced by analyze().

validateCSCLayerGeometry()

void CSCGeometryValidate::validateCSCLayerGeometry ( )

private

Definition at line 164 of file CSCGeometryValidate.cc.

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

References CSCRadialStripTopology::angularWidth(), angularWidth_, CSCRadialStripTopology::centreToIntersection(), centreToIntersection_, clearData2(), cscGeometry_, fwGeometry_, FWGeometry::getParameters(), 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().

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_

edm::ESHandle<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 CSCGeometryValidate(), 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.

Referenced by CSCGeometryValidate().

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(), CSCGeometryValidate(), endJob(), makeHistograms(), and makeHistograms2().

outfileName_

string CSCGeometryValidate::outfileName_

private

Definition at line 116 of file CSCGeometryValidate.cc.

Referenced by CSCGeometryValidate(), and 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().

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