ME0NumberingScheme.cc

Go to the documentation of this file.

#include "Geometry/MuonNumbering/interface/ME0NumberingScheme.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonDDDConstants.h"
#include "DataFormats/MuonDetId/interface/ME0DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"

#include <iostream>

//#define LOCAL_DEBUG

ME0NumberingScheme::ME0NumberingScheme( const MuonDDDConstants& muonConstants ) {
  initMe(muonConstants);
}

ME0NumberingScheme::ME0NumberingScheme( const DDCompactView& cpv ){
  MuonDDDConstants muonConstants(cpv);
  initMe(muonConstants);
}

void ME0NumberingScheme::initMe ( const MuonDDDConstants& muonConstants ) {
  int theLevelPart= muonConstants.getValue("level");
  theRegionLevel  = muonConstants.getValue("m0_region")/theLevelPart;
  theLayerLevel   = muonConstants.getValue("m0_layer")/theLevelPart;
  theSectorLevel  = muonConstants.getValue("m0_sector")/theLevelPart;
  theRollLevel    = muonConstants.getValue("m0_roll")/theLevelPart;
  theNEtaPart     = muonConstants.getValue("m0_nroll");

  // Debug using LOCAL_DEBUG
  #ifdef LOCAL_DEBUG
    std::cout << "Initialize ME0NumberingScheme" <<std::endl;
    std::cout << "theRegionLevel " << theRegionLevel <<std::endl;
    std::cout << "theLayerLevel "  << theLayerLevel   <<std::endl;
    std::cout << "theSectorLevel " << theSectorLevel <<std::endl;
    std::cout << "theRollLevel "   << theRollLevel   <<std::endl;
    std::cout << "theNEtaPart  "   << theNEtaPart    <<std::endl;
  #endif
  // -----------------------

  // Debug using LogDebug
  std::stringstream DebugStringStream;
  DebugStringStream << "Initialize ME0NumberingScheme" <<std::endl;
  DebugStringStream << "theRegionLevel " << theRegionLevel <<std::endl;
  DebugStringStream << "theLayerLevel "  << theLayerLevel   <<std::endl;
  DebugStringStream << "theSectorLevel " << theSectorLevel <<std::endl;
  DebugStringStream << "theRollLevel "   << theRollLevel   <<std::endl;
  DebugStringStream << "theNEtaPart  "   << theNEtaPart    <<std::endl;
  std::string DebugString = DebugStringStream.str();
  edm::LogVerbatim("ME0NumberingScheme")<<DebugString;
  // --------------------
}

int ME0NumberingScheme::baseNumberToUnitNumber(const MuonBaseNumber& num) {

    edm::LogVerbatim("ME0NumberingScheme")<<"ME0NumberingScheme::baseNumberToUnitNumber BEGIN "<<std::endl;
  // Debug using LOCAL_DEBUG
  #ifdef LOCAL_DEBUG
    std::cout << "ME0Numbering "<<num.getLevels()<<std::endl;
    for (int level=1;level<=num.getLevels();level++) {
      std::cout << "level "<<level << " " << num.getSuperNo(level)
          << " " << num.getBaseNo(level) << std::endl;
    }
  #endif
  // -----------------------

  // Debug using LogDebug
  std::stringstream DebugStringStream;
  DebugStringStream << "ME0Numbering :: number of levels = "<<num.getLevels()<<std::endl;
  DebugStringStream << "Level \t SuperNo \t BaseNo"<<std::endl;
  for (int level=1;level<=num.getLevels();level++) {
    DebugStringStream <<level << " \t " << num.getSuperNo(level)
              << " \t " << num.getBaseNo(level) << std::endl;
  }
  std::string DebugString = DebugStringStream.str();
  edm::LogVerbatim("ME0NumberingScheme")<<DebugString;
  // -----------------------



  int maxLevel = 0;
  if(theNEtaPart==1)  maxLevel = theLayerLevel;
  else                maxLevel = theRollLevel;
  if (num.getLevels()!=maxLevel) {
    throw cms::Exception("MuonNumbering") << "MuonME0NS::BNToUN "
          << "BaseNumber has " << num.getLevels() << " levels,"
          << "need "<<maxLevel<<std::endl;
    return 0;
  }

  int region(0), layer(0), chamber(0), roll(0);

  //decode significant ME0 levels
  
  if (num.getBaseNo(theRegionLevel) == 0) 
    region = 1;
  else                                    
    region =-1;
  layer   = num.getBaseNo(theLayerLevel)+1;
  chamber = num.getBaseNo(theSectorLevel) + 1;

  if(theNEtaPart==1)  roll = 0;
  else                roll = num.getBaseNo(theRollLevel)+1;

  // collect all info
  
  // Debug using LOCAL_DEBUG
  #ifdef LOCAL_DEBUG
    std::cout << "ME0NumberingScheme: Region " << region 
          << " Layer " << layer
          << " Chamber " << chamber << " Roll " << roll << std::endl;
  #endif
  // -----------------------

  // Debug using LogDebug
  edm::LogVerbatim("ME0NumberingScheme") << "ME0NumberingScheme: Region " << region 
                     << " Layer " << layer
                     << " Chamber " << chamber << " Roll " << roll << std::endl;
  // -----------------------                                  

  // Build the actual numbering
  ME0DetId id(region,layer,chamber,roll);
  
  // Debug using LOCAL_DEBUG  
  #ifdef LOCAL_DEBUG
    std::cout << " DetId " << id << std::endl;
  #endif
  // ---------------------
    
  // Debug using LogDebug
  edm::LogVerbatim("ME0NumberingScheme")<< " DetId " << id << std::endl;
  // -----------------------                                    
    
  return id.rawId();
}