250 std::cout <<
"L1TCompare: analyze...." << std::endl;
274 edm::LogInfo(
"DataNotFound") <<
"can't find L1CaloEmCollection with label "
284 <<
"L1CaloRegionCollection with label "
294 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
297 if (!gctIsoEmCands.
isValid()) {
298 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
301 if (!gctNonIsoEmCands.
isValid()) {
302 std::cerr <<
"L1TGCT: could not find one of the classes?" << std::endl;
309 for ( L1GctEmCandCollection::const_iterator iem = gctIsoEmCands->begin();
310 iem != gctIsoEmCands->end(); ++iem) {
314 << iem->etaIndex() <<
"("
329 for (L1CaloEmCollection::const_iterator iem = em->begin();
330 iem != em->end(); ++iem) {
337 RctObject
h(
id.rctEta(),
id.rctPhi(), iem->rank());
338 if ( !iem->isolated() )
339 rcj_non_iso.push_back(
h);
341 rcj_iso.push_back(
h);
345 std::sort(rcj.begin(), rcj.end(), RctObjectComp());
346 std::sort(rcj_non_iso.begin(), rcj_non_iso.end(), RctObjectComp());
347 std::sort(rcj_iso.begin(), rcj_iso.end(), RctObjectComp());
349 for (RctObjectCollection::reverse_iterator ij = rcj_iso.rbegin();
350 ij != rcj_iso.rend() && ij != rcj_iso.rbegin()+8; ++ij) {
352 << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_
356 L1GctEmCandCollection::const_iterator lead_em = gctIsoEmCands->begin();
357 if ( !lead_em->empty() ) {
365 for ( L1GctEmCandCollection::const_iterator iem
366 = gctNonIsoEmCands->begin(); iem != gctNonIsoEmCands->end();
368 if ( ! iem->empty() )
369 std::cout <<
"GCT EM non: " << iem->rank()
380 for (RctObjectCollection::reverse_iterator ij = rcj_non_iso.rbegin();
381 ij != rcj_non_iso.rend() && ij != rcj_non_iso.rbegin()+8; ++ij) {
383 << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_
387 lead_em = gctNonIsoEmCands->begin();
388 if ( !lead_em->empty() ) {
390 rcj_non_iso.rbegin()->eta_);
392 rcj_non_iso.rbegin()->phi_);
394 rcj_non_iso.rbegin()->rank_);
403 <<
"can't find EcalTrigPrimCollection with label "
409 ieTP != eTP->end(); ieTP++) {
410 ecalobs.push_back(RctObject(ieTP->id().ieta(),
412 ieTP->compressedEt()));
414 std::sort(ecalobs.begin(), ecalobs.end(), RctObjectComp());
416 for (RctObjectCollection::reverse_iterator ij = ecalobs.rbegin();
417 ij != ecalobs.rend() && ij != ecalobs.rbegin()+8; ++ij) {
419 << ij->rank_ <<
", " << ij->eta_ <<
", " << ij->phi_
424 if ( rcj.rbegin()->rank_ > 4 ) {
426 ecalobs.rbegin()->eta_);
427 int e2 = (rcj.rbegin()->eta_&0x7UL)* ((rcj.rbegin()->eta_&0x8UL)?1:-1);
431 ecalobs.rbegin()->rank_);
434 int seta = rcj.rbegin()->eta_;
435 seta = (seta&0x7UL)*(seta&0x8?-1:1);
437 << rcj.rbegin()->eta_ <<
" "
438 << rcj.rbegin()->phi_ <<
" "
439 << rcj.rbegin()->rank_ <<
" "
440 << (++rcj.rbegin())->rank_<<
" "
441 << ecalobs.rbegin()->eta_ <<
" "
442 << ecalobs.rbegin()->phi_ <<
" "
443 << ecalobs.rbegin()->rank_ <<
" "
444 << (++ecalobs.rbegin())->rank_<<
" "
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmCandsToken_
edm::EDGetTokenT< L1CaloRegionCollection > rctSourceRctEmRgn_token_
MonitorElement * rctGctLeadingIsoEmRank_
MonitorElement * rctGctLeadingIsoEmEta_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * ecalTpgRctLeadingEmEta2_
std::vector< EcalTriggerPrimitiveDigi >::const_iterator const_iterator
MonitorElement * ecalTpgRctLeadingEmRank_
MonitorElement * rctGctLeadingNonIsoEmEta_
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
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_
MonitorElement * ecalTpgRctLeadingEmPhi_