DTSequentialCellNumber.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  $Date: 2010/05/06 16:23:00 $
 *  $Revision: 1.1 $
 *  \author Paolo Ronchese INFN Padova
 *
 */

//----------------------
// This Class' Header --
//----------------------
#include "CondFormats/DTObjects/interface/DTSequentialCellNumber.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------

//---------------
// C++ Headers --
//---------------
namespace {
  constexpr int cellsInTheta = 58;
  constexpr int cellsPerLayer[] = {0, 50, 60, 72, 96};
  constexpr int cellsInMB1 = (cellsPerLayer[1] * 8) + (cellsInTheta * 4);
  constexpr int cellsInMB2 = (cellsPerLayer[2] * 8) + (cellsInTheta * 4);
  constexpr int cellsInMB3 = (cellsPerLayer[3] * 8) + (cellsInTheta * 4);
  constexpr int cellsInMB4 = cellsPerLayer[4] * 8;

  constexpr int offsetChamber[] = {0, 0, cellsInMB1, cellsInMB1 + cellsInMB2, cellsInMB1 + cellsInMB2 + cellsInMB3};
  constexpr int cellsPerSector = cellsInMB1 + cellsInMB2 + cellsInMB3 + cellsInMB4;
  constexpr int cellsIn13Sectors = (cellsPerSector * 12) + cellsInMB4;
  constexpr int cellsPerWheel = cellsIn13Sectors + cellsInMB4;
}  // namespace

//--------------
// Operations --
//--------------
int DTSequentialCellNumber::id(int wheel, int station, int sector, int superlayer, int layer, int cell) {
  wheel += 3;
  if (wheel <= 0)
    return -1;
  if (station <= 0)
    return -2;
  if (sector <= 0)
    return -3;
  if (superlayer <= 0)
    return -4;
  if (layer <= 0)
    return -5;
  if (cell <= 0)
    return -6;

  int seqWireNum = 0;

  if (wheel > 5)
    return -1;
  seqWireNum += (wheel - 1) * cellsPerWheel;

  if (sector > 14)
    return -2;
  if (sector > 12 && station < 4)
    return -2;
  if (sector > 13)
    seqWireNum += cellsIn13Sectors;
  else
    seqWireNum += (sector - 1) * cellsPerSector;

  if (station > 4)
    return -3;
  if (sector < 13)
    seqWireNum += offsetChamber[station];

  if (superlayer > 3)
    return -4;
  if (layer > 4)
    return -5;
  if (superlayer != 2) {
    if (cell > cellsPerLayer[station])
      return -6;
    if (superlayer == 3)
      layer += 4;
    seqWireNum += (layer - 1) * cellsPerLayer[station];
  } else {
    if (station == 4)
      return -4;
    if (cell > cellsInTheta)
      return -6;
    seqWireNum += (8 * cellsPerLayer[station]) + ((layer - 1) * cellsInTheta);
  }

  return seqWireNum + cell;
}

int DTSequentialCellNumber::max() { return 5 * cellsPerWheel; }