HGCalCalibrationCell Class Reference

#include <HGCalCalibrationCell.h>

Public Member Functions
int	findCell (int zside, int layer, int waferU, int waferV, int cellUV, const std::pair< double, double > &xy) const
	HGCalCalibrationCell (const HGCalDDDConstants *cons)
	HGCalCalibrationCell ()

Private Attributes
std::vector< int >	cells_ [4]
const HGCalDDDConstants *	cons_
double	radius_ [2]
std::unique_ptr< HGCalCell >	wafer_

Detailed Description

Definition at line 10 of file HGCalCalibrationCell.h.

Constructor & Destructor Documentation

HGCalCalibrationCell() [1/2]

HGCalCalibrationCell::HGCalCalibrationCell

(

const HGCalDDDConstants *

cons

)

Definition at line 10 of file HGCalCalibrationCell.cc.

References HGCalParameters::calibCellFullHD_, HGCalParameters::calibCellFullLD_, HGCalParameters::calibCellPartHD_, HGCalParameters::calibCellPartLD_, HGCalParameters::calibCellRHD_, HGCalParameters::calibCellRLD_, cells_, cons_, mps_fire::end, HGCalDDDConstants::getParameter(), isotrackApplyRegressor::k, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, radius_, mathSSE::sqrt(), wafer_, and HGCalParameters::waferSize_.

                                                                        : cons_(cons) {
  wafer_ = std::make_unique<HGCalCell>(
      cons_->getParameter()->waferSize_, cons_->getParameter()->nCellsFine_, cons_->getParameter()->nCellsCoarse_);

  radius_[0] = cons_->getParameter()->calibCellRHD_;
  cells_[0].insert(
      cells_[0].end(), cons_->getParameter()->calibCellFullHD_.begin(), cons_->getParameter()->calibCellFullHD_.end());
  cells_[1].insert(
      cells_[1].end(), cons_->getParameter()->calibCellPartHD_.begin(), cons_->getParameter()->calibCellPartHD_.end());
  radius_[1] = cons_->getParameter()->calibCellRLD_;
  cells_[2].insert(
      cells_[2].end(), cons_->getParameter()->calibCellFullLD_.begin(), cons_->getParameter()->calibCellFullLD_.end());
  cells_[3].insert(
      cells_[3].end(), cons_->getParameter()->calibCellPartLD_.begin(), cons_->getParameter()->calibCellPartLD_.end());
  for (int k = 0; k < 2; ++k)
    radius_[k] *= radius_[k];

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalCalibrationCell: " << cells_[0].size() << " HD calibration cells of radius "
                                << std::sqrt(radius_[0]);
  for (unsigned int k = 0; k < cells_[0].size(); ++k)
    edm::LogVerbatim("HGCalGeom") << " [" << k << "] " << cells_[0][k] << ":" << cells_[1][k];
  edm::LogVerbatim("HGCalGeom") << "HGCalCalibrationCell: " << cells_[2].size() << " LD calibration cells of radius "
                                << std::sqrt(radius_[1]);
  for (unsigned int k = 0; k < cells_[2].size(); ++k)
    edm::LogVerbatim("HGCalGeom") << " [" << k << "] " << cells_[2][k] << ":" << cells_[3][k];
#endif
}

HGCalCalibrationCell() [2/2]

HGCalCalibrationCell::HGCalCalibrationCell ( )

inline

Definition at line 13 of file HGCalCalibrationCell.h.

{}

Member Function Documentation

findCell()

int HGCalCalibrationCell::findCell	(	int	zside,
		int	layer,
		int	waferU,
		int	waferV,
		int	cellUV,
		const std::pair< double, double > &	xy
	)		const

Definition at line 39 of file HGCalCalibrationCell.cc.

References SplitLinear::begin, HGCalCell::cellPlacementIndex(), cells_, cons_, mps_fire::end, spr::find(), info(), nano_mu_digi_cff::layer, HGCalDDDConstants::layerType(), convertSQLiteXML::ok, diffTwoXMLs::r2, radius_, wafer_, HGCalTypes::WaferCenterB, HGCalTypes::WaferFull, HGCalTypes::WaferHD120, HGCalTypes::WaferHD200, HGCalDDDConstants::waferInfo(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, geometryCSVtoXML::yy, and ecaldqm::zside().

                                                                                                       {
  const auto& info = cons_->waferInfo(layer, waferU, waferV);
  int ld = ((info.type == HGCalTypes::WaferHD120) || (info.type == HGCalTypes::WaferHD200)) ? 1 : 0;
  int part = (info.part == HGCalTypes::WaferFull) ? 1 : 0;
  int indx = 2 * ld + part;
#ifdef EDM_ML_DEBUG
  std::ostringstream st1;
  st1 << "HGCalCalibrationCell::findCell::input " << layer << ":" << waferU << ":" << waferV << ":" << cellUV << ":"
      << xy.first << ":" << xy.second << " Type:Part " << info.type << ":" << info.part << ":" << ld << ":" << part
      << ":" << indx;
#endif
  bool ok = (std::find(cells_[indx].begin(), cells_[indx].end(), cellUV) != cells_[indx].end());
  int retval(-1);
  if (ok) {
    int layertype = (cons_->layerType(layer) == HGCalTypes::WaferCenterB) ? 1 : 0;
    int place = HGCalCell::cellPlacementIndex(zside, layertype, info.orient);
    int cellU = (cellUV / 100) % 100;
    int cellV = cellUV % 100;
    auto xyc = wafer_->cellUV2XY1(cellU, cellV, place, (1 - ld));
    double xx = xy.first - xyc.first;
    double yy = xy.second - xyc.second;
    double r2 = (xx * xx + yy * yy);
    retval = (r2 > radius_[ld]) ? 0 : 1;
#ifdef EDM_ML_DEBUG
    st1 << " layertype:place:cellU:cellV " << layertype << ":" << place << ":" << cellU << ":" << cellV << " xx "
        << xy.first << ":" << xyc.first << ":" << xx << " yy " << xy.second << ":" << xyc.second << ":" << yy << " R2 "
        << r2 << ":" << radius_[ld];
#endif
  }
#ifdef EDM_ML_DEBUG
  st1 << " Return Value " << ok << ":" << retval;
  edm::LogVerbatim("HGCalGeom") << st1.str();
#endif
  return retval;
}

Member Data Documentation

cells_

std::vector<int> HGCalCalibrationCell::cells_[4]

private

Definition at line 21 of file HGCalCalibrationCell.h.

Referenced by findCell(), and HGCalCalibrationCell().

cons_

const HGCalDDDConstants* HGCalCalibrationCell::cons_

private

Definition at line 18 of file HGCalCalibrationCell.h.

Referenced by findCell(), and HGCalCalibrationCell().

radius_

double HGCalCalibrationCell::radius_[2]

private

Definition at line 20 of file HGCalCalibrationCell.h.

Referenced by findCell(), and HGCalCalibrationCell().

wafer_

std::unique_ptr<HGCalCell> HGCalCalibrationCell::wafer_

private

Definition at line 19 of file HGCalCalibrationCell.h.

Referenced by findCell(), and HGCalCalibrationCell().