EEDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_EEDETID_H
#define ECALDETID_EEDETID_H

#include <iosfwd>
#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) {}
  
  // fast  
  EEDetId(int crystal_ix, int crystal_iy, int iz) : 
    DetId( Ecal, EcalEndcap ) {
    id_|=(crystal_iy&0x7f)|((crystal_ix&0x7f)<<7)|((iz>0)?(0x4000):(0));
  }
  // slow
  EEDetId(int i, int j, int iz, int mode);
  
  EEDetId(const DetId& id) : DetId(id){}
  
  EEDetId& operator=(const DetId& id) {id_ = id.rawId(); return *this;}
  
  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 (  ( positiveZ() ? kEEhalf : 0) + 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) {
    return 
      crystal_ix >= IX_MIN && crystal_ix <= IX_MAX &&
      crystal_iy >= IY_MIN && crystal_iy <= IY_MAX &&  
      std::abs(iz)==1 && 
      ( fastValidDetId(crystal_ix,crystal_iy) ||
    slowValidDetId(crystal_ix,crystal_iy) );
  }
  static bool slowValidDetId(int crystal_ix, int crystal_iy);

  static bool fastValidDetId(int crystal_ix, int crystal_iy) {
    float x =  crystal_ix; float y =  crystal_iy;
    float r = (x - 50.5f) * (x - 50.5f) + (y - 50.5f) * (y - 50.5f);
    return r > 12.f * 12.f && r < 48.f * 48.f;
  }

  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