CSCNumberingScheme.cc

Go to the documentation of this file.

#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "Geometry/MuonNumbering/interface/CSCNumberingScheme.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonDDDConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>

//#define LOCAL_DEBUG

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

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

void CSCNumberingScheme::initMe (  const MuonDDDConstants& muonConstants ) {
  int theLevelPart=muonConstants.getValue("level");
  theRegionLevel=muonConstants.getValue("me_region")/theLevelPart;
  theStationLevel=muonConstants.getValue("me_station")/theLevelPart;
  theSubringLevel=muonConstants.getValue("me_subring")/theLevelPart;
  theSectorLevel=muonConstants.getValue("me_sector")/theLevelPart;
  theRingLevel=muonConstants.getValue("me_ring")/theLevelPart;
  theLayerLevel=muonConstants.getValue("me_layer")/theLevelPart;
#ifdef LOCAL_DEBUG
  std::cout << "Initialize CSCNumberingScheme" << std::endl;
  std::cout << "theRegionLevel " << theRegionLevel <<std::endl;
  std::cout << "theStationLevel " << theStationLevel <<std::endl;
  std::cout << "theSubringLevel " << theSubringLevel <<std::endl;
  std::cout << "theSectorLevel " << theSectorLevel <<std::endl;
  std::cout << "theRingLevel " << theRingLevel <<std::endl;
  std::cout << "theLayerLevel " << theLayerLevel <<std::endl;
#endif
}

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

#ifdef LOCAL_DEBUG
  std::cout << "CSCNumbering "<<num.getLevels()<<std::endl;
  for (int level=1;level<=num.getLevels();level++) {
    std::cout << level << " " << num.getSuperNo(level)
     << " " << num.getBaseNo(level) << std::endl;
  }
#endif

  int fwbw_id=0;
  int station_id=0;
  int ring_id=0;
  int subring_id=0;
  int sector_id=0;
  int layer_id=0;

  // Decode endcap levels
  // We should be able to work with 6 (layer-level) or 5 (chamber-level)
  
    for (int level=1;level<=num.getLevels();level++) {

    if (level==theRegionLevel) {
      const int copyno=num.getBaseNo(level);
      fwbw_id=copyno+1;
      LogDebug("CSCNumbering") << "endcap=" << fwbw_id ;

    } else if (level==theStationLevel) {
      const int station_tag = num.getSuperNo(level);      
      station_id=station_tag;
      LogDebug("CSCNumbering")<< "station=" << station_id ;
      
    } else if (level==theSubringLevel) {
      const int copyno=num.getBaseNo(level);
      subring_id=copyno+1;
      LogDebug("CSCNumbering")<<"subring=" << subring_id ;

    } else if (level==theSectorLevel) {
      const int copyno=num.getBaseNo(level);
      sector_id=copyno+1;
      LogDebug("CSCNumbering")<< "sector=" << sector_id ;

    } else if (level==theLayerLevel) {
      const int copyno=num.getBaseNo(level);
      layer_id=copyno+1;
      LogDebug("CSCNumbering") << "layer=" << layer_id ;

    } else if (level==theRingLevel) {
      const int ring_tag = num.getSuperNo(level);      
      ring_id=ring_tag;
      LogDebug("CSCNumbering") << "ring=" << ring_id ;
    }
  }      

  // check validity 
    
  if ((fwbw_id < 1) || (fwbw_id > 2)) {
    edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "forward/backward id out of range:" << fwbw_id;
  }
    
  if ((station_id < 1) || (station_id > 4)) {
    edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "station id out of range:" << station_id;
  }
    
   if ((ring_id < 1) || (ring_id > 4)) {
     edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "ring id out of range:" << ring_id ;
   }
    
  if ((subring_id < 1) || (subring_id > 2)) {
    edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "subring id out of range:" << subring_id ;
  }
    
  if ((sector_id < 1) || (sector_id > 36)) {
    edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "sector id out of range:" << sector_id ;
  }

  // Allow id=0 since that means a chamber    
  if ((layer_id < 0) || (layer_id > 6)) {
    edm::LogError("CSCNumbering") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber" << "layer id out of range" << layer_id ;
  }
    
  // find appropriate chamber label
    
   int chamber_id=chamberIndex(station_id, ring_id,
                  subring_id, sector_id);
    
  // convert into raw id of appropriate DetId
    
  int intIndex=CSCDetId::rawIdMaker(fwbw_id, station_id, ring_id,
             chamber_id, layer_id);

#ifdef LOCAL_DEBUG
  std::cout << "CSCNumberingScheme : ";
  std::cout << " fw/bw " <<  fwbw_id;
  std::cout << " station " <<  station_id;
  std::cout << " ring " <<  ring_id;
  std::cout << " subring " <<  subring_id;
  std::cout << " chamber " <<  chamber_id;
  std::cout << " sector " <<  sector_id;
  std::cout << " layer " <<  layer_id;
  std::cout << std::endl;
#endif

  return intIndex;
}

int CSCNumberingScheme::chamberIndex(int station_id,  
           int ring_id, int subring_id, int sector_id) const {

  int chamber_id=0;

  // chamber label is related to sector_id but we need to
  // adjust to real hardware labelling
  // Tim confirms this works properly according to CMS IN 2000/004 Version 2.5 March 2007.

  if (ring_id == 3) {
    chamber_id=sector_id;
  } else {
    if (subring_id == 1) {
      chamber_id=2*sector_id-1;
    } else { 
      chamber_id=2*sector_id;
    }
  }

  return chamber_id;

}