EBDetId.cc

Go to the documentation of this file.

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

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

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

 
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);
}
  


//Following TB 2004  numbering scheme 
int EBDetId::ic() const {
  int ie = ietaAbs() -1;
  return  (ie * kCrystalsInPhi)
    +  ( positiveZ() ?
     ( kCrystalsInPhi - ( (iphi() -1 ) % kCrystalsInPhi ) )
     : ( ( iphi() -1 ) % kCrystalsInPhi  + 1)  
     );
}


//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 || (kModuleBoundaries + 4)!=std::find( kModuleBoundaries, kModuleBoundaries + 4, ieta );
}

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) {
  int PI = 180;
  int  result = a.iphi() - b.iphi();
  
  while  (result > PI)    result -= 2*PI;
  while  (result <= -PI)  result += 2*PI;
  return abs(result);

}

float EBDetId::approxEta( const DetId id ) {
  if( id.subdetId() == EcalBarrel ) {
    EBDetId ebId( id );
    return ebId.approxEta();
  } else {
    return 0;
  }
}
  

#include<ostream>
std::ostream& operator<<(std::ostream& s,const EBDetId& id) {
  return s << "(EB ieta " << id.ieta() << ", iphi " << id.iphi() 
       << " ; ism " << id.ism() << " , ic " << id.ic()  << ')';
}