CSCIndexerStartup.h

Go to the documentation of this file.

#ifndef CSCIndexerStartup_H
#define CSCIndexerStartup_H

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

class CSCIndexerStartup : public CSCIndexerBase {
public:
  ~CSCIndexerStartup() override;

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


  LongIndexType maxStripChannelIndex() const override { return 252288; }
  IndexType maxChipIndex() const override { return 15768; }
  IndexType maxGasGainIndex() const override { return 55944; }



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

  IndexType stripChannelsPerOfflineLayer(IndexType is, IndexType ir) const override {
    const IndexType nSC[16] = {64, 80, 64, 16, 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 override {
    const IndexType nSC[16] = {80, 80, 64, 80, 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 override {
    const IndexType nCinL[16] = {5, 5, 4, 5, 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 override {
    const IndexType nSCinC[16] = {80, 80, 64, 80, 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 override {
    // 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[32] = {1, 17281,  34561,  1,      48385, 57025,  0,      0,      74305,  82945,  0,
                                      0, 100225, 217729, 0,      0,     108865, 126145, 143425, 108865, 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 override {
    const IndexType nCinL[16] = {5, 5, 4, 5, 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 override {
    // 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[32] = {1,    1081, 2161, 1,    3025, 3565,  0, 0, 4645,  5185,  0, 0, 6265,  13609, 0, 0,
                                  6805, 7885, 8965, 6805, 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 override {
    // 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[32] = {1,     1081,  4321,
                                  1,  // ME+1/1,ME+1/2,ME+1/3,ME+1/4
                                  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,ME+4/3 (note, ME+4/2 index follows ME-4/1...)
                                  22573, 23653, 26893,
                                  22573,  // ME-1/1,ME-1/2,ME-1/3, ME-1/4
                                  29485, 31105, 0,
                                  0,  // ME-2/1,ME-2/2,ME-2/3
                                  36505, 38125, 0,
                                  0,  // ME-3/1,ME-3/2,ME-3/3
                                  43525, 50545, 0,
                                  0};  // ME-4/1,ME-4/2,ME-4/3 (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 override;
  std::pair<CSCDetId, IndexType> detIdFromChipIndex(IndexType ichi) const override;
  GasGainIndexType detIdFromGasGainIndex(IndexType igg) const override;

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

#endif