DTGeometryValidate Class Reference

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 39 of file DTGeometryValidate.cc.

Constructor & Destructor Documentation

DTGeometryValidate::DTGeometryValidate ( const ParameterSet &  iConfig )

explicit

Definition at line 83 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 42 of file DTGeometryValidate.cc.

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

{}

Member Function Documentation

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

overrideprivate

Definition at line 92 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 303 of file DTGeometryValidate.cc.

References outFile_.

{ outFile_->cd(); }

void DTGeometryValidate::clearData ( )

inlineprivate

Definition at line 61 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 205 of file DTGeometryValidate.cc.

References bottomWidths_, Surface::bounds(), Bounds::length(), lengths_, GeomDet::surface(), Bounds::thickness(), Calorimetry_cff::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 195 of file DTGeometryValidate.cc.

References HLT_2018_cff::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 305 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 260 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 255 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 288 of file DTGeometryValidate.cc.

References begin, end, compare::hist, and HLT_2018_cff::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 267 of file DTGeometryValidate.cc.

References bottomWidths_, ztail::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 105 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 134 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(), Calorimetry_cff::thickness, 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");
}

Member Data Documentation

vector<float> DTGeometryValidate::bottomWidths_

private

Definition at line 75 of file DTGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

edm::ESHandle<DTGeometry> DTGeometryValidate::dtGeometry_

private

Definition at line 70 of file DTGeometryValidate.cc.

Referenced by analyze(), validateDTChamberGeometry(), and validateDTLayerGeometry().

edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTGeometryValidate::dtGeometryToken_

private

Definition at line 69 of file DTGeometryValidate.cc.

Referenced by analyze().

FWGeometry DTGeometryValidate::fwGeometry_

private

Definition at line 71 of file DTGeometryValidate.cc.

Referenced by DTGeometryValidate(), validateDTChamberGeometry(), and validateDTLayerGeometry().

vector<float> DTGeometryValidate::globalDistances_

private

Definition at line 73 of file DTGeometryValidate.cc.

Referenced by compareTransform(), and makeHistograms().

string DTGeometryValidate::infileName_

private

Definition at line 78 of file DTGeometryValidate.cc.

Referenced by DTGeometryValidate().

vector<float> DTGeometryValidate::lengths_

private

Definition at line 76 of file DTGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

TFile* DTGeometryValidate::outFile_

private

Definition at line 72 of file DTGeometryValidate.cc.

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

string DTGeometryValidate::outfileName_

private

Definition at line 79 of file DTGeometryValidate.cc.

Referenced by DTGeometryValidate(), and endJob().

vector<float> DTGeometryValidate::thicknesses_

private

Definition at line 77 of file DTGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

int DTGeometryValidate::tolerance_

private

Definition at line 80 of file DTGeometryValidate.cc.

Referenced by DTGeometryValidate(), and getDiff().

vector<float> DTGeometryValidate::topWidths_

private

Definition at line 74 of file DTGeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().