|
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
|
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ESTransientHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "DetectorDescription/DDCMS/interface/DDDetector.h"
#include "Geometry/Records/interface/DDVectorRegistryRcd.h"
#include "DetectorDescription/DDCMS/interface/DDVectorRegistry.h"
#include "DD4hep/Detector.h"
#include <memory>
#include <string>
using namespace std;
using namespace cms;
using namespace edm;
using namespace dd4hep;
class DDCMSDetector : public one::EDAnalyzer<> {
public:
explicit DDCMSDetector(const ParameterSet& p);
void beginJob() override {}
void analyze(Event const& iEvent, EventSetup const&) override;
void endJob() override;
private:
const ESInputTag m_tag;
const ESGetToken<DDDetector, IdealGeometryRecord> m_detectorToken;
const ESGetToken<DDVectorRegistry, DDVectorRegistryRcd> m_registryToken;
};
DDCMSDetector::DDCMSDetector(const ParameterSet& iConfig)
: m_tag(iConfig.getParameter<ESInputTag>("DDDetector")),
m_detectorToken(esConsumes(m_tag)),
m_registryToken(esConsumes(m_tag)) {}
void DDCMSDetector::analyze(const Event&, const EventSetup& iEventSetup) {
ESTransientHandle<DDDetector> det = iEventSetup.getTransientHandle(m_detectorToken);
LogVerbatim("Geometry") << "Iterate over the detectors:\n";
LogVerbatim("Geometry").log([&](auto& log) {
for (auto const& it : det->description()->detectors()) {
dd4hep::DetElement det(it.second);
log << it.first << ": " << det.path();
}
});
LogVerbatim("Geometry") << "..done!";
ESTransientHandle<DDVectorRegistry> registry = iEventSetup.getTransientHandle(m_registryToken);
LogVerbatim("Geometry") << "DD Vector Registry size: " << registry->vectors.size();
LogVerbatim("Geometry").log([&](auto& log) {
for (const auto& p : registry->vectors) {
log << " " << p.first << " => ";
for (const auto& i : p.second)
log << i << ", ";
log << '\n';
}
});
}
void DDCMSDetector::endJob() {}
DEFINE_FWK_MODULE(DDCMSDetector);
|