194 std::cout <<
"L1TCompare: analyze...." << std::endl;
234 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
237 if (!gctIsoEmCands.
isValid()) {
238 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
241 if (!gctNonIsoEmCands.
isValid()) {
242 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
248 for (L1GctEmCandCollection::const_iterator iem = gctIsoEmCands->begin(); iem != gctIsoEmCands->end(); ++iem) {
250 std::cout <<
"GCT EM: " << iem->rank() <<
", " << iem->etaIndex()
253 <<
"), " << iem->phiIndex() << std::endl;
263 for (L1CaloEmCollection::const_iterator iem = em->begin(); iem != em->end(); ++iem) {
269 RctObject
h(
id.rctEta(),
id.rctPhi(), iem->rank());
270 if (!iem->isolated())
271 rcj_non_iso.push_back(
h);
273 rcj_iso.push_back(
h);
278 std::sort(rcj_non_iso.begin(), rcj_non_iso.end(),
rctObjectComp);
281 for (RctObjectCollection::reverse_iterator ij = rcj_iso.rbegin();
282 ij != rcj_iso.rend() && ij != rcj_iso.rbegin() + 8;
284 std::cout <<
"RCT cj: " << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_ << std::endl;
287 L1GctEmCandCollection::const_iterator lead_em = gctIsoEmCands->begin();
288 if (!lead_em->empty()) {
296 for (L1GctEmCandCollection::const_iterator iem = gctNonIsoEmCands->begin(); iem != gctNonIsoEmCands->end(); ++iem) {
298 std::cout <<
"GCT EM non: " << iem->rank() <<
", "
302 <<
", " << iem->phiIndex() << std::endl;
306 for (RctObjectCollection::reverse_iterator ij = rcj_non_iso.rbegin();
307 ij != rcj_non_iso.rend() && ij != rcj_non_iso.rbegin() + 8;
309 std::cout <<
"RCT cj non: " << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_ << std::endl;
312 lead_em = gctNonIsoEmCands->begin();
313 if (!lead_em->empty()) {
329 ecalobs.push_back(RctObject(ieTP->id().ieta(), ieTP->id().iphi(), ieTP->compressedEt()));
333 for (RctObjectCollection::reverse_iterator ij = ecalobs.rbegin();
334 ij != ecalobs.rend() && ij != ecalobs.rbegin() + 8;
336 std::cout <<
"ECAL cj : " << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_ << std::endl;
340 if (rcj.rbegin()->rank_ > 4) {
342 int e2 = (rcj.rbegin()->eta_ & 0x7UL) * ((rcj.rbegin()->eta_ & 0x8UL) ? 1 : -1);
348 int seta = rcj.rbegin()->eta_;
349 seta = (seta & 0x7UL) * (seta & 0x8 ? -1 : 1);
350 std::cout <<
"ZZ: " << rcj.rbegin()->eta_ <<
" " << rcj.rbegin()->phi_ <<
" " << rcj.rbegin()->rank_ <<
" "
351 << (++rcj.rbegin())->rank_ <<
" " << ecalobs.rbegin()->eta_ <<
" " << ecalobs.rbegin()->phi_ <<
" "
352 << ecalobs.rbegin()->rank_ <<
" " << (++ecalobs.rbegin())->rank_ <<
" " << seta <<
" " << std::endl;
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmCandsToken_
edm::EDGetTokenT< L1CaloRegionCollection > rctSourceRctEmRgn_token_
MonitorElement * rctGctLeadingIsoEmRank_
uint16_t *__restrict__ id
MonitorElement * rctGctLeadingIsoEmEta_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * ecalTpgRctLeadingEmEta2_
std::vector< T >::const_iterator const_iterator
static bool rctObjectComp(const RctObject &a, const RctObject &b)
MonitorElement * ecalTpgRctLeadingEmRank_
MonitorElement * rctGctLeadingNonIsoEmEta_
Log< level::Info, false > LogInfo
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsToken_
edm::EDGetTokenT< L1CaloEmCollection > rctSourceEm_token_
MonitorElement * rctGctLeadingIsoEmPhi_
std::vector< L1TCompare::RctObject > RctObjectCollection
MonitorElement * ecalTpgRctLeadingEmEta_
edm::InputTag ecalTpgSource_
MonitorElement * rctGctLeadingNonIsoEmPhi_
edm::EDGetTokenT< L1GctEmCandCollection > gctIsoEmCandsToken_
edm::EDGetTokenT< EcalTrigPrimDigiCollection > ecalTpgSource_token_
MonitorElement * rctGctLeadingNonIsoEmRank_
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
MonitorElement * ecalTpgRctLeadingEmPhi_