29 for (
unsigned int i = 0;
i < 4;
i++) {
50 sprintf(
title,
"DT Station %i Number of Segments (No Arbitration)",
station + 1);
88 sprintf(
title,
"CSC Station %i Number of Segments (No Arbitration)",
station + 1);
100 sprintf(
title,
"CSC Station %i Delta DxDz",
station + 1);
104 sprintf(
title,
"CSC Station %i Pull DxDz",
station + 1);
116 sprintf(
title,
"CSC Station %i Delta DyDz",
station + 1);
120 sprintf(
title,
"CSC Station %i Pull DyDz",
station + 1);
128 using namespace reco;
137 for (
unsigned int i = 0;
i < 4;
i++) {
194 LocalPoint segmentLocalPosition = segment->localPosition();
195 LocalVector segmentLocalDirection = segment->localDirection();
196 LocalError segmentLocalPositionError = segment->localPositionError();
197 LocalError segmentLocalDirectionError = segment->localDirectionError();
198 bool segmentFound =
false;
201 if (!
muon->isMatchesValid())
204 for (std::vector<MuonChamberMatch>::const_iterator chamberMatch =
muon->matches().begin();
205 chamberMatch !=
muon->matches().end();
207 for (std::vector<MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
208 segmentMatch != chamberMatch->segmentMatches.end();
210 if (fabs(segmentMatch->x - segmentLocalPosition.
x()) < 1E-6 &&
211 fabs(segmentMatch->y - segmentLocalPosition.
y()) < 1E-6 &&
212 fabs(segmentMatch->dXdZ - segmentLocalDirection.
x() / segmentLocalDirection.
z()) < 1E-6 &&
213 fabs(segmentMatch->dYdZ - segmentLocalDirection.
y() / segmentLocalDirection.
z()) < 1E-6 &&
214 fabs(segmentMatch->xErr -
sqrt(segmentLocalPositionError.
xx())) < 1E-6 &&
215 fabs(segmentMatch->yErr -
sqrt(segmentLocalPositionError.
yy())) < 1E-6 &&
216 fabs(segmentMatch->dXdZErr -
sqrt(segmentLocalDirectionError.
xx())) < 1E-6 &&
217 fabs(segmentMatch->dYdZErr -
sqrt(segmentLocalDirectionError.
yy())) < 1E-6) {
238 LocalPoint segmentLocalPosition = segment->localPosition();
239 LocalVector segmentLocalDirection = segment->localDirection();
240 LocalError segmentLocalPositionError = segment->localPositionError();
241 LocalError segmentLocalDirectionError = segment->localDirectionError();
242 bool segmentFound =
false;
245 if (!
muon->isMatchesValid())
248 for (std::vector<MuonChamberMatch>::const_iterator chamberMatch =
muon->matches().begin();
249 chamberMatch !=
muon->matches().end();
251 for (std::vector<MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
252 segmentMatch != chamberMatch->segmentMatches.end();
254 if (fabs(segmentMatch->x - segmentLocalPosition.
x()) < 1E-6 &&
255 fabs(segmentMatch->y - segmentLocalPosition.
y()) < 1E-6 &&
256 fabs(segmentMatch->dXdZ - segmentLocalDirection.
x() / segmentLocalDirection.
z()) < 1E-6 &&
257 fabs(segmentMatch->dYdZ - segmentLocalDirection.
y() / segmentLocalDirection.
z()) < 1E-6 &&
258 fabs(segmentMatch->xErr -
sqrt(segmentLocalPositionError.
xx())) < 1E-6 &&
259 fabs(segmentMatch->yErr -
sqrt(segmentLocalPositionError.
yy())) < 1E-6 &&
260 fabs(segmentMatch->dXdZErr -
sqrt(segmentLocalDirectionError.
xx())) < 1E-6 &&
261 fabs(segmentMatch->dYdZErr -
sqrt(segmentLocalDirectionError.
yy())) < 1E-6) {
281 if (fabs(
f) > 900000)