|
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
|
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "Geometry/HcalTowerAlgo/interface/HcalGeometry.h"
#include <iostream>
class HcalGeometryDetIdTester : public edm::one::EDAnalyzer<> {
public:
explicit HcalGeometryDetIdTester(const edm::ParameterSet&);
static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
void beginJob() override {}
void analyze(edm::Event const& iEvent, edm::EventSetup const&) override;
void endJob() override {}
private:
static constexpr int ndetMax_ = 4;
int detMin_, detMax_;
edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> tok_htopo_;
edm::ESGetToken<CaloGeometry, CaloGeometryRecord> tok_geom_;
};
HcalGeometryDetIdTester::HcalGeometryDetIdTester(const edm::ParameterSet& iConfig) {
detMin_ = std::min(ndetMax_, std::max(iConfig.getParameter<int>("DetectorMin"), 1));
detMax_ = std::min(ndetMax_, std::max(iConfig.getParameter<int>("DetectorMax"), 1));
if (detMin_ > detMax_)
detMin_ = detMax_;
tok_htopo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
edm::LogVerbatim("HCalGeom") << "Study DetIds for SubDetId in the range " << detMin_ << ":" << detMax_;
}
void HcalGeometryDetIdTester::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
desc.add<int>("DetectorMin", 1);
desc.add<int>("DetectorMax", 4);
descriptions.add("hcalGeometryDetIdTester", desc);
}
void HcalGeometryDetIdTester::analyze(const edm::Event& /*iEvent*/, const edm::EventSetup& iSetup) {
const HcalTopology topology = iSetup.getData(tok_htopo_);
const CaloGeometry* geo = &iSetup.getData(tok_geom_);
const HcalGeometry* hcalGeom = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
std::string subdets[ndetMax_] = {"HB", "HE", "HO", "HF"};
HcalSubdetector subdetd[ndetMax_] = {HcalBarrel, HcalEndcap, HcalOuter, HcalForward};
int ietaMin[ndetMax_] = {1, 16, 1, 29};
int ietaMax[ndetMax_] = {16, 29, 15, 41};
int depthMin[ndetMax_] = {1, 1, 4, 1};
int depthMax[ndetMax_] = {4, 7, 4, 4};
for (int subd = (detMin_ - 1); subd < detMax_; ++subd) {
edm::LogVerbatim("HCalGeom") << "\n\nStudy Detector = Hcal SubDetector = " << subdets[subd]
<< "\n======================================\n";
int nall(0), nbad(0);
const std::vector<DetId>& ids = hcalGeom->getValidDetIds(DetId::Hcal, subdetd[subd]);
for (auto id : ids) {
++nall;
if (!(topology.valid(id))) {
++nbad;
edm::LogVerbatim("HCalGeom") << "Check " << HcalDetId(id) << " *****";
}
}
edm::LogVerbatim("HCalGeom") << "\n"
<< nbad << " bad out of " << nall
<< " detIds\n========================\n\nNow List All IDs\n================";
int k(0);
for (auto id : ids) {
edm::LogVerbatim("HCalGeom") << "[ " << std::setw(4) << k << "] " << HcalDetId(id);
++k;
}
int n(0);
edm::LogVerbatim("HCalGeom")
<< "\nNow List all IDs declared valid by Topology\n===========================================\n";
for (int ieta = ietaMin[subd]; ieta <= ietaMax[subd]; ++ieta) {
for (int depth = depthMin[subd]; depth <= depthMax[subd]; ++depth) {
HcalDetId id(subdetd[subd], ieta, 1, depth);
if (topology.validHcal(id)) {
edm::LogVerbatim("HCalGeom") << "[ " << std::setw(2) << n << "] " << id;
++n;
}
}
}
edm::LogVerbatim("HCalGeom") << "\nFinds a total of " << n << " IDs";
}
}
DEFINE_FWK_MODULE(HcalGeometryDetIdTester);
|