DTNumberingScheme.cc

Go to the documentation of this file.

#include "Geometry/MuonNumbering/interface/DTNumberingScheme.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonDDDConstants.h"
#include <iostream>

//#define LOCAL_DEBUG

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

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

void DTNumberingScheme::initMe ( const MuonDDDConstants& muonConstants ) {
  int theLevelPart=muonConstants.getValue("level");
  theRegionLevel=muonConstants.getValue("mb_region")/theLevelPart;
  theWheelLevel=muonConstants.getValue("mb_wheel")/theLevelPart;
  theStationLevel=muonConstants.getValue("mb_station")/theLevelPart;
  theSuperLayerLevel=muonConstants.getValue("mb_superlayer")/theLevelPart;
  theLayerLevel=muonConstants.getValue("mb_layer")/theLevelPart;
  theWireLevel=muonConstants.getValue("mb_wire")/theLevelPart;
 #ifdef LOCAL_DEBUG
   std::cout << "theRegionLevel " << theRegionLevel <<std::endl;
   std::cout << "theWheelLevel " << theWheelLevel <<std::endl;
   std::cout << "theStationLevel " << theStationLevel <<std::endl;
   std::cout << "theSuperLayerLevel " << theSuperLayerLevel <<std::endl;
   std::cout << "theLayerLevel " << theLayerLevel <<std::endl;
   std::cout << "theWireLevel " << theWireLevel <<std::endl;
 #endif

}

int DTNumberingScheme::baseNumberToUnitNumber(const MuonBaseNumber num){
  
#ifdef LOCAL_DEBUG
  std::cout << "DTNumbering "<<num.getLevels()<<std::endl;
  for (int level=1;level<=num.getLevels();level++) {
    std::cout << level << " " << num.getSuperNo(level)
       << " " << num.getBaseNo(level) << std::endl;
  }
#endif
  if (num.getLevels()!=theWireLevel) {
    std::cout << "DTNS::BNToUN "
       << "BaseNumber has " << num.getLevels() << " levels,"
       << "need "<<theWireLevel<<std::endl;
    return 0;
  }
  

//   // Meaningful ranges are enforced by DTWireId, (which
//   // however allows for 0 in wire, layer, superlayer!!!)
// 
//   if ((wire_id < 1) || (wire_id > 100)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "wire id out of range: ";
//     std::cout << wire_id <<std::endl;
//   }
    
//   if ((layer_id < 1) || (layer_id > 4)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "layer id out of range: ";
//     std::cout << layer_id <<std::endl;
//   }
    
//   if ((superlayer_id < 1) || (superlayer_id > 3)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "super-layer id out of range: ";
//     std::cout << superlayer_id <<std::endl;
//   }


  return getDetId(num);
}

int DTNumberingScheme::getDetId(const MuonBaseNumber num) const {
  
  int wire_id=0;
  int layer_id=0;
  int superlayer_id=0;
  int sector_id=0;
  int station_id=0;
  int wheel_id=0;

  //decode significant barrel levels
  decode(num,
         wire_id,
         layer_id,
         superlayer_id,
         sector_id,
         station_id,
         wheel_id);
  
// These ranges are enforced by DTWireId
//   if ((sector_id < 1) || (sector_id > 14)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "sector id out of range: ";
//     std::cout << sector_id <<std::endl;
//   }
    
//   if ((station_id < 1) || (station_id > 4)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "station id out of range: ";
//     std::cout << station_id <<std::endl;
//   }
    
//   if ((wheel_id < -2) || (wheel_id > 2)) {
//     std::cout << "DTNumberingScheme: ";
//     std::cout << "wheel id out of range: ";
//     std::cout << wheel_id <<std::endl;
//   }
    
  DTWireId id(wheel_id,station_id,sector_id,superlayer_id,layer_id,wire_id);
  
#ifdef LOCAL_DEBUG
  std::cout << "DTNumberingScheme: " << id << std::endl;
#endif
  
  return id.rawId();
}

void DTNumberingScheme::decode(const MuonBaseNumber& num,
                                       int& wire_id,
                                       int& layer_id,
                                       int& superlayer_id,
                                       int& sector_id,
                                       int& station_id,
                                       int& wheel_id) const {
  for (int level=1;level<=num.getLevels();level++) {

    //decode
    if (level==theWheelLevel) {
      const int copyno=num.getBaseNo(level);      
      wheel_id=copyno-2;

    } else if (level==theStationLevel) {
      const int station_tag = num.getSuperNo(level);
      const int copyno = num.getBaseNo(level);
      station_id=station_tag;
      sector_id=copyno+1;   

    } else if (level==theSuperLayerLevel) {
      const int copyno = num.getBaseNo(level);
      superlayer_id = copyno + 1;

    } else if (level==theLayerLevel) {
      const int copyno = num.getBaseNo(level);
      layer_id=copyno+1;

    } else if (level==theWireLevel) {
      const int copyno = num.getBaseNo(level);
      wire_id = copyno+1;
    }
  }
}