HGCalMouseBite Class Reference

#include <HGCalMouseBite.h>

Public Member Functions
template<class T >
bool	exclude (const T &id)
	HGCalMouseBite (const HGCalDDDConstants &hgc, bool waferRotate=false)

Private Attributes
std::vector< int >	rejectCoarse_
std::vector< int >	rejectFine_

Detailed Description

Definition at line 13 of file HGCalMouseBite.h.

Constructor & Destructor Documentation

HGCalMouseBite::HGCalMouseBite	(	const HGCalDDDConstants &	hgc,
		bool	waferRotate = false
	)

Definition at line 7 of file HGCalMouseBite.cc.

References funct::abs(), angle(), angle_units::operators::convertRadToDeg(), funct::cos(), GOODCOLL_filter_cfg::cut, HGCalDDDConstants::getParameter(), HGCSiliconDetId::HGCalCoarseN, HGCSiliconDetId::HGCalFineN, isotrackApplyRegressor::k, HGCalDDDConstants::mouseBite(), hcal_runs::rf, funct::sin(), mathSSE::sqrt(), funct::tan(), findQualityFiles::v, HGCalDDDConstants::waferHexagon8(), HGCalDDDConstants::waferSize(), and HGCalParameters::waferSize_.

                                                                           {
  if (hgc.waferHexagon8()) {
    const std::vector<double> angle = {90._deg, 30._deg};
    std::vector<std::pair<double, double> > projXY;
    projXY.reserve(angle.size());
    for (auto ang : angle)
      projXY.emplace_back(std::make_pair(cos(ang), sin(ang)));
    const double mousebite(hgc.mouseBite(true));
    const double wafersize(hgc.waferSize(true));
    double cut = wafersize * tan(angle[1]) - mousebite;
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Creating HGCMosueBite with cut at " << cut << " along " << angle.size()
                                  << " axes in wafers of size " << wafersize;
    for (unsigned int k = 0; k < angle.size(); ++k)
      edm::LogVerbatim("HGCalGeom") << "Axis[" << k << "] " << convertRadToDeg(angle[k]) << " with projections "
                                    << projXY[k].first << ":" << projXY[k].second;
#endif
    static const double sqrt3 = std::sqrt(3.0);
    int nf(HGCSiliconDetId::HGCalFineN);
    int nf2 = nf / 2;
    double Rf = wafersize / (3.0 * nf);
    double rf = 0.5 * Rf * sqrt3;
    for (int u = 0; u < 2 * nf; ++u) {
      for (int v = 0; v < 2 * nf; ++v) {
        if (((v - u) < nf) && ((u - v) <= nf)) {
          double yp = std::abs((u - 0.5 * v - nf2) * 2 * rf);
          double xp = std::abs((1.5 * (v - nf) + 1.0) * Rf);
          for (auto proj : projXY) {
            double dist = (rot ? (yp * proj.first + xp * proj.second) : (xp * proj.first + yp * proj.second));
            if (dist > cut) {
              rejectFine_.emplace_back(100 * u + v);
              break;
            }
          }
        }
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalMouseBite:: " << rejectFine_.size()
                                  << " masked u, v's among the fine grain wafers:";
    for (unsigned int k = 0; k < rejectFine_.size(); ++k)
      edm::LogVerbatim("HGCalGeom") << "[" << k << "] = (" << rejectFine_[k] / 100 << ", " << rejectFine_[k] % 100
                                    << ")";
#endif
    int nc(HGCSiliconDetId::HGCalCoarseN);
    int nc2 = nc / 2;
    double Rc = hgc.getParameter()->waferSize_ / (3.0 * nc);
    double rc = 0.5 * Rc * sqrt3;
    for (int u = 0; u < 2 * nc; ++u) {
      for (int v = 0; v < 2 * nc; ++v) {
        if (((v - u) < nc) && ((u - v) <= nc)) {
          double yp = (u - 0.5 * v - nc2) * 2 * rc;
          double xp = (1.5 * (v - nc) + 1.0) * Rc;
          for (auto proj : projXY) {
            double dist = (rot ? (yp * proj.first + xp * proj.second) : (xp * proj.first + yp * proj.second));
            if (dist > cut)
              rejectCoarse_.emplace_back(100 * u + v);
            break;
          }
        }
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalMouseBite:: " << rejectCoarse_.size()
                                  << " masked u, v's among the coarse grain wafers:";
    for (unsigned int k = 0; k < rejectCoarse_.size(); ++k)
      edm::LogVerbatim("HGCalGeom") << "[" << k << "] = (" << rejectCoarse_[k] / 100 << ", " << rejectCoarse_[k] % 100
                                    << ")";
#endif
  }
}

Member Function Documentation

template<class T >

bool HGCalMouseBite::exclude ( const T &  id )

inline

Definition at line 17 of file HGCalMouseBite.h.

References trackerTree::check(), rejectCoarse_, and rejectFine_.

                            {
    int iuv = (100 * id.cellU() + id.cellV());
    bool check = ((id.type() == 0) ? (std::binary_search(rejectFine_.begin(), rejectFine_.end(), iuv))
                                   : (std::binary_search(rejectCoarse_.begin(), rejectCoarse_.end(), iuv)));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalMouseBite:: DetId " << id
                                  << " is checked to be in the list of masked ID's with flag " << check;
#endif
    return check;
  }

Member Data Documentation

std::vector<int> HGCalMouseBite::rejectCoarse_

private

Definition at line 29 of file HGCalMouseBite.h.

Referenced by exclude().

std::vector<int> HGCalMouseBite::rejectFine_

private

Definition at line 29 of file HGCalMouseBite.h.

Referenced by exclude().