|
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
|
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "Geometry/CSCGeometry/interface/CSCGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCLayer.h"
#include "Geometry/CSCGeometry/interface/CSCLayerGeometry.h"
#include <string>
#include <iomanip> // for setw() etc.
#include <vector>
class CSCGeometryOfStrips : public edm::one::EDAnalyzer<> {
public:
explicit CSCGeometryOfStrips(const edm::ParameterSet&);
~CSCGeometryOfStrips() override = default;
void beginJob() override {}
void analyze(edm::Event const&, edm::EventSetup const&) override;
void endJob() override {}
const std::string& myName() { return myName_; }
private:
const int dashedLineWidth_;
const std::string dashedLine_;
const std::string myName_;
const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> tokGeom_;
};
CSCGeometryOfStrips::CSCGeometryOfStrips(const edm::ParameterSet& iConfig)
: dashedLineWidth_(101),
dashedLine_(std::string(dashedLineWidth_, '-')),
myName_("CSCGeometryOfStrips"),
tokGeom_(esConsumes()) {}
void CSCGeometryOfStrips::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
// const double dPi = Geom::pi();
// const double radToDeg = 180. / dPi;
std::cout << myName() << ": Analyzer..." << std::endl;
std::cout << "start " << dashedLine_ << std::endl;
const edm::ESHandle<CSCGeometry>& pDD = iSetup.getHandle(tokGeom_);
std::cout << " Geometry node for CSCGeometry is " << &(*pDD) << std::endl;
std::cout << " I have " << pDD->detTypes().size() << " detTypes" << std::endl;
std::cout << " I have " << pDD->detUnits().size() << " detUnits" << std::endl;
std::cout << " I have " << pDD->dets().size() << " dets" << std::endl;
std::cout << " I have " << pDD->layers().size() << " layers" << std::endl;
std::cout << " I have " << pDD->chambers().size() << " chambers" << std::endl;
std::cout << myName() << ": Begin iteration over geometry..." << std::endl;
std::cout << "iter " << dashedLine_ << std::endl;
int icount = 0;
// Check the DetUnits
for (auto it : pDD->detUnits()) {
// Do we really have a CSC layer?
auto layer = dynamic_cast<CSCLayer const*>(it);
if (layer) {
++icount;
DetId detId = layer->geographicalId();
int id = detId(); // or detId.rawId()
std::cout << "\n"
<< "Testing CSC layer# " << icount << " id= " << id << " = " << std::oct << id << std::dec
<< " (octal) "
<< " E" << CSCDetId::endcap(id) << " S" << CSCDetId::station(id) << " R" << CSCDetId::ring(id) << " C"
<< CSCDetId::chamber(id) << " L" << CSCDetId::layer(id) << std::endl;
const CSCLayerGeometry* geom = layer->geometry();
// std::cout << *geom;
const CSCStripTopology* topol = geom->topology();
// std::cout << "\n" << *topol;
int nStrips = geom->numberOfStrips();
float mps = 30.;
if (nStrips < static_cast<int>(mps))
mps = 8.;
MeasurementPoint mp(mps, 0.3);
MeasurementError merr(0.25, 0, 0.25);
LocalPoint lp = topol->localPosition(mp);
LocalError lerr = topol->localError(mp, merr);
MeasurementPoint mp2 = topol->measurementPosition(lp);
MeasurementError merr2 = topol->measurementError(lp, lerr);
float strip = 30.;
if (nStrips < static_cast<int>(strip))
strip = 8.;
LocalPoint lps = topol->localPosition(strip);
float stripLP = topol->strip(lps);
std::cout << "Strip check: strip in=" << strip << ", strip out=" << stripLP << std::endl;
std::cout << "mp.x in=" << mp.x() << ", mp.x out=" << mp2.x() << std::endl;
std::cout << "mp.y in=" << mp.y() << ", mp.y out=" << mp2.y() << std::endl;
const float kEps = 1.E-04;
if (fabs(strip - stripLP) > kEps) {
std::cout << "...strip mismatch! " << std::endl;
}
if (fabs(mp.x() - mp2.x()) > kEps || fabs(mp.y() - mp2.y()) > kEps) {
std::cout << "...measurement point mismatch!" << std::endl;
}
std::cout << "merr.uu in=" << merr.uu() << ", merr.uu out=" << merr2.uu() << std::endl;
std::cout << "merr.uv in=" << merr.uv() << ", merr.uv out=" << merr2.uv() << std::endl;
std::cout << "merr.vv in=" << merr.vv() << ", merr.vv out=" << merr2.vv() << std::endl;
if (fabs(merr.uu() - merr2.uu()) > kEps || fabs(merr.vv() - merr2.vv()) > kEps ||
fabs(merr.uv() - merr2.uv()) > kEps) {
std::cout << "...measurement point error mismatch!" << std::endl;
}
} else {
std::cout << "Could not dynamic_cast to a CSCLayer* " << std::endl;
}
}
std::cout << dashedLine_ << " end" << std::endl;
std::cout << " UPDATED 28.10.2012" << std::endl;
}
//define this as a plug-in
DEFINE_FWK_MODULE(CSCGeometryOfStrips);
|