CSCIndexer Class Reference

#include <CSCIndexer.h>

Public Types
typedef uint16_t	IndexType
typedef uint32_t	LongIndexType

Public Member Functions
IndexType	chamberIndex (const CSCDetId &id) const
IndexType	chamberIndex (IndexType ie, IndexType is, IndexType ir, IndexType ic) const
IndexType	chamberLabelFromChamberIndex (IndexType) const
IndexType	chipIndex (IndexType ie, IndexType is, IndexType ir, IndexType ic, IndexType il, IndexType ichip) const
IndexType	chipIndex (const CSCDetId &id, IndexType ichip) const
IndexType	chipIndex (IndexType istrip) const
IndexType	chipsPerLayer (IndexType is, IndexType ir) const
IndexType	chipStart (IndexType ie, IndexType is, IndexType ir) const
	CSCIndexer ()
int	dbIndex (const CSCDetId &id, int &channel)
CSCDetId	detIdFromChamberIndex (IndexType ici) const
CSCDetId	detIdFromChamberIndex_OLD (IndexType ici) const
CSCDetId	detIdFromChamberLabel (IndexType ie, IndexType icl) const
std::pair< CSCDetId, IndexType >	detIdFromChipIndex (IndexType ichi) const
CSCDetId	detIdFromLayerIndex (IndexType ili) const
std::pair< CSCDetId, IndexType >	detIdFromStripChannelIndex (LongIndexType ichi) const
IndexType	gasGainIndex (const CSCDetId &id, IndexType istrip, IndexType iwire) const
IndexType	gasGainIndex (IndexType ie, IndexType is, IndexType ir, IndexType ic, IndexType il, IndexType istrip, IndexType iwire) const
IndexType	hvSegmentIndex (IndexType is, IndexType ir, IndexType iwire) const
IndexType	hvSegmentsPerLayer (IndexType is, IndexType ir) const
IndexType	layerIndex (const CSCDetId &id) const
IndexType	layerIndex (IndexType ie, IndexType is, IndexType ir, IndexType ic, IndexType il) const
IndexType	sectorsPerLayer (IndexType is, IndexType ir) const
IndexType	sectorStart (IndexType ie, IndexType is, IndexType ir) const
IndexType	startChamberIndexInEndcap (IndexType ie, IndexType is, IndexType ir) const
LongIndexType	stripChannelIndex (IndexType ie, IndexType is, IndexType ir, IndexType ic, IndexType il, IndexType istrip) const
LongIndexType	stripChannelIndex (const CSCDetId &id, IndexType istrip) const
IndexType	stripChannelsPerLayer (IndexType is, IndexType ir) const
LongIndexType	stripChannelStart (IndexType ie, IndexType is, IndexType ir) const
	~CSCIndexer ()

Static Public Member Functions
static IndexType	chambersInRingOfStation (IndexType is, IndexType ir)
static IndexType	ringsInStation (IndexType is)

Static Private Member Functions
static std::vector< IndexType > const &	chamberLabel ()
static std::vector< IndexType >	fillChamberLabel ()

Detailed Description

Creates a linear index for various sublevels of the endcap muon CSC system.

It supplies a linear index for:

1. Each chamber of the CSC system: range 1-468 (CSC system as installed 2008) 469-540 for ME42
   
2. Each layer of the CSC system: range 1-2808 (CSCs 2008) 2809-3240 for ME42
   
3. Each strip channel of the CSC system: range 1-217728 (CSCs 2008) 217729-252288 for ME42
   
4. Each Buckeye chip (1 for each layer in each CFEB)of the CSC system: range 1-13608 (CSCs 2008) 13609-15768 for ME42
   
5. Each Gas Gain Sector (1 for each [CFEB*HV segment] combination in each layer): range 1-45144 (CSCs 2008) 45145-55944 for ME42
   

The chamber and layer may be specified by CSCDetId or labels for endcap, station, ring, chamber, layer. The strip channel is a value 1-80 (or 64: ME13 chambers have only 64 channels.) The chip number is a value 1-30 (or 24: ME13 chambers have only 24 chips.)

The main user interface is the set of functions
chamberIndex(.)
layerIndex(.)
stripChannelIndex(.)
chipIndex(.)
But the other functions are public since they may be useful in contexts other than for Conditions Data for which the above functions are intended.

Warning

This class is hard-wired for the CSC system at start-up of CMS in 2008. with rings ME11, ME12, ME13, ME21, ME22, ME31, ME32, ME41 totalling 234 chambers per endcap. But ME42 is appended (to permit simulation studies), so the chamber order is
+z ME11, ME12, ME13, ME21, ME22, ME31, ME32, ME41,
-z ME11, ME12, ME13, ME21, ME22, ME31, ME32, ME41,
+z ME42, -z ME42

This uses magic numbers galore!!

EVERY LABEL COUNTS FROM ONE NOT ZERO.

Definition at line 44 of file CSCIndexer.h.

Member Typedef Documentation

IndexType

typedef uint16_t CSCIndexer::IndexType

Definition at line 48 of file CSCIndexer.h.

LongIndexType

typedef uint32_t CSCIndexer::LongIndexType

Definition at line 49 of file CSCIndexer.h.

Constructor & Destructor Documentation

CSCIndexer()

CSCIndexer::CSCIndexer ( )

inline

Definition at line 51 of file CSCIndexer.h.

{};

~CSCIndexer()

CSCIndexer::~CSCIndexer ( )

inline

Definition at line 52 of file CSCIndexer.h.

{};

Member Function Documentation

chamberIndex() [1/2]

IndexType CSCIndexer::chamberIndex ( const CSCDetId &  id ) const

inline

Linear index to label each CSC in CSC system. Argument is the CSCDetId of some CSCChamber.

Output is 1-468 (CSCs 2008) 469-540 (ME42)

WARNING: Do not input ME1a values (i.e. ring '4'): does not consider ME1a and ME1b to be separate, No sanity checking on input value: if you supply an ME1a CSCDetId then you'll get nonsense.

Definition at line 64 of file CSCIndexer.h.

References relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, relativeConstraints::ring, and relativeConstraints::station.

Referenced by layerIndex(), and CSCChamberTimeCorrectionsValues::prefill().

                                                   {
    return chamberIndex(id.endcap(), id.station(), id.ring(), id.chamber());
  }

chamberIndex() [2/2]

IndexType CSCIndexer::chamberIndex ( IndexType  ie, IndexType  is, IndexType  ir, IndexType  ic  ) const

inline

Linear index for chamber 'ic' in ring 'ir' of station 'is' in endcap 'ie', in range 1-468 (CSCs 2008) or 469-540 (ME42)

Definition at line 96 of file CSCIndexer.h.

References startChamberIndexInEndcap().

                                                                                       {
    return startChamberIndexInEndcap(ie, is, ir) + ic - 1;  // -1 so start index _is_ ic=1
  }

chamberLabel()

const std::vector< CSCIndexer::IndexType > & CSCIndexer::chamberLabel ( )

staticprivate

Definition at line 27 of file CSCIndexer.cc.

References fillChamberLabel().

Referenced by chamberLabelFromChamberIndex(), and detIdFromChamberIndex().

                                                               {
  static const auto s_chamberLabel = fillChamberLabel();
  return s_chamberLabel;
}

chamberLabelFromChamberIndex()

CSCIndexer::IndexType CSCIndexer::chamberLabelFromChamberIndex

(

IndexType

ici

)

const

Definition at line 91 of file CSCIndexer.cc.

References chamberLabel().

                                                                                {
  // 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];
}

chambersInRingOfStation()

static IndexType CSCIndexer::chambersInRingOfStation ( IndexType  is, IndexType  ir  )

inlinestatic

How many chambers are there in ring ir of station is?

Works for ME1a (ring 4 of ME1) too.

Definition at line 124 of file CSCIndexer.h.

Referenced by fillChamberLabel().

                                                                       {
    IndexType nc = 36;  // most rings have 36 chambers
    if (is > 1 && ir < 2)
      nc = 18;  // but 21, 31, 41 have 18
    return nc;
  }

chipIndex() [1/3]

IndexType CSCIndexer::chipIndex ( IndexType  ie, IndexType  is, IndexType  ir, IndexType  ic, IndexType  il, IndexType  ichip  ) const

inline

Linear index for Buckeye chip 'ichip' in layer 'il' of chamber 'ic' of ring 'ir' in station 'is' of endcap 'ie'.

Output is 1-13608 (CSCs 2008) or 13609-15768 (ME42).

WARNING: Use at your own risk! You must input labels within hardware ranges. No trapping on out-of-range values!

Definition at line 234 of file CSCIndexer.h.

References chipsPerLayer(), and chipStart().

Referenced by chipIndex(), gasGainIndex(), and CSCChipSpeedCorrectionDBConditions::prefillDBChipSpeedCorrection().

                                                                                                                   {
    //printf("ME%d/%d/%d/%d layer %d  chip %d chipindex %d\n",ie,is,ir,ic,il,ichip,chipStart(ie,is,ir)+( (ic-1)*6 + il - 1 )*chipsPerLayer(is,ir) + (ichip-1));
    return chipStart(ie, is, ir) + ((ic - 1) * 6 + il - 1) * chipsPerLayer(is, ir) + (ichip - 1);
  }

chipIndex() [2/3]

IndexType CSCIndexer::chipIndex ( const CSCDetId &  id, IndexType  ichip  ) const

inline

Linear index for Buckeye chip 'ichip' in layer labelled by CSCDetId 'id'.

Output is 1-13608 (CSCs 2008) or 13609-15768 (ME42).

WARNING: Use at your own risk! The supplied CSCDetId must be a layer id. No trapping on out-of-range values!

Definition at line 248 of file CSCIndexer.h.

References relativeConstraints::chamber, chipIndex(), makeMuonMisalignmentScenario::endcap, phase1PixelTopology::layer, relativeConstraints::ring, and relativeConstraints::station.

                                                                 {
    return chipIndex(id.endcap(), id.station(), id.ring(), id.chamber(), id.layer(), ichip);
  }

chipIndex() [3/3]

IndexType CSCIndexer::chipIndex ( IndexType  istrip ) const

inline

Linear index for Buckeye chip processing strip 'istrip'.

Output is 1-5.

WARNING: Use at your own risk! The supplied CSCDetId must be a strip id 1-80 ME1/1a strips must be 65-80 No trapping on out-of-range values!

Definition at line 262 of file CSCIndexer.h.

{ return (istrip - 1) / 16 + 1; }

chipsPerLayer()

IndexType CSCIndexer::chipsPerLayer ( IndexType  is, IndexType  ir  ) const

inline

Number of Buckeye chips per layer in a chamber in ring ir of station is.

Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 1 of the 5 total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.
Considers ME42 as standard 5 chip per layer chambers.

Definition at line 200 of file CSCIndexer.h.

Referenced by chipIndex(), gasGainIndex(), and sectorsPerLayer().

                                                            {
    const IndexType nCinL[12] = {5, 5, 4, 5, 5, 0, 5, 5, 0, 5, 5, 0};
    return nCinL[(is - 1) * 3 + ir - 1];
  }

chipStart()

IndexType CSCIndexer::chipStart ( IndexType  ie, IndexType  is, IndexType  ir  ) const

inline

Linear index for 1st Buckey chip in ring 'ir' of station 'is' in endcap 'ie'.

Endcap label range 1-2, Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 1 of the 5 chips total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.

Definition at line 213 of file CSCIndexer.h.

Referenced by chipIndex().

                                                                      {
    // These are in the ranges 1-13608 (CSCs 2008) and 13609-15768 (ME42).
    // There are 1-6804 chips per endcap (CSCs 2008) and 1080 channels per endcap (ME42).
    // Start of -z channels (CSCs 2008) is 6804 + 1 = 6805
    // Start of +z (ME42) is 13608 + 1 = 13609
    // Start of -z (ME42) is 13608 + 1 + 1080 = 14689
    const IndexType nStart[24] = {1,    1081, 2161, 3025, 3565,  0, 4645,  5185,  0, 6265,  13609, 0,
                                  6805, 7885, 8965, 9829, 10369, 0, 11449, 11989, 0, 13069, 14689, 0};

    return nStart[(ie - 1) * 12 + (is - 1) * 3 + ir - 1];
  }

dbIndex()

int CSCIndexer::dbIndex	(	const CSCDetId &	id,
		int &	channel
	)

Build index used internally in online CSC conditions databases (the 'Igor Index')

This is the decimal integer ie*100000 + is*10000 + ir*1000 + ic*10 + il
(ie=1-2, is=1-4, ir=1-4, ic=1-36, il=1-6)
Channels 1-16 in ME1A (is=1, ir=4) are reset to channels 65-80 of ME11.

Definition at line 235 of file CSCIndexer.cc.

Referenced by CSCNoiseMatrix::item(), CSCPedestals::item(), CSCGains::item(), and CSCcrosstalk::item().

                                                        {
  int ec = id.endcap();
  int st = id.station();
  int rg = id.ring();
  int ch = id.chamber();
  int la = id.layer();

  // The channels of ME1A are channels 65-80 of ME11
  if (st == 1 && rg == 4) {
    rg = 1;
    if (channel <= 16)
      channel += 64;  // no trapping for any bizarreness
  }
  return ec * 100000 + st * 10000 + rg * 1000 + ch * 10 + la;
}

detIdFromChamberIndex()

CSCDetId CSCIndexer::detIdFromChamberIndex

(

IndexType

ici

)

const

Definition at line 69 of file CSCIndexer.cc.

References chamberLabel(), detIdFromChamberLabel(), and label.

Referenced by detIdFromLayerIndex(), and CSCChipSpeedCorrectionDBConditions::prefillDBChipSpeedCorrection().

                                                              {
  // 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);
}

detIdFromChamberIndex_OLD()

CSCDetId CSCIndexer::detIdFromChamberIndex_OLD

(

IndexType

ici

)

const

Definition at line 32 of file CSCIndexer.cc.

References mps_fire::i, relativeConstraints::ring, startChamberIndexInEndcap(), and relativeConstraints::station.

                                                                  {
  // Will not work as is for ME42
  // ============================

  IndexType ie = 1;
  if (ici > 234) {
    ie = 2;
    ici -= 234;
  }
  // Now ici is in range 1-234 (assuming valid input in range 1-468)

  // MEij pairs...
  const IndexType station[] = {0, 1, 1, 1, 2, 2, 3, 3, 4};
  const IndexType ring[] = {0, 1, 2, 3, 1, 2, 1, 2, 1};

  // MEij preceding a given MEij matching linear index above
  const IndexType prevs[] = {0, 0, 1, 1, 1, 2, 2, 3, 3};
  const IndexType prevr[] = {0, 0, 1, 2, 3, 1, 2, 1, 2};

  IndexType is = 4;
  IndexType ir = 1;
  for (IndexType i = 2; i <= 8; ++i) {
    IndexType js = station[i];
    IndexType jr = ring[i];
    // if it's before start of MEjs/jr then it's in the previous MEis/ir
    if (ici < startChamberIndexInEndcap(ie, js, jr)) {
      is = prevs[i];
      ir = prevr[i];
      break;
    }
    // otherwise it's in ME41
  }
  IndexType ic = ici - startChamberIndexInEndcap(ie, is, ir) + 1;

  return CSCDetId(ie, is, ir, ic);
}

detIdFromChamberLabel()

CSCDetId CSCIndexer::detIdFromChamberLabel	(	IndexType	ie,
		IndexType	icl
	)		const

Definition at line 111 of file CSCIndexer.cc.

References label.

Referenced by detIdFromChamberIndex().

                                                                              {
  IndexType is = label / 1000;
  label -= is * 1000;
  IndexType ir = label / 100;
  label -= ir * 100;
  IndexType ic = label;

  return CSCDetId(ie, is, ir, ic);
}

detIdFromChipIndex()

std::pair< CSCDetId, CSCIndexer::IndexType > CSCIndexer::detIdFromChipIndex

(

IndexType

ichi

)

const

Definition at line 182 of file CSCIndexer.cc.

References detIdFromLayerIndex().

                                                                                         {
  const LongIndexType lastnonme42 = 13608;      // chips in 2008 installed chambers
  const LongIndexType lastplusznonme42 = 6804;  // = 13608/2
  const LongIndexType firstme13 = 2161;         // First channel of ME13
  const LongIndexType lastme13 = 3024;          // Last channel of ME13

  const IndexType lastnonme42layer = 2808;
  const IndexType lastplusznonme42layer = 1404;  // = 2808/2
  const IndexType firstme13layer = 433;          // = 72*6 + 1 (ME13 chambers are 72-108 in range 1-234)
  const IndexType lastme13layer = 648;           // = 108*6

  // All chambers but ME13 have 5 chips/layer
  IndexType nchipPerLayer = 5;

  // Set endcap to +z. This should work for ME42 channels too, since we don't need to calculate its endcap explicitly.
  IndexType ie = 1;

  LongIndexType istart = 0;
  IndexType layerOffset = 0;

  if (ici <= lastnonme42) {
    // Chambers as of 2008 Installation

    if (ici > lastplusznonme42) {
      ie = 2;
      ici -= lastplusznonme42;
    }

    if (ici > lastme13) {  // after ME13
      istart = lastme13;
      layerOffset = lastme13layer;
    } else if (ici >= firstme13) {  // ME13
      istart = firstme13 - 1;
      layerOffset = firstme13layer - 1;
      nchipPerLayer = 4;
    }
  } else {
    // ME42 chambers

    istart = lastnonme42;
    layerOffset = lastnonme42layer;
  }

  ici -= istart;  // remove earlier group(s)
  IndexType ichip = (ici - 1) % nchipPerLayer + 1;
  IndexType ili = (ici - 1) / nchipPerLayer + 1;
  ili += layerOffset;  // add appropriate offset for earlier group(s)
  if (ie != 1)
    ili += lastplusznonme42layer;  // add offset to -z endcap; ME42 doesn't need this.

  return std::pair<CSCDetId, IndexType>(detIdFromLayerIndex(ili), ichip);
}

detIdFromLayerIndex()

CSCDetId CSCIndexer::detIdFromLayerIndex

(

IndexType

ili

)

const

Decode CSCDetId from various indexes and labels

Definition at line 121 of file CSCIndexer.cc.

References relativeConstraints::chamber, detIdFromChamberIndex(), makeMuonMisalignmentScenario::endcap, relativeConstraints::ring, and relativeConstraints::station.

Referenced by detIdFromChipIndex(), and detIdFromStripChannelIndex().

                                                            {
  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);
}

detIdFromStripChannelIndex()

std::pair< CSCDetId, CSCIndexer::IndexType > CSCIndexer::detIdFromStripChannelIndex

(

LongIndexType

ichi

)

const

Definition at line 129 of file CSCIndexer.cc.

References detIdFromLayerIndex().

                                                                                                     {
  const LongIndexType lastnonme42 = 217728;       // channels in 2008 installed chambers
  const LongIndexType lastplusznonme42 = 108864;  // = 217728/2
  const LongIndexType firstme13 = 34561;          // First channel of ME13
  const LongIndexType lastme13 = 48384;           // Last channel of ME13

  const IndexType lastnonme42layer = 2808;
  const IndexType lastplusznonme42layer = 1404;  // = 2808/2
  const IndexType firstme13layer = 433;          // = 72*6 + 1 (ME13 chambers are 72-108 in range 1-234)
  const IndexType lastme13layer = 648;           // = 108*6

  // All chambers but ME13 have 80 channels
  IndexType nchan = 80;

  // Set endcap to +z. This should work for ME42 channels too, since we don't need to calculate its endcap explicitly.
  IndexType ie = 1;

  LongIndexType istart = 0;
  IndexType layerOffset = 0;

  if (isi <= lastnonme42) {
    // Chambers as of 2008 Installation

    if (isi > lastplusznonme42) {
      ie = 2;
      isi -= lastplusznonme42;
    }

    if (isi > lastme13) {  // after ME13
      istart = lastme13;
      layerOffset = lastme13layer;
    } else if (isi >= firstme13) {  // ME13
      istart = firstme13 - 1;
      layerOffset = firstme13layer - 1;
      nchan = 64;
    }
  } else {
    // ME42 chambers

    istart = lastnonme42;
    layerOffset = lastnonme42layer;
  }

  isi -= istart;  // remove earlier group(s)
  IndexType ichan = (isi - 1) % nchan + 1;
  IndexType ili = (isi - 1) / nchan + 1;
  ili += layerOffset;  // add appropriate offset for earlier group(s)
  if (ie != 1)
    ili += lastplusznonme42layer;  // add offset to -z endcap; ME42 doesn't need this.

  return std::pair<CSCDetId, IndexType>(detIdFromLayerIndex(ili), ichan);
}

fillChamberLabel()

std::vector< CSCIndexer::IndexType > CSCIndexer::fillChamberLabel ( )

staticprivate

Definition at line 4 of file CSCIndexer.cc.

References chambersInRingOfStation(), submitPVResolutionJobs::count, and ringsInStation().

Referenced by chamberLabel().

                                                            {
  // Fill the member vector which permits decoding of the linear chamber index
  // Logically const since initializes cache only,
  // Beware that the ME42 indices 235-270 within this vector do NOT correspond to
  // their 'real' linear indices (which are 469-504 for +z)
  std::vector<IndexType> tChamberLabel;

  tChamberLabel.resize(271);  // one more than #chambers per endcap. Includes ME42.
  IndexType count = 0;
  tChamberLabel[count] = 0;

  for (IndexType is = 1; is != 5; ++is) {
    IndexType irmax = ringsInStation(is);
    for (IndexType ir = 1; ir != irmax + 1; ++ir) {
      IndexType icmax = chambersInRingOfStation(is, ir);
      for (IndexType ic = 1; ic != icmax + 1; ++ic) {
        tChamberLabel[++count] = is * 1000 + ir * 100 + ic;
      }
    }
  }
  return tChamberLabel;
}

gasGainIndex() [1/2]

IndexType CSCIndexer::gasGainIndex ( const CSCDetId &  id, IndexType  istrip, IndexType  iwire  ) const

inline

Linear index for Gas gain sector, based on CSCDetId 'id', cathode strip 'istrip' and anode wire 'iwire'

Output is 1-45144 (CSCs 2008) and 45145-55944 (ME42).

WARNING: Use at your own risk! You must input labels within hardware ranges (e.g., 'id' must correspond to a specific layer 1-6). No trapping on out-of-range values!

Definition at line 379 of file CSCIndexer.h.

References relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, phase1PixelTopology::layer, relativeConstraints::ring, and relativeConstraints::station.

                                                                                      {
    return gasGainIndex(id.endcap(), id.station(), id.ring(), id.chamber(), id.layer(), istrip, iwire);
  }

gasGainIndex() [2/2]

IndexType CSCIndexer::gasGainIndex ( IndexType  ie, IndexType  is, IndexType  ir, IndexType  ic, IndexType  il, IndexType  istrip, IndexType  iwire  ) const

inline

Linear index for Gas gain sector, based on the cathode strip 'istrip' and anode wire 'iwire' located in layer 'il' of chamber 'ic' of ring 'ir' in station 'is' of endcap 'ie'.

Output is 1-45144 (CSCs 2008) and 45145-55944 (ME42).

WARNING: Use at your own risk! You must input labels within hardware ranges. No trapping on out-of-range values!

Definition at line 392 of file CSCIndexer.h.

References chipIndex(), chipsPerLayer(), hvSegmentIndex(), sectorsPerLayer(), and sectorStart().

                                                                                                                     {
    IndexType ichip = this->chipIndex(istrip);
    IndexType ihvsegment = this->hvSegmentIndex(is, ir, iwire);
    return sectorStart(ie, is, ir) + ((ic - 1) * 6 + il - 1) * sectorsPerLayer(is, ir) +
           (ihvsegment - 1) * chipsPerLayer(is, ir) + (ichip - 1);
  }

hvSegmentIndex()

IndexType CSCIndexer::hvSegmentIndex ( IndexType  is, IndexType  ir, IndexType  iwire  ) const

inline

Linear index for HV segment

Output is 1-5.

WARNING: Use at your own risk! The supplied CSCDetId must be chamber station, ring, and wire. No trapping on out-of-range values!

Definition at line 297 of file CSCIndexer.h.

Referenced by gasGainIndex().

                                                                              {
    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;
  }

hvSegmentsPerLayer()

IndexType CSCIndexer::hvSegmentsPerLayer ( IndexType  is, IndexType  ir  ) const

inline

Number of HV segments per layer in a chamber in ring ir of station is.

Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 1 of the 5 total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.

Definition at line 272 of file CSCIndexer.h.

Referenced by sectorsPerLayer().

                                                                 {
    const IndexType nSinL[12] = {1, 3, 3, 3, 5, 0, 3, 5, 0, 3, 5, 0};
    return nSinL[(is - 1) * 3 + ir - 1];
  }

layerIndex() [1/2]

IndexType CSCIndexer::layerIndex ( const CSCDetId &  id ) const

inline

Linear index to label each hardware layer in CSC system. Argument is the CSCDetId of some CSCLayer.

Output is 1-2808 (CSCs 2008) 2809-3240 (ME42)

WARNING: Do not input ME1a values (i.e. ring '4'): does not consider ME1a and ME1b to be separate, No sanity checking on input value: if you supply an ME1a CSCDetId then you'll get nonsense.

Definition at line 78 of file CSCIndexer.h.

References relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, phase1PixelTopology::layer, relativeConstraints::ring, and relativeConstraints::station.

                                                 {
    return layerIndex(id.endcap(), id.station(), id.ring(), id.chamber(), id.layer());
  }

layerIndex() [2/2]

IndexType CSCIndexer::layerIndex ( IndexType  ie, IndexType  is, IndexType  ir, IndexType  ic, IndexType  il  ) const

inline

Linear index for layer 'il' of chamber 'ic' in ring 'ir' of station 'is' in endcap 'ie', in range 1-2808 (CSCs 2008) or 2809-3240 (ME42).

Definition at line 104 of file CSCIndexer.h.

References chamberIndex().

                                                                                                   {
    const IndexType layersInChamber = 6;
    return (chamberIndex(ie, is, ir, ic) - 1) * layersInChamber + il;
  }

ringsInStation()

static IndexType CSCIndexer::ringsInStation ( IndexType  is )

inlinestatic

How many rings are there in station is=1, 2, 3, 4 ?

BEWARE! Includes ME42 so claims 2 rings in station 4. There is only 1 at CSC installation 2008.

Definition at line 114 of file CSCIndexer.h.

Referenced by fillChamberLabel().

                                                {
    const IndexType nrins[5] = {0, 3, 2, 2, 2};  // rings per station
    return nrins[is];
  }

sectorsPerLayer()

IndexType CSCIndexer::sectorsPerLayer ( IndexType  is, IndexType  ir  ) const

inline

Number of Gas Gain sectors per layer in a chamber in ring ir of station is.

Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 1 of the 5 total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.

Definition at line 285 of file CSCIndexer.h.

References chipsPerLayer(), and hvSegmentsPerLayer().

Referenced by gasGainIndex().

                                                              {
    return chipsPerLayer(is, ir) * hvSegmentsPerLayer(is, ir);
  }

sectorStart()

IndexType CSCIndexer::sectorStart ( IndexType  ie, IndexType  is, IndexType  ir  ) const

inline

Linear index for 1st Gas gain sector in ring 'ir' of station 'is' in endcap 'ie'.

Endcap label range 1-2, Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 1 of the 5 chips total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.

Definition at line 351 of file CSCIndexer.h.

Referenced by gasGainIndex().

                                                                        {
    // There are 36 chambers * 6 layers * 5 CFEB's * 1 HV segment = 1080 gas-gain sectors in ME1/1
    // There are 36*6*5*3 = 3240 gas-gain sectors in ME1/2
    // There are 36*6*4*3 = 2592 gas-gain sectors in ME1/3
    // There are 18*6*5*3 = 1620 gas-gain sectors in ME[2-4]/1
    // There are 36*6*5*5 = 5400 gas-gain sectors in ME[2-4]/2
    // Start of -z channels (CSCs 2008) is 22572 + 1 = 22573
    // Start of +z (ME42) is 45144 + 1 = 45145
    // Start of -z (ME42) is 45144 + 1 + 5400 = 50545
    const IndexType nStart[24] = {1,     1081,  4321,   //ME+1/1,ME+1/2,ME+1/3
                                  6913,  8533,  0,      //ME+2/1,ME+2/2,ME+2/3
                                  13933, 15553, 0,      //ME+3/1,ME+3/2,ME+3/3
                                  20953, 45145, 0,      //ME+4/1,ME+4/2,ME+4/3 (note, ME+4/2 index follows ME-4/1...)
                                  22573, 23653, 26893,  //ME-1/1,ME-1/2,ME-1/3
                                  29485, 31105, 0,      //ME-2/1,ME-2/2,ME-2/3
                                  36505, 38125, 0,      //ME-3/1,ME-3/2,ME-3/3
                                  43525, 50545, 0};     //ME-4/1,ME-4/2,ME-4/3 (note, ME-4/2 index follows ME+4/2...)
    return nStart[(ie - 1) * 12 + (is - 1) * 3 + ir - 1];
  }

startChamberIndexInEndcap()

IndexType CSCIndexer::startChamberIndexInEndcap ( IndexType  ie, IndexType  is, IndexType  ir  ) const

inline

Starting index for first chamber in ring 'ir' of station 'is' in endcap 'ie', in range 1-468 (CSCs 2008) or 469-540 (ME42).

Definition at line 86 of file CSCIndexer.h.

Referenced by chamberIndex(), and detIdFromChamberIndex_OLD().

                                                                                      {
    const IndexType nschin[24] = {1,   37,  73,  109, 127, 0, 163, 181, 0, 217, 469, 0,
                                  235, 271, 307, 343, 361, 0, 397, 415, 0, 451, 505, 0};
    return nschin[(ie - 1) * 12 + (is - 1) * 3 + ir - 1];
  }

stripChannelIndex() [1/2]

LongIndexType CSCIndexer::stripChannelIndex ( IndexType  ie, IndexType  is, IndexType  ir, IndexType  ic, IndexType  il, IndexType  istrip  ) const

inline

Linear index for strip channel istrip in layer 'il' of chamber 'ic' of ring 'ir' in station 'is' of endcap 'ie'.

Output is 1-217728 (CSCs 2008) or 217729-252288 (ME42).

WARNING: Use at your own risk! You must input labels within hardware ranges. No trapping on out-of-range values!

Definition at line 173 of file CSCIndexer.h.

References stripChannelsPerLayer(), and stripChannelStart().

Referenced by stripChannelIndex().

                                                                                                    {
    return stripChannelStart(ie, is, ir) + ((ic - 1) * 6 + il - 1) * stripChannelsPerLayer(is, ir) + (istrip - 1);
  }

stripChannelIndex() [2/2]

LongIndexType CSCIndexer::stripChannelIndex ( const CSCDetId &  id, IndexType  istrip  ) const

inline

Linear index for strip channel 'istrip' in layer labelled by CSCDetId 'id'.

Output is 1-217728 (CSCs 2008) or 217729-252288 (ME42).

WARNING: Use at your own risk! The supplied CSCDetId must be a layer id. No trapping on out-of-range values!

Definition at line 187 of file CSCIndexer.h.

References relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, phase1PixelTopology::layer, relativeConstraints::ring, relativeConstraints::station, and stripChannelIndex().

                                                                              {
    return stripChannelIndex(id.endcap(), id.station(), id.ring(), id.chamber(), id.layer(), istrip);
  }

stripChannelsPerLayer()

IndexType CSCIndexer::stripChannelsPerLayer ( IndexType  is, IndexType  ir  ) const

inline

Number of strip channels per layer in a chamber in ring ir of station is.

Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 16 of the 80 total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.
Considers ME42 as standard 80-strip per layer chambers.

Definition at line 140 of file CSCIndexer.h.

Referenced by stripChannelIndex().

                                                                    {
    const IndexType nSCinC[12] = {80, 80, 64, 80, 80, 0, 80, 80, 0, 80, 80, 0};
    return nSCinC[(is - 1) * 3 + ir - 1];
  }

stripChannelStart()

LongIndexType CSCIndexer::stripChannelStart ( IndexType  ie, IndexType  is, IndexType  ir  ) const

inline

Linear index for 1st strip channel in ring 'ir' of station 'is' in endcap 'ie'.

Endcap label range 1-2, Station label range 1-4, Ring label range 1-3.

WARNING: ME1a channels are the last 16 of the 80 total in each layer of an ME11 chamber, and an input ir=4 is invalid and will give nonsense.

Definition at line 153 of file CSCIndexer.h.

Referenced by stripChannelIndex().

                                                                                  {
    // These are in the ranges 1-217728 (CSCs 2008) and 217729-252288 (ME42).
    // There are 1-108884 channels per endcap (CSCs 2008) and 17280 channels per endcap (ME42).
    // Start of -z channels (CSCs 2008) is 108864 + 1 = 108865
    // Start of +z (ME42) is 217728 + 1 = 217729
    // Start of -z (ME42) is 217728 + 1 + 17280 = 235009
    const LongIndexType nStart[24] = {1,      17281,  34561,  48385,  57025,  0, 74305,  82945,  0, 100225, 217729, 0,
                                      108865, 126145, 143425, 157249, 165889, 0, 183169, 191809, 0, 209089, 235009, 0};
    return nStart[(ie - 1) * 12 + (is - 1) * 3 + ir - 1];
  }