27 : m_rctEmCands(consumes(iConfig.getParameter<
edm::
InputTag>(
"rctEmCandsLabel"))),
28 m_rctRegions(consumes(iConfig.getParameter<
edm::
InputTag>(
"rctRegionsLabel"))) {
33 h_emRank =
fs->make<TH1F>(
"emRank",
"emRank", 64, 0., 64.);
34 h_emIeta =
fs->make<TH1F>(
"emOccupancyIeta",
"emOccupancyIeta", 22, 0., 22.);
35 h_emIphi =
fs->make<TH1F>(
"emOccupancyIphi",
"emOccupancyIphi", 18, 0., 18.);
36 h_emIsoOccIetaIphi =
fs->make<TH2F>(
"emIsoOccupancy2D",
"emIsoOccupancy2D", 22, 0., 22., 18, 0., 18.);
39 h_regionSum =
fs->make<TH1F>(
"regionSum",
"regionSum", 100, 0., 100.);
40 h_regionSumIetaIphi =
fs->make<TH2F>(
"regionSumEtWeighted2D",
"regionSumEtWeighted2D", 22, 0., 22., 18, 0., 18.);
41 h_regionOccIetaIphi =
fs->make<TH2F>(
"regionOccupancy2D",
"regionOccupancy2D", 22, 0., 22., 18, 0., 18.);
61 for (
auto em = rctEmCands->begin(); em != rctEmCands->end(); em++) {
62 if ((*em).rank() > 0) {
64 h_emIeta->Fill((*em).regionId().ieta());
65 h_emIphi->Fill((*em).regionId().iphi());
66 if ((*em).isolated()) {
74 for (
auto rgn = rctRegions->begin(); rgn != rctRegions->end(); rgn++) {
75 if ((*rgn).et() > 0) {
static const std::string kSharedResource
void analyze(const edm::Event &, const edm::EventSetup &) override
TH2F * h_regionOccIetaIphi
~L1RCTRelValAnalyzer() override
edm::EDGetTokenT< L1CaloRegionCollection > m_rctRegions
L1RCTRelValAnalyzer(const edm::ParameterSet &)
TH2F * h_emNonIsoOccIetaIphi
edm::EDGetTokenT< L1CaloEmCollection > m_rctEmCands
TH2F * h_emIsoOccIetaIphi
TH2F * h_regionSumIetaIphi