/data/refman/pasoursint/CMSSW_5_3_3/src/DataFormats/HcalDetId/src/HcalCastorDetId.cc

Go to the documentation of this file.

#include "DataFormats/HcalDetId/interface/HcalCastorDetId.h"
#include "DataFormats/HcalDetId/interface/HcalGenericDetId.h"
#include "FWCore/Utilities/interface/Exception.h"

HcalCastorDetId::HcalCastorDetId() : DetId() {}

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

void
HcalCastorDetId::buildMe( Section section, 
                          bool    true_for_positive_eta, 
                          int     sector, 
                          int     module          ) 
{
    sector -= 1; // we count sector from 1-16 instead of 0-15
    id_ |= ( true_for_positive_eta << 8 ) | ( sector << 4 ) | module;
}

HcalCastorDetId::HcalCastorDetId( Section section, 
                                  bool    true_for_positive_eta, 
                                  int     sector, 
                                  int     module          ) 
   : DetId(DetId::Calo, SubdetectorId) 
{
   buildMe( section, true_for_positive_eta, sector, module ) ;
}


HcalCastorDetId::HcalCastorDetId( bool true_for_positive_eta, 
                                  int sector, 
                                  int module ) 
   : DetId( DetId::Calo, SubdetectorId ) 
{
   buildMe( Section(Unknown), true_for_positive_eta, sector, module ) ;
}

HcalCastorDetId::HcalCastorDetId( const DetId& gen ) 
{
    if( !gen.null() &&
        ( gen.det()      != DetId::Calo   ||
          gen.subdetId() != SubdetectorId    ) ) 
    {
       throw cms::Exception("Invalid DetId") << "Cannot initialize CASTORDetId from " << std::hex << gen.rawId() << std::dec; 
    }  
    id_= gen.rawId();
}

HcalCastorDetId& 
HcalCastorDetId::operator=( const DetId& gen ) 
{
   if( !gen.null() &&
       ( gen.det()      != DetId::Calo    ||
         gen.subdetId() != SubdetectorId    ) ) 
    {
       throw cms::Exception("Invalid DetId") << "Cannot assign Castor DetId from " << std::hex << gen.rawId() << std::dec; 
    }
  
   id_ = gen.rawId();
    
   return *this;
}

/*
int HcalCastorDetId::channel() const {
  int channelid = 16*(sector-1)+module;
  return channelid;
}
*/

HcalCastorDetId::Section 
HcalCastorDetId::section() const
{
   const int mod = module();

   Section sect ;
   if ( mod <= 2 )
   {
      sect = HcalCastorDetId::EM ; 
   }
   else
   {
      if( mod > 2 && mod <= 14 )
      {
         sect = HcalCastorDetId::HAD ;
      }
      else
      {
         sect = HcalCastorDetId::Unknown;
      }
   }
   return sect ;
}

uint32_t 
HcalCastorDetId::denseIndex() const
{
   return ( kNumberCellsPerEnd*( zside() + 1 )/2 +
            kNumberSectorsPerEnd*( module() - 1 ) + sector() - 1 ) ;
}

bool 
HcalCastorDetId::validDetId( Section iSection,
                             bool    posEta,
                             int     iSector,
                             int     iModule  )
{
   return ( 0                   <   iSector &&
            kNumberSectorsPerEnd >= iSector && 
            0                   <   iModule &&
            kNumberModulesPerEnd >= iModule    ) ;
}

HcalCastorDetId 
HcalCastorDetId::detIdFromDenseIndex( uint32_t di ) 
{
   return HcalCastorDetId( kNumberCellsPerEnd <= di ,
                           di%kNumberSectorsPerEnd + 1 ,
                           (di%kNumberCellsPerEnd)/kNumberSectorsPerEnd + 1) ;
}

std::ostream& operator<<(std::ostream& s,const HcalCastorDetId& id) 
{
    s << "(CASTOR" << ((id.zside()==1)?("+"):("-"));
    
    switch (id.section()) 
    {
    case(HcalCastorDetId::EM) : s << " EM "; break;
    case(HcalCastorDetId::HAD) : s << " HAD "; break;
    default : s <<" UNKNOWN ";
    }
  
    return s << id.sector() << ',' << id.module() << ',' << ')';
}