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	chambersInRingOfStation (IndexType is, IndexType ir) 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	ringsInStation (IndexType is) 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 ()

Private Member Functions
void	fillChamberLabel () const

Private Attributes
std::vector< IndexType >	chamberLabel

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

typedef uint16_t CSCIndexer::IndexType

Definition at line 50 of file CSCIndexer.h.

typedef uint32_t CSCIndexer::LongIndexType

Definition at line 51 of file CSCIndexer.h.

Constructor & Destructor Documentation

CSCIndexer::CSCIndexer ( )

inline

Definition at line 53 of file CSCIndexer.h.

{};

CSCIndexer::~CSCIndexer ( )

inline

Definition at line 54 of file CSCIndexer.h.

{};

Member Function Documentation

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 66 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() );
   }

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 99 of file CSCIndexer.h.

References startChamberIndexInEndcap().

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

CSCIndexer::IndexType CSCIndexer::chamberLabelFromChamberIndex

(

IndexType

ici

)

const

Definition at line 86 of file CSCIndexer.cc.

References chamberLabel, and fillChamberLabel().

Referenced by gasGainIndex().

                                                                                  {
  // 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
    }
  }
  if (chamberLabel.empty()) fillChamberLabel();
  return chamberLabel[ici];  

}

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

inline

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

Works for ME1a (ring 4 of ME1) too.

Definition at line 127 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;
   }

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 237 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);

  }

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 252 of file CSCIndexer.h.

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

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

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 266 of file CSCIndexer.h.

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

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 202 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];
   }

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 215 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];
   }

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 240 of file CSCIndexer.cc.

References RecoEcal_EventContent_cff::ec.

Referenced by CSCValidation::doADCTiming(), CSCValidation::doAFEBTiming(), CSCValidation::doCompTiming(), CSCValidation::doGasGain(), gasGainIndex(), CSCPedestals::item(), CSCGains::item(), CSCNoiseMatrix::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;
}

CSCDetId CSCIndexer::detIdFromChamberIndex

(

IndexType

ici

)

const

Definition at line 62 of file CSCIndexer.cc.

References chamberLabel, detIdFromChamberLabel(), fillChamberLabel(), and diffTwoXMLs::label.

Referenced by detIdFromLayerIndex(), gasGainIndex(), 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
    }
  }
  if (chamberLabel.empty()) fillChamberLabel();
  IndexType label = chamberLabel[ici];    
  return detIdFromChamberLabel( ie, label );
}

CSCDetId CSCIndexer::detIdFromChamberIndex_OLD

(

IndexType

ici

)

const

Definition at line 24 of file CSCIndexer.cc.

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

Referenced by gasGainIndex().

                                                                    {

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

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

Definition at line 109 of file CSCIndexer.cc.

References diffTwoXMLs::label.

Referenced by detIdFromChamberIndex(), and gasGainIndex().

                                                                                {

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

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

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

(

IndexType

ichi

)

const

Definition at line 185 of file CSCIndexer.cc.

References detIdFromLayerIndex().

Referenced by gasGainIndex().

                                                                                            {

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

CSCDetId CSCIndexer::detIdFromLayerIndex

(

IndexType

ili

)

const

Decode CSCDetId from various indexes and labels

Definition at line 120 of file CSCIndexer.cc.

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

Referenced by detIdFromChipIndex(), detIdFromStripChannelIndex(), and gasGainIndex().

                                                              {

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

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

(

LongIndexType

ichi

)

const

Definition at line 129 of file CSCIndexer.cc.

References detIdFromLayerIndex(), and nchan.

Referenced by gasGainIndex().

                                                                                                        {

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

void CSCIndexer::fillChamberLabel ( ) const

private

Definition at line 4 of file CSCIndexer.cc.

References chamberLabel, chambersInRingOfStation(), KineDebug3::count(), and ringsInStation().

Referenced by chamberLabelFromChamberIndex(), detIdFromChamberIndex(), and gasGainIndex().

                                        {
  // 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)
   chamberLabel.resize( 271 ); // one more than #chambers per endcap. Includes ME42.
   IndexType count = 0;
   chamberLabel[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 ) {
       chamberLabel[ ++count ] = is*1000 + ir*100 + ic ;
         }
      } 
   }
}

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 370 of file CSCIndexer.h.

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

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

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 383 of file CSCIndexer.h.

References chamberLabelFromChamberIndex(), chipIndex(), chipsPerLayer(), dbIndex(), detIdFromChamberIndex(), detIdFromChamberIndex_OLD(), detIdFromChamberLabel(), detIdFromChipIndex(), detIdFromLayerIndex(), detIdFromStripChannelIndex(), fillChamberLabel(), 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);
  }

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

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 278 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];
  }

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 80 of file CSCIndexer.h.

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

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

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 107 of file CSCIndexer.h.

References chamberIndex().

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

IndexType CSCIndexer::ringsInStation ( IndexType  is ) const

inline

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 117 of file CSCIndexer.h.

Referenced by fillChamberLabel().

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

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 291 of file CSCIndexer.h.

References chipsPerLayer(), and hvSegmentsPerLayer().

Referenced by gasGainIndex().

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

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 342 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];
  }

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 88 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];

   }

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

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 189 of file CSCIndexer.h.

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

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

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 142 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];
   }

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 155 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];
   }

Member Data Documentation

std::vector<IndexType> CSCIndexer::chamberLabel

mutableprivate

Definition at line 413 of file CSCIndexer.h.

Referenced by chamberLabelFromChamberIndex(), detIdFromChamberIndex(), and fillChamberLabel().