CrystalPad Class Reference

#include <CrystalPad.h>

Classes
class	padEqual
	equality operator More...

Public Types
typedef ROOT::Math::Transform3DPJ::Point	Point
typedef ROOT::Math::Transform3DPJ	Transform3D
typedef math::XYZVector	XYZPoint
typedef math::XYZVector	XYZVector

Public Member Functions
const CLHEP::Hep2Vector &	center () const
	get the center More...
	CrystalPad ()
	CrystalPad (unsigned number, const std::vector< CLHEP::Hep2Vector > &corners)
	Order matters. 1234 2341 3412 4123 are ok but not 1324 .... More...
	CrystalPad (unsigned number, int onEcal, const std::vector< XYZPoint > &corners, const XYZPoint &origin, const XYZVector &vec1, const XYZVector &vec2)
	CrystalPad (unsigned number, const std::vector< XYZPoint > &corners, const Transform3D &, double scaf=1., bool bothdirections=false)
	CrystalPad (const CrystalPad &right)
CLHEP::Hep2Vector &	edge (unsigned iside, int n)
	access to the corners in direction iside; n=0,1 More...
CLHEP::Hep2Vector &	edge (CaloDirection)
	access to one corner (NE,NW,SE,SW) More...
void	extrems (double &xmin, double &xmax, double &ymin, double &ymax) const
	xmin xmax, ymin ymax of the quad More...
const std::vector< CLHEP::Hep2Vector > &	getCorners () const
	Check that the point (in the global frame) is inside the crystal. More...
void	getDrawingCoordinates (std::vector< float > &x, std::vector< float > &y) const
	for graphic debugging More...
unsigned	getNumber () const
	access to the number More...
bool	inside (const CLHEP::Hep2Vector &point, bool debug=false) const
	Check that the point (in the local frame) is inside the crystal. More...
bool	operator< (const CrystalPad &quad) const
CrystalPad &	operator= (const CrystalPad &rhs)
bool	operator== (const CrystalPad &quad) const
	get the coordinates in the original frame More...
void	print () const
	print More...
void	resetCorners ()
	Rescale the Quad to allow for some inaccuracy ... More...
void	setSurvivalProbability (double val)
double	survivalProbability () const
	access methods to the survivalProbability More...
	~CrystalPad ()

Private Attributes
CLHEP::Hep2Vector	center_
std::vector< CLHEP::Hep2Vector >	corners_
std::vector< CLHEP::Hep2Vector >	dir_
bool	dummy_
double	epsilon_
unsigned	number_
ROOT::Math::Rotation3D	rotation_
double	survivalProbability_
Transform3D	trans_
XYZVector	translation_
double	yscalefactor_

Static Private Attributes
static std::vector< CLHEP::Hep2Vector >	aVector

Detailed Description

Definition at line 13 of file CrystalPad.h.

Member Typedef Documentation

Point

typedef ROOT::Math::Transform3DPJ::Point CrystalPad::Point

Definition at line 18 of file CrystalPad.h.

Transform3D

typedef ROOT::Math::Transform3DPJ CrystalPad::Transform3D

Definition at line 17 of file CrystalPad.h.

XYZPoint

typedef math::XYZVector CrystalPad::XYZPoint

Definition at line 16 of file CrystalPad.h.

XYZVector

typedef math::XYZVector CrystalPad::XYZVector

Definition at line 15 of file CrystalPad.h.

Constructor & Destructor Documentation

CrystalPad() [1/5]

CrystalPad::CrystalPad ( )

inline

Definition at line 20 of file CrystalPad.h.

References dummy_.

{ dummy_ = true; };

CrystalPad() [2/5]

CrystalPad::CrystalPad	(	unsigned	number,
		const std::vector< CLHEP::Hep2Vector > &	corners
	)

Order matters. 1234 2341 3412 4123 are ok but not 1324 ....

Definition at line 30 of file CrystalPad.cc.

References center_, corners_, gather_cfg::cout, dir_, and dummy_.

    : corners_(corners), dir_(aVector), number_(number), survivalProbability_(1.), center_(0., 0.), epsilon_(0.001) {
  //  std::cout << " Hello " << std::endl;
  if (corners.size() != 4) {
    std::cout << " Try to construct a quadrilateral with " << corners.size() << " points ! " << std::endl;
    dummy_ = true;
  } else {
    dummy_ = false;
    // Set explicity the z to 0 !
    for (unsigned ic = 0; ic < 4; ++ic) {
      dir_[ic] = (corners[(ic + 1) % 4] - corners[ic]).unit();
      center_ += corners_[ic];
    }
    center_ *= 0.25;
  }
  //  std::cout << " End of 1 constructor " << std::endl;
  //  std::cout << " Ncorners " << corners_.size() << std::endl;
  //  std::cout << " Ndirs " << dir_.size() << std::endl;
}

CrystalPad() [3/5]

CrystalPad::CrystalPad	(	unsigned	number,
		int	onEcal,
		const std::vector< XYZPoint > &	corners,
		const XYZPoint &	origin,
		const XYZVector &	vec1,
		const XYZVector &	vec2
	)

Constructor from space points, with the description of the local frame (origin,vec1,vec2) where vec1 is normal to the plane and vec2 in the plane

Definition at line 50 of file CrystalPad.cc.

References center_, TCMET_cfi::corner, corners_, gather_cfg::cout, dir_, dummy_, ROOT::Math::Transform3DPJ::GetDecomposition(), rotation_, Validation_hcalonly_cfi::sign, trans_, and translation_.

    : corners_(aVector), dir_(aVector), number_(number), survivalProbability_(1.), center_(0., 0.), epsilon_(0.001) {
  //  std::cout << " We are in the 2nd constructor " << std::endl;
  if (corners.size() != 4) {
    std::cout << " Try to construct a quadrilateral with " << corners.size() << " points ! " << std::endl;
    dummy_ = true;
  } else {
    dummy_ = false;
    double sign = (onEcal == 1) ? -1. : 1.;

    // the good one in the central
    trans_ = Transform3D((Point)origin,
                         (Point)(origin + vec1),
                         (Point)(origin + vec2),
                         Point(0., 0., 0.),
                         Point(0., 0., sign),
                         Point(0., 1., 0.));
    trans_.GetDecomposition(rotation_, translation_);
    //      std::cout << " Constructor 2; input corners "  << std::endl;
    for (unsigned ic = 0; ic < 4; ++ic) {
      //      std::cout << corners[ic]<< " " ;
      XYZPoint corner = rotation_(corners[ic]) + translation_;
      //      std::cout << corner << std::endl ;
      corners_[ic] = CLHEP::Hep2Vector(corner.X(), corner.Y());
      center_ += corners_[ic];
    }
    for (unsigned ic = 0; ic < 4; ++ic) {
      dir_[ic] = (corners_[(ic + 1) % 4] - corners_[ic]).unit();
    }
    center_ *= 0.25;
  }
  //  std::cout << " End of 2 constructor " << std::endl;
  //  std::cout << " Corners(constructor) " ;
  //  std::cout << corners_[0] << std::endl;
  //  std::cout << corners_[1] << std::endl;
  //  std::cout << corners_[2] << std::endl;
  //  std::cout << corners_[3] << std::endl;
}

CrystalPad() [4/5]

CrystalPad::CrystalPad	(	unsigned	number,
		const std::vector< XYZPoint > &	corners,
		const Transform3D &	trans,
		double	scaf = 1.,
		bool	bothdirections = false
	)

Definition at line 93 of file CrystalPad.cc.

References center_, TCMET_cfi::corner, corners_, gather_cfg::cout, dir_, dummy_, ROOT::Math::Transform3DPJ::GetDecomposition(), rotation_, trans_, translation_, and yscalefactor_.

    : corners_(aVector),
      dir_(aVector),
      number_(number),
      survivalProbability_(1.),
      center_(0., 0.),
      epsilon_(0.001),
      yscalefactor_(scaf) {
  //  std::cout << " We are in the 2nd constructor " << std::endl;
  if (corners.size() != 4) {
    std::cout << " Try to construct a quadrilateral with " << corners.size() << " points ! " << std::endl;
    dummy_ = true;
  } else {
    dummy_ = false;

    // the good one in the central
    trans_ = trans;
    //      std::cout << " Constructor 2; input corners "  << std::endl;
    trans_.GetDecomposition(rotation_, translation_);
    for (unsigned ic = 0; ic < 4; ++ic) {
      XYZPoint corner = rotation_(corners[ic]) + translation_;
      //      std::cout << corner << std::endl ;
      double xscalefactor = (bothdirections) ? yscalefactor_ : 1.;
      corners_[ic] = CLHEP::Hep2Vector(corner.X() * xscalefactor, corner.Y() * yscalefactor_);
      center_ += corners_[ic];
    }
    for (unsigned ic = 0; ic < 4; ++ic) {
      dir_[ic] = (corners_[(ic + 1) % 4] - corners_[ic]).unit();
    }
    center_ *= 0.25;
  }
}

CrystalPad() [5/5]

CrystalPad::CrystalPad

(

const CrystalPad &

right

)

Definition at line 7 of file CrystalPad.cc.

References center_, corners_, dir_, dummy_, epsilon_, number_, and survivalProbability_.

                                              {
  corners_ = right.corners_;
  dir_ = right.dir_;
  number_ = right.number_;
  survivalProbability_ = right.survivalProbability_;
  center_ = right.center_;
  epsilon_ = right.epsilon_;
  dummy_ = right.dummy_;
}

~CrystalPad()

CrystalPad::~CrystalPad ( )

inline

Definition at line 43 of file CrystalPad.h.

{ ; };

Member Function Documentation

center()

const CLHEP::Hep2Vector& CrystalPad::center ( ) const

inline

get the center

Definition at line 95 of file CrystalPad.h.

References center_.

Referenced by operator<().

{ return center_; }

edge() [1/2]

CLHEP::Hep2Vector & CrystalPad::edge	(	unsigned	iside,
		int	n
	)

access to the corners in direction iside; n=0,1

Definition at line 222 of file CrystalPad.cc.

References corners_, and create_idmaps::n.

Referenced by EcalHitMaker::cracksPads(), EcalHitMaker::gapsLifting(), and operator<<().

{ return corners_[(iside + n) % 4]; }

edge() [2/2]

CLHEP::Hep2Vector & CrystalPad::edge

(

CaloDirection

dir

)

access to one corner (NE,NW,SE,SW)

Definition at line 224 of file CrystalPad.cc.

References corners_, gather_cfg::cout, DeadROC_duringRun::dir, NORTHEAST, NORTHWEST, SOUTHEAST, and SOUTHWEST.

                                                   {
  switch (dir) {
    case NORTHWEST:
      return corners_[0];
      break;
    case NORTHEAST:
      return corners_[1];
      break;
    case SOUTHEAST:
      return corners_[2];
      break;
    case SOUTHWEST:
      return corners_[3];
      break;
    default: {
      std::cout << " Serious problem in CrystalPad ! " << dir << std::endl;
      return corners_[0];
    }
  }
  return corners_[0];
}

extrems()

void CrystalPad::extrems	(	double &	xmin,
		double &	xmax,
		double &	ymin,
		double &	ymax
	)		const

xmin xmax, ymin ymax of the quad

Definition at line 246 of file CrystalPad.cc.

References corners_, x, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, y, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

                                                                                     {
  xmin = ymin = 999;
  xmax = ymax = -999;
  for (unsigned ic = 0; ic < 4; ++ic) {
    if (corners_[ic].x() < xmin)
      xmin = corners_[ic].x();
    if (corners_[ic].x() > xmax)
      xmax = corners_[ic].x();
    if (corners_[ic].y() < ymin)
      ymin = corners_[ic].y();
    if (corners_[ic].y() > ymax)
      ymax = corners_[ic].y();
  }
}

getCorners()

const std::vector<CLHEP::Hep2Vector>& CrystalPad::getCorners ( ) const

inline

Check that the point (in the global frame) is inside the crystal.

coordinates of the point in the local frame get the corners

Definition at line 54 of file CrystalPad.h.

References corners_.

{ return corners_; }

getDrawingCoordinates()

void CrystalPad::getDrawingCoordinates	(	std::vector< float > &	x,
		std::vector< float > &	y
	)		const

for graphic debugging

Definition at line 284 of file CrystalPad.cc.

References corners_, x, and y.

                                                                                     {
  x.clear();
  y.clear();
  x.push_back(corners_[0].x());
  x.push_back(corners_[1].x());
  x.push_back(corners_[2].x());
  x.push_back(corners_[3].x());
  x.push_back(corners_[0].x());
  y.push_back(corners_[0].y());
  y.push_back(corners_[1].y());
  y.push_back(corners_[2].y());
  y.push_back(corners_[3].y());
  y.push_back(corners_[0].y());
}

getNumber()

unsigned CrystalPad::getNumber ( ) const

inline

access to the number

Definition at line 73 of file CrystalPad.h.

References number_.

Referenced by CrystalPad::padEqual::operator()(), operator<<(), and operator==().

{ return number_; };

inside()

bool CrystalPad::inside	(	const CLHEP::Hep2Vector &	point,
		bool	debug = false
	)		const

Check that the point (in the local frame) is inside the crystal.

Definition at line 127 of file CrystalPad.cc.

References corners_, dir_, diffTwoXMLs::r1, and diffTwoXMLs::r2.

                                                                       {
  //  std::cout << "Inside " << ppoint <<std::endl;
  //  std::cout << "Corners " << corners_.size() << std::endl;
  //  std::cout << corners_[0] << std::endl;
  //  std::cout << corners_[1] << std::endl;
  //  std::cout << corners_[2] << std::endl;
  //  std::cout << corners_[3] << std::endl;
  //  std::cout << " Got the 2D point " << std::endl;
  CLHEP::Hep2Vector pv0(ppoint - corners_[0]);
  CLHEP::Hep2Vector pv2(ppoint - corners_[2]);
  CLHEP::Hep2Vector n1(pv0 - (pv0 * dir_[0]) * dir_[0]);
  CLHEP::Hep2Vector n2(pv2 - (pv2 * dir_[2]) * dir_[2]);

  //  double N1(n1.mag());
  //  double N2(n2.mag());
  double r1(n1 * n2);
  bool inside1(r1 <= 0.);

  if (!inside1)
    return false;

  //  if(debug)
  //    {
  //      std::cout << n1 << std::endl;
  //      std::cout << n2 << std::endl;
  //      std::cout << r1 << std::endl;
  //      std::cout << inside1 << std::endl;
  //    }

  //  bool close1=(N1<epsilon_||N2<epsilon_);
  //
  //  if(!close1&&!inside1) return false;
  //  std::cout << " First calculation " << std::endl;
  CLHEP::Hep2Vector pv1(ppoint - corners_[1]);
  CLHEP::Hep2Vector pv3(ppoint - corners_[3]);
  CLHEP::Hep2Vector n3(pv1 - (pv1 * dir_[1]) * dir_[1]);
  CLHEP::Hep2Vector n4(pv3 - (pv3 * dir_[3]) * dir_[3]);
  //  double N3(n3.mag());
  //  double N4(n4.mag());
  //  bool close2=(N3<epsilon_||N4<epsilon_);
  double r2(n3 * n4);
  bool inside2(r2 <= 0.);
  //  //  std::cout << " pv1 & pv3 " << pv1.mag() << " " << pv3.mag() << std::endl;
  //  //  double tmp=(pv1-(pv1*dir_[1])*dir_[1])*(pv3-(pv3*dir_[3])*dir_[3]);
  //  //  std::cout << " Computed tmp " << tmp << std::endl;
  //  if(debug)
  //    {
  //      std::cout << n3 << std::endl;
  //      std::cout << n4 << std::endl;
  //      std::cout << r2 << std::endl;
  //      std::cout << inside2 << std::endl;
  //    }
  //  if(!close2&&!inside2) return false;
  //  std::cout << " Second calculation " << std::endl;
  //  std::cout << " True " << std::endl;
  //    return (!close1&&!close2||(close2&&inside1||close1&&inside2));

  return inside2;
}

operator<()

bool CrystalPad::operator< ( const CrystalPad &  quad ) const

inline

Definition at line 89 of file CrystalPad.h.

References center(), and center_.

{ return (center_.mag() < quad.center().mag()); }

operator=()

CrystalPad & CrystalPad::operator=

(

const CrystalPad &

rhs

)

Definition at line 17 of file CrystalPad.cc.

References center_, corners_, dir_, dummy_, epsilon_, number_, and survivalProbability_.

                                                         {
  if (this != &right) {  // don't copy into yourself
    corners_ = right.corners_;
    dir_ = right.dir_;
    number_ = right.number_;
    survivalProbability_ = right.survivalProbability_;
    center_ = right.center_;
    epsilon_ = right.epsilon_;
    dummy_ = right.dummy_;
  }
  return *this;
}

operator==()

bool CrystalPad::operator== ( const CrystalPad &  quad ) const

inline

get the coordinates in the original frame

Definition at line 84 of file CrystalPad.h.

References getNumber().

                                                       {
    //      std::cout << " First " << quad.getCellID() << " Second " << this->getCellID() << std::endl;
    return quad.getNumber() == this->getNumber();
  }

print()

void CrystalPad::print

(

void

)

const

print

Definition at line 205 of file CrystalPad.cc.

References corners_, and gather_cfg::cout.

                             {
  std::cout << " Corners " << std::endl;
  std::cout << corners_[0] << std::endl;
  std::cout << corners_[1] << std::endl;
  std::cout << corners_[2] << std::endl;
  std::cout << corners_[3] << std::endl;
}

resetCorners()

void CrystalPad::resetCorners

(

)

Rescale the Quad to allow for some inaccuracy ...

Definition at line 261 of file CrystalPad.cc.

References center_, and corners_.

                              {
  // Find the centre-of-gravity of the Quad (after re-organization)
  center_ = CLHEP::Hep2Vector(0., 0.);
  for (unsigned ic = 0; ic < 4; ++ic)
    center_ += corners_[ic];
  center_ *= 0.25;

  // Rescale the corners to allow for some inaccuracies in
  // in the inside test
  for (unsigned ic = 0; ic < 4; ++ic)
    corners_[ic] += 0.001 * (corners_[ic] - center_);
}

setSurvivalProbability()

void CrystalPad::setSurvivalProbability ( double  val )

inline

Definition at line 64 of file CrystalPad.h.

References survivalProbability_, and heppy_batch::val.

Referenced by EcalHitMaker::cracksPads().

{ survivalProbability_ = val; };

survivalProbability()

double CrystalPad::survivalProbability ( ) const

inline

access methods to the survivalProbability

Definition at line 63 of file CrystalPad.h.

References survivalProbability_.

{ return survivalProbability_; };

Member Data Documentation

aVector

std::vector< CLHEP::Hep2Vector > CrystalPad::aVector

staticprivate

Definition at line 101 of file CrystalPad.h.

center_

CLHEP::Hep2Vector CrystalPad::center_

private

Definition at line 110 of file CrystalPad.h.

Referenced by center(), CrystalPad(), operator<(), operator=(), and resetCorners().

corners_

std::vector<CLHEP::Hep2Vector> CrystalPad::corners_

private

Definition at line 103 of file CrystalPad.h.

Referenced by CrystalPad(), edge(), extrems(), getCorners(), getDrawingCoordinates(), inside(), operator=(), print(), and resetCorners().

dir_

std::vector<CLHEP::Hep2Vector> CrystalPad::dir_

private

Definition at line 104 of file CrystalPad.h.

Referenced by CrystalPad(), inside(), and operator=().

dummy_

bool CrystalPad::dummy_

private

Definition at line 112 of file CrystalPad.h.

Referenced by CrystalPad(), and operator=().

epsilon_

double CrystalPad::epsilon_

private

Definition at line 111 of file CrystalPad.h.

Referenced by CrystalPad(), and operator=().

number_

unsigned CrystalPad::number_

private

Definition at line 105 of file CrystalPad.h.

Referenced by CrystalPad(), getNumber(), and operator=().

rotation_

ROOT::Math::Rotation3D CrystalPad::rotation_

private

Definition at line 107 of file CrystalPad.h.

Referenced by CrystalPad().

survivalProbability_

double CrystalPad::survivalProbability_

private

Definition at line 109 of file CrystalPad.h.

Referenced by CrystalPad(), operator=(), setSurvivalProbability(), and survivalProbability().

trans_

Transform3D CrystalPad::trans_

private

Definition at line 106 of file CrystalPad.h.

Referenced by CrystalPad().

translation_

XYZVector CrystalPad::translation_

private

Definition at line 108 of file CrystalPad.h.

Referenced by CrystalPad().

yscalefactor_

double CrystalPad::yscalefactor_

private

Definition at line 113 of file CrystalPad.h.

Referenced by CrystalPad().