EEDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_EEDETID_H
#define ECALDETID_EEDETID_H

#include <ostream>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalScDetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"


class EEDetId : public DetId {
   public:
      enum {
     Subdet=EcalEndcap
      };
      
      EEDetId() {}
      
      EEDetId(uint32_t rawid) : DetId(rawid) {}
      
      EEDetId(int i, int j, int iz, int mode = XYMODE);
      
      EEDetId(const DetId& id);
      
      EEDetId& operator=(const DetId& id);
      
      static EcalSubdetector subdet() { return EcalEndcap;}
      
      int zside() const { return (id_&0x4000)?(1):(-1); }
      
      int ix() const { return (id_>>7)&0x7F; }
      
      int iy() const { return id_&0x7F; }
      
      EcalScDetId sc() const {
     const int scEdge = 5;
     return EcalScDetId(1+(ix()-1)/scEdge, 1+(iy()-1)/scEdge, zside());
      }
      
      int isc() const ;
      
      int ic() const;
      
      int iquadrant() const ;
      
      bool positiveZ() const { return id_&0x4000;}
      
      int iPhiOuterRing() const ; // 1-360 else==0 if not on outer ring!

      static EEDetId idOuterRing( int iPhi , int zEnd ) ;

      int hashedIndex() const 
      {
     const uint32_t jx ( ix() ) ;
     const uint32_t jd ( 2*( iy() - 1 ) + ( jx - 1 )/50 ) ;
     return (  ( zside()<0 ? 0 : kEEhalf ) + kdi[jd] + jx - kxf[jd] ) ;
      }

      uint32_t denseIndex() const { return hashedIndex() ; }

      EEDetId offsetBy( int nrStepsX, int nrStepsY ) const;

      EEDetId switchZSide() const;
 
      static DetId offsetBy( const DetId startId, int nrStepsX, int nrStepsY );
      static DetId switchZSide( const DetId startId );

      static bool validDenseIndex( uint32_t din ) { return validHashIndex( din ) ; }

      static EEDetId detIdFromDenseIndex( uint32_t din ) { return unhashIndex( din ) ; }

      static bool isNextToBoundary(     EEDetId id ) ;
      
      static bool isNextToDBoundary(    EEDetId id ) ;

      static bool isNextToRingBoundary( EEDetId id ) ;

      static EEDetId unhashIndex( int hi ) ;

      static bool validHashIndex( int i ) { return ( i < kSizeForDenseIndexing ) ; }

      static bool validDetId(int crystal_ix, int crystal_iy, int iz);

      static int distanceX(const EEDetId& a,const EEDetId& b);
      
      static int distanceY(const EEDetId& a,const EEDetId& b); 


      static int isc( int iscCol,   // output is 1-316
              int iscRow ) ; // 

      static const int IX_MIN =1;

      static const int IY_MIN =1;
      
      static const int IX_MAX =100;

      static const int IY_MAX =100;

      static const int ISC_MIN=1;

      static const int ICR_MIN=1;

      static const int ISC_MAX=316;

      static const int ICR_MAX=25;
      
      enum {
     kEEhalf = 7324 ,
     kSizeForDenseIndexing = 2*kEEhalf
      };

      static const int XYMODE        = 0;
      static const int SCCRYSTALMODE = 1;
   private:

      bool        isOuterRing() const ;

      static bool isOuterRingXY( int ax, int ay ) ;

      //Functions from B. Kennedy to retrieve ix and iy from SC and Crystal number

      static const int nCols = 10;
      static const int nCrys = 5; /* Number of crystals per row in SC */
      static const int QuadColLimits[nCols+1];
      static const int iYoffset[nCols+1];

      static const unsigned short kxf[2*IY_MAX] ;
      static const unsigned short kdi[2*IY_MAX] ;
  
      int ix( int iSC, int iCrys ) const;
      int iy( int iSC, int iCrys ) const;
      int ixQuadrantOne() const;
      int iyQuadrantOne() const;
};


std::ostream& operator<<(std::ostream& s,const EEDetId& id);

#endif