EBDetId.cc

Go to the documentation of this file.

#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "FWCore/Utilities/interface/Exception.h"

//#include <iostream>
#include <algorithm>
const int EBDetId::kModuleBoundaries[4] = { 25, 45, 65, 85 };

// pi / 180.
const float EBDetId::crystalUnitToEta = 0.017453292519943295;

EBDetId
EBDetId::unhashIndex( int hi )  
{
   const int pseudo_eta = hi/MAX_IPHI - MAX_IETA;
   return ( validHashIndex( hi ) ?
        EBDetId(pseudo_eta<0 ? pseudo_eta :  pseudo_eta+1, hi%MAX_IPHI+1) :
        EBDetId() ) ;
}
  
EBDetId::EBDetId(int index1, int index2, int mode) 
  : DetId(Ecal,EcalBarrel)
{
  int crystal_ieta;
  int crystal_iphi;
  if (mode == ETAPHIMODE) {
    crystal_ieta = index1;
    crystal_iphi = index2;  
  } else if (mode == SMCRYSTALMODE) {
    int SM = index1;
    int crystal = index2;
    int i = (int)  floor((crystal-1) / kCrystalsInPhi);
    int j = ((crystal-1) - (kCrystalsInPhi*i));
    if (SM <= 18) {
      crystal_ieta = i + 1;
      crystal_iphi = ((SM-1) * kCrystalsInPhi) + (kCrystalsInPhi-j);
    } else {
      crystal_ieta = -(i+1);
      crystal_iphi = ((SM-19) * kCrystalsInPhi) + j+1;
    }
  } else {
    throw cms::Exception("InvalidDetId") << "EBDetId:  Cannot create object.  Unknown mode for (int, int) constructor."; 
  }

  if ( !validDetId(crystal_ieta, crystal_iphi) ) {
    //    std::cout << "crystal_eta " << crystal_ieta << "crystal_phi " << crystal_iphi << std::endl;
    throw cms::Exception("InvalidDetId") << "EBDetId:  Cannot create object.  Indexes out of bounds \n" 
                                         << "eta = " << crystal_ieta << " phi = " << crystal_iphi;
  }
  id_|=((crystal_ieta>0)?(0x10000|(crystal_ieta<<9)):((-crystal_ieta)<<9))|(crystal_iphi&0x1FF);
}
  
EBDetId::EBDetId(const DetId& gen) {
  if (!gen.null() && ( gen.det()!=Ecal || gen.subdetId()!=EcalBarrel )) {
    throw cms::Exception("InvalidDetId");
  }
  id_=gen.rawId();
}

bool EBDetId::validDetId(int i, int j) {

  bool valid = true;
  if (i < -MAX_IETA || i == 0 || i > MAX_IETA ||
      j < MIN_IPHI || j > MAX_IPHI) {
    valid = false;
  }  
  return valid;

}
  
EBDetId& EBDetId::operator=(const DetId& gen) {
  if (!gen.null() && ( gen.det()!=Ecal || gen.subdetId()!=EcalBarrel )) {
    throw cms::Exception("InvalidDetId");
  }
  id_=gen.rawId();
  return *this;
}

//New SM numbering scheme. Avoids discontinuity in phi crossing \eta=0  
int EBDetId::ism() const {
  int id = ( iphi() - 1 ) / kCrystalsInPhi + 1;
  if ( zside() < 0 ) id += 18;
  return id;
}

int EBDetId::im() const {
  for (int i=0; i < 4 ; i++)
    if ( ietaAbs() <= kModuleBoundaries[i] )
      return i+1;
  //Shold never be reached!
  return -1;
}

//Following TB 2004  numbering scheme 
int EBDetId::ic() const {
  int ie = ietaAbs() -1;
  int ip;
  if (zside() < 0) {
    ip = (( iphi() -1 ) % kCrystalsInPhi ) + 1;
  } else {
    ip = kCrystalsInPhi - ((iphi() -1 ) % kCrystalsInPhi );
  }

  return (ie * kCrystalsInPhi) + ip;
}

//Maintains SM crystals in bunch of 1700 indices
int EBDetId::numberBySM() const {
  return (ism()-1) * kCrystalsPerSM + ic() -1;
}

EBDetId EBDetId::offsetBy(int nrStepsEta, int nrStepsPhi ) const
{
        int newEta = ieta()+nrStepsEta;
        if( newEta*ieta() <= 0 ) {
                if( ieta() < 0 ) {
                        newEta++;
                } else if ( ieta() > 0 ) {
                        newEta--;
                }
        }
        int newPhi = iphi() + nrStepsPhi;
        while ( newPhi>360 ) newPhi -= 360;
        while ( newPhi<=0  ) newPhi += 360;

        if( validDetId( newEta, newPhi ) ) {
                return EBDetId( newEta, newPhi);
        } else {
                return EBDetId(0);
        }
}

EBDetId EBDetId::switchZSide() const
{
        int newEta = ieta()*-1;
        if( validDetId( newEta, iphi() ) ) {
                return EBDetId( newEta, iphi() );
        } else {
                return EBDetId(0);
        }
}


DetId EBDetId::offsetBy(const DetId startId, int nrStepsEta, int nrStepsPhi )
{
        if( startId.det() == DetId::Ecal && startId.subdetId() == EcalBarrel ) {
                EBDetId ebStartId(startId);
                return ebStartId.offsetBy( nrStepsEta, nrStepsPhi ).rawId();
        } else {
                return DetId(0);
        }
}

DetId EBDetId::switchZSide( const DetId startId )
{
        if( startId.det() == DetId::Ecal && startId.subdetId() == EcalBarrel ) {
                EBDetId ebStartId(startId);
                return ebStartId.switchZSide().rawId();
        } else {
                return DetId(0);
        }
}

//corrects for HB/EB differing iphi=1
int EBDetId::tower_iphi() const { 
  int iphi_simple=((iphi()-1)/5)+1; 
  iphi_simple-=2;
  return ((iphi_simple<=0)?(iphi_simple+72):(iphi_simple));
}


bool EBDetId::isNextToBoundary(EBDetId id) {
        return isNextToEtaBoundary( id ) || isNextToPhiBoundary( id );
}

bool EBDetId::isNextToEtaBoundary(EBDetId id) {
        int ieta = id.ietaSM();
        return ieta == 1 || std::binary_search( kModuleBoundaries, kModuleBoundaries + 4, ieta ) || ieta == 85;
}

bool EBDetId::isNextToPhiBoundary(EBDetId id) {
        int iphi = id.iphiSM();
        return iphi == 1 || iphi == 20;
}

int EBDetId::distanceEta(const EBDetId& a,const EBDetId& b)
{
  if (a.ieta() * b.ieta() > 0)
    return abs(a.ieta()-b.ieta());
  else
    return abs(a.ieta()-b.ieta())-1;
}

int EBDetId::distancePhi(const EBDetId& a,const EBDetId& b)
{
  if (abs(a.iphi() -b.iphi()) > 180)
    return abs(a.iphi() - b.iphi()) - 180;
  else
    return abs(a.iphi()-b.iphi());
}

float EBDetId::approxEta( const DetId id ) {
        if( id.subdetId() == EcalBarrel ) {
                EBDetId ebId( id );
                return ebId.approxEta();
        } else {
                return 0;
        }
}
  
std::ostream& operator<<(std::ostream& s,const EBDetId& id) {
  return s << "(EB ieta " << id.ieta() << ", iphi" << id.iphi() 
       << " ; ism " << id.ism() << " , ic " << id.ic()  << ')';
}