CSCIndexerPostls1.h

Go to the documentation of this file.

#ifndef CSCIndexerPostls1_H
#define CSCIndexerPostls1_H

#include <CalibMuon/CSCCalibration/interface/CSCIndexerBase.h>

class CSCIndexerPostls1 : public CSCIndexerBase
{
public:

  ~CSCIndexerPostls1();

  std::string name() const { return "CSCIndexerPostls1"; }



  LongIndexType maxStripChannelIndex() const { return 273024; }
  IndexType maxChipIndex() const             { return 17064; }
  IndexType maxGasGainIndex() const          { return 57240; }




  IndexType onlineRingsInStation( IndexType is ) const
  {
    const IndexType nrings[5] = { 0, 4, 2, 2, 2 };
    return nrings[is];
  }

  IndexType stripChannelsPerOfflineLayer( IndexType is, IndexType ir ) const
  {
    const IndexType nSC[16] = { 64,80,64,48,  80,80,0,0,  80,80,0,0,  80,80,0,0 };
    return nSC[(is-1)*4 + ir - 1];
  }

  IndexType stripChannelsPerOnlineLayer( IndexType is, IndexType ir ) const
  {
    const IndexType nSC[16] = { 64,80,64,48,  80,80,0,0,  80,80,0,0,  80,80,0,0 };
    return nSC[(is-1)*4 + ir - 1];
  }

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





  IndexType stripChannelsPerLayer( IndexType is, IndexType ir ) const
  {
    const IndexType nSCinC[16] = { 80,80,64,48,  80,80,0,0,  80,80,0,0,  80,80,0,0 };
    return nSCinC[(is - 1)*4 + ir - 1];
  }

  LongIndexType stripChannelStart( IndexType ie, IndexType is, IndexType ir ) const
  {
    // These are in the ranges 1-217728 (CSCs 2008), 217729-252288 (ME42), and 252289-273024 (unganged ME1a)
    // There are 1-108884 channels per endcap (CSCs 2008), 17280 channels per endcap (ME42),
    // and 10368 channels per endcap (unganged ME1a)
    // 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
    // Start of +z (unganged ME1a) is 252288 + 1         = 252289
    // Start of -z (unganged ME1a) is 252288 + 1 + 10368 = 262657
    const LongIndexType nStart[32] =
      {      1, 17281, 34561,252289,   48385, 57025,0,0,   74305, 82945,0,0,  100225,217729,0,0,
        108865,126145,143425,262657,  157249,165889,0,0,  183169,191809,0,0,  209089,235009,0,0 };
    return  nStart[(ie - 1)*16 + (is - 1)*4 + ir - 1];
  }




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

  IndexType chipStart( IndexType ie, IndexType is, IndexType ir ) const
  {
    // These are in the ranges 1-13608 (CSCs 2008) and 13609-15768 (ME42) and 15769-17064 (ME1a).
    // There are 1-6804 chips per endcap (CSCs 2008) and 1080 chips per endcap (ME42) and 648 chips per endcap (ME1a).
    // 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
    // Start of +z (ME1a) is 15768 + 1 = 15769
    // Start of -z (ME1a) is 15768 + 1 + 648 = 16417
    const IndexType nStart[32] =
      {1,   1081, 2161, 15769,   3025, 3565, 0,0,  4645, 5185, 0,0,  6265, 13609,0,0,
       6805,7885, 8965, 16417,   9829, 10369,0,0,  11449,11989,0,0,  13069,14689,0,0 };
    return  nStart[(ie - 1)*16 + (is - 1)*4 + ir - 1];
  }




  IndexType sectorStart( IndexType ie, IndexType is, IndexType ir ) const
  {
    // There are 36 chambers * 6 layers * 5 CFEB's * 1 HV segment = 1080 gas-gain sectors in ME1/1 (non-upgraded)
    // There are 36 chambers * 6 layers * 3 CFEB's * 1 HV segment = 648 gas-gain sectors in ME1/1a (upgraded)
    // 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
    // Start of +z (ME1a) is 45144 + 1 + 2*5400 = 55945
    // Start of -z (ME42) is 45144 + 1 + 2*5400 + 648 = 56593
    const IndexType nStart[32] =
      {1    , 1081 ,  4321, 55945,  //ME+1/1,ME+1/2,ME+1/3,ME+1/a
       6913 , 8533 ,     0,     0,  //ME+2/1,ME+2/2
       13933, 15553,     0,     0,  //ME+3/1,ME+3/2
       20953, 45145,     0,     0,  //ME+4/1,ME+4/2 (note, ME+4/2 index follows ME-4/1...)
       22573, 23653, 26893, 56593,  //ME-1/1,ME-1/2,ME-1/3,ME+1/a
       29485, 31105,     0,     0,  //ME-2/1,ME-2/2
       36505, 38125,     0,     0,  //ME-3/1,ME-3/2
       43525, 50545,     0,     0 };//ME-4/1,ME-4/2 (note, ME-4/2 index follows ME+4/2...)
    return  nStart[(ie-1)*16 + (is-1)*4 + ir - 1];
  }

  std::pair<CSCDetId, IndexType> detIdFromStripChannelIndex( LongIndexType ichi ) const;
  std::pair<CSCDetId, IndexType> detIdFromChipIndex( IndexType ichi ) const;
  CSCIndexerBase::GasGainIndexType detIdFromGasGainIndex( IndexType igg ) const;

  int dbIndex(const CSCDetId & id, int & channel) const;

};

#endif