BaseCrystal Class Reference

#include <BaseCrystal.h>

Public Types
typedef ROOT::Math::Plane3D	Plane3D
typedef math::XYZVector	XYZPoint
typedef math::XYZVector	XYZVector

Public Member Functions
	BaseCrystal ()
	Empty constructor. More...
	BaseCrystal (const DetId &cell)
	constructor from DetId More...
const XYZVector &	exitingNormal (const CaloDirection &side) const
	normal exiting vector for the surface More...
const XYZVector &	getAxis () const
	get crystal axis More...
const XYZPoint &	getBackCenter () const
	get front center More...
const Plane3D &	getBackPlane () const
	back plane More...
void	getBackSide (XYZPoint &a, XYZPoint &b, XYZPoint &c, XYZPoint &d) const
	Coordinates of the back side. More...
void	getBackSide (std::vector< XYZPoint > &corners) const
const XYZPoint &	getCenter () const
	get 1/8*(Sum of corners) More...
const XYZPoint &	getCorner (unsigned i) const
	get the i-th corner More...
const DetId &	getDetId () const
	get the DetId More...
void	getDrawingCoordinates (std::vector< float > &x, std::vector< float > &y, std::vector< float > &z) const
	for debugging. More...
const XYZVector &	getFifthEdge () const
	Direction of the fifth edge. More...
const XYZVector &	getFirstEdge () const
	Direction of the first edge. More...
const XYZPoint &	getFrontCenter () const
	get front center More...
const Plane3D &	getFrontPlane () const
	front plane More...
void	getFrontSide (XYZPoint &a, XYZPoint &b, XYZPoint &c, XYZPoint &d) const
	coordinates of the front side More...
void	getFrontSide (std::vector< XYZPoint > &corners) const
const XYZVector &	getLateralEdge (unsigned i) const
	lateral directions More...
void	getLateralEdges (unsigned i, XYZPoint &, XYZPoint &) const
	get the lateral edges More...
const Plane3D &	getLateralPlane (unsigned i) const
	lateral planes More...
void	getLateralSide (unsigned i, XYZPoint &a, XYZPoint &b, XYZPoint &c, XYZPoint &d) const
	Coordinates of the i=th lateral side. More...
void	getLateralSide (unsigned i, std::vector< XYZPoint > &corners) const
const Plane3D &	getPlane (const CaloDirection &side) const
	generic access More...
void	getSide (const CaloDirection &side, XYZPoint &a, XYZPoint &b, XYZPoint &c, XYZPoint &d) const
	generic access More...
void	getSide (const CaloDirection &side, std::vector< XYZPoint > &corners) const
const int	getSubdetNumber () const
	get the subdector More...
void	print () const
void	setCorners (const CaloCellGeometry::CornersVec &vec, const GlobalPoint &pos)
	~BaseCrystal ()
	Copy constructor. More...

Static Public Member Functions
static unsigned	oppositeDirection (unsigned iside)

Private Member Functions
void	computeBasicProperties ()

Private Attributes
XYZPoint	backcenter_
DetId	cellid_
XYZPoint	center_
XYZPoint	corners_ [8]
XYZVector	crystalaxis_
XYZVector	exitingNormal_ [6]
XYZVector	fifthedgedirection_
XYZVector	firstedgedirection_
XYZPoint	frontcenter_
XYZVector	lateraldirection_ [4]
Plane3D	lateralPlane_ [6]
int	subdetn_

Detailed Description

Definition at line 19 of file BaseCrystal.h.

Member Typedef Documentation

typedef ROOT::Math::Plane3D BaseCrystal::Plane3D

Definition at line 26 of file BaseCrystal.h.

typedef math::XYZVector BaseCrystal::XYZPoint

Definition at line 25 of file BaseCrystal.h.

typedef math::XYZVector BaseCrystal::XYZVector

Definition at line 24 of file BaseCrystal.h.

Constructor & Destructor Documentation

BaseCrystal::BaseCrystal ( )

inline

Empty constructor.

Definition at line 31 of file BaseCrystal.h.

{;};

BaseCrystal::BaseCrystal

(

const DetId &

cell

)

constructor from DetId

Definition at line 9 of file BaseCrystal.cc.

References cellid_, DetId::subdetId(), and subdetn_.

                                          :cellid_(cell)
{
  subdetn_ = cellid_.subdetId();
}

BaseCrystal::~BaseCrystal ( )

inline

Copy constructor.

Definition at line 42 of file BaseCrystal.h.

References setCorners().

{;}

Member Function Documentation

void BaseCrystal::computeBasicProperties ( )

private

Definition at line 86 of file BaseCrystal.cc.

References backcenter_, center_, corners_, crystalaxis_, exitingNormal_, fifthedgedirection_, firstedgedirection_, frontcenter_, mps_fire::i, lateraldirection_, and lateralPlane_.

Referenced by getAxis(), and setCorners().

{
  //if(corners_.size()==0) return;
  center_=XYZPoint(0.,0.,0.);  
  for(unsigned ic=0;ic<8;++ic)
    {
      center_+=corners_[ic];
    }
  
  center_*=0.125;

  //  std::cout << " Ncorners ? " << corners_.size() << std::endl;
  frontcenter_ = 0.25*(corners_[0]+corners_[1]+corners_[2]+corners_[3]);
  backcenter_ = 0.25*(corners_[4]+corners_[5]+corners_[6]+corners_[7]);
  crystalaxis_ = backcenter_-frontcenter_;
  firstedgedirection_=-(corners_[1]-corners_[0]).Unit();
  fifthedgedirection_=-(corners_[5]-corners_[4]).Unit();
  //  std::cout << " Direction laterales " << std::endl;
  for(unsigned il=0;il<4;++il)
    {
      lateraldirection_[il]=-(corners_[(il+1)%4]-corners_[il]).Unit();
    }
  
  Plane3D frontPlane((Point)corners_[0],(Point)corners_[1],(Point)corners_[2]);
  Plane3D backPlane ((Point)corners_[4],(Point)corners_[5],(Point)corners_[6]);
  for(unsigned i=0;i<4;++i)
    {
      lateralPlane_[i]=
    Plane3D((Point)corners_[i],(Point)corners_[(i+1)%4],(Point)corners_[i+4]);
    }
  // Front plane i=4 (UP)
  lateralPlane_[4] = frontPlane;
  // Back plane i =5 (DOWN)
  lateralPlane_[5] = backPlane;

  for(unsigned i=0;i<6;++i)
    {
      exitingNormal_[i] = 
    (lateralPlane_[i].Distance(Point(center_.X(),center_.Y(),center_.Z())) < 0.) ? 
    lateralPlane_[i].Normal().Unit() : -lateralPlane_[i].Normal().Unit();
    }
}

const XYZVector& BaseCrystal::exitingNormal ( const CaloDirection &  side ) const

inline

normal exiting vector for the surface

Definition at line 95 of file BaseCrystal.h.

References exitingNormal_, getDrawingCoordinates(), oppositeDirection(), CaloDirectionOperations::Side(), x, y, and z.

Referenced by Crystal::exitingNormal().

{return exitingNormal_[CaloDirectionOperations::Side(side)];};

const XYZVector& BaseCrystal::getAxis ( ) const

inline

get crystal axis

Definition at line 103 of file BaseCrystal.h.

References computeBasicProperties(), crystalaxis_, and print().

Referenced by Crystal::getAxis().

{ return crystalaxis_;}

const XYZPoint& BaseCrystal::getBackCenter ( ) const

inline

get front center

Definition at line 55 of file BaseCrystal.h.

References backcenter_.

Referenced by Crystal::getBackCenter().

{return backcenter_;}

const Plane3D& BaseCrystal::getBackPlane ( ) const

inline

back plane

Definition at line 85 of file BaseCrystal.h.

References lateralPlane_.

Referenced by Crystal::getBackPlane().

{return lateralPlane_[5];}

void BaseCrystal::getBackSide	(	XYZPoint &	a,
		XYZPoint &	b,
		XYZPoint &	c,
		XYZPoint &	d
	)		const

Coordinates of the back side.

Definition at line 157 of file BaseCrystal.cc.

References corners_.

Referenced by Crystal::getBackSide(), getSide(), and getSubdetNumber().

{
  a=corners_[4];
  b=corners_[5];
  c=corners_[6];
  d=corners_[7];
}

void BaseCrystal::getBackSide

(

std::vector< XYZPoint > &

corners

)

const

Definition at line 165 of file BaseCrystal.cc.

References corners_.

{
  if(corners.size()==4)
    {
      corners[0]=corners_[4];
      corners[1]=corners_[5];
      corners[2]=corners_[6];
      corners[3]=corners_[7];
    }
}

const XYZPoint& BaseCrystal::getCenter ( ) const

inline

get 1/8*(Sum of corners)

Definition at line 51 of file BaseCrystal.h.

References center_.

Referenced by Crystal::getCenter().

{return center_;};

const XYZPoint& BaseCrystal::getCorner ( unsigned  i ) const

inline

get the i-th corner

Definition at line 49 of file BaseCrystal.h.

References corners_, and mps_fire::i.

Referenced by Crystal::getCorner().

{ return corners_[i];};

const DetId& BaseCrystal::getDetId ( void  ) const

inline

get the DetId

Definition at line 61 of file BaseCrystal.h.

References cellid_.

{return cellid_;};

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

for debugging.

Definition at line 196 of file BaseCrystal.cc.

References corners_, X, DOFs::Y, and DOFs::Z.

Referenced by exitingNormal(), and Crystal::getDrawingCoordinates().

{
  x.clear();
  y.clear();
  z.clear();

  x.push_back(corners_[0].X());
  x.push_back(corners_[3].X());
  x.push_back(corners_[2].X());
  x.push_back(corners_[1].X());
  x.push_back(corners_[5].X());
  x.push_back(corners_[6].X());
  x.push_back(corners_[7].X());
  x.push_back(corners_[4].X());
  x.push_back(corners_[0].X());
  x.push_back(corners_[1].X());
  x.push_back(corners_[2].X());
  x.push_back(corners_[6].X());
  x.push_back(corners_[5].X());
  x.push_back(corners_[4].X());
  x.push_back(corners_[7].X());
  x.push_back(corners_[3].X());

  y.push_back(corners_[0].Y());
  y.push_back(corners_[3].Y());
  y.push_back(corners_[2].Y());
  y.push_back(corners_[1].Y());
  y.push_back(corners_[5].Y());
  y.push_back(corners_[6].Y());
  y.push_back(corners_[7].Y());
  y.push_back(corners_[4].Y());
  y.push_back(corners_[0].Y());
  y.push_back(corners_[1].Y());
  y.push_back(corners_[2].Y());
  y.push_back(corners_[6].Y());
  y.push_back(corners_[5].Y());
  y.push_back(corners_[4].Y());
  y.push_back(corners_[7].Y());
  y.push_back(corners_[3].Y());

  z.push_back(corners_[0].Z());
  z.push_back(corners_[3].Z());
  z.push_back(corners_[2].Z());
  z.push_back(corners_[1].Z());
  z.push_back(corners_[5].Z());
  z.push_back(corners_[6].Z());
  z.push_back(corners_[7].Z());
  z.push_back(corners_[4].Z());
  z.push_back(corners_[0].Z());
  z.push_back(corners_[1].Z());
  z.push_back(corners_[2].Z());
  z.push_back(corners_[6].Z());
  z.push_back(corners_[5].Z());
  z.push_back(corners_[4].Z());
  z.push_back(corners_[7].Z());
  z.push_back(corners_[3].Z());
}

const XYZVector& BaseCrystal::getFifthEdge ( ) const

inline

Direction of the fifth edge.

Definition at line 59 of file BaseCrystal.h.

References fifthedgedirection_.

Referenced by Crystal::getFifthEdge().

{return  fifthedgedirection_;}

const XYZVector& BaseCrystal::getFirstEdge ( ) const

inline

Direction of the first edge.

Definition at line 57 of file BaseCrystal.h.

References firstedgedirection_.

Referenced by Crystal::getFirstEdge().

{return  firstedgedirection_;}

const XYZPoint& BaseCrystal::getFrontCenter ( ) const

inline

get front center

Definition at line 53 of file BaseCrystal.h.

References frontcenter_.

Referenced by Crystal::getFrontCenter().

{return frontcenter_;};

const Plane3D& BaseCrystal::getFrontPlane ( ) const

inline

front plane

Definition at line 83 of file BaseCrystal.h.

References lateralPlane_.

Referenced by Crystal::getFrontPlane().

{return lateralPlane_[4];}

void BaseCrystal::getFrontSide	(	XYZPoint &	a,
		XYZPoint &	b,
		XYZPoint &	c,
		XYZPoint &	d
	)		const

coordinates of the front side

Definition at line 138 of file BaseCrystal.cc.

References corners_.

Referenced by Crystal::getFrontSide(), getSide(), and getSubdetNumber().

{
  a=corners_[0];
  b=corners_[1];
  c=corners_[2];
  d=corners_[3];
}

void BaseCrystal::getFrontSide

(

std::vector< XYZPoint > &

corners

)

const

Definition at line 146 of file BaseCrystal.cc.

References corners_.

{
  if(corners.size()==4)
    {
      corners[0]=corners_[0];
      corners[1]=corners_[1];
      corners[2]=corners_[2];
      corners[3]=corners_[3];
    }
}

const XYZVector& BaseCrystal::getLateralEdge ( unsigned  i ) const

inline

lateral directions

Definition at line 92 of file BaseCrystal.h.

References mps_fire::i, and lateraldirection_.

Referenced by Crystal::getLateralEdge().

{return lateraldirection_[i];};

void BaseCrystal::getLateralEdges	(	unsigned	i,
		XYZPoint &	a,
		XYZPoint &	b
	)		const

get the lateral edges

Definition at line 129 of file BaseCrystal.cc.

References corners_, mps_fire::i, and mitigatedMETSequence_cff::U.

Referenced by Crystal::getLateralEdges(), getLateralSide(), and getSubdetNumber().

{
  if(i<4U) // i >= 0, since i is unsigned
    {
      a=corners_[i];
      b=corners_[i+4]; 
    }
}

const Plane3D& BaseCrystal::getLateralPlane ( unsigned  i ) const

inline

lateral planes

Definition at line 87 of file BaseCrystal.h.

References mps_fire::i, and lateralPlane_.

Referenced by Crystal::getLateralPlane().

{return lateralPlane_[i];};

void BaseCrystal::getLateralSide	(	unsigned	i,
		XYZPoint &	a,
		XYZPoint &	b,
		XYZPoint &	c,
		XYZPoint &	d
	)		const

Coordinates of the i=th lateral side.

Definition at line 176 of file BaseCrystal.cc.

References getLateralEdges(), and mitigatedMETSequence_cff::U.

Referenced by Crystal::getLateralSide(), getSide(), and getSubdetNumber().

{
  if(i<4U) // i >= 0, since i is unsigned
    {
      getLateralEdges(i,a,b);
      getLateralEdges((i+1)%4,c,d);
    }
}

void BaseCrystal::getLateralSide	(	unsigned	i,
		std::vector< XYZPoint > &	corners
	)		const

Definition at line 185 of file BaseCrystal.cc.

References corners_, mps_fire::i, and mitigatedMETSequence_cff::U.

{
  if(corners.size()==4&&i<4U) // i >= 0, since i is unsigned
    {
      corners[0]=corners_[i];
      corners[1]=corners_[i+4];
      corners[2]=corners_[4+(i+1)%4];
      corners[3]=corners_[(i+1)%4];
    }
}

const Plane3D& BaseCrystal::getPlane ( const CaloDirection &  side ) const

inline

generic access

Definition at line 89 of file BaseCrystal.h.

References lateralPlane_, and CaloDirectionOperations::Side().

Referenced by Crystal::getPlane().

{return lateralPlane_[CaloDirectionOperations::Side(side)];}

void BaseCrystal::getSide	(	const CaloDirection &	side,
		XYZPoint &	a,
		XYZPoint &	b,
		XYZPoint &	c,
		XYZPoint &	d
	)		const

generic access

Definition at line 257 of file BaseCrystal.cc.

References DOWN, getBackSide(), getFrontSide(), getLateralSide(), CaloDirectionOperations::Side(), and UP.

Referenced by Crystal::getSide(), and getSubdetNumber().

{
  switch (side)
    {
    case UP:
      getFrontSide(a,b,c,d);
      break;
    case DOWN:
      getBackSide(a,b,c,d);
      break;
    default:
      getLateralSide(CaloDirectionOperations::Side(side),a,b,c,d);
    }
}

void BaseCrystal::getSide	(	const CaloDirection &	side,
		std::vector< XYZPoint > &	corners
	)		const

Definition at line 285 of file BaseCrystal.cc.

References DOWN, getBackSide(), getFrontSide(), getLateralSide(), CaloDirectionOperations::Side(), and UP.

{
  switch (side)
    {
    case UP:
      getFrontSide(corners);
      break;
    case DOWN:
      getBackSide(corners);
      break;
    default:
      getLateralSide(CaloDirectionOperations::Side(side),corners);
    }
}

const int BaseCrystal::getSubdetNumber ( ) const

inline

get the subdector

Definition at line 63 of file BaseCrystal.h.

References a, b, EnergyCorrector::c, edmIntegrityCheck::d, getBackSide(), getFrontSide(), getLateralEdges(), getLateralSide(), getSide(), mps_fire::i, and subdetn_.

Referenced by Crystal::getSubdetNumber().

{return subdetn_;}

static unsigned BaseCrystal::oppositeDirection ( unsigned  iside )

static

Referenced by exitingNormal().

void BaseCrystal::print

(

void

)

const

Definition at line 272 of file BaseCrystal.cc.

References backcenter_, cellid_, center_, corners_, gather_cfg::cout, exitingNormal_, frontcenter_, hcalTTPDigis_cfi::id, and DetId::rawId().

Referenced by getAxis(), and Crystal::print().

                             {
  std::cout << "CellID " << cellid_.rawId() << std::endl;
  std::cout << " Corners " << std::endl;
  for(unsigned ic=0;ic<8;++ic)
    std::cout << corners_[ic] << std::endl;
  std::cout << " Center " << center_ << std::endl;
  std::cout << " Front Center " << frontcenter_ << std::endl;
  std::cout << " Back Center " << backcenter_ << std::endl;
  std::cout << " Normales sortantes " << std::endl;
  for(unsigned id=0;id<6;++id)
    std::cout << exitingNormal_[id] << std::endl;
}

void BaseCrystal::setCorners	(	const CaloCellGeometry::CornersVec &	vec,
		const GlobalPoint &	pos
	)

Definition at line 14 of file BaseCrystal.cc.

References computeBasicProperties(), corners_, EcalBarrel, EcalEndcap, PFRecoTauDiscriminationByIsolation_cfi::offset, EZArrayFL< T >::size(), subdetn_, x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by CaloGeometryHelper::buildCrystalArray(), and ~BaseCrystal().

{
  unsigned ncorners= vec.size();
  if(ncorners!=8) return;

  // This is really a pity to have to make the conversion GlobalPoint to XYZPoint, but the latter have many
  // useful properties (transformations, conversions....) that aren't implemented
  // for the GlobalPoints and GlobalVectors. 

  if(subdetn_==EcalBarrel)
    {
      if(pos.z()>0.)
    {
      for(unsigned ic=0;ic<8;++ic)
        {
          corners_[ic]=XYZPoint(vec[ic].x(),vec[ic].y(),vec[ic].z());
        }
    }
      else
    {
      corners_[0]=XYZPoint(vec[2].x(),vec[2].y(),vec[2].z());
      corners_[1]=XYZPoint(vec[3].x(),vec[3].y(),vec[3].z());
      corners_[2]=XYZPoint(vec[0].x(),vec[0].y(),vec[0].z());
      corners_[3]=XYZPoint(vec[1].x(),vec[1].y(),vec[1].z());
      corners_[4]=XYZPoint(vec[6].x(),vec[6].y(),vec[6].z());
      corners_[5]=XYZPoint(vec[7].x(),vec[7].y(),vec[7].z());
      corners_[6]=XYZPoint(vec[4].x(),vec[4].y(),vec[4].z());
      corners_[7]=XYZPoint(vec[5].x(),vec[5].y(),vec[5].z());
    }
    }
  else if(subdetn_==EcalEndcap)
    {
      double x=pos.x();
      double y=pos.y();
      double z=pos.z();
      unsigned offset=0;
      int zsign=1;
      if(z>0) 
    {
      if(x>0&&y>0) 
        offset=1;
      else  if(x<0&&y>0) 
        offset=2;
      else if(x>0&&y<0) 
        offset=0;
      else if (x<0&&y<0) 
        offset=3;
      zsign=1;
    }
      else
    {
      if(x>0&&y>0) 
        offset=3;
      else if(x<0&&y>0) 
        offset=2;
      else if(x>0&&y<0) 
        offset=0;
      else if(x<0&&y<0) 
        offset=1;
      zsign=-1;
    }
      for(unsigned ic=0;ic<4;++ic)
    {
      unsigned i1=(unsigned)((zsign*ic+offset)%4);
      unsigned i2=i1+4;
      corners_[ic]=XYZPoint(vec[i1].x(),vec[i1].y(),vec[i1].z());
      corners_[4+ic]=XYZPoint(vec[i2].x(),vec[i2].y(),vec[i2].z());
    }
    }
  computeBasicProperties();
}

Member Data Documentation

XYZPoint BaseCrystal::backcenter_

private

Definition at line 116 of file BaseCrystal.h.

Referenced by computeBasicProperties(), getBackCenter(), and print().

DetId BaseCrystal::cellid_

private

Definition at line 112 of file BaseCrystal.h.

Referenced by BaseCrystal(), getDetId(), and print().

XYZPoint BaseCrystal::center_

private

Definition at line 114 of file BaseCrystal.h.

Referenced by computeBasicProperties(), getCenter(), and print().

XYZPoint BaseCrystal::corners_[8]

private

Definition at line 111 of file BaseCrystal.h.

Referenced by computeBasicProperties(), getBackSide(), getCorner(), getDrawingCoordinates(), getFrontSide(), getLateralEdges(), getLateralSide(), print(), and setCorners().

XYZVector BaseCrystal::crystalaxis_

private

Definition at line 119 of file BaseCrystal.h.

Referenced by computeBasicProperties(), and getAxis().

XYZVector BaseCrystal::exitingNormal_[6]

private

Definition at line 122 of file BaseCrystal.h.

Referenced by computeBasicProperties(), exitingNormal(), and print().

XYZVector BaseCrystal::fifthedgedirection_

private

Definition at line 118 of file BaseCrystal.h.

Referenced by computeBasicProperties(), and getFifthEdge().

XYZVector BaseCrystal::firstedgedirection_

private

Definition at line 117 of file BaseCrystal.h.

Referenced by computeBasicProperties(), and getFirstEdge().

XYZPoint BaseCrystal::frontcenter_

private

Definition at line 115 of file BaseCrystal.h.

Referenced by computeBasicProperties(), getFrontCenter(), and print().

XYZVector BaseCrystal::lateraldirection_[4]

private

Definition at line 120 of file BaseCrystal.h.

Referenced by computeBasicProperties(), and getLateralEdge().

Plane3D BaseCrystal::lateralPlane_[6]

private

Definition at line 121 of file BaseCrystal.h.

Referenced by computeBasicProperties(), getBackPlane(), getFrontPlane(), getLateralPlane(), and getPlane().

int BaseCrystal::subdetn_

private

Definition at line 113 of file BaseCrystal.h.

Referenced by BaseCrystal(), getSubdetNumber(), and setCorners().