BTLNumberingScheme.cc

Go to the documentation of this file.

#include "Geometry/MTDCommonData/interface/BTLNumberingScheme.h"
#include "DataFormats/ForwardDetId/interface/BTLDetId.h"

#include <iostream>
#include <cstring>

//#define EDM_ML_DEBUG

BTLNumberingScheme::BTLNumberingScheme() : 
  MTDNumberingScheme() {
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "Creating BTLNumberingScheme";
#endif  
}

BTLNumberingScheme::~BTLNumberingScheme() {
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "Deleting BTLNumberingScheme";
#endif  
}

uint32_t BTLNumberingScheme::getUnitID(const MTDBaseNumber& baseNumber) const {

  const uint32_t nLevels ( baseNumber.getLevels() ) ;

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "BTLNumberingScheme geometry levels = " << nLevels;
#endif  

  if( 9 > nLevels )
  {   
     edm::LogWarning("MTDGeom") << "BTLNumberingScheme::getUnitID(): "
                                << "Not enough levels found in MTDBaseNumber ( "
                                <<  nLevels 
                                << ") Returning 0" ;
     return 0;
  }

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << baseNumber.getLevelName(0) << ", "
                          << baseNumber.getLevelName(1) << ", "
                          << baseNumber.getLevelName(2) << ", "
                          << baseNumber.getLevelName(3) << ", "
                          << baseNumber.getLevelName(4) << ", "
                          << baseNumber.getLevelName(5) << ", "
                          << baseNumber.getLevelName(6) << ", "
                          << baseNumber.getLevelName(7) << ", "
                          << baseNumber.getLevelName(8);
#endif

  const uint32_t crystal ( baseNumber.getCopyNumber( 0 ) ) ;

  const uint32_t modCopy ( baseNumber.getCopyNumber( 2 ) ) ;
  const uint32_t rodCopy ( baseNumber.getCopyNumber( 3 ) ) ;

  const std::string& modName ( baseNumber.getLevelName( 2 ) ) ; // name of module volume
  uint32_t pos = modName.find("Positive");

  const uint32_t zside = ( pos <= strlen(modName.c_str()) ? 1 : 0 );
  std::string baseName = modName.substr(modName.find(":")+1);

  // trick to accomodate both 54, 42 and 48 modules designs
  const int modtyp ( ::atoi( (baseName.substr(8,1)).c_str() ) == 9 || 
                     ::atoi( (baseName.substr(8,1)).c_str() ) == 5 ||
                     (::atoi( (baseName.substr(8,1)).c_str() ) == 7 && ::atoi( (baseName.substr(7,1)).c_str() ) == 1 ) ?
                     ::atoi( (baseName.substr(7,1)).c_str() ) + 1 : ::atoi( (baseName.substr(7,1)).c_str() ) ) ;
  
  // error checking

  if( 1 > crystal ||
    64 < crystal ) 
  {
    edm::LogWarning("MTDGeom") << "BTLNumberingScheme::getUnitID(): "
                               << "****************** Bad crystal number = " 
                               << crystal 
                               << ", Volume Number = " 
                               << baseNumber.getCopyNumber(0);
    return 0 ;
  }
  
  if( 1  > modtyp ||
      3 < modtyp    )
  {
    edm::LogWarning("MTDGeom") << "BTLNumberingScheme::getUnitID(): "
                               << "****************** Bad module name = " 
                               << modName 
                               << ", Volume Name = " 
                               << baseNumber.getLevelName(2)              ;
    return 0 ;
  }
  
  if( 1 > modCopy ||
      54 < modCopy    )
    {
      edm::LogWarning("MTDGeom") << "BTLNumberingScheme::getUnitID(): "
                                 << "****************** Bad module copy = " 
                                 << modCopy 
                                 << ", Volume Number = " 
                                 << baseNumber.getCopyNumber(2)              ;
      return 0 ;
    }
  
  if( 1 > rodCopy ||
      36 < rodCopy    )
    {
      edm::LogWarning("MTDGeom") << "BTLNumberingScheme::getUnitID(): "
                                 << "****************** Bad rod copy = " 
                                 << rodCopy  
                                 << ", Volume Number = " 
                                 << baseNumber.getCopyNumber(4)              ;
      return 0 ;
    }
  
  // all inputs are fine. Go ahead and decode
  
  BTLDetId thisBTLdetid( zside, rodCopy, modCopy, modtyp, crystal );
  const int32_t intindex = thisBTLdetid.rawId() ;

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "BTL Numbering scheme: " 
                          << " rod = " << rodCopy
                          << " zside = " << zside
                          << " module = " << modCopy
                          << " modtyp = " << modtyp
                          << " crystal = " << crystal
                          << " Raw Id = " << intindex
                          << thisBTLdetid;
#endif

  return intindex ;  
}