/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/DataFormats/CaloTowers/src/CaloTowerDetId.cc

Go to the documentation of this file.

#include "DataFormats/CaloTowers/interface/CaloTowerDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <iostream>

CaloTowerDetId::CaloTowerDetId() : DetId() {
}
  
CaloTowerDetId::CaloTowerDetId(uint32_t rawid) : DetId(rawid&0xFFF0FFFFu) {
  
}
  
CaloTowerDetId::CaloTowerDetId(int ieta, int iphi) : DetId(Calo,SubdetId) {
  id_|= 
    ((ieta>0)?(0x2000|((ieta&0x3F)<<7)):(((-ieta)&0x3f)<<7)) |
    (iphi&0x7F);
}
  
CaloTowerDetId::CaloTowerDetId(const DetId& gen) {
  if (!gen.null() && (gen.det()!=Calo || gen.subdetId()!=SubdetId)) {
    throw cms::Exception("Invalid DetId") << "Cannot initialize CaloTowerDetId from " << std::hex << gen.rawId() << std::dec;
  }
  id_=gen.rawId(); 
}
  
CaloTowerDetId& CaloTowerDetId::operator=(const DetId& gen) {
  if (!gen.null() && (gen.det()!=Calo || gen.subdetId()!=SubdetId)) {
    throw cms::Exception("Invalid DetId") << "Cannot assign CaloTowerDetId from " << std::hex << gen.rawId() << std::dec;
  }
  id_=gen.rawId();
  return *this;
}

int CaloTowerDetId::iphi() const {
  int retval=id_&0x7F;
  return retval;
}  

bool 
CaloTowerDetId::validDetId( int ie , int ip ) 
{
   const int ia ( abs( ie ) ) ;
   return ( ( ia >= 1 )        &&
            ( ip >= 1 )        &&
            ( ia <= kMaxIEta ) &&
            ( ip <= kMaxIPhi ) &&
            ( ( ( ia <= kBarIEta ) &&
                ( ip <= kBarNPhi )    ) ||
             ( ( ia >  kBarIEta ) &&
                ( ia <= kEndIEta ) &&
                ( (ip-1)%2 == 0  )    ) ||
              ( ( ia >  kEndIEta ) &&
                ( ia <= kForIEta ) &&
                ( (ip-3)%4 == 0 )    )    )  ) ;
}

uint32_t 
CaloTowerDetId::denseIndex() const 
{
   const uint32_t ie ( ietaAbs()     ) ;
   const uint32_t ip ( iphi()    - 1 ) ;
   
   return ( ( 0 > zside() ? 0 : kAllNTot ) +
            ( ( kBarIEta >= ie ? ( ie - 1 )*kBarNPhi + ip :
                ( kEndIEta >= ie ?  kBarNTot + ( ie - 1 - kBarIEta )*kEndNPhi + ip/2 :
                  kBarNTot + kEndNTot + ( ie - 1 - kEndIEta )*kForNPhi + ip/4 ) ) ) ) ;
}

CaloTowerDetId 
CaloTowerDetId::detIdFromDenseIndex( uint32_t din ) 
{
   const int iz ( din < kAllNTot ? -1 : 1 ) ;
   din %= kAllNTot ;
   const uint32_t ie ( ( kBarNTot + kEndNTot ) - 1 < din ?
                       kEndIEta + ( din - kBarNTot - kEndNTot )/kForNPhi + 1 :
                       ( kBarNTot - 1 < din ?
                         kBarIEta + ( din - kBarNTot )/kEndNPhi + 1 :
                         din/kBarNPhi + 1 ) ) ;

   const uint32_t ip ( ( kBarNTot + kEndNTot ) - 1 < din ?
                       ( ( din - kBarNTot - kEndNTot )%kForNPhi )*4 + 3 :
                       ( kBarNTot - 1 < din ?
                         ( ( din - kBarNTot )%kEndNPhi )*2 + 1 :
                         din%kBarNPhi + 1 ) ) ;

   return ( validDenseIndex( din ) ? CaloTowerDetId( iz*ie, ip ) : CaloTowerDetId() ) ;
}

std::ostream& operator<<(std::ostream& s, const CaloTowerDetId& id) {
  return s << "Tower (" << id.ieta() << "," << id.iphi() << ")";
}