#include <HGCMouseBite.h>

Public Member Functions
bool	exclude (G4ThreeVector &point, int zside, int waferU, int waferV)

	HGCMouseBite (const HGCalDDDConstants &hgc, const std::vector< double > &angle, double maxLength, bool waferRotate)

Private Attributes
double	cut_

const HGCalDDDConstants &	hgcons_

bool	modeUV_

std::vector< std::pair< double, double > >	projXY_

bool	rot_

Detailed Description

Definition at line 9 of file HGCMouseBite.h.

Constructor & Destructor Documentation

HGCMouseBite::HGCMouseBite	(	const HGCalDDDConstants &	hgc,
		const std::vector< double > &	angle,
		double	maxLength,
		bool	waferRotate
	)

Definition at line 9 of file HGCMouseBite.cc.

References funct::cos(), cut_, HGCalDDDConstants::geomMode(), HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, hgcons_, gen::k, modeUV_, projXY_, and funct::sin().

                          : hgcons_(hgc), cut_(maxL), rot_(rot) {
 
   modeUV_ = ((hgcons_.geomMode() == HGCalGeometryMode::Hexagon8) ||
          (hgcons_.geomMode() == HGCalGeometryMode::Hexagon8Full));
   for (auto ang : angle) {
     projXY_.push_back(std::pair<double,double>(cos(ang*CLHEP::deg),sin(ang*CLHEP::deg)));
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "Creating HGCMosueBite with cut at " << cut_ 
                  << " with mode " << modeUV_ << " along " 
                  << angle.size() << " axes";
   for (unsigned int k=0; k<angle.size(); ++k) 
     edm::LogVerbatim("HGCSim") << "Axis[" << k << "] " << angle[k] 
                    << " with projections " << projXY_[k].first
                    << ":" << projXY_[k].second;
 #endif
 }

Member Function Documentation

bool HGCMouseBite::exclude	(	G4ThreeVector &	point,
		int	zside,
		int	waferU,
		int	waferV
	)

Definition at line 29 of file HGCMouseBite.cc.

References funct::abs(), trackerTree::check(), cut_, PVValHelper::dx, PVValHelper::dy, hgcons_, modeUV_, projXY_, rot_, HGCalDDDConstants::waferPosition(), geometryCSVtoXML::xx, and geometryCSVtoXML::xy.

                            {
   bool check(false);
   std::pair<double,double> xy = (modeUV_ ? 
                  hgcons_.waferPosition(waferU,waferV,false) :
                  hgcons_.waferPosition(waferU,false));
   double xx = (zside > 0) ? xy.first : -xy.first;
   double dx(0), dy(0);
   if (rot_) {
     dx = std::abs(point.y() - xy.second);
     dy = std::abs(point.x() - xx);
   } else {
     dx = std::abs(point.x() - xx);
     dy = std::abs(point.y() - xy.second);
   }
   for (auto proj : projXY_) {
     double dist = dx*proj.first + dy*proj.second;
     if (dist > cut_) {check = true; break;}
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "HGCMouseBite:: Point " << point << " zside "
                  << zside << " wafer " << waferU << ":" << waferV 
                  << " position " << xy.first << ":" << xx << ":" 
                  << xy.second << " dxy " << dx << ":" << dy 
                  << " check " << check;
 #endif
   return check;
 }