EcalScDetId.h

Go to the documentation of this file.

// -*- Mode: C++; c-basic-offset: 2; indent-tabs-mode: t; tab-width: 8; -*-
//
// \author Philippe Gras (CEA/Saclay). Code adapted from EEDetId.
//
#ifndef EcalDetId_EcalScDetId_h
#define EcalDetId_EcalScDetId_h

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

class EcalScDetId : public DetId {
 public:

  EcalScDetId();
  
  EcalScDetId(uint32_t rawid);

  EcalScDetId(int ix, int iy, int iz);
  
  EcalScDetId(const DetId& id);

  EcalScDetId& operator=(const DetId& id);

  EcalSubdetector subdet() const { return EcalSubdetector(subdetId()); }
  
  int zside() const { return (id_&0x4000)?(1):(-1); }
  
  int ix() const { return (id_>>7)&0x7F; }
  
  int iy() const { return id_&0x7F; }
  
  int iquadrant() const ;
  
  
  int hashedIndex() const{
    checkHashedIndexMap();
    if(!validDetId(ix(),iy(),zside())) return -1;
    return xyz2HashedIndex[ix()-IX_MIN][iy()-IY_MIN][zside()>0?1:0];
  }

  static EcalScDetId unhashIndex(int hi){
    checkHashedIndexMap();
    if(hi < 0 || hi >= kSizeForDenseIndexing) return EcalScDetId();
    return hashedIndex2DetId[hi];
  }
  
  uint32_t denseIndex() const { return hashedIndex() ; }

  static bool validDenseIndex(uint32_t din) { return din < kSizeForDenseIndexing; }

  
  static bool validHashIndex(int hi) { return validDenseIndex(hi) ; }

  static const int SC_PER_EE_CNT = 312;
  
  static const int IX_MIN=1;

  static const int IY_MIN=1;

  static const int IX_MAX=20;

  static const int IY_MAX=20;

  static const int IHASHED_MIN = 0;

  static const int IHASHED_MAX = SC_PER_EE_CNT*2 - 1;
  
  static bool validDetId(int ix, int iy, int iz) ;

private:
  static void checkHashedIndexMap();
  

  //fields
public:
  enum {
    kSizeForDenseIndexing = SC_PER_EE_CNT * 2
  };

private:
  static const int nEndcaps = 2;
  
  CMS_THREAD_SAFE static short xyz2HashedIndex[IX_MAX][IY_MAX][nEndcaps];
  
  CMS_THREAD_SAFE static EcalScDetId hashedIndex2DetId[kSizeForDenseIndexing];
  
  /*The two arrays are thread safe since they are filled safely using std::call_once and
    then only read and never modified.
   */
};


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


#endif //EcalDetId_EcalScDetId_h not defined