|
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
|
#include <cstdlib>
#include <iostream>
#include <memory>
#include <set>
#include <string>
#include <utility>
#include "DetectorDescription/Core/interface/DDRotationMatrix.h"
#include "DetectorDescription/Core/interface/DDTranslation.h"
#include "DetectorDescription/Core/interface/DDExpandedNode.h"
#include "DetectorDescription/Core/interface/DDExpandedView.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/DDPosData.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
using namespace std;
void DDExpandedViewDump(ostream& os, DDExpandedView& ex, size_t skip, size_t sto) {
typedef set<string>::iterator s_iter;
set<string> result;
bool go(true);
int count(0);
bool dotrans(true);
if (std::getenv("DDNOTRANS"))
dotrans = false;
++skip;
while (go) {
if ((count % 5000) == 0)
cout << count << ' ' << flush;
if (sto > 0)
if ((count % sto) == 0)
break;
++count;
if ((count % skip) == 0) {
stringstream s;
s.setf(ios_base::fixed, ios_base::floatfield);
s.precision(2);
s << ex.logicalPart().name() << ' ' << ex.copyno() << ' ' << ex.geoHistory()
<< " r=" << ex.geoHistory().back().posdata()->ddrot().name() << "\n";
DDRotationMatrix rm = ex.rotation();
{
double v[9];
rm.GetComponents(v, v + 9);
s << "R=(";
s << v[0] << ' ' << v[1] << ' ' << v[2] << endl;
s << v[3] << ' ' << v[4] << ' ' << v[5] << endl;
s << v[6] << ' ' << v[7] << ' ' << v[7] << endl;
}
rm = ex.geoHistory().back().posdata()->rotation();
{
double v[9];
rm.GetComponents(v, v + 9);
s << "r=(";
s << v[0] << ' ' << v[1] << ' ' << v[2] << endl;
s << v[3] << ' ' << v[4] << ' ' << v[5] << endl;
s << v[6] << ' ' << v[7] << ' ' << v[7] << endl;
}
if (dotrans) {
s << "T=(" << ex.translation().x() << ',' << ex.translation().y() << ',' << ex.translation().z() << ") ";
}
pair<s_iter, bool> ins = result.insert(s.str());
if (!ins.second) {
cout << "DDExpandedViewDump:ERROR: duplicated=" << s.str() << endl;
}
}
go = ex.next();
}
s_iter it(result.begin()), ed(result.end());
for (; it != ed; ++it) {
os << *it << endl;
}
}
|