47 : g4Label_(ps.getParameter<
std::
string>(
"moduleLabel")),
48 caloHitSource_(ps.getParameter<
std::
string>(
"caloHitSource")),
49 nameSense_(ps.getParameter<
std::
string>(
"nameSense")),
50 nameDetector_(ps.getParameter<
std::
string>(
"nameDevice")),
51 verbosity_(ps.getParameter<
int>(
"Verbosity")),
64 desc.add<
int>(
"Verbosity", 0);
65 descriptions.
add(
"hgcalHitIdCheckEE",
desc);
70 edm::LogVerbatim(
"HitStudy") <<
"Run = " <<
e.id().run() <<
" Event = " <<
e.id().event();
74 const std::vector<DetId>& validIds =
geom->getValidDetIds();
78 bool getHits = (hitsCalo.
isValid());
79 uint32_t
nhits = (getHits) ? hitsCalo->size() : 0;
80 uint32_t good(0),
all(0);
82 edm::LogVerbatim(
"HitStudy") <<
"HGcalHitIdCheck: Input flags Hits " << getHits <<
" with " <<
nhits <<
" hits";
85 std::vector<PCaloHit>
hits;
86 hits.insert(
hits.end(), hitsCalo->begin(), hitsCalo->end());
91 if (
std::find(validIds.begin(), validIds.end(),
id) != validIds.end()) {
110 bool ok = (
std::find(validIds.begin(), validIds.end(),
DetId(hid2)) != validIds.end());
111 edm::LogVerbatim(
"HitStudy") <<
"Hit[" <<
all <<
":" << good <<
"]" << hid1 <<
" not valid ***** but "
112 << hid2 <<
" in list " <<
ok;
115 <<
"Hit[" <<
all <<
":" << good <<
"]" << std::hex <<
id.rawId() <<
std::dec <<
" not valid *****";
123 <<
" Invalid DetIds = " << (
all - good);