CSCIndexerBase.cc

Go to the documentation of this file.

#include "CalibMuon/CSCCalibration/interface/CSCIndexerBase.h"
 
CSCIndexerBase::CSCIndexerBase()
    : chamberLabel_(271)  // # of physical chambers per endcap + 1. Includes ME42.
{
  // Fill the member vector which permits decoding of the linear chamber index.
  // Beware that the ME42 indices 235-270 within this vector do NOT correspond
  // to their 'real' linear indices (which are 469-504 for +z)
  IndexType count = 0;
  chamberLabel_[count] = 0;
 
  for (IndexType is = 1; is <= 4; ++is) {
    IndexType irmax = ringsInStation(is);
    for (IndexType ir = 1; ir <= irmax; ++ir) {
      IndexType icmax = chambersInRingOfStation(is, ir);
      for (IndexType ic = 1; ic <= icmax; ++ic) {
        chamberLabel_[++count] = is * 1000 + ir * 100 + ic;
      }
    }
  }
}
 
CSCIndexerBase::~CSCIndexerBase() {}
 
CSCIndexerBase::IndexType CSCIndexerBase::chamberLabelFromChamberIndex(IndexType ici) const {
  // This is just for cross-checking
 
  // Expected range of input range argument is 1-540.
  // 1-468 for CSCs installed at 2008 start-up. 469-540 for ME42.
 
  if (ici > 468) {
    // ME42
    ici -= 234;     // now in range 235-306
    if (ici > 270)  // -z
    {
      ici -= 36;  // now in range 235-270
    }
  } else  // in range 1-468
  {
    if (ici > 234)  // -z
    {
      ici -= 234;  // now in range 1-234
    }
  }
  return chamberLabel_[ici];
}
 
CSCIndexerBase::IndexType CSCIndexerBase::hvSegmentIndex(IndexType is, IndexType ir, IndexType iwire) const {
  IndexType hvSegment = 1;  // There is only one HV segment in ME1/1
 
  if (is > 2 && ir == 1)  // HV segments are the same in ME3/1 and ME4/1
  {
    if (iwire >= 33 && iwire <= 64) {
      hvSegment = 2;
    } else if (iwire >= 65 && iwire <= 96) {
      hvSegment = 3;
    }
  } else if (is > 1 && ir == 2)  // HV segments are the same in ME2/2, ME3/2, and ME4/2
  {
    if (iwire >= 17 && iwire <= 28) {
      hvSegment = 2;
    } else if (iwire >= 29 && iwire <= 40) {
      hvSegment = 3;
    } else if (iwire >= 41 && iwire <= 52) {
      hvSegment = 4;
    } else if (iwire >= 53 && iwire <= 64) {
      hvSegment = 5;
    }
  } else if (is == 1 && ir == 2) {
    if (iwire >= 25 && iwire <= 48) {
      hvSegment = 2;
    } else if (iwire >= 49 && iwire <= 64) {
      hvSegment = 3;
    }
  } else if (is == 1 && ir == 3) {
    if (iwire >= 13 && iwire <= 22) {
      hvSegment = 2;
    } else if (iwire >= 23 && iwire <= 32) {
      hvSegment = 3;
    }
  } else if (is == 2 && ir == 1) {
    if (iwire >= 45 && iwire <= 80) {
      hvSegment = 2;
    } else if (iwire >= 81 && iwire <= 112) {
      hvSegment = 3;
    }
  }
  return hvSegment;
}
 
CSCDetId CSCIndexerBase::detIdFromChamberLabel(IndexType ie, IndexType label) const {
  IndexType is = label / 1000;
  label -= is * 1000;
  IndexType ir = label / 100;
  label -= ir * 100;
  IndexType ic = label;
 
  return CSCDetId(ie, is, ir, ic);
}
 
CSCDetId CSCIndexerBase::detIdFromChamberIndex(IndexType ici) const {
  // Expected range of input range argument is 1-540.
  // 1-468 for CSCs installed at 2008 start-up. 469-540 for ME42.
 
  IndexType ie = 1;
  if (ici > 468) {
    // ME42
    ici -= 234;     // now in range 235-306
    if (ici > 270)  // -z
    {
      ie = 2;
      ici -= 36;  // now in range 235-270
    }
  } else  // in range 1-468
  {
    if (ici > 234)  // -z
    {
      ie = 2;
      ici -= 234;  // now in range 1-234
    }
  }
 
  IndexType label = chamberLabel_[ici];
  return detIdFromChamberLabel(ie, label);
}
 
CSCDetId CSCIndexerBase::detIdFromLayerIndex(IndexType ili) const {
  IndexType il = (ili - 1) % 6 + 1;
  IndexType ici = (ili - 1) / 6 + 1;
  CSCDetId id = detIdFromChamberIndex(ici);
 
  return CSCDetId(id.endcap(), id.station(), id.ring(), id.chamber(), il);
}