RPCDetId.cc

Go to the documentation of this file.

#include <DataFormats/MuonDetId/interface/RPCDetId.h>
#include <DataFormats/MuonDetId/interface/MuonSubdetId.h>

#include <iostream>

RPCDetId::RPCDetId() : DetId(DetId::Muon, MuonSubdetId::RPC), trind(0) {}

RPCDetId::RPCDetId(uint32_t id) : DetId(id), trind(0) {
  //  std::cout<<" constructor of the RPCDetId" <<std::endl;
  if (det() != DetId::Muon || subdetId() != MuonSubdetId::RPC) {
    throw cms::Exception("InvalidDetId") << "RPCDetId ctor:"
                                         << " det: " << det() << " subdet: " << subdetId() << " is not a valid RPC id";
  }
}
RPCDetId::RPCDetId(DetId id) : DetId(id), trind(0) {
  //  std::cout<<" constructor of the RPCDetId" <<std::endl;
  if (det() != DetId::Muon || subdetId() != MuonSubdetId::RPC) {
    throw cms::Exception("InvalidDetId") << "RPCDetId ctor:"
                                         << " det: " << det() << " subdet: " << subdetId() << " is not a valid RPC id";
  }
}

RPCDetId::RPCDetId(int region, int ring, int station, int sector, int layer, int subsector, int roll)
    : DetId(DetId::Muon, MuonSubdetId::RPC), trind(0) {
  this->init(region, ring, station, sector, layer, subsector, roll);
}

void RPCDetId::buildfromDB(int region,
                           int ring,
                           int trlayer,
                           int sector,
                           const std::string& subs,
                           const std::string& roll,
                           const std::string& dbname) {
  bool barrel = (region == 0);
  //STATION
  int station = -1;
  if (barrel) {
    if (trlayer == 1 || trlayer == 2)
      station = 1;
    else if (trlayer == 3 || trlayer == 4)
      station = 2;
    else
      station = trlayer - 2;
  } else {
    station = abs(ring);
  }

  //LAYER
  //int layer = 1;
  //if (barrel && station==1) layer = trlayer;
  //if (barrel && station==2) layer = trlayer-2;

  //SUBSECTOR
  int subsector = 1;

  if (barrel) {
    if (station == 3 && subs == "+")
      subsector = 2;
    if (station == 4 &&
        (sector == 1 || sector == 2 || sector == 3 || sector == 5 || sector == 6 || sector == 7 || sector == 8 ||
         sector == 10 || sector == 12) &&
        (subs == "+")) {
      subsector = 2;
    }

    if (station == 4 && sector == 4) {
      if (subs == "--")
        subsector = 1;
      if (subs == "-")
        subsector = 2;
      if (subs == "+")
        subsector = 3;
      if (subs == "++")
        subsector = 4;
    }
  }

  // ROLL
  int iroll = 0;

  if (roll == "Backward" || roll == "A")
    iroll = 1;
  else if (roll == "Central" || roll == "B")
    iroll = 2;
  else if (roll == "Forward" || roll == "C")
    iroll = 3;
  else if (roll == "D")
    iroll = 4;
  else {
    std::cout << "** RPC: DBSpecToDetUnit, how to assigne roll to: " << roll << " ???" << std::endl;
  }

  int trIndex = 0;
  if (barrel) {
    //cout <<" BARREL: " << endl;
    int eta_id = 6 + ring;
    int plane_id = station;
    if (trlayer == 2)
      plane_id = 5;
    if (trlayer == 4)
      plane_id = 6;
    int sector_id = sector * 3;
    int copy_id = subsector;
    int roll_id = iroll;
    trIndex = (eta_id * 10000 + plane_id * 1000 + sector_id * 10 + copy_id) * 10 + roll_id;
  } else {
    //    cout << "ENDCAP : " << endl;
    int eta_id = trlayer;
    if (ring > 0)
      eta_id = 12 - trlayer;
    int plane_id = abs(ring);
    int sector_id = sector;

    if (region < 0) {
      if (sector_id < 20) {
        sector_id = 19 + 1 - sector_id;
      } else {
        sector_id = 36 + 20 - sector_id;
      }
    }
    sector_id -= 1;

    //
    int copy_id = 1;
    int roll_id = iroll;
    trIndex = (eta_id * 10000 + plane_id * 1000 + sector_id * 10 + copy_id) * 10 + roll_id;
  }
  this->buildfromTrIndex(trIndex);
}

void RPCDetId::buildfromTrIndex(int trIndex) {
  trind = trIndex;
  int eta_id = trIndex / 100000;
  int region = 0;
  int ring = 0;
  if (eta_id <= 3) {
    region = -1;
    ring = eta_id;
  } else if (eta_id >= 9) {
    region = 1;
    ring = 12 - eta_id;
  } else {
    region = 0;
    ring = eta_id - 6;
  }
  trIndex = trIndex % 100000;
  int plane_id = trIndex / 10000;
  int station = 0;
  int layer = 0;
  if (plane_id <= 4) {
    station = plane_id;
    layer = 1;
  } else {
    station = plane_id - 4;
    layer = 2;
  }
  trIndex = trIndex % 10000;
  int sector_id = trIndex / 100;

  // RE+1/1 :: the chamber at x=0 (phi=0) start as CH02 instead of CH01, which is not desired
  // while for other chambers: RE+1/(2,3) and RE+2,3,4/(2,3) the rotation seems to be arbitrary
  // I will leave the code for the existing chambers as it is, but will remove RE+1/1 from selection
  // These lines are programmed very asymmetric between Pos \& Neg endcap:
  // - it affects the whole Negative Endcap
  // - it affects the whole RE+/-1 Station
  // - it affects all RE+(1,2,3,4)/(2,3)
  // ==> why does it act differently on RE-(2,3,4)/1 and RE+(2,3,4)/1 ???
  if (region != 0) {
    if (!(ring == 1 && station > 1 && region == 1)) {
      // skip RE+1/1 (ri=1, st=1, re=-1,+1)
      if (!(ring == 1 && station == 1 && region != 0)) {
        sector_id += 1;
        if (sector_id == 37)
          sector_id = 1;
      }
    }
  }

  if (region == -1) {
    if (sector_id < 20) {
      sector_id = 19 + 1 - sector_id;
    } else {
      sector_id = 36 + 20 - sector_id;
    }
  }
  trIndex = trIndex % 100;
  int copy_id = trIndex / 10;
  int sector = (sector_id - 1) / 3 + 1;
  if (region != 0) {
    sector = (sector + 1) / 2;
  }
  int subsector = 0;
  if (region == 0) {
    subsector = copy_id;
  } else {
    if (ring == 1 && station > 1) {
      // 20 degree chambers
      subsector = ((sector_id + 1) / 2 - 1) % 3 + 1;
    } else {
      // 10 degree chambers
      subsector = (sector_id - 1) % 6 + 1;
    }
    //     std::cout <<" RE"<<station*region<<"/"<<ring<<" sector_id "<<sector_id
    //            << " sector "<<sector <<" sub "<<subsector<<std::endl;
  }

  int roll = trIndex % 10;
  this->init(region, ring, station, sector, layer, subsector, roll);
}

void RPCDetId::init(int region, int ring, int station, int sector, int layer, int subsector, int roll) {
  int minRing = 0;
  int maxRing = RPCDetId::maxRingForwardId;
  if (!region) {
    minRing = RPCDetId::minRingBarrelId;
    maxRing = RPCDetId::maxRingBarrelId;
  }

  if (region < minRegionId || region > maxRegionId || ring < minRing || ring > maxRing || station < minStationId ||
      station > maxStationId || sector < minSectorId || sector > maxSectorId || layer < minLayerId ||
      layer > maxLayerId || subsector < minSubSectorId || subsector > maxSubSectorId || roll < minRollId ||
      roll > maxRollId) {
    throw cms::Exception("InvalidDetId") << "RPCDetId ctor:"
                                         << " Invalid parameters: "
                                         << " region " << region << " ring " << ring << " station " << station
                                         << " sector " << sector << " layer " << layer << " subsector " << subsector
                                         << " roll " << roll << std::endl;
  }

  int regionInBits = region - minRegionId;
  int ringInBits = 0;
  if (region != 0)
    ringInBits = ring - minRingForwardId;
  if (!region)
    ringInBits = ring + RingBarrelOffSet - minRingBarrelId;

  int stationInBits = station - minStationId;
  int sectorInBits = sector - (minSectorId + 1);
  int layerInBits = layer - minLayerId;
  int subSectorInBits = subsector - (minSubSectorId + 1);
  int rollInBits = roll;

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

std::ostream& operator<<(std::ostream& os, const RPCDetId& id) {
  os << " Re " << id.region() << " Ri " << id.ring() << " St " << id.station() << " Se " << id.sector() << " La "
     << id.layer() << " Su " << id.subsector() << " Ro " << id.roll() << " Tr " << id.trIndex() << " ";

  return os;
}