|
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
|
#include "CondFormats/DataRecord/interface/SiStripNoisesRcd.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/SiStripDetId/interface/SiStripDetId.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include <iostream>
#include <exception>
class SiStripNoisesGetAllChecker : public edm::one::EDAnalyzer<> {
public:
explicit SiStripNoisesGetAllChecker(const edm::ParameterSet&);
~SiStripNoisesGetAllChecker() override = default;
void analyze(const edm::Event&, const edm::EventSetup&) override;
private:
const edm::ESGetToken<SiStripNoises, SiStripNoisesRcd> noisesToken_;
const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
void checkModuleNoise(const SiStripNoises&, const uint32_t detID, uint16_t maxNStrips);
};
SiStripNoisesGetAllChecker::SiStripNoisesGetAllChecker(const edm::ParameterSet&)
: noisesToken_(esConsumes()), tkGeomToken_(esConsumes()) {}
void SiStripNoisesGetAllChecker::analyze(const edm::Event&, const edm::EventSetup& iSetup) {
const auto& siStripNoises = iSetup.getData(noisesToken_);
const auto& tkGeom = &iSetup.getData(tkGeomToken_);
const auto& tkDets = tkGeom->dets();
edm::LogInfo("SiStripNoisesGetAllChecker") << "Starting to loop over all SiStrip modules...";
// Get all DetIDs associated with SiStripNoises
std::vector<uint32_t> detIDs;
siStripNoises.getDetIds(detIDs);
size_t exceptionCounts{0};
for (const auto& detID : detIDs) {
uint16_t maxNStrips{0};
auto det = std::find_if(tkDets.begin(), tkDets.end(), [detID](auto& elem) -> bool {
return (elem->geographicalId().rawId() == detID);
});
const StripTopology& p = dynamic_cast<const StripGeomDetUnit*>(*det)->specificTopology();
maxNStrips = p.nstrips();
try {
checkModuleNoise(siStripNoises, detID, maxNStrips);
} catch (const std::exception& e) {
// Increment the exception counter if checkModuleNoise itself throws an exception
edm::LogError("SiStripNoisesGetAllChecker")
<< "Exception in checkModuleNoise for detID " << detID << ": " << e.what();
++exceptionCounts;
} catch (...) {
edm::LogError("SiStripNoisesGetAllChecker") << "Unknown exception in checkModuleNoise for detID " << detID;
++exceptionCounts;
}
}
std::ostringstream message;
// Define the box width
const int boxWidth = 50;
message << "\n"
<< std::string(boxWidth, '*') << "\n"
<< "* " << std::setw(boxWidth - 4) << std::left << "SiStripNoisesGetAllChecker Summary"
<< " *\n"
<< std::string(boxWidth, '*') << "\n"
<< "* " << std::setw(boxWidth - 4) << std::left
<< ("Completed loop over " + std::to_string(detIDs.size()) + " SiStrip modules.") << " *\n"
<< "* " << std::setw(boxWidth - 4) << std::left
<< ("Encountered " + std::to_string(exceptionCounts) + " exceptions.") << " *\n"
<< std::string(boxWidth, '*');
edm::LogSystem("SiStripNoisesGetAllChecker") << message.str();
}
void SiStripNoisesGetAllChecker::checkModuleNoise(const SiStripNoises& siStripNoises,
const uint32_t detID,
uint16_t maxNStrips) {
try {
SiStripNoises::Range detNoiseRange = siStripNoises.getRange(detID);
std::vector<float> noises;
noises.resize(maxNStrips);
siStripNoises.allNoises(noises, detNoiseRange);
edm::LogInfo("SiStripNoisesGetAllChecker") << "Successfully processed detID: " << detID;
} catch (const std::exception& e) {
edm::LogError("SiStripNoisesGetAllChecker") << "Exception caught for detID " << detID << ": " << e.what();
throw; // Re-throw the exception to allow the outer loop to handle it
} catch (...) {
edm::LogError("SiStripNoisesGetAllChecker") << "Unknown exception caught for detID " << detID;
throw; // Re-throw the unknown exception
}
}
// Define this as a plug-in
DEFINE_FWK_MODULE(SiStripNoisesGetAllChecker);
|