ME0GeometryValidate Class Reference

Inheritance diagram for ME0GeometryValidate:

Public Member Functions
	ME0GeometryValidate (const ParameterSet &)
	~ME0GeometryValidate () override

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	clearData ()
void	clearData2 ()
void	compareShape (const GeomDet *, const float *)
void	compareTransform (const GlobalPoint &, const TGeoMatrix *)
void	endJob () override
float	getDiff (const float, const float)
float	getDistance (const GlobalPoint &, const GlobalPoint &)
void	makeHistogram (const string &, vector< float > &)
void	makeHistograms (const char *)
void	makeHistograms2 (const char *)
void	validateME0ChamberGeometry ()
void	validateME0EtaPartitionGeometry ()
void	validateME0EtaPartitionGeometry2 ()

Private Attributes
vector< float >	bottomWidths_
FWGeometry	fwGeometry_
vector< float >	globalDistances_
string	infileName_
vector< float >	lengths_
const ME0Geometry *	me0Geometry_
vector< float >	nstrips_
TFile *	outFile_
string	outfileName_
vector< float >	pitch_
vector< float >	stripslen_
vector< float >	thicknesses_
const edm::ESGetToken< ME0Geometry, MuonGeometryRecord >	tokGeom_
int	tolerance_
vector< float >	topWidths_

Detailed Description

Definition at line 49 of file ME0GeometryValidate.cc.

Constructor & Destructor Documentation

ME0GeometryValidate()

ME0GeometryValidate::ME0GeometryValidate ( const ParameterSet &  iConfig )

explicit

Definition at line 104 of file ME0GeometryValidate.cc.

    : tokGeom_{esConsumes<ME0Geometry, MuonGeometryRecord>(edm::ESInputTag{})},
      infileName_(iConfig.getUntrackedParameter<string>("infileName", "cmsRecoGeom-2026.root")),
      outfileName_(iConfig.getUntrackedParameter<string>("outfileName", "validateME0Geometry.root")),
      tolerance_(iConfig.getUntrackedParameter<int>("tolerance", 6)) {
  fwGeometry_.loadMap(infileName_.c_str());
  outFile_ = TFile::Open(outfileName_.c_str(), "RECREATE");
}

~ME0GeometryValidate()

ME0GeometryValidate::~ME0GeometryValidate ( )

inlineoverride

Definition at line 52 of file ME0GeometryValidate.cc.

{}

Member Function Documentation

analyze()

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

overrideprivate

Definition at line 113 of file ME0GeometryValidate.cc.

References options_cfi::eventSetup, LogTrace, me0Geometry_, tokGeom_, validateME0ChamberGeometry(), validateME0EtaPartitionGeometry(), and validateME0EtaPartitionGeometry2().

                                                                                        {
  me0Geometry_ = &eventSetup.getData(tokGeom_);
  LogTrace("ME0Geometry") << "Validating ME0 chamber geometry";
  validateME0ChamberGeometry();
  validateME0EtaPartitionGeometry2();
  validateME0EtaPartitionGeometry();
}

beginJob()

void ME0GeometryValidate::beginJob ( void  )

overrideprivate

Definition at line 329 of file ME0GeometryValidate.cc.

References outFile_.

{ outFile_->cd(); }

clearData()

void ME0GeometryValidate::clearData ( )

inlineprivate

Definition at line 73 of file ME0GeometryValidate.cc.

Referenced by validateME0ChamberGeometry(), and validateME0EtaPartitionGeometry2().

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

clearData2()

void ME0GeometryValidate::clearData2 ( )

inlineprivate

Definition at line 81 of file ME0GeometryValidate.cc.

Referenced by validateME0EtaPartitionGeometry().

                    {
    nstrips_.clear();
    pitch_.clear();
    stripslen_.clear();
  }

compareShape()

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

private

Definition at line 213 of file ME0GeometryValidate.cc.

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

Referenced by validateME0ChamberGeometry(), and validateME0EtaPartitionGeometry2().

                                                                             {
  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 {
    LogTrace("ME0Geometry") << "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 {
    LogTrace("ME0Geometry") << "Failed to get bounds";

    return;
  }
  topWidths_.push_back(std::abs(shapeTopWidth - topWidth));
  bottomWidths_.push_back(std::abs(shapeBottomWidth - bottomWidth));
  lengths_.push_back(std::abs(shapeLength - length));
  thicknesses_.push_back(std::abs(shapeThickness - thickness));
}

compareTransform()

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

private

Definition at line 201 of file ME0GeometryValidate.cc.

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

Referenced by validateME0ChamberGeometry(), and validateME0EtaPartitionGeometry2().

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

endJob()

void ME0GeometryValidate::endJob ( void  )

overrideprivate

Definition at line 331 of file ME0GeometryValidate.cc.

References LogTrace, outFile_, and outfileName_.

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

getDiff()

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

private

Definition at line 271 of file ME0GeometryValidate.cc.

References almost_equal(), and tolerance_.

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

getDistance()

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

private

Definition at line 266 of file ME0GeometryValidate.cc.

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

Referenced by compareTransform().

                                                                                   {
  return sqrt((p1.x() - p2.x()) * (p1.x() - p2.x()) + (p1.y() - p2.y()) * (p1.y() - p2.y()) +
              (p1.z() - p2.z()) * (p1.z() - p2.z()));
}

makeHistogram()

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

private

Definition at line 314 of file ME0GeometryValidate.cc.

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

Referenced by makeHistograms(), and makeHistograms2().

                                                                               {
  if (data.empty())
    return;

  const auto [minE, maxE] = minmax_element(begin(data), end(data));

  TH1D hist(name.c_str(), name.c_str(), 100, *minE * (1 + 0.10), *maxE * (1 + 0.10));

  for (auto const& it : data)
    hist.Fill(it);

  hist.GetXaxis()->SetTitle("[cm]");
  hist.Write();
}

makeHistograms()

void ME0GeometryValidate::makeHistograms ( const char *  detector )

private

Definition at line 278 of file ME0GeometryValidate.cc.

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

Referenced by validateME0ChamberGeometry(), and validateME0EtaPartitionGeometry2().

                                                             {
  outFile_->cd();

  string d(detector);

  string gdn = d + ": distance between points in global coordinates";
  makeHistogram(gdn, globalDistances_);

  string twn = d + ": absolute difference between top widths (along X)";
  makeHistogram(twn, topWidths_);

  string bwn = d + ": absolute difference between bottom widths (along X)";
  makeHistogram(bwn, bottomWidths_);

  string ln = d + ": absolute difference between lengths (along Y)";
  makeHistogram(ln, lengths_);

  string tn = d + ": absolute difference between thicknesses (along Z)";
  makeHistogram(tn, thicknesses_);
}

makeHistograms2()

void ME0GeometryValidate::makeHistograms2 ( const char *  detector )

private

Definition at line 299 of file ME0GeometryValidate.cc.

References ztail::d, hgcalTestNeighbor_cfi::detector, makeHistogram(), nstrips_, outFile_, pi, pitch_, and stripslen_.

Referenced by validateME0EtaPartitionGeometry().

                                                              {
  outFile_->cd();

  string d(detector);

  string ns = d + ": absolute difference between nStrips";
  makeHistogram(ns, nstrips_);

  string pi = d + ": absolute difference between Strips Pitch";
  makeHistogram(pi, pitch_);

  string pl = d + ": absolute difference between Strips Length";
  makeHistogram(pl, stripslen_);
}

validateME0ChamberGeometry()

void ME0GeometryValidate::validateME0ChamberGeometry ( )

private

Definition at line 121 of file ME0GeometryValidate.cc.

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

Referenced by analyze().

                                                     {
  clearData();

  for (auto const& it : me0Geometry_->chambers()) {
    ME0DetId chId = it->id();
    GlobalPoint gp = it->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
    const TGeoMatrix* matrix = fwGeometry_.getMatrix(chId.rawId());

    if (!matrix) {
      LogTrace("ME0Geometry") << "Failed to get matrix of ME0 chamber with detid: " << chId.rawId();
      continue;
    }
    compareTransform(gp, matrix);

    auto const& shape = fwGeometry_.getShapePars(chId.rawId());

    if (!shape) {
      LogTrace("ME0Geometry") << "Failed to get shape of ME0 chamber with detid: " << chId.rawId();
      continue;
    }
    compareShape(it, shape);
  }
  makeHistograms("ME0 Chamber");
}

validateME0EtaPartitionGeometry()

void ME0GeometryValidate::validateME0EtaPartitionGeometry ( )

private

Definition at line 171 of file ME0GeometryValidate.cc.

References funct::abs(), clearData2(), MillePedeFileConverter_cfg::e, ME0Geometry::etaPartitions(), f, fwGeometry_, FWGeometry::getParameters(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeHistograms2(), me0Geometry_, nstrips_, pitch_, DetId::rawId(), StripTopology::stripLength(), and stripslen_.

Referenced by analyze().

                                                          {
  clearData2();

  for (auto const& it : me0Geometry_->etaPartitions()) {
    ME0DetId chId = it->id();
    const int n_strips = it->nstrips();
    const float n_pitch = it->pitch();
    const StripTopology& topo = it->specificTopology();
    const float stripLen = topo.stripLength();
    const float* parameters = fwGeometry_.getParameters(chId.rawId());
    nstrips_.push_back(std::abs(n_strips - parameters[0]));

    if (n_strips) {
      for (int istrips = 1; istrips <= n_strips; istrips++) {
        if (std::abs(n_pitch - parameters[2]) < 1.0e-05)
          pitch_.push_back(0.0f);
        else
          pitch_.push_back(std::abs(n_pitch - parameters[2]));
        if (std::abs(stripLen - parameters[1]) < 1.0e-05)
          pitch_.push_back(0.0f);
        else
          stripslen_.push_back(std::abs(stripLen - parameters[1]));
      }
    } else {
      LogWarning("ME0Geometry") << "ATTENTION! nStrips == 0";
    }
  }
  makeHistograms2("ME0 Eta Partition");
}

validateME0EtaPartitionGeometry2()

void ME0GeometryValidate::validateME0EtaPartitionGeometry2 ( )

private

Definition at line 146 of file ME0GeometryValidate.cc.

References clearData(), compareShape(), compareTransform(), ME0Geometry::etaPartitions(), fwGeometry_, FWGeometry::getMatrix(), FWGeometry::getShapePars(), runTauDisplay::gp, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogTrace, makeHistograms(), makeMuonMisalignmentScenario::matrix, me0Geometry_, DetId::rawId(), and l1trig_cff::shape.

Referenced by analyze().

                                                           {
  clearData();

  for (auto const& it : me0Geometry_->etaPartitions()) {
    ME0DetId chId = it->id();
    GlobalPoint gp = it->surface().toGlobal(LocalPoint(0.0, 0.0, 0.0));
    const TGeoMatrix* matrix = fwGeometry_.getMatrix(chId.rawId());

    if (!matrix) {
      LogTrace("ME0Geometry") << "Failed to get matrix of ME0 eta partition 2 with detid: " << chId.rawId();
      continue;
    }
    compareTransform(gp, matrix);

    auto const& shape = fwGeometry_.getShapePars(chId.rawId());

    if (!shape) {
      LogTrace("ME0Geometry") << "Failed to get shape of ME0 eta partition 2 with detid: " << chId.rawId();
      continue;
    }
    compareShape(it, shape);
  }
  makeHistograms("ME0 Eta Partition");
}

Member Data Documentation

bottomWidths_

vector<float> ME0GeometryValidate::bottomWidths_

private

Definition at line 93 of file ME0GeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

fwGeometry_

FWGeometry ME0GeometryValidate::fwGeometry_

private

Definition at line 89 of file ME0GeometryValidate.cc.

Referenced by validateME0ChamberGeometry(), validateME0EtaPartitionGeometry(), and validateME0EtaPartitionGeometry2().

globalDistances_

vector<float> ME0GeometryValidate::globalDistances_

private

Definition at line 91 of file ME0GeometryValidate.cc.

Referenced by compareTransform(), and makeHistograms().

infileName_

string ME0GeometryValidate::infileName_

private

Definition at line 99 of file ME0GeometryValidate.cc.

lengths_

vector<float> ME0GeometryValidate::lengths_

private

Definition at line 94 of file ME0GeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

me0Geometry_

const ME0Geometry* ME0GeometryValidate::me0Geometry_

private

Definition at line 88 of file ME0GeometryValidate.cc.

Referenced by analyze(), validateME0ChamberGeometry(), validateME0EtaPartitionGeometry(), and validateME0EtaPartitionGeometry2().

nstrips_

vector<float> ME0GeometryValidate::nstrips_

private

Definition at line 96 of file ME0GeometryValidate.cc.

Referenced by makeHistograms2(), and validateME0EtaPartitionGeometry().

outFile_

TFile* ME0GeometryValidate::outFile_

private

Definition at line 90 of file ME0GeometryValidate.cc.

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

outfileName_

string ME0GeometryValidate::outfileName_

private

Definition at line 100 of file ME0GeometryValidate.cc.

Referenced by endJob().

pitch_

vector<float> ME0GeometryValidate::pitch_

private

Definition at line 97 of file ME0GeometryValidate.cc.

Referenced by makeHistograms2(), and validateME0EtaPartitionGeometry().

stripslen_

vector<float> ME0GeometryValidate::stripslen_

private

Definition at line 98 of file ME0GeometryValidate.cc.

Referenced by makeHistograms2(), and validateME0EtaPartitionGeometry().

thicknesses_

vector<float> ME0GeometryValidate::thicknesses_

private

Definition at line 95 of file ME0GeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().

tokGeom_

const edm::ESGetToken<ME0Geometry, MuonGeometryRecord> ME0GeometryValidate::tokGeom_

private

Definition at line 87 of file ME0GeometryValidate.cc.

Referenced by analyze().

tolerance_

int ME0GeometryValidate::tolerance_

private

Definition at line 101 of file ME0GeometryValidate.cc.

Referenced by getDiff().

topWidths_

vector<float> ME0GeometryValidate::topWidths_

private

Definition at line 92 of file ME0GeometryValidate.cc.

Referenced by compareShape(), and makeHistograms().