CMSSW_4_4_3_patch1/src/DataFormats/HcalDetId/src/HcalDetId.cc

Go to the documentation of this file.

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <ostream>

const HcalDetId HcalDetId::Undefined(HcalEmpty,0,0,0);

HcalDetId::HcalDetId() : DetId() {
}

HcalDetId::HcalDetId(uint32_t rawid) : DetId(rawid) {
}

HcalDetId::HcalDetId(HcalSubdetector subdet, int tower_ieta, int tower_iphi, int depth) : DetId(Hcal,subdet) {
  // (no checking at this point!)
  id_ |= ((depth&0x1F)<<14) |
    ((tower_ieta>0)?(0x2000|(tower_ieta<<7)):((-tower_ieta)<<7)) |
    (tower_iphi&0x7F);
}

HcalDetId::HcalDetId(const DetId& gen) {
  if (!gen.null()) {
    HcalSubdetector subdet=(HcalSubdetector(gen.subdetId()));
    if (gen.det()!=Hcal || 
        (subdet!=HcalBarrel && subdet!=HcalEndcap && 
         subdet!=HcalOuter && subdet!=HcalForward ))
      {
        throw cms::Exception("Invalid DetId") << "Cannot initialize HcalDetId from " << std::hex << gen.rawId() << std::dec; 
      }  
  }
  id_=gen.rawId();
}

HcalDetId& HcalDetId::operator=(const DetId& gen) {
  if (!gen.null()) {
    HcalSubdetector subdet=(HcalSubdetector(gen.subdetId()));
    if (gen.det()!=Hcal || 
        (subdet!=HcalBarrel && subdet!=HcalEndcap && 
         subdet!=HcalOuter && subdet!=HcalForward ))
      {
        throw cms::Exception("Invalid DetId") << "Cannot assign HcalDetId from " << std::hex << gen.rawId() << std::dec; 
      }  
  }
  id_=gen.rawId();
  return (*this);
}

int HcalDetId::crystal_iphi_low() const { 
  int simple_iphi=((iphi()-1)*5)+1; 
  simple_iphi+=10;
  return ((simple_iphi>360)?(simple_iphi-360):(simple_iphi));
}

int HcalDetId::crystal_iphi_high() const { 
  int simple_iphi=((iphi()-1)*5)+5; 
  simple_iphi+=10;
  return ((simple_iphi>360)?(simple_iphi-360):(simple_iphi));
}

std::ostream& operator<<(std::ostream& s,const HcalDetId& id) {
  switch (id.subdet()) {
  case(HcalBarrel) : return s << "(HB " << id.ieta() << ',' << id.iphi() << ',' << id.depth() << ')';
  case(HcalEndcap) : return s << "(HE " << id.ieta() << ',' << id.iphi() << ',' << id.depth() << ')';
  case(HcalForward) : return s << "(HF " << id.ieta() << ',' << id.iphi() << ',' << id.depth() << ')';
  case(HcalOuter) : return s << "(HO " << id.ieta() << ',' << id.iphi() << ')';
  default : return s << id.rawId();
  }
}

bool 
HcalDetId::validDetId( HcalSubdetector sd,
                       int             ies,
                       int             ip,
                       int             dp      ) 
{
   const int ie ( abs( ies ) ) ;

   return ( ( ip >=  1         ) &&
            ( ip <= 72         ) &&
            ( dp >=  1         ) &&
            ( ie >=  1         ) &&
            ( ( ( sd == HcalBarrel ) &&
                ( ( ( ie <= 14         ) &&
                    ( dp ==  1         )    ) ||
                  ( ( ( ie == 15 ) || ( ie == 16 ) ) && 
                    ( dp <= 2          )                ) ) ) ||
              (  ( sd == HcalEndcap ) &&
                 ( ( ( ie == 16 ) &&
                     ( dp ==  3 )          ) ||
                   ( ( ie == 17 ) &&
                     ( dp ==  1 )          ) ||
                   ( ( ie >= 18 ) &&
                     ( ie <= 20 ) &&
                     ( dp <=  2 )          ) ||
                   ( ( ie >= 21 ) &&
                     ( ie <= 26 ) &&
                     ( dp <=  2 ) &&
                     ( ip%2 == 1 )         ) ||
                   ( ( ie >= 27 ) &&
                     ( ie <= 28 ) &&
                     ( dp <=  3 ) &&
                     ( ip%2 == 1 )         ) ||
                   ( ( ie == 29 ) &&
                     ( dp <=  2 ) &&
                     ( ip%2 == 1 )         )          )      ) ||
              (  ( sd == HcalOuter ) &&
                 ( ie <= 15 ) &&
                 ( dp ==  4 )           ) ||
              (  ( sd == HcalForward ) &&
                 ( dp <=  2 )          &&
                 ( ( ( ie >= 29 ) &&
                     ( ie <= 39 ) &&
                     ( ip%2 == 1 )    ) ||
                   ( ( ie >= 40 ) &&
                     ( ie <= 41 ) &&
                     ( ip%4 == 3 )         )  ) ) ) ) ;
}

int 
HcalDetId::hashed_index() const 
{
   const HcalSubdetector sd ( subdet()  ) ;
   const int             ip ( iphi()    ) ;
   const int             ie ( ietaAbs() ) ;
   const int             dp ( depth()   ) ;
   const int             zn ( zside() < 0 ? 1 : 0 ) ;

   // HB valid DetIds: phi=1-72,eta=1-14,depth=1; phi=1-72,eta=15-16,depth=1-2
  
   // HE valid DetIds: phi=1-72,eta=16-17,depth=1; phi=1-72,eta=18-20,depth=1-2; 
   //                  phi=1-71(in steps of 2),eta=21-26,depth=1-2; phi=1-71(in steps of 2),eta=27-28,depth=1-3
   //                  phi=1-71(in steps of 2),eta=29,depth=1-2

   // HO valid DetIds: phi=1-72,eta=1-15,depth=4!

   // HF valid DetIds: phi=1-71(in steps of 2),eta=29-39,depth=1-2; phi=3-71(in steps of 4),eta=40-41,depth=1-2

   return ( ( sd == HcalBarrel ) ?
            ( ip - 1 )*18 + dp - 1 + ie - ( ie<16 ? 1 : 0 ) + zn*kHBhalf :
            ( ( sd == HcalEndcap ) ?
              2*kHBhalf + ( ip - 1 )*8 + ( ip/2 )*20 +
              ( ( ie==16 || ie==17 ) ? ie - 16 :
                ( ( ie>=18 && ie<=20 ) ? 2 + 2*( ie - 18 ) + dp - 1 :
                  ( ( ie>=21 && ie<=26 ) ? 8 + 2*( ie - 21 ) + dp - 1 :
                    ( ( ie>=27 && ie<=28 ) ? 20 + 3*( ie - 27 ) + dp - 1 :
                      26 + 2*( ie - 29 ) + dp - 1 ) ) ) ) + zn*kHEhalf :
              ( ( sd == HcalOuter ) ?
                2*kHBhalf + 2*kHEhalf + ( ip - 1 )*15 + ( ie - 1 ) + zn*kHOhalf :
                ( ( sd == HcalForward ) ?
                  2*kHBhalf + 2*kHEhalf + 2*kHOhalf + 
                  ( ( ip - 1 )/4 )*4 + ( ( ip - 1 )/2 )*22 + 
                  2*( ie - 29 ) + ( dp - 1 ) + zn*kHFhalf : -1 ) ) ) ) ; 
}

HcalDetId 
HcalDetId::detIdFromDenseIndex( uint32_t di )
{
    if( validDenseIndex( di ) )
    {
      HcalSubdetector sd ( HcalBarrel ) ;
      int ie ( 0 ) ;
      int ip ( 0 ) ;
      int dp ( 0 ) ;
      int in ( di ) ;
      int iz ( 1 ) ;
      if( in > 2*( kHBhalf + kHEhalf + kHOhalf ) - 1 )
      {
         sd  = HcalForward ;
         in -= 2*( kHBhalf + kHEhalf + kHOhalf ) ; 
         iz  = ( in<kHFhalf ? 1 : -1 ) ;
         in %= kHFhalf ; 
         ip  = 4*( in/48 ) ;
         in %= 48 ;
         ip += 1 + ( in>21 ? 2 : 0 ) ;
         if( 3 == ip%4 ) in -= 22 ;
         ie  = 29 + in/2 ;
         dp  = 1 + in%2 ;
      }
      else
      {
         if( in > 2*( kHBhalf + kHEhalf ) - 1 )
         {
            sd  = HcalOuter ;
            in -= 2*( kHBhalf + kHEhalf ) ; 
            iz  = ( in<kHOhalf ? 1 : -1 ) ;
            in %= kHOhalf ; 
            dp  = 4 ;
            ip  = 1 + in/15 ;
            ie  = 1 + ( in - 15*( ip - 1 ) ) ;
         }
         else
         {
            if( in > 2*kHBhalf - 1 )
            {
               sd  = HcalEndcap ;
               in -= 2*kHBhalf ;
               iz  = ( in<kHEhalf ? 1 : -1 ) ;
               in %= kHEhalf ; 
               ip  = 2*( in/36 ) ;
               in %= 36 ;
               ip += 1 + in/28 ;
               if( 0 == ip%2 ) in %= 28 ;
               ie  = 15 + ( in<2 ? 1 + in : 2 + 
                            ( in<20 ? 1 + ( in - 2 )/2 : 9 +
                              ( in<26 ? 1 + ( in - 20 )/3 : 3 ) ) ) ;
               dp  = ( in<1 ? 3 :
                       ( in<2 ? 1 : 
                         ( in<20 ? 1 + ( in - 2 )%2 : 
                           ( in<26 ? 1 + ( in - 20 )%3 : 
                             ( 1 + ( in - 26 )%2 ) ) ) ) ) ;
            }
            else // barrel
            {
               iz  = ( di<kHBhalf ? 1 : -1 ) ;
               in %= kHBhalf ; 
               ip = in/18 + 1 ;
               in %= 18 ;
               if( in < 14 )
               {
                  dp = 1 ;
                  ie = in + 1 ;
               }
               else
               {
                  in %= 14 ;
                  dp =  1 + in%2 ;
                  ie = 15 + in/2 ;
               }
            }
         }
      }
      return HcalDetId( sd, iz*int(ie), ip, dp ) ;
   }
   else
   {
      return HcalDetId() ;
   }
}