RPCDetId.h

Go to the documentation of this file.

#ifndef MuonDetId_RPCDetId_h
#define MuonDetId_RPCDetId_h

#include <DataFormats/DetId/interface/DetId.h>
#include <FWCore/Utilities/interface/Exception.h>

#include <iosfwd>

class RPCDetId : public DetId {
public:
  RPCDetId();

  RPCDetId(uint32_t id);
  RPCDetId(DetId id);

  RPCDetId(int region, int ring, int station, int sector, int layer, int subsector, int roll);

  bool operator<(const RPCDetId& r) const {
    if (r.station() == this->station()) {
      if (this->layer() == r.layer()) {
        return this->rawId() < r.rawId();
      } else {
        return (this->layer() < r.layer());
      }
    } else {
      return this->station() < r.station();
    }
  }

  void buildfromDB(int region,
                   int ring,
                   int layer,
                   int sector,
                   const std::string& subsector,
                   const std::string& roll,
                   const std::string& dbname);

  void buildfromTrIndex(int trIndex);

  int region() const { return int((id_ >> RegionStartBit_) & RegionMask_) + minRegionId; }

  int ring() const {
    int ring_ = (id_ >> RingStartBit_) & RingMask_;

    if (ring_ < RingBarrelOffSet) {
      if (this->region() == 0) {
        throw cms::Exception("InvalidDetId") << "RPCDetId ctor:"
                                             << " Ring - Region Inconsistency, "
                                             << " region " << this->region() << " ring " << ring_ << std::endl;
      }

      return int(ring_ + minRingForwardId);

    } else {  // if(ring_ >= RingBarrelOffSet)
      return int(ring_ - RingBarrelOffSet + minRingBarrelId);
    }
  }

  int station() const { return int((id_ >> StationStartBit_) & StationMask_) + minStationId; }

  int sector() const { return int((id_ >> SectorStartBit_) & SectorMask_) + (minSectorId + 1); }

  int layer() const { return int((id_ >> LayerStartBit_) & LayerMask_) + minLayerId; }

  int subsector() const { return int((id_ >> SubSectorStartBit_) & SubSectorMask_) + (minSubSectorId + 1); }

  int roll() const {
    return int((id_ >> RollStartBit_) & RollMask_);  // value 0 is used as wild card
  }

  int trIndex() const { return trind; }

  RPCDetId chamberId() const { return RPCDetId(id_ & chamberIdMask_); }

  static const int minRegionId = -1;
  static const int maxRegionId = 1;

  static const int minRingForwardId = 1;
  static const int maxRingForwardId = 3;
  static const int minRingBarrelId = -2;
  static const int maxRingBarrelId = 2;
  static const int RingBarrelOffSet = 3;

  static const int minStationId = 1;
  static const int maxStationId = 4;

  static const int minSectorId = 0;
  static const int maxSectorId = 12;
  static const int minSectorBarrelId = 1;
  static const int maxSectorBarrelId = 12;
  static const int minSectorForwardId = 1;
  static const int maxSectorForwardId = 6;

  static const int minLayerId = 1;
  static const int maxLayerId = 2;

  static const int minSubSectorId = 0;
  static const int maxSubSectorId = 6;
  static const int minSubSectorBarrelId = 1;
  static const int maxSubSectorBarrelId = 4;
  static const int minSubSectorForwardId = 1;
  static const int maxSubSectorForwardId = 6;

  static const int minRollId = 0;
  static const int maxRollId = 5;  // used to be 4 ... need 5 for upgrade
                                   // if we decide to divide the upgrade RE3/1 and RE4/1 up to eta=2.4
                                   // and we want basically to have 0.10 eta / roll ...

private:
  static const int RegionNumBits_ = 2;
  static const int RegionStartBit_ = 0;
  static const int RegionMask_ = 0X3;

  static const int RingNumBits_ = 3;
  static const int RingStartBit_ = RegionStartBit_ + RegionNumBits_;
  static const unsigned int RingMask_ = 0X7;

  static const int StationNumBits_ = 2;
  static const int StationStartBit_ = RingStartBit_ + RingNumBits_;
  static const unsigned int StationMask_ = 0X3;

  static const int SectorNumBits_ = 4;
  static const int SectorStartBit_ = StationStartBit_ + StationNumBits_;
  static const unsigned int SectorMask_ = 0XF;

  static const int LayerNumBits_ = 1;
  static const int LayerStartBit_ = SectorStartBit_ + SectorNumBits_;
  static const unsigned int LayerMask_ = 0X1;

  static const int SubSectorNumBits_ = 3;
  static const int SubSectorStartBit_ = LayerStartBit_ + LayerNumBits_;
  static const unsigned int SubSectorMask_ = 0X7;

  static const int RollNumBits_ = 3;
  static const int RollStartBit_ = SubSectorStartBit_ + SubSectorNumBits_;
  static const unsigned int RollMask_ = 0X7;

  static const uint32_t chamberIdMask_ = ~(RollMask_ << RollStartBit_);

private:
  void init(int region, int ring, int station, int sector, int layer, int subsector, int roll);

  int trind;
};  // RPCDetId

std::ostream& operator<<(std::ostream& os, const RPCDetId& id);

#endif