|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
/*
* 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; }
|