|
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
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
|
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "Geometry/MuonNumbering/interface/CSCNumberingScheme.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonGeometryConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
//#define EDM_ML_DEBUG
CSCNumberingScheme::CSCNumberingScheme(const MuonGeometryConstants& muonConstants) { initMe(muonConstants); }
void CSCNumberingScheme::initMe(const MuonGeometryConstants& muonConstants) {
int theLevelPart = muonConstants.getValue("level");
theRegionLevel = muonConstants.getValue("me_region") / theLevelPart;
theStationLevel = muonConstants.getValue("me_station") / theLevelPart;
theSubringLevel = muonConstants.getValue("me_subring") / theLevelPart;
theSectorLevel = muonConstants.getValue("me_sector") / theLevelPart;
theRingLevel = muonConstants.getValue("me_ring") / theLevelPart;
theLayerLevel = muonConstants.getValue("me_layer") / theLevelPart;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "Initialize CSCNumberingScheme"
<< "\ntheRegionLevel " << theRegionLevel << "\ntheStationLevel " << theStationLevel
<< "\ntheSubringLevel " << theSubringLevel << "\ntheSectorLevel " << theSectorLevel
<< "\ntheRingLevel " << theRingLevel << "\ntheLayerLevel " << theLayerLevel;
#endif
}
int CSCNumberingScheme::baseNumberToUnitNumber(const MuonBaseNumber& num) const {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "CSCNumbering " << num.getLevels();
for (int level = 1; level <= num.getLevels(); level++) {
edm::LogVerbatim("MuonGeom") << level << " " << num.getSuperNo(level) << " " << num.getBaseNo(level);
}
#endif
int fwbw_id = 0;
int station_id = 0;
int ring_id = 0;
int subring_id = 0;
int sector_id = 0;
int layer_id = 0;
// Decode endcap levels
// We should be able to work with 6 (layer-level) or 5 (chamber-level)
for (int level = 1; level <= num.getLevels(); level++) {
if (level == theRegionLevel) {
const int copyno = num.getBaseNo(level);
fwbw_id = copyno + 1;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "endcap=" << fwbw_id;
#endif
} else if (level == theStationLevel) {
const int station_tag = num.getSuperNo(level);
station_id = station_tag;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "station=" << station_id;
#endif
} else if (level == theSubringLevel) {
const int copyno = num.getBaseNo(level);
subring_id = copyno + 1;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "subring=" << subring_id;
#endif
} else if (level == theSectorLevel) {
const int copyno = num.getBaseNo(level);
sector_id = copyno + 1;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "sector=" << sector_id;
#endif
} else if (level == theLayerLevel) {
const int copyno = num.getBaseNo(level);
layer_id = copyno + 1;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "layer=" << layer_id;
#endif
} else if (level == theRingLevel) {
const int ring_tag = num.getSuperNo(level);
ring_id = ring_tag;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "ring=" << ring_id;
#endif
}
}
// check validity
if ((fwbw_id < 1) || (fwbw_id > 2)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "forward/backward id out of range:" << fwbw_id;
}
if ((station_id < 1) || (station_id > 4)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "station id out of range:" << station_id;
}
if ((ring_id < 1) || (ring_id > 4)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "ring id out of range:" << ring_id;
}
if ((subring_id < 1) || (subring_id > 2)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "subring id out of range:" << subring_id;
}
if ((sector_id < 1) || (sector_id > 36)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "sector id out of range:" << sector_id;
}
// Allow id=0 since that means a chamber
if ((layer_id < 0) || (layer_id > 6)) {
edm::LogError("MuonGeom") << "@SUB=CSCNumberingScheme::baseNumberToUnitNumber"
<< "layer id out of range" << layer_id;
}
// find appropriate chamber label
int chamber_id = chamberIndex(station_id, ring_id, subring_id, sector_id);
// convert into raw id of appropriate DetId
int intIndex = CSCDetId::rawIdMaker(fwbw_id, station_id, ring_id, chamber_id, layer_id);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "CSCNumberingScheme : fw/bw " << fwbw_id << " station " << station_id << " ring "
<< ring_id << " subring " << subring_id << " chamber " << chamber_id << " sector "
<< sector_id << " layer " << layer_id;
#endif
return intIndex;
}
int CSCNumberingScheme::chamberIndex(int station_id, int ring_id, int subring_id, int sector_id) const {
int chamber_id = 0;
// chamber label is related to sector_id but we need to
// adjust to real hardware labelling
// Tim confirms this works properly according to CMS IN 2000/004 Version 2.5 March 2007.
if (ring_id == 3) {
chamber_id = sector_id;
} else {
if (subring_id == 1) {
chamber_id = 2 * sector_id - 1;
} else {
chamber_id = 2 * sector_id;
}
}
return chamber_id;
}
|