CSCDetId.cc

Go to the documentation of this file.

#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include <FWCore/Utilities/interface/Exception.h>
#include <iostream>

CSCDetId::CSCDetId():DetId(DetId::Muon, MuonSubdetId::CSC){}


CSCDetId::CSCDetId(uint32_t id):DetId(id) {
  if (det()!=DetId::Muon || subdetId()!=MuonSubdetId::CSC) {
    throw cms::Exception("InvalidDetId") << "CSCDetId ctor:"
                     << " det: " << det()
                     << " subdet: " << subdetId()
                     << " is not a valid CSC id";  
  }
}

CSCDetId::CSCDetId(DetId id):DetId(id) {
  if (det()!=DetId::Muon || subdetId()!=MuonSubdetId::CSC) {
    throw cms::Exception("InvalidDetId") << "CSCDetId ctor:"
                     << " det: " << det()
                     << " subdet: " << subdetId()
                     << " is not a valid CSC id";  
  }
}

CSCDetId::CSCDetId( int iendcap, int istation, int iring, int ichamber, 
            int ilayer ) : 
  DetId(DetId::Muon, MuonSubdetId::CSC) 
{    
  if (iendcap  < 0 || iendcap  > MAX_ENDCAP ||
      istation < 0 || istation > MAX_STATION ||
      iring    < 0 || iring    > MAX_RING ||
      ichamber < 0 || ichamber > MAX_CHAMBER ||
      ilayer   < 0 || ilayer   > MAX_LAYER ) {
    throw cms::Exception("InvalidDetId") << "CSCDetId ctor:" 
                     << " Invalid parameters: " 
                     << " E:"<< iendcap
                     << " S:"<< istation
                     << " R:"<< iring
                     << " C:"<< ichamber
                     << " L:"<< ilayer   
                     << std::endl;
  }
  id_ |= init(iendcap, istation, iring, ichamber, ilayer);
}

int CSCDetId::triggerSector() const
{
  // UPDATED TO OCT 2005 - LGRAY Feb 2006

  int result;
  int ring    = this->ring();
  int station = this->station();
  int chamber = this->chamber();

    if(station > 1 && ring > 1 ) {
      result = ((static_cast<unsigned>(chamber-3) & 0x7f) / 6) + 1; // ch 3-8->1, 9-14->2, ... 1,2 -> 6
    }
    else {
      result =  (station != 1) ? ((static_cast<unsigned>(chamber-2) & 0x1f) / 3) + 1 : // ch 2-4-> 1, 5-7->2, ...
                             ((static_cast<unsigned>(chamber-3) & 0x7f) / 6) + 1;
    }

  return (result <= 6) ? result : 6; // max sector is 6, some calculations give a value greater than six but this is expected.
}

int CSCDetId::triggerCscId() const
{
  // UPDATED TO OCT 2005 - LGRAY Feb 2006

  int result;
  int ring    = this->ring();
  int station = this->station();
  int chamber = this->chamber();

  if( station == 1 ) {
    result = (chamber) % 3 + 1; // 1,2,3
    switch (ring) {
    case 1:
      break;
    case 2:
      result += 3; // 4,5,6
      break;
    case 3:
      result += 6; // 7,8,9
      break;
    }
  }
  else {
    if( ring == 1 ) {
      result = (chamber+1) % 3 + 1; // 1,2,3
    }
    else {
      result = (chamber+3) % 6 + 4; // 4,5,6,7,8,9
    }
  }
  return result;
}

unsigned short CSCDetId::iChamberType( unsigned short istation, unsigned short iring ) {
  int i = 2 * istation + iring; // i=2S+R ok for S=2, 3, 4
  if ( istation == 1 ) {
    --i;                       // ring 1R -> i=1+R (2S+R-1=1+R for S=1)
    if ( i > 4 ) i = 1;        // But ring 1A (R=4) -> i=1
  }   
  return i;
}


std::ostream& operator<<( std::ostream& os, const CSCDetId& id )
{
  // Note that there is no endl to end the output

   os << " E:" << id.endcap()
      << " S:" << id.station()
      << " R:" << id.ring()
      << " C:" << id.chamber()
      << " L:" << id.layer();
   return os;
}