RPCCompDetId.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     src
// Class  :     RPCCompDetId
// 
// Implementation:
//     [Notes on implementation]
//
// Author:      Marcello Maggi
// Created:     Wed Nov  2 12:09:10 CET 2011
#include <iostream>
#include <sstream>
#include <iomanip>
#include "DataFormats/MuonDetId/interface/RPCCompDetId.h"
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h" 

RPCCompDetId::RPCCompDetId():DetId(DetId::Muon, MuonSubdetId::RPC),_dbname(""),_type(0){}

RPCCompDetId::RPCCompDetId(uint32_t id):DetId(id),_dbname(""),_type(0) {
  if (det()!=DetId::Muon || subdetId()!=MuonSubdetId::RPC) {
    throw cms::Exception("InvalidDetId") << "RPCCompDetId ctor:"
                     << " det: " << det()
                     << " subdet: " << subdetId()
                         << " is not a valid RPC id";  
  }
}



RPCCompDetId::RPCCompDetId(DetId id):DetId(id),_dbname(""),_type(0) {
  if (det()!=DetId::Muon || subdetId()!=MuonSubdetId::RPC) {
    throw cms::Exception("InvalidDetId") << "RPCCompDetId ctor:"
                     << " det: " << det()
                     << " subdet: " << subdetId()
                     << " is not a valid RPC id";
  }
}



RPCCompDetId::RPCCompDetId(int region,
               int ring,
               int station,
               int sector,
               int layer,
               int subsector,
               int type):
  DetId(DetId::Muon, MuonSubdetId::RPC),_dbname(""),_type(type)
{
  this->init(region,ring,station,sector,layer,subsector); 
}

RPCCompDetId::RPCCompDetId(const std::string& name, int type):
  DetId(DetId::Muon, MuonSubdetId::RPC),_dbname(name),_type(type)
{
  this->init();
}

bool 
RPCCompDetId::operator < (const RPCCompDetId& r) const{
  return this->dbname()<r.dbname();
}

int 
RPCCompDetId::region() const{
  return int((id_>>RegionStartBit_) & RegionMask_) + allRegionId;
}

int
RPCCompDetId::ring() const{
  return int((id_>>RingStartBit_) & RingMask_) + allRingId;
}


int
RPCCompDetId::wheel() const{
  int w=allRingId;
  if (this->region()==0)
    w=this->ring();
  return w;
}

int
RPCCompDetId::station() const{
  return int((id_>>StationStartBit_) & StationMask_) + allStationId;
}

int
RPCCompDetId::disk() const{
  int d=allStationId;
  if (this->region()!=0)
    d=this->station();
  return d;
}


int
RPCCompDetId::sector() const{
  return int((id_>>SectorStartBit_) & SectorMask_) + allSectorId;
}


int
RPCCompDetId::layer() const{
  return int((id_>>LayerStartBit_) & LayerMask_) + allLayerId;
}


int
RPCCompDetId::subsector() const{
  return int((id_>>SubSectorStartBit_) & SubSectorMask_) + allSubSectorId;
}


int
RPCCompDetId::type() const{
  return _type;
}

std::string
RPCCompDetId::dbname() const{
  std::string a=_dbname;
  if (a.empty()){
    if(this->type() == 0){
      a=this->gasDBname();
    }
  }
  return a;
}

void
RPCCompDetId::init(int region,
           int ring,
           int station,
           int sector,
           int layer,
           int subsector)
{
  int maxRing=maxRingForwardId;
  if (!region)
    {
      maxRing=maxRingBarrelId;
    }
  
  if ( region     < allRegionId    || region    > maxRegionId ||
       ring       < allRingId      || ring      > maxRing ||
       station    < allStationId   || station   > maxStationId ||
       sector     < allSectorId    || sector    > maxSectorId ||
       layer      < allLayerId     || layer     > maxLayerId ||
       subsector  < allSubSectorId || subsector > maxSubSectorId ){
    
    throw cms::Exception("InvalidDetId") << "RPCDetId ctor:"
                                         << " Invalid parameters: "
                                         << " region "<<region
                                         << " ring "<<ring
                                         << " station "<<station
                                         << " sector "<<sector
                                         << " layer "<<layer
                                         << " subsector "<<subsector
                                         << std::endl;
  }
  int regionInBits    = region  -  allRegionId;
  int ringInBits      = ring     - allRingId;
  int stationInBits   = station  - allStationId;
  int sectorInBits    = sector   - allSectorId;
  int layerInBits     = layer    - allLayerId;
  int subSectorInBits = subsector- allSubSectorId;

  id_ |= ( regionInBits    & RegionMask_)    << RegionStartBit_    |
         ( ringInBits      & RingMask_)      << RingStartBit_      |
         ( stationInBits   & StationMask_)   << StationStartBit_   |
         ( sectorInBits    & SectorMask_)    << SectorStartBit_    |
         ( layerInBits     & LayerMask_)     << LayerStartBit_     |
         ( subSectorInBits & SubSectorMask_) << SubSectorStartBit_ ;
}

void 
RPCCompDetId::init()
{
  if (this->type()==0){
    this->initGas();
  }
}

void
RPCCompDetId::initGas()
{
  std::string buf(this->dbname()); 
  // check if the name contains the dcs namespace
  if (buf.find(':')!=buf.npos){
    buf = buf.substr(buf.find(':')+1,buf.npos);
  }
  _dbname=buf;
  // Check if endcap o barrel
  int region=0;
  if(buf.substr(0,1)=="W"){
    region=0;
  }else if(buf.substr(0,2)=="EP"){
    region=1;
  }else if(buf.substr(0,2)=="EM"){
    region=-1;
  }else{
    throw cms::Exception("InvalidDBName")<<" RPCCompDetId: "<<this->dbname()
                     <<" is not a valid DB Name for RPCCompDetId"
                                         << " det: " << det()
                     << " subdet: " << subdetId();
  }
  int ring=allRingId;
  int station = allStationId;
  int sector=allSectorId;
  int layer=allLayerId;
  int subsector=allSubSectorId;
    //Barrel
  if (region==0) {
    // Extract the Wheel (named ring)
    {
      std::stringstream os;
      os<<buf.substr(2,1);
      os>>ring;
      if (buf.substr(1,1)=="M"){
    ring *= -1;
      }
    }
    //Extract the station
    {
      std::stringstream os;
      os<<buf.substr(buf.find("RB")+2,1);
      os>>station;
    }
    //Extract the sector
    {
      std::stringstream os;
      os<<buf.substr(buf.find("S")+1,2);
      os>>sector;
    }
    //Extract subsector of sectors 4 and 10
    {
      if (buf.find("4L")!=buf.npos)
    subsector=1;
      if (buf.find("4R")!=buf.npos)
    subsector=2;
    }
  }else{
  // Extract the Ring 
    {
      std::stringstream os;
      os<<buf.substr(buf.find("_R0")+3,1);
      os>>ring;
    }
  //Extract the disk (named station)
    {
      std::stringstream os;
      os<<buf.substr(2,1);
      os>>station;
    }
    //Extract the sector or chamber
    {
      std::stringstream os;
      os<<buf.substr(buf.find("_C")+2,2);
      os>>sector;
    }
    //Extract layer
    {
      if (buf.find("UP")!=buf.npos)
    layer=1;
      if (buf.find("DW")!=buf.npos)
    layer=2;
    }
  }
  this->init(region,ring,station,sector,layer,subsector); 
}

std::string
RPCCompDetId::gasDBname() const{
  std::stringstream os;
  if(this->region()==0){
    // Barrel
    std::string wsign="P";
    if (this->wheel()<0)wsign= "M";
    std::string lr="";
    if (this->subsector()==1) lr="L";
    if (this->subsector()==2) lr="R";
    os<<"W"<<wsign<<abs(this->wheel())<<"_S"<<std::setw(2)<<std::setfill('0')<<this->sector()<<"_RB"<<this->station()<<lr;
  } else {
    // Endcap
    std::string esign="P";
    if (this->region()<0)
      esign="M";

    os<<"E"<<esign<<this->disk();


    if (this->disk()==1){
      os<<"_R"<<std::setw(2)<<std::setfill('0')<<this->ring()
    <<"_C"<<std::setw(2)<<std::setfill('0')<<this->sector()
    <<"_C"<<std::setw(2)<<std::setfill('0')<<this->sector()+5;
    }else{
      os<<"_R"<<std::setw(2)<<std::setfill('0')<<this->ring()
    <<"_R"<<std::setw(2)<<std::setfill('0')<<this->ring()+1
    <<"_C"<<std::setw(2)<<std::setfill('0')<<this->sector()
    <<"_C"<<std::setw(2)<<std::setfill('0')<<this->sector()+2;
     
    }
    std::string lay="";
    if(this->layer()==1)
      lay="UP";
    else if (this->layer()==2)
      lay="DW";

    os<<"_"<<lay;
      
 }
  return os.str();
}

std::ostream& operator<<( std::ostream& os, const RPCCompDetId& id ){

  os <<id.dbname();

  return os;
}