|
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
|
/**
\file
test file for CSCDetId
\author Stefano ARGIRO
\date 27 Jul 2005
*/
#include <cppunit/extensions/HelperMacros.h>
#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include <FWCore/Utilities/interface/Exception.h>
#include <iostream>
using namespace std;
class testCSCDetId : public CppUnit::TestFixture {
CPPUNIT_TEST_SUITE(testCSCDetId);
CPPUNIT_TEST(testOne);
CPPUNIT_TEST(testFail);
CPPUNIT_TEST(testStatic);
CPPUNIT_TEST_SUITE_END();
public:
void setUp() {}
void tearDown() {}
void testOne();
void testFail();
void testStatic();
private:
void testValid(CSCDetId&);
};
///registration of the test so that the runner can find it
CPPUNIT_TEST_SUITE_REGISTRATION(testCSCDetId);
void testCSCDetId::testOne() {
/*
std::cout << "\ntestCSCDetId: testOne starting... " << std::endl;
std::cout << "min endcap = " << CSCDetId::minEndcapId() << std::endl;
std::cout << "max endcap = " << CSCDetId::maxEndcapId() << std::endl;
std::cout << "min station = " << CSCDetId::minStationId() << std::endl;
std::cout << "max station = " << CSCDetId::maxStationId() << std::endl;
*/
unsigned short rings[] = {4, 2, 2, 2};
unsigned short chambers[] = {36, 36, 36, 36, 18, 36, 18, 36, 18, 36};
for (int endcap = CSCDetId::minEndcapId(); endcap <= CSCDetId::maxEndcapId(); ++endcap)
for (int station = CSCDetId::minStationId(); station <= CSCDetId::maxStationId(); ++station)
for (int ring = CSCDetId::minRingId(); ring <= rings[station - 1]; ++ring) {
unsigned short itype = CSCDetId::iChamberType(station, ring);
unsigned short nchambers = chambers[itype - 1];
for (int chamber = CSCDetId::minChamberId(); chamber <= nchambers; ++chamber)
for (int layer = CSCDetId::minLayerId(); layer <= CSCDetId::maxLayerId(); ++layer) {
CSCDetId detid(endcap, station, ring, chamber, layer);
// std::cout << "detid = " << detid.rawId() << " " << hex << detid.rawId() <<
// " " << oct << detid.rawId() << dec << std::endl;
// std::cout << "\ndetid.endcap()= " << detid.endcap() << " endcap = " << endcap << std::endl;
CPPUNIT_ASSERT(detid.endcap() == endcap);
// std::cout << "\ndetid.station()= " << detid.station() << " station = " << station << std::endl;
CPPUNIT_ASSERT(detid.station() == station);
CPPUNIT_ASSERT(detid.ring() == ring);
CPPUNIT_ASSERT(detid.chamber() == chamber);
CPPUNIT_ASSERT(detid.layer() == layer);
// check chamber type number
// std::cout << "E" << endcap << " S" << station << " R" << ring << " C" << chamber <<
// " L" << layer << " chamber type=" << detid.iChamberType() << std::endl;
CPPUNIT_ASSERT(detid.iChamberType() == CSCDetId::iChamberType(station, ring));
// test constructor from id
int myId = detid.rawId();
CSCDetId anotherId(myId);
CPPUNIT_ASSERT(detid == anotherId);
}
}
}
void testCSCDetId::testFail() {
// construct using an invalid input index
// Invalid layer
std::cout << "\nConstruct CSCDetId using an invalid input index\n";
CSCDetId detid1(3, 1, 1, 1, 7);
testValid(detid1);
// contruct using an invalid input id
std::cout << "\nConstruct CSCDetId using an invalid input id\n";
CSCDetId detid2(3211);
testValid(detid2);
}
void testCSCDetId::testStatic() {
int ie = 2;
int is = 3;
int ir = 2;
int ic = 10;
int il = 3;
int id1 = CSCDetId::rawIdMaker(2, 3, 2, 10, 3);
int id2 = CSCDetId::rawIdMaker(ie, is, ir, ic, il);
int id3 = CSCDetId::rawIdMaker(ie, is, ir, ic, 0); // all layers i.e. chamber id
// std::cout << "\nE" << ie << " S" << is << " R" << ir << " C" << ic
// << " L" << il << " has rawId = " << id2 << " (dec) = "
// << hex << id2 << " (hex) " << oct << id2 << " (oct)" << dec << std::endl;
//
// std::cout << "\nE" << ie << " S" << is << " R" << ir << " C" << ic
// << " L0" << " has rawId = " << id3 << " = "
// << hex << id3 << " (hex) " << oct << id3 << " (oct)" << dec << std::endl;
CPPUNIT_ASSERT(id1 == id2);
CPPUNIT_ASSERT(CSCDetId::endcap(id2) == ie);
CPPUNIT_ASSERT(CSCDetId::station(id2) == is);
CPPUNIT_ASSERT(CSCDetId::ring(id2) == ir);
CPPUNIT_ASSERT(CSCDetId::chamber(id2) == ic);
CPPUNIT_ASSERT(CSCDetId::layer(id2) == il);
CPPUNIT_ASSERT(CSCDetId::chamber(id3) == ic);
}
void testCSCDetId::testValid(CSCDetId& detId) {
if (detId.det() != DetId::Muon || detId.subdetId() != MuonSubdetId::CSC) {
std::cout << "Invalid CSCDetId:"
<< " det: " << detId.det() << " subdet: " << detId.subdetId() << " is not a valid CSC id\n";
}
int iendcap = detId.endcap();
int istation = detId.station();
int iring = detId.ring();
int ichamber = detId.chamber();
int ilayer = detId.layer();
if (iendcap < detId.minEndcapId() || iendcap > detId.maxEndcapId() || istation < detId.minStationId() ||
istation > detId.maxStationId() || iring < detId.minRingId() || iring > detId.maxRingId() ||
ichamber < detId.minChamberId() || ichamber > detId.maxChamberId() || ilayer < detId.minLayerId() ||
ilayer > detId.maxLayerId()) {
std::cout << "Invalid CSCDetId:"
<< " Invalid parameters: "
<< " E:" << iendcap << "(" << detId.minEndcapId() << ":" << detId.maxEndcapId() << ")"
<< " S:" << istation << "(" << detId.minStationId() << ":" << detId.maxStationId() << ")"
<< " R:" << iring << "(" << detId.minRingId() << ":" << detId.maxRingId() << ")"
<< " C:" << ichamber << "(" << detId.minChamberId() << ":" << detId.maxChamberId() << ")"
<< " L:" << ilayer << "(" << detId.minLayerId() << ":" << detId.maxLayerId() << ")"
<< "\n";
}
}
|