CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EcalEndcapGeometry.cc
Go to the documentation of this file.
1 #include "Geometry/CaloGeometry/interface/CaloGenericDetId.h"
2 #include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"
3 #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
4 #include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
5 #include "FWCore/Utilities/interface/Exception.h"
6 #include <CLHEP/Geometry/Point3D.h>
7 #include <CLHEP/Geometry/Plane3D.h>
8 
9 typedef CaloCellGeometry::CCGFloat CCGFloat ;
10 typedef CaloCellGeometry::Pt3D Pt3D ;
11 typedef CaloCellGeometry::Pt3DVec Pt3DVec ;
12 typedef HepGeom::Plane3D<CCGFloat> Pl3D ;
13 
14 EcalEndcapGeometry::EcalEndcapGeometry( void )
15  : _nnmods( 316 ),
16  _nncrys( 25 ),
17  zeP( 0. ),
18  zeN( 0. ),
19  m_wref( 0. ),
20  m_del( 0. ),
21  m_nref( 0 ),
22  m_borderMgr( nullptr ),
23  m_borderPtrVec( nullptr ),
24  m_avgZ( -1 ),
25  m_check( nullptr ),
26  m_cellVec( k_NumberOfCellsForCorners )
27 {
28  m_xlo[0] = 999.;
29  m_xlo[1] = 999.;
30  m_xhi[0] = -999.;
31  m_xhi[1] = -999.;
32  m_ylo[0] = 999.;
33  m_ylo[1] = 999.;
34  m_yhi[0] = -999.;
35  m_yhi[1] = -999.;
36  m_xoff[0] = 0.;
37  m_xoff[1] = 0.;
38  m_yoff[0] = 0.;
39  m_yoff[1] = 0.;
40 }
41 
42 EcalEndcapGeometry::~EcalEndcapGeometry()
43 {
44  if(m_borderPtrVec)
45  {
46  auto ptr = m_borderPtrVec.load(std::memory_order_acquire);
47  for(auto& v: (*ptr)) {
48  if(v) delete v;
49  v = nullptr;
50  }
51  delete m_borderPtrVec ;
52  }
53  delete m_borderMgr ;
54 }
55 
56 unsigned int
57 EcalEndcapGeometry::alignmentTransformIndexLocal( const DetId& id )
58 {
59  const CaloGenericDetId gid ( id ) ;
60 
61  assert( gid.isEE() ) ;
62  unsigned int index ( EEDetId(id).ix()/51 + ( EEDetId(id).zside()<0 ? 0 : 2 ) ) ;
63 
64  return index ;
65 }
66 
67 DetId
68 EcalEndcapGeometry::detIdFromLocalAlignmentIndex( unsigned int iLoc )
69 {
70  return EEDetId( 20 + 50*( iLoc%2 ), 50, 2*( iLoc/2 ) - 1 ) ;
71 }
72 
73 unsigned int
74 EcalEndcapGeometry::alignmentTransformIndexGlobal( const DetId& /*id*/ )
75 {
76  return (unsigned int)DetId::Ecal - 1 ;
77 }
78 
79 void
80 EcalEndcapGeometry::initializeParms()
81 {
82  zeP=0.;
83  zeN=0.;
84  unsigned nP=0;
85  unsigned nN=0;
86  m_nref = 0 ;
87 
88  for( uint32_t i ( 0 ) ; i != m_cellVec.size() ; ++i )
89  {
90  const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
91  if( 0 != cell )
92  {
93  const CCGFloat z ( cell->getPosition().z() ) ;
94  if(z>0.)
95  {
96  zeP+=z;
97  ++nP;
98  }
99  else
100  {
101  zeN+=z;
102  ++nN;
103  }
104  const EEDetId myId ( EEDetId::detIdFromDenseIndex(i) ) ;
105  const unsigned int ix ( myId.ix() ) ;
106  const unsigned int iy ( myId.iy() ) ;
107  if( ix > m_nref ) m_nref = ix ;
108  if( iy > m_nref ) m_nref = iy ;
109  }
110  }
111  if( 0 < nP ) zeP/=(CCGFloat)nP;
112  if( 0 < nN ) zeN/=(CCGFloat)nN;
113 
114  m_xlo[0] = 999 ;
115  m_xhi[0] = -999 ;
116  m_ylo[0] = 999 ;
117  m_yhi[0] = -999 ;
118  m_xlo[1] = 999 ;
119  m_xhi[1] = -999 ;
120  m_ylo[1] = 999 ;
121  m_yhi[1] = -999 ;
122  for( uint32_t i ( 0 ) ; i != m_cellVec.size() ; ++i )
123  {
124  const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
125  if( 0 != cell )
126  {
127  const GlobalPoint p ( cell->getPosition() ) ;
128  const CCGFloat z ( p.z() ) ;
129  const CCGFloat zz ( 0 > z ? zeN : zeP ) ;
130  const CCGFloat x ( p.x()*zz/z ) ;
131  const CCGFloat y ( p.y()*zz/z ) ;
132 
133  if( 0 > z && x < m_xlo[0] ) m_xlo[0] = x ;
134  if( 0 < z && x < m_xlo[1] ) m_xlo[1] = x ;
135  if( 0 > z && y < m_ylo[0] ) m_ylo[0] = y ;
136  if( 0 < z && y < m_ylo[1] ) m_ylo[1] = y ;
137 
138  if( 0 > z && x > m_xhi[0] ) m_xhi[0] = x ;
139  if( 0 < z && x > m_xhi[1] ) m_xhi[1] = x ;
140  if( 0 > z && y > m_yhi[0] ) m_yhi[0] = y ;
141  if( 0 < z && y > m_yhi[1] ) m_yhi[1] = y ;
142  }
143  }
144 
145  m_xoff[0] = ( m_xhi[0] + m_xlo[0] )/2. ;
146  m_xoff[1] = ( m_xhi[1] + m_xlo[1] )/2. ;
147  m_yoff[0] = ( m_yhi[0] + m_ylo[0] )/2. ;
148  m_yoff[1] = ( m_yhi[1] + m_ylo[1] )/2. ;
149 
150  m_del = ( m_xhi[0] - m_xlo[0] + m_xhi[1] - m_xlo[1] +
151  m_yhi[0] - m_ylo[0] + m_yhi[1] - m_ylo[1] ) ;
152 
153  if( 1 != m_nref ) m_wref = m_del/(4.*(m_nref-1)) ;
154 
155  m_xlo[0] -= m_wref/2 ;
156  m_xlo[1] -= m_wref/2 ;
157  m_xhi[0] += m_wref/2 ;
158  m_xhi[1] += m_wref/2 ;
159 
160  m_ylo[0] -= m_wref/2 ;
161  m_ylo[1] -= m_wref/2 ;
162  m_yhi[0] += m_wref/2 ;
163  m_yhi[1] += m_wref/2 ;
164 
165  m_del += m_wref ;
166 /*
167  std::cout<<"zeP="<<zeP<<", zeN="<<zeN<<", nP="<<nP<<", nN="<<nN<<std::endl ;
168 
169  std::cout<<"xlo[0]="<<m_xlo[0]<<", xlo[1]="<<m_xlo[1]<<", xhi[0]="<<m_xhi[0]<<", xhi[1]="<<m_xhi[1]
170  <<"\nylo[0]="<<m_ylo[0]<<", ylo[1]="<<m_ylo[1]<<", yhi[0]="<<m_yhi[0]<<", yhi[1]="<<m_yhi[1]<<std::endl ;
171 
172  std::cout<<"xoff[0]="<<m_xoff[0]<<", xoff[1]"<<m_xoff[1]<<", yoff[0]="<<m_yoff[0]<<", yoff[1]"<<m_yoff[1]<<std::endl ;
173 
174  std::cout<<"nref="<<m_nref<<", m_wref="<<m_wref<<std::endl ;
175 */
176 }
177 
178 
179 unsigned int
180 EcalEndcapGeometry::xindex( CCGFloat x,
181  CCGFloat z ) const
182 {
183  const CCGFloat xlo ( 0 > z ? m_xlo[0] : m_xlo[1] ) ;
184  const int i ( 1 + int( ( x - xlo )/m_wref ) ) ;
185 
186  return ( 1 > i ? 1 :
187  ( m_nref < (unsigned int) i ? m_nref : (unsigned int) i ) ) ;
188 
189 }
190 
191 unsigned int
192 EcalEndcapGeometry::yindex( CCGFloat y,
193  CCGFloat z ) const
194 {
195  const CCGFloat ylo ( 0 > z ? m_ylo[0] : m_ylo[1] ) ;
196  const int i ( 1 + int( ( y - ylo )/m_wref ) ) ;
197 
198  return ( 1 > i ? 1 :
199  ( m_nref < (unsigned int) i ? m_nref : (unsigned int) i ) ) ;
200 }
201 
202 EEDetId
203 EcalEndcapGeometry::gId( float x,
204  float y,
205  float z ) const
206 {
207  const CCGFloat fac ( fabs( ( 0 > z ? zeN : zeP )/z ) ) ;
208  const unsigned int ix ( xindex( x*fac, z ) ) ;
209  const unsigned int iy ( yindex( y*fac, z ) ) ;
210  const unsigned int iz ( z>0 ? 1 : -1 ) ;
211 
212  if( EEDetId::validDetId( ix, iy, iz ) )
213  {
214  return EEDetId( ix, iy, iz ) ; // first try is on target
215  }
216  else // try nearby coordinates, spiraling out from center
217  {
218  for( unsigned int i ( 1 ) ; i != 6 ; ++i )
219  {
220  for( unsigned int k ( 0 ) ; k != 8 ; ++k )
221  {
222  const int jx ( 0 == k || 4 == k || 5 == k ? +i :
223  ( 1 == k || 5 < k ? -i : 0 ) ) ;
224  const int jy ( 2 == k || 4 == k || 6 == k ? +i :
225  ( 3 == k || 5 == k || 7 == k ? -i : 0 ) ) ;
226  if( EEDetId::validDetId( ix + jx, iy + jy, iz ) )
227  {
228  return EEDetId( ix + jx, iy + jy, iz ) ;
229  }
230  }
231  }
232  }
233  return EEDetId() ; // nowhere near any crystal
234 }
235 
236 
237 // Get closest cell, etc...
238 DetId
239 EcalEndcapGeometry::getClosestCell( const GlobalPoint& r ) const
240 {
241  try
242  {
243  EEDetId mycellID ( gId( r.x(), r.y(), r.z() ) ) ; // educated guess
244 
245  if( EEDetId::validDetId( mycellID.ix(),
246  mycellID.iy(),
247  mycellID.zside() ) )
248  {
249  // now get points in convenient ordering
250 
251  Pt3D A;
252  Pt3D B;
253  Pt3D C;
254  Pt3D point(r.x(),r.y(),r.z());
255  // D.K. : equation of plane : AA*x+BB*y+CC*z+DD=0;
256  // finding equation for each edge
257 
258  // ================================================================
259  CCGFloat x,y,z;
260  unsigned offset=0;
261  int zsign=1;
262  //================================================================
263 
264  // compute the distance of the point with respect of the 4 crystal lateral planes
265 
266 
267  if( 0 != getGeometry(mycellID) )
268  {
269  const GlobalPoint& myPosition=getGeometry(mycellID)->getPosition();
270 
271  x=myPosition.x();
272  y=myPosition.y();
273  z=myPosition.z();
274 
275  offset=0;
276  // This will disappear when Andre has applied his fix
277  zsign=1;
278 
279  if(z>0)
280  {
281  if(x>0&&y>0)
282  offset=1;
283  else if(x<0&&y>0)
284  offset=2;
285  else if(x>0&&y<0)
286  offset=0;
287  else if (x<0&&y<0)
288  offset=3;
289  zsign=1;
290  }
291  else
292  {
293  if(x>0&&y>0)
294  offset=3;
295  else if(x<0&&y>0)
296  offset=2;
297  else if(x>0&&y<0)
298  offset=0;
299  else if(x<0&&y<0)
300  offset=1;
301  zsign=-1;
302  }
303  GlobalPoint corners[8];
304  for(unsigned ic=0;ic<4;++ic)
305  {
306  corners[ic]=getGeometry(mycellID)->getCorners()[(unsigned)((zsign*ic+offset)%4)];
307  corners[4+ic]=getGeometry(mycellID)->getCorners()[(unsigned)(4+(zsign*ic+offset)%4)];
308  }
309 
310  CCGFloat SS[4];
311  for (short i=0; i < 4 ; ++i)
312  {
313  A = Pt3D(corners[i%4].x(),corners[i%4].y(),corners[i%4].z());
314  B = Pt3D(corners[(i+1)%4].x(),corners[(i+1)%4].y(),corners[(i+1)%4].z());
315  C = Pt3D(corners[4+(i+1)%4].x(),corners[4+(i+1)%4].y(),corners[4+(i+1)%4].z());
316  Pl3D plane(A,B,C);
317  plane.normalize();
318  CCGFloat distance = plane.distance(point);
319  if (corners[0].z()<0.) distance=-distance;
320  SS[i] = distance;
321  }
322 
323  // Only one move in necessary direction
324 
325  const bool yout ( 0 > SS[0]*SS[2] ) ;
326  const bool xout ( 0 > SS[1]*SS[3] ) ;
327 
328  if( yout || xout )
329  {
330  const int ydel ( !yout ? 0 : ( 0 < SS[0] ? -1 : 1 ) ) ;
331  const int xdel ( !xout ? 0 : ( 0 < SS[1] ? -1 : 1 ) ) ;
332  const unsigned int ix ( mycellID.ix() + xdel ) ;
333  const unsigned int iy ( mycellID.iy() + ydel ) ;
334  const unsigned int iz ( mycellID.zside() ) ;
335  if( EEDetId::validDetId( ix, iy, iz ) )
336  mycellID = EEDetId( ix, iy, iz ) ;
337  }
338 
339  return mycellID;
340  }
341  return DetId(0);
342  }
343  }
344  catch ( cms::Exception &e )
345  {
346  return DetId(0);
347  }
348  return DetId(0);
349 }
350 
351 CaloSubdetectorGeometry::DetIdSet
352 EcalEndcapGeometry::getCells( const GlobalPoint& r,
353  double dR ) const
354 {
355  CaloSubdetectorGeometry::DetIdSet dis ; // return object
356  if( 0.000001 < dR )
357  {
358  if( dR > M_PI/2. ) // this code assumes small dR
359  {
360  dis = CaloSubdetectorGeometry::getCells( r, dR ) ; // base class version
361  }
362  else
363  {
364  const float dR2 ( dR*dR ) ;
365  const float reta ( r.eta() ) ;
366  const float rphi ( r.phi() ) ;
367  const float rx ( r.x() ) ;
368  const float ry ( r.y() ) ;
369  const float rz ( r.z() ) ;
370  const float fac ( std::abs( zeP/rz ) ) ;
371  const float xx ( rx*fac ) ; // xyz at endcap z
372  const float yy ( ry*fac ) ;
373  const float zz ( rz*fac ) ;
374 
375  const float xang ( std::atan( xx/zz ) ) ;
376  const float lowX ( zz>0 ? zz*std::tan( xang - dR ) : zz*std::tan( xang + dR ) ) ;
377  const float highX ( zz>0 ? zz*std::tan( xang + dR ) : zz*std::tan( xang - dR ) ) ;
378  const float yang ( std::atan( yy/zz ) ) ;
379  const float lowY ( zz>0 ? zz*std::tan( yang - dR ) : zz*std::tan( yang + dR ) ) ;
380  const float highY ( zz>0 ? zz*std::tan( yang + dR ) : zz*std::tan( yang - dR ) ) ;
381 
382  const float refxlo ( 0 > rz ? m_xlo[0] : m_xlo[1] ) ;
383  const float refxhi ( 0 > rz ? m_xhi[0] : m_xhi[1] ) ;
384  const float refylo ( 0 > rz ? m_ylo[0] : m_ylo[1] ) ;
385  const float refyhi ( 0 > rz ? m_yhi[0] : m_yhi[1] ) ;
386 
387  if( lowX < refxhi && // proceed if any possible overlap with the endcap
388  lowY < refyhi &&
389  highX > refxlo &&
390  highY > refylo )
391  {
392  const int ix_ctr ( xindex( xx, rz ) ) ;
393  const int iy_ctr ( yindex( yy, rz ) ) ;
394  const int iz ( rz>0 ? 1 : -1 ) ;
395 
396  const int ix_hi ( ix_ctr + int( ( highX - xx )/m_wref ) + 2 ) ;
397  const int ix_lo ( ix_ctr - int( ( xx - lowX )/m_wref ) - 2 ) ;
398 
399  const int iy_hi ( iy_ctr + int( ( highY - yy )/m_wref ) + 2 ) ;
400  const int iy_lo ( iy_ctr - int( ( yy - lowY )/m_wref ) - 2 ) ;
401 
402  for( int kx ( ix_lo ) ; kx <= ix_hi ; ++kx )
403  {
404  if( kx > 0 &&
405  kx <= (int) m_nref )
406  {
407  for( int ky ( iy_lo ) ; ky <= iy_hi ; ++ky )
408  {
409  if( ky > 0 &&
410  ky <= (int) m_nref )
411  {
412  if( EEDetId::validDetId( kx, ky, iz ) ) // reject invalid ids
413  {
414  const EEDetId id ( kx, ky, iz ) ;
415  const CaloCellGeometry* cell = &m_cellVec[ id.denseIndex()];
416  const float eta (cell->etaPos() ) ;
417  const float phi (cell->phiPos() ) ;
418  if( reco::deltaR2( eta, phi, reta, rphi ) < dR2 ) dis.insert( id ) ;
419  }
420  }
421  }
422  }
423  }
424  }
425  }
426  }
427  return dis;
428 }
429 
430 const EcalEndcapGeometry::OrderedListOfEBDetId*
431 EcalEndcapGeometry::getClosestBarrelCells( EEDetId id ) const
432 {
433  OrderedListOfEBDetId* ptr ( nullptr ) ;
434  auto ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
435  if(!ptrVec) {
436  if(0 != id.rawId() && 0 != getGeometry(id)) {
437  const float phi(370.+getGeometry(id)->getPosition().phi().degrees());
438  const int iPhi(1+int(phi)%360) ;
439  const int iz(id.zside()) ;
440  if (!m_borderMgr.load(std::memory_order_acquire)) {
441  EZMgrFL<EBDetId>* expect = nullptr;
442  auto ptrMgr = new EZMgrFL<EBDetId>( 720*9, 9 );
443  bool exchanged = m_borderMgr.compare_exchange_strong(expect, ptrMgr, std::memory_order_acq_rel);
444  if(!exchanged) delete ptrMgr;
445  }
446  VecOrdListEBDetIdPtr* expect = nullptr;
447  auto ptrVec = new VecOrdListEBDetIdPtr();
448  ptrVec->reserve(720);
449  for( unsigned int i ( 0 ) ; i != 720 ; ++i )
450  {
451  const int kz ( 360>i ? -1 : 1 ) ;
452  const int iEta ( kz*85 ) ;
453  const int iEtam1 ( kz*84 ) ;
454  const int iEtam2 ( kz*83 ) ;
455  const int jPhi ( i%360 + 1 ) ;
456  OrderedListOfEBDetId& olist ( *new OrderedListOfEBDetId( m_borderMgr.load(std::memory_order_acquire) ) );
457  olist[0]=EBDetId( iEta , jPhi ) ;
458  olist[1]=EBDetId( iEta , myPhi( jPhi+1 ) ) ;
459  olist[2]=EBDetId( iEta , myPhi( jPhi-1 ) ) ;
460  olist[3]=EBDetId( iEtam1, jPhi ) ;
461  olist[4]=EBDetId( iEtam1, myPhi( jPhi+1 ) ) ;
462  olist[5]=EBDetId( iEtam1, myPhi( jPhi-1 ) ) ;
463  olist[6]=EBDetId( iEta , myPhi( jPhi+2 ) ) ;
464  olist[7]=EBDetId( iEta , myPhi( jPhi-2 ) ) ;
465  olist[8]=EBDetId( iEtam2, jPhi ) ;
466  ptrVec->push_back( &olist ) ;
467  }
468  bool exchanged = m_borderPtrVec.compare_exchange_strong(expect, ptrVec, std::memory_order_acq_rel);
469  if(!exchanged) delete ptrVec;
470  ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
471  ptr = (*ptrVec)[ ( iPhi - 1 ) + ( 0>iz ? 0 : 360 ) ] ;
472  }
473  }
474  return ptr;
475 }
476 
477 void
478 EcalEndcapGeometry::localCorners( Pt3DVec& lc ,
479  const CCGFloat* pv ,
480  unsigned int /*i*/,
481  Pt3D& ref )
482 {
483  TruncatedPyramid::localCorners( lc, pv, ref ) ;
484 }
485 
486 void
487 EcalEndcapGeometry::newCell( const GlobalPoint& f1 ,
488  const GlobalPoint& f2 ,
489  const GlobalPoint& f3 ,
490  const CCGFloat* parm ,
491  const DetId& detId )
492 {
493  const unsigned int cellIndex ( EEDetId( detId ).denseIndex() ) ;
494  m_cellVec[ cellIndex ] =
495  TruncatedPyramid( cornersMgr(), f1, f2, f3, parm ) ;
496  addValidID( detId ) ;
497 }
498 
499 
500 CCGFloat
501 EcalEndcapGeometry::avgAbsZFrontFaceCenter() const
502 {
503  if(!m_check.load(std::memory_order_acquire))
504  {
505  CCGFloat sum ( 0 ) ;
506  for( unsigned int i ( 0 ) ; i != m_cellVec.size() ; ++i )
507  {
508  const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
509  if( 0 != cell )
510  {
511  sum += fabs( cell->getPosition().z() ) ;
512  }
513  }
514  m_avgZ = sum/m_cellVec.size();
515  m_check.store(true, std::memory_order_release);
516  }
517  return m_avgZ;
518 }
519 
520 const CaloCellGeometry*
521 EcalEndcapGeometry::cellGeomPtr( uint32_t index ) const
522 {
523  const CaloCellGeometry* cell ( &m_cellVec[ index ] ) ;
524  return ( m_cellVec.size() < index ||
525  0 == cell->param() ? 0 : cell ) ;
526 }
EEDetId::detIdFromDenseIndex
static EEDetId detIdFromDenseIndex(uint32_t din)
Definition: EEDetId.h:220
degrees
double degrees(double radiants)
def degrees
Definition: EcalEleCalibLooper.cc:411
i
int i
Definition: DBlmapReader.cc:9
B
double_binary B
Definition: DDStreamer.cc:234
CaloCellGeometry::phiPos
float phiPos() const
Definition: CaloCellGeometry.h:78
nullptr
#define nullptr
Definition: GCC11Compatibility.h:38
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
CaloCellGeometry::Pt3DVec
std::vector< Pt3D > Pt3DVec
Definition: CaloCellGeometry.h:57
eta
T eta() const
Definition: Basic3DVectorLD.h:179
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
Pl3D
HepGeom::Plane3D< CCGFloat > Pl3D
Definition: EcalBarrelGeometry.cc:18
CaloSubdetectorGeometry::getCells
virtual DetIdSet getCells(const GlobalPoint &r, double dR) const
Get a list of all cells within a dR of the given cell.
Definition: CaloSubdetectorGeometry.cc:85
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
CCGFloat
CaloCellGeometry::CCGFloat CCGFloat
Definition: CaloCellGeometry.cc:4
PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.dR
tuple dR
Definition: PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.py:19
evf::evtn::offset
unsigned int offset(bool)
Definition: GlobalEventNumber.cc:60
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:81
Pt3D
CaloCellGeometry::Pt3D Pt3D
Definition: CaloCellGeometry.cc:5
DetId
Definition: DetId.h:18
M_PI
#define M_PI
Definition: BFit3D.cc:3
EEDetId::validDetId
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:248
Pt3DVec
CaloCellGeometry::Pt3DVec Pt3DVec
Definition: CaloCellGeometry.cc:6
CaloCellGeometry::Pt3D
HepGeom::Point3D< CCGFloat > Pt3D
Definition: CaloCellGeometry.h:56
EZMgrFL
Definition: EZMgrFL.h:8
reco::deltaR2
T1 deltaR2(T1 eta1, T2 phi1, T3 eta2, T4 phi2)
Definition: deltaR.h:58
TruncatedPyramid
A base class to handle the particular shape of Ecal Xtals. Taken from ORCA Calorimetry Code...
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
CaloCellGeometry::etaPos
float etaPos() const
Definition: CaloCellGeometry.h:77
x
Definition: DDAxes.h:10
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
point
*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
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5