ESDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_ESDETID_H
#define ECALDETID_ESDETID_H

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

class ESDetId : public DetId {
 public:

  enum { Subdet = EcalPreshower } ;
  ESDetId(){}
  ESDetId( uint32_t rawid ) : DetId( rawid ) {}


  ESDetId(int strip, int ixs, int iys, int plane, int iz, bool doverify=false) :
   DetId( Ecal, EcalPreshower ) {
   id_ |=  (strip&0x3F) |
         ((ixs&0x3F)<<6) |
         ((iys&0x3F)<<12) |
         (((plane-1)&0x1)<<18) |
         ((iz>0)?(1<<19):(0))
       ;
   if (doverify) verify(strip, ixs, iys, plane, iz); 
  }

  ESDetId(const DetId& id);
  ESDetId& operator=(const DetId& id);

  EcalSubdetector subdet() const { return EcalSubdetector(subdetId()); }
  int zside() const { return (id_&0x80000)?(1):(-1); }
  int plane() const { return ((id_>>18)&0x1)+1; }
  int six() const { return (id_>>6)&0x3F; }
  int siy() const { return (id_>>12)&0x3F; }
  int strip() const { return (id_&0x3F); }
  int hashedIndex() const ;

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

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

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

  static ESDetId unhashIndex(    int hi ) ;
  static bool    validHashIndex( int hi ) { return ( hi < kSizeForDenseIndexing ) ; }
  static bool validDetId(int istrip, int ixs, int iys, int iplane, int iz) ;
  static void verify(int istrip, int ixs, int iys, int iplane, int iz) ;


  static const int IX_MIN=1;
  static const int IY_MIN=1;
  static const int IX_MAX=40;
  static const int IY_MAX=40;
  static const int ISTRIP_MIN=1;
  static const int ISTRIP_MAX=32;
  static const int PLANE_MIN=1;
  static const int PLANE_MAX=2;
  static const int IZ_NUM=2;

   private :

      enum { kXYMAX=1072,
         kXYMIN=   1,
         kXMAX =  40,
         kYMAX =  40,
         kXMIN =   1,
         kYMIN =   1,
         // now normalize to A-D notation for ease of use
         kNa   =IZ_NUM,
         kNb   =PLANE_MAX - PLANE_MIN + 1,
         kNc   =kXYMAX - kXYMIN + 1,
         kNd   =ISTRIP_MAX - ISTRIP_MIN + 1,
         kLd   =kNd,
         kLc   =kLd*kNc,
         kLb   =kLc*kNb,
         kLa   =kLb*kNa } ;
      
      static const unsigned short hxy1[ kXMAX ][ kYMAX ] ;
      static const unsigned short hx1[ kXYMAX ] ;
      static const unsigned short hy1[ kXYMAX ] ;
      
      static const unsigned short hxy2[ kXMAX ][ kYMAX ] ;
      static const unsigned short hx2[ kXYMAX ] ;
      static const unsigned short hy2[ kXYMAX ] ;

   public:

      enum { kSizeForDenseIndexing = kLa } ;

};

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

inline
ESDetId::ESDetId( const DetId& gen ) : DetId(gen) 
{
#ifdef EDM_ML_DEBUG
   if( !gen.null() &&
       ( gen.det()  != Ecal          ||
         gen.subdetId() != EcalPreshower    ) )
   {
      throw cms::Exception("InvalidDetId");
   }
#endif
}

inline
ESDetId&
ESDetId::operator=( const DetId& gen )
{
#ifdef EDM_ML_DEBUG
   if (!gen.null() &&
       ( gen.det()      != Ecal          ||                              
         gen.subdetId() != EcalPreshower    ) )
   {
      throw cms::Exception("InvalidDetId");
   }
#endif
   id_=gen.rawId();
   return *this;
}


#endif