|
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
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
|
// -*- C++ -*-
//
// Package: TestTrackerHierarchy
// Class: TestTrackerHierarchy
//
//
// Description: Module to test the Alignment software
//
//
// Original Author: Frederic Ronga
// Created: March 16, 2006
// $Id: TestTrackerHierarchy.cpp,v 1.6 2012/06/30 08:59:35 eulisse Exp $
// system include files
#include <sstream>
#include <iomanip>
#include <memory>
// user include files
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "CondFormats/Alignment/interface/AlignTransform.h"
#include "CondFormats/Alignment/interface/Alignments.h"
#include "Alignment/TrackerAlignment/interface/AlignableTracker.h"
#include "Alignment/CommonAlignment/interface/AlignableNavigator.h"
#include "Alignment/CommonAlignment/interface/AlignableObjectId.h"
static const int kLEAD_WIDTH = 40; // First field width
//
//
// class declaration
//
class TestTrackerHierarchy : public edm::one::EDAnalyzer<> {
public:
explicit TestTrackerHierarchy(const edm::ParameterSet& pSet)
: tTopoToken_(esConsumes()),
tkGeomToken_(esConsumes()),
aliToken_(esConsumes()),
aliErrToken_(esConsumes()),
dumpAlignments_(pSet.getUntrackedParameter<bool>("dumpAlignments")) {}
virtual void analyze(const edm::Event&, const edm::EventSetup&);
private:
// ----------member data ---------------------------
const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tTopoToken_;
const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
const edm::ESGetToken<Alignments, TrackerAlignmentRcd> aliToken_;
const edm::ESGetToken<AlignmentErrorsExtended, TrackerAlignmentErrorExtendedRcd> aliErrToken_;
void dumpAlignable(const Alignable*, unsigned int, unsigned int);
void printInfo(const Alignable*, unsigned int);
void dumpAlignments(const edm::EventSetup& setup, AlignableTracker* aliTracker) const;
std::string leaders_, blank_, filled_;
const bool dumpAlignments_;
std::unique_ptr<AlignableTracker> alignableTracker_;
};
void TestTrackerHierarchy::analyze(const edm::Event&, const edm::EventSetup& setup) {
//Retrieve tracker topology from geometry
const TrackerTopology* const tTopo = &setup.getData(tTopoToken_);
edm::LogInfo("TrackerHierarchy") << "Starting!";
const TrackerGeometry* trackerGeometry = &setup.getData(tkGeomToken_);
alignableTracker_ = std::make_unique<AlignableTracker>(&(*trackerGeometry), tTopo);
leaders_ = "";
blank_ = " "; // These two...
filled_ = "| "; // ... must have the same length
// Now dump mother of each alignable
//const Alignable* alignable = (&(*theAlignableTracker))->pixelHalfBarrels()[0];
this->dumpAlignable(alignableTracker_.get(), 1, 1);
edm::LogInfo("TrackerHierarchy") << "Done!";
if (dumpAlignments_) {
this->dumpAlignments(setup, alignableTracker_.get());
}
}
//__________________________________________________________________________________________________
// Recursive loop on alignable hierarchy
void TestTrackerHierarchy::dumpAlignable(const Alignable* alignable, unsigned int idau, unsigned int ndau) {
printInfo(alignable, idau);
if (ndau != idau)
leaders_ += filled_;
else
leaders_ += blank_;
const align::Alignables& comps = alignable->components();
if (unsigned int ndau = comps.size()) {
unsigned int idau = 0;
for (align::Alignables::const_iterator iter = comps.begin(); iter != comps.end(); ++iter)
dumpAlignable(*iter, ++idau, ndau);
}
leaders_ = leaders_.substr(0, leaders_.length() - blank_.length());
}
//__________________________________________________________________________________________________
// Do the actual printout
void TestTrackerHierarchy::printInfo(const Alignable* alignable, unsigned int idau) {
int width = kLEAD_WIDTH - leaders_.length();
std::ostringstream name, pos, rot;
name << alignableTracker_->objectIdProvider().idToString(alignable->alignableObjectId()) << idau;
// Position
pos.setf(std::ios::fixed);
pos << "(" << std::right << std::setw(8) << std::setprecision(4) << alignable->globalPosition().x() << ","
<< std::setw(8) << std::setprecision(4) << alignable->globalPosition().y() << "," << std::setw(8)
<< std::setprecision(4) << alignable->globalPosition().z() << ")";
edm::LogVerbatim("DumpAlignable") << leaders_ << "+-> " << std::setw(width) << std::left << name.str() << " | "
<< std::setw(3) << std::left << alignable->components().size() << " | "
<< std::setw(11) << std::left << alignable->id() << " | " << pos.str();
}
//__________________________________________________________________________________________________
void TestTrackerHierarchy::dumpAlignments(const edm::EventSetup& setup, AlignableTracker* aliTracker) const {
const Alignments* alignments = &setup.getData(aliToken_);
if (alignments->empty()) {
edm::LogWarning("TrackerAlignment") << "@SUB=dumpAlignments"
<< "No TrackerAlignmentRcd.";
} else {
AlignableNavigator navi(aliTracker);
edm::LogInfo("TrackerAlignment") << "@SUB=dumpAlignments"
<< "Start dumping alignments.";
unsigned int nProblems = 0;
for (std::vector<AlignTransform>::const_iterator iAlign = alignments->m_align.begin(),
iEnd = alignments->m_align.end();
iAlign != iEnd;
++iAlign) {
const align::ID id = (*iAlign).rawId();
const AlignTransform::Translation pos((*iAlign).translation());
edm::LogVerbatim("DumpAlignable") << (*iAlign).rawId() << " | " << pos;
AlignableDetOrUnitPtr aliPtr = navi.alignableFromDetId(id);
if (!aliPtr.isNull()) {
const Alignable::PositionType& aliPos = aliPtr->globalPosition();
double dR = aliPos.perp() - pos.perp();
double dRphi = (aliPos.phi() - pos.phi()) * pos.perp();
double dZ = aliPos.z() - pos.z();
if (dR * dR + dRphi * dRphi + dZ * dZ) {
++nProblems;
edm::LogWarning("Alignment") << "@SUB=analyze"
<< "Delta r,rphi,z: " << dR << " " << dRphi << " " << dZ
<< "\nPos r,phi,z: " << pos.perp() << " " << pos.phi() << " " << pos.z();
}
} else {
++nProblems;
edm::LogWarning("Alignment") << "@SUB=dumpAlignments"
<< "No Alignable for Id " << id;
}
} // ending loop
if (nProblems) {
edm::LogWarning("TrackerAlignment") << "@SUB=dumpAlignments"
<< "Ending: " << nProblems << " Alignments with problems.";
} else {
edm::LogInfo("TrackerAlignment") << "@SUB=dumpAlignments"
<< "Ending without problem.";
}
}
}
//define this as a plug-in
DEFINE_FWK_MODULE(TestTrackerHierarchy);
|