EcalBarrelNumberingScheme.cc

Go to the documentation of this file.

// File: EcalBarrelNumberingScheme.cc
// Description: Numbering scheme for barrel electromagnetic calorimeter
#include "Geometry/EcalCommonData/interface/EcalBarrelNumberingScheme.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"

#include <iostream>

EcalBarrelNumberingScheme::EcalBarrelNumberingScheme() : 
  EcalNumberingScheme() {
  edm::LogInfo("EcalGeom") << "Creating EcalBarrelNumberingScheme";
}

EcalBarrelNumberingScheme::~EcalBarrelNumberingScheme() {
  edm::LogInfo("EcalGeom") << "Deleting EcalBarrelNumberingScheme";
}

uint32_t EcalBarrelNumberingScheme::getUnitID(const EcalBaseNumber& baseNumber) const {

  const uint32_t nLevels ( baseNumber.getLevels() ) ;

//  LogDebug("EcalGeom") << "ECalBarrelNumberingScheme geometry levels = " << nLevels;
  

  if( 12 > nLevels )
  {   
     edm::LogWarning("EcalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "Not enough levels found in EcalBaseNumber ( "
                 <<  nLevels 
                 << ") Returning 0" ;
     return 0;
  }

  const std::string& cryName ( baseNumber.getLevelName( 0 ) ) ; // name of crystal volume

  const int cryType ( ::atoi( cryName.c_str() + 5 ) ) ;

  const int off ( 12 < nLevels ? 3 : 0 ) ;

  const uint32_t wallCopy ( baseNumber.getCopyNumber( 3 + off ) ) ;
  const uint32_t hawCopy  ( baseNumber.getCopyNumber( 4 + off ) ) ;
  const uint32_t fawCopy  ( baseNumber.getCopyNumber( 5 + off ) ) ;
  const uint32_t supmCopy ( baseNumber.getCopyNumber( 6 + off ) ) ;

//      LogDebug("EcalGeom") << baseNumber.getLevelName(0) << ", "
//                           << baseNumber.getLevelName(1) << ", "
//                           << baseNumber.getLevelName(2) << ", "
//                           << baseNumber.getLevelName(3) << ", "
//                           << baseNumber.getLevelName(4) << ", "
//                           << baseNumber.getLevelName(5) << ", "
//                           << baseNumber.getLevelName(6) << ", "
//                           << baseNumber.getLevelName(7)         ;

      // error checking

  if( 1  > cryType ||
      17 < cryType    )
  {
     edm::LogWarning("EdalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "****************** Bad crystal name = " 
                 << cryName 
                 << ", Volume Name = " 
                 << baseNumber.getLevelName(0)              ;
     return 0 ;
  }

  if( 1 > wallCopy ||
      5 < wallCopy    )
  {
     edm::LogWarning("EcalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "****************** Bad wall copy = " 
                 << wallCopy 
                 << ", Volume Name = " 
                 << baseNumber.getLevelName(3)              ;
     return 0 ;
  }

  if( 1 > hawCopy ||
      2 < hawCopy    )
  {
     edm::LogWarning("EcalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "****************** Bad haw copy = " 
                 << hawCopy  
                 << ", Volume Name = " 
                 << baseNumber.getLevelName(4)              ;
     return 0 ;
  }

  if( 1  > fawCopy ||
      10 < fawCopy    )
  {
     edm::LogWarning("EcalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "****************** Bad faw copy = " 
                 << fawCopy  
                 << ", Volume Name = " 
                 << baseNumber.getLevelName(5)              ;
     return 0 ;
  }

  if( 1  > supmCopy ||
      36 < supmCopy    )
  {
     edm::LogWarning("EcalGeom") << "ECalBarrelNumberingScheme::getUnitID(): "
                 << "****************** Bad supermodule copy = " 
                 << supmCopy 
                 << ", Volume Name = " 
                 << baseNumber.getLevelName(6)              ;
     return 0 ;
  }

  // all inputs are fine. Go ahead and decode

  const int32_t zsign    ( 18 < supmCopy ? -1 : 1 ) ;

  const int32_t eta      ( 5*( cryType - 1 ) + wallCopy ) ;

  const int32_t phi      ( 18 < supmCopy ?
               20*( supmCopy - 19 ) + 2*( 10 - fawCopy ) + 3 - hawCopy :
               20*( supmCopy -  1 ) + 2*( fawCopy - 1  ) +     hawCopy   ) ;

  const int32_t intindex ( EBDetId( zsign*eta, phi ).rawId() ) ;

/*
  static int count ( 1 ) ;
  if( 0==count%1000 )
  {
     std::cout<<"************************** NLEVELS="<<nLevels
          <<", eta="<<eta<<", phi="<<phi<<", zsign="<<zsign<<std::endl;
  }
  ++count;
*/
//      LogDebug("EcalGeom") << "EcalBarrelNumberingScheme: "
//                           << "supmCopy = " << supmCopy
//                           << ", fawCopy = " << fawCopy
//                           << ", hawCopy = " << hawCopy
//                           << ", wallCopy = " << wallCopy
//                           << ", cryType = " << cryType
//                           << "\n           zsign = "  << zsign 
//                           << ", eta = " << eta 
//                           << ", phi = " << phi 
//                           << ", packed index = 0x" << std::hex << intindex << std::dec ;

  return intindex ;  
}