HGCMouseBite Class Reference

#include <HGCMouseBite.h>

Public Member Functions
bool	exclude (G4ThreeVector &point, int zside, int wafer)
	HGCMouseBite (const HGCalDDDConstants &hgc, const std::vector< double > &angle, double maxLength, bool waferRotate)

Private Attributes
double	cut_
const HGCalDDDConstants &	hgcons_
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(), gather_cfg::cout, cut_, plotBeamSpotDB::first, gen::k, projXY_, edm::second(), and funct::sin().

                         : hgcons_(hgc), cut_(maxL), rot_(rot) {

  for (auto ang : angle) {
    projXY_.push_back(std::pair<double,double>(cos(ang*CLHEP::deg),sin(ang*CLHEP::deg)));
  }
#ifdef EDM_ML_DEBUG
  std::cout << "Creating HGCMosueBite with cut at " << cut_ << " along " 
        << angle.size() << " axes" << std::endl;
  for (unsigned int k=0; k<angle.size(); ++k) 
    std::cout << "Axis[" << k << "] " << angle[k] << " with projections "
          << projXY_[k].first << ":" << projXY_[k].second << std::endl;
#endif
}

Member Function Documentation

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

Definition at line 25 of file HGCMouseBite.cc.

References funct::abs(), trackerTree::check(), gather_cfg::cout, cut_, hgcons_, projXY_, rot_, HGCalDDDConstants::waferPosition(), geometryCSVtoXML::xx, and geometryCSVtoXML::xy.

Referenced by HGCSD::setDetUnitId().

                                                                     {
  bool check(false);
  std::pair<double,double> xy = hgcons_.waferPosition(wafer,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
  std::cout << "HGCMouseBite:: Point " << point << " zside " << zside
        << " wafer " << wafer << " position " << xy.first << ":" << xx
        << ":" << xy.second << " dxy " << dx << ":" << dy << " check "
        << check << std::endl;
#endif
  return check;
}

Member Data Documentation

double HGCMouseBite::cut_

private

Definition at line 20 of file HGCMouseBite.h.

Referenced by exclude(), and HGCMouseBite().

const HGCalDDDConstants& HGCMouseBite::hgcons_

private

Definition at line 19 of file HGCMouseBite.h.

Referenced by exclude().

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

private

Definition at line 22 of file HGCMouseBite.h.

Referenced by exclude(), and HGCMouseBite().

bool HGCMouseBite::rot_

private

Definition at line 21 of file HGCMouseBite.h.

Referenced by exclude().