CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
CaloCellGeometry.cc
Go to the documentation of this file.
2 #include <CLHEP/Geometry/Plane3D.h>
3 
8 
9 typedef HepGeom::Vector3D<CCGFloat> FVec3D ;
10 typedef HepGeom::Plane3D<CCGFloat> Plane3D ;
11 
12 typedef HepGeom::Vector3D<double> DVec3D ;
13 typedef HepGeom::Plane3D<double> DPlane3D ;
14 typedef HepGeom::Point3D<double> DPt3D ;
15 typedef std::vector<DPt3D> DPt3DVec ;
16 
17 const float CaloCellGeometry::k_ScaleFromDDDtoGeant ( 0.1 ) ;
18 
20  m_refPoint ( 0., 0., 0. ),
21  m_corners ( ) ,
22  m_parms ( (CCGFloat*) 0 ), m_eta(0), m_phi(0)
23 {}
24 
25 
27 {}
28 
29 
30 
31 
32 
34  const CornersMgr* mgr,
35  const CCGFloat* par ) :
36  m_refPoint ( gp ),
37  m_corners ( mgr ),
38  m_parms ( par ),
39  m_eta(gp.eta()), m_phi(gp.phi())
40 {}
41 
43  const CCGFloat* par ) :
44  m_refPoint ( 0.25*( cv[0].x() + cv[1].x() + cv[2].x() + cv[3].x() ),
45  0.25*( cv[0].y() + cv[1].y() + cv[2].y() + cv[3].y() ),
46  0.25*( cv[0].z() + cv[1].z() + cv[2].z() + cv[3].z() ) ),
47  m_corners ( cv ),
48  m_parms ( par ),
49  m_eta( m_refPoint.eta()), m_phi(m_refPoint.phi())
50 {}
51 
54 {
55  return m_corners ;
56 }
57 
60 {
61  return m_corners ;
62 }
63 
64 std::ostream& operator<<( std::ostream& s, const CaloCellGeometry& cell )
65 {
66  s << ", Center: " << cell.getPosition() << std::endl;
67 
68  if( cell.emptyCorners() )
69  {
70  s << "Corners vector is empty." << std::endl ;
71  }
72  else
73  {
74  const CaloCellGeometry::CornersVec& corners ( cell.getCorners() ) ;
75  for ( unsigned int i ( 0 ) ; i != corners.size() ; ++i )
76  {
77  s << "Corner: " << corners[ i ] << std::endl;
78  }
79  }
80  return s ;
81 }
82 
83 void
85 {
87  const Pt3D gFront ( p.x(), p.y(), p.z() ) ;
88  const DPt3D dgFront ( p.x(), p.y(), p.z() ) ;
89 
90  Pt3D lFront ;
91  assert( 0 != param() ) ;
92  Pt3DVec lc ( 8, Pt3D(0,0,0) ) ;
93  vocalCorners( lc, param(), lFront ) ;
94 
95  DPt3D dlFront ( lFront.x(), lFront.y(), lFront.z() ) ;
96 
97  const Pt3D lBack ( 0.25*(lc[4]+lc[5]+lc[6]+lc[7]) ) ;
98  const DPt3D dlBack ( lBack.x(), lBack.y(), lBack.z() ) ;
99 
100  const Pt3D dlOne ( lc[0].x(), lc[0].y(), lc[0].z() ) ;
101 
102  const CornersVec& cor ( getCorners() ) ;
103  DPt3DVec dkor ( 8, DPt3D(0,0,0) ) ;
104  for( unsigned int i ( 0 ) ; i != 8 ; ++i )
105  {
106  dkor[i] = DPt3D ( cor[i].x(), cor[i].y(), cor[i].z() ) ;
107  }
108 
109  DPt3D dgBack ( 0.25*( dkor[4]+dkor[5]+dkor[6]+dkor[7] ) ) ;
110 
111  const DVec3D dgAxis ( (dgBack-dgFront).unit() ) ;
112 
113  dgBack = ( dgFront + (dlBack-dlFront).mag()*dgAxis ) ;
114  const DPt3D dgOneT ( dgFront + ( dlOne - dlFront ).mag()*( dkor[0] - dgFront ).unit() ) ;
115 
116  const double dlangle ( ( dlBack - dlFront).angle( dlOne - dlFront ) ) ;
117  const double dgangle ( ( dgBack - dgFront).angle( dgOneT- dgFront ) ) ;
118  const double ddangle ( dlangle - dgangle ) ;
119 
120  const DPlane3D dgPl ( dgFront, dgOneT, dgBack ) ;
121  const DPt3D dp2 ( dgFront + dgPl.normal().unit() ) ;
122 
123  const DPt3D dgOne ( dgFront + HepGeom::Rotate3D( -ddangle, dgFront, dp2 )*
124  DVec3D ( dgOneT - dgFront ) ) ;
125 
126  tr = Tr3D( dlFront , dlBack , dlOne ,
127  dgFront , dgBack , dgOne ) ;
128 
129  if( 0 != lptr ) (*lptr) = lc ;
130 }
131 
132 const float*
134  const std::vector<float>& vv ,
136 {
137  const float* pP ( 0 ) ;
138 
139  for( unsigned int ii ( 0 ) ; ii != pvv.size() ; ++ii )
140  {
141  const ParVec& v ( pvv[ii] ) ;
142  assert( v.size() == vv.size() ) ;
143 
144  bool same ( true ) ;
145  for( unsigned int j ( 0 ) ; j != vv.size() ; ++j )
146  {
147  same = same && ( fabs( vv[j] - v[j] )<1.e-6 ) ;
148  if( !same ) break ;
149  }
150  if( same )
151  {
152  pP = &(*v.begin()) ;
153  break ;
154  }
155  }
156  return pP ;
157 }
158 
159 const float*
161  const std::vector<float>& vv ,
164 {
165  const float* pP ( checkParmPtr( vv, pvv ) ) ;
166 
167  if( 0 == pP )
168  {
169  pvv.push_back( ParVec( mgr ) ) ;
170  ParVec& back ( pvv.back() ) ;
171  for( unsigned int i ( 0 ) ; i != vv.size() ; ++i )
172  {
173  back[i] = vv[i] ;
174  }
175  pP = &(*back.begin()) ;
176  }
177  return pP ;
178 }
179 
180 bool
182 {
183  bool ans ( false ) ;
184  const Pt3D p ( point.x(), point.y(), point.z() ) ;
185  const CornersVec& cog ( getCorners() ) ;
186  Pt3D co[8] ;
187  for( unsigned int i ( 0 ) ; i != 8 ; ++i )
188  {
189  co[i] = Pt3D ( cog[i].x(), cog[i].y(), cog[i].z() ) ;
190  }
191 
192  const Plane3D AA ( co[0], co[1], co[2] ) ; // z<0
193  const Plane3D BB ( co[6], co[5], co[4] ) ; // z>0
194 
195  if( AA.distance(p)*BB.distance(p) >= 0 )
196  {
197  const Plane3D CC ( co[0], co[4], co[5] ) ; // x<0
198  const Plane3D DD ( co[2], co[6], co[7] ) ; // x>0
199  if( CC.distance(p)*DD.distance(p) >= 0 )
200  {
201  const Plane3D EE ( co[3], co[7], co[4] ) ; // y<0
202  const Plane3D FF ( co[1], co[5], co[6] ) ; // y>0
203  if( EE.distance(p)*FF.distance(p) >= 0 )
204  {
205  ans = true ;
206  }
207  }
208  }
209  return ans ;
210 }
int i
Definition: DBlmapReader.cc:9
bool emptyCorners() const
EZArrayFL< CCGFloat > ParVec
std::vector< ParVec > ParVecVec
HepGeom::Point3D< double > DPt3D
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
T y() const
Definition: PV3DBase.h:63
const_iterator begin() const
Definition: EZArrayFL.h:63
HepGeom::Transform3D Tr3D
std::vector< Pt3D > Pt3DVec
T eta() const
std::ostream & operator<<(std::ostream &out, const ALILine &li)
Definition: ALILine.cc:187
int ii
Definition: cuy.py:588
float float float z
MgrType::const_reference const_reference
Definition: EZArrayFL.h:29
const CCGFloat * param() const
virtual void getTransform(Tr3D &tr, Pt3DVec *lptr) const
--------— only needed by specific utility; overloaded when needed -—
virtual ~CaloCellGeometry()
string unit
Definition: csvLumiCalc.py:46
T z() const
Definition: PV3DBase.h:64
int j
Definition: DBlmapReader.cc:9
HepGeom::Vector3D< double > DVec3D
CaloCellGeometry::CCGFloat CCGFloat
static const CCGFloat * getParmPtr(const std::vector< CCGFloat > &vd, ParMgr *mgr, ParVecVec &pvv)
static const CCGFloat k_ScaleFromDDDtoGeant
virtual void vocalCorners(Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const =0
dictionary cv
Definition: cuy.py:362
CaloCellGeometry::Pt3D Pt3D
bool inside(const GlobalPoint &point) const
Returns true if the specified point is inside this cell.
size_type size() const
Definition: EZArrayFL.h:81
HepGeom::Plane3D< double > DPlane3D
CaloCellGeometry::Pt3DVec Pt3DVec
HepGeom::Point3D< CCGFloat > Pt3D
Definition: EZMgrFL.h:8
CaloCellGeometry::Tr3D Tr3D
CornersVec & setCorners() const
std::vector< DPt3D > DPt3DVec
Definition: DDAxes.h:10
const GlobalPoint & getPosition() const
Returns the position of reference for this cell.
T x() const
Definition: PV3DBase.h:62
HepGeom::Vector3D< CCGFloat > FVec3D
ROOT::Math::Plane3D Plane3D
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
virtual const CornersVec & getCorners() const =0
Returns the corner points of this cell&#39;s volume.
static const CCGFloat * checkParmPtr(const std::vector< CCGFloat > &vd, ParVecVec &pvv)
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
Definition: DDAxes.h:10