#include <DQM/L1TMonitor/interface/L1TGCT.h>
Definition at line 98 of file L1TGCT.h.
L1TGCT::L1TGCT | ( | const edm::ParameterSet & | ps | ) |
Definition at line 161 of file L1TGCT.cc.
References dbe, lat::endl(), edm::ParameterSet::getUntrackedParameter(), NULL, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), and verbose_.
00161 : 00162 gctCenJetsSource_(ps.getParameter<edm::InputTag>("gctCentralJetsSource")), 00163 gctForJetsSource_(ps.getParameter<edm::InputTag>("gctForwardJetsSource")), 00164 gctTauJetsSource_(ps.getParameter<edm::InputTag>("gctTauJetsSource")), 00165 gctEnergySumsSource_(ps.getParameter<edm::InputTag>("gctEnergySumsSource")), 00166 gctIsoEmSource_(ps.getParameter<edm::InputTag>("gctIsoEmSource")), 00167 gctNonIsoEmSource_(ps.getParameter<edm::InputTag>("gctNonIsoEmSource")) 00168 { 00169 00170 // verbosity switch 00171 verbose_ = ps.getUntrackedParameter < bool > ("verbose", false); 00172 00173 if (verbose_) 00174 edm::LogInfo("L1TGCT") << "L1TGCT: constructor...." << std::endl; 00175 00176 00177 dbe = NULL; 00178 if (ps.getUntrackedParameter < bool > ("DQMStore", false)) { 00179 dbe = edm::Service < DQMStore > ().operator->(); 00180 dbe->setVerbose(0); 00181 } 00182 00183 outputFile_ = ps.getUntrackedParameter < std::string > ("outputFile", ""); 00184 if (outputFile_.size() != 0) { 00185 edm::LogInfo("L1TGCT") << "L1T Monitoring histograms will be saved to " 00186 << outputFile_ << std::endl; 00187 } 00188 00189 bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false); 00190 if(disable){ 00191 outputFile_=""; 00192 } 00193 00194 00195 if (dbe != NULL) { 00196 dbe->setCurrentFolder("L1T/L1TGCT"); 00197 } 00198 00199 }
void L1TGCT::analyze | ( | const edm::Event & | e, | |
const edm::EventSetup & | c | |||
) | [protected, virtual] |
Implements edm::EDAnalyzer.
Definition at line 373 of file L1TGCT.cc.
References lat::endl(), MonitorElement::Fill(), gctCenJetsSource_, gctEnergySumsSource_, gctForJetsSource_, gctIsoEmSource_, gctNonIsoEmSource_, gctTauJetsSource_, edm::Event::getByLabel(), edm::Handle< T >::isValid(), l1GctCenJetsEtEtaPhi_, l1GctCenJetsOccEta_, l1GctCenJetsOccEtaPhi_, l1GctCenJetsOccPhi_, l1GctCenJetsRank_, l1GctCenJetsRankCand0_, l1GctCenJetsRankCand1_, l1GctCenJetsRankCand2_, l1GctCenJetsRankCand3_, l1GctCenJetsRankDiff01_, l1GctCenJetsRankDiff12_, l1GctCenJetsRankDiff23_, l1GctEtHad_, l1GctEtMiss_, l1GctEtMissPhi_, l1GctEtTotal_, l1GctForJetsEtEtaPhi_, l1GctForJetsOccEta_, l1GctForJetsOccEtaPhi_, l1GctForJetsOccPhi_, l1GctForJetsRank_, l1GctForJetsRankCand0_, l1GctForJetsRankCand1_, l1GctForJetsRankCand2_, l1GctForJetsRankCand3_, l1GctForJetsRankDiff01_, l1GctForJetsRankDiff12_, l1GctForJetsRankDiff23_, l1GctHFRing0ETSumNegEta_, l1GctHFRing0ETSumPosEta_, l1GctHFRing0PosEtaNegEta_, l1GctHFRing0TowerCountNegEta_, l1GctHFRing0TowerCountPosEta_, l1GctHFRing0TowerCountPosEtaNegEta_, l1GctHFRing1ETSumNegEta_, l1GctHFRing1ETSumPosEta_, l1GctHFRing1PosEtaNegEta_, l1GctHFRing1TowerCountNegEta_, l1GctHFRing1TowerCountPosEta_, l1GctHFRing1TowerCountPosEtaNegEta_, l1GctHFRingRatioNegEta_, l1GctHFRingRatioPosEta_, l1GctIsoEmOccEta_, l1GctIsoEmOccEtaPhi_, l1GctIsoEmOccPhi_, l1GctIsoEmRank_, l1GctIsoEmRankCand0_, l1GctIsoEmRankCand1_, l1GctIsoEmRankCand2_, l1GctIsoEmRankCand3_, l1GctIsoEmRankDiff01_, l1GctIsoEmRankDiff12_, l1GctIsoEmRankDiff23_, l1GctIsoEmRankEtaPhi_, l1GctNonIsoEmOccEta_, l1GctNonIsoEmOccEtaPhi_, l1GctNonIsoEmOccPhi_, l1GctNonIsoEmRank_, l1GctNonIsoEmRankCand0_, l1GctNonIsoEmRankCand1_, l1GctNonIsoEmRankCand2_, l1GctNonIsoEmRankCand3_, l1GctNonIsoEmRankDiff01_, l1GctNonIsoEmRankDiff12_, l1GctNonIsoEmRankDiff23_, l1GctNonIsoEmRankEtaPhi_, l1GctTauJetsEtEtaPhi_, l1GctTauJetsOccEta_, l1GctTauJetsOccEtaPhi_, l1GctTauJetsOccPhi_, l1GctTauJetsRank_, l1GctTauJetsRankCand0_, l1GctTauJetsRankCand1_, l1GctTauJetsRankCand2_, l1GctTauJetsRankCand3_, l1GctTauJetsRankDiff01_, l1GctTauJetsRankDiff12_, l1GctTauJetsRankDiff23_, nev_, and verbose_.
00374 { 00375 nev_++; 00376 if (verbose_) { 00377 edm::LogInfo("L1TGCT") << "L1TGCT: analyze...." << std::endl; 00378 } 00379 00380 // Get all the collections 00381 edm::Handle < L1GctEmCandCollection > l1IsoEm; 00382 edm::Handle < L1GctEmCandCollection > l1NonIsoEm; 00383 edm::Handle < L1GctJetCandCollection > l1CenJets; 00384 edm::Handle < L1GctJetCandCollection > l1ForJets; 00385 edm::Handle < L1GctJetCandCollection > l1TauJets; 00386 edm::Handle < L1GctHFRingEtSumsCollection > l1HFSums; 00387 edm::Handle < L1GctHFBitCountsCollection > l1HFCounts; 00388 edm::Handle < L1GctEtMissCollection > l1EtMiss; 00389 edm::Handle < L1GctEtHadCollection > l1EtHad; 00390 edm::Handle < L1GctEtTotalCollection > l1EtTotal; 00391 00392 e.getByLabel(gctIsoEmSource_, l1IsoEm); 00393 e.getByLabel(gctNonIsoEmSource_, l1NonIsoEm); 00394 e.getByLabel(gctCenJetsSource_, l1CenJets); 00395 e.getByLabel(gctForJetsSource_, l1ForJets); 00396 e.getByLabel(gctTauJetsSource_, l1TauJets); 00397 e.getByLabel(gctEnergySumsSource_, l1HFSums); 00398 e.getByLabel(gctEnergySumsSource_, l1HFCounts); 00399 e.getByLabel(gctEnergySumsSource_, l1EtMiss); 00400 e.getByLabel(gctEnergySumsSource_, l1EtHad); 00401 e.getByLabel(gctEnergySumsSource_, l1EtTotal); 00402 00403 // Fill histograms 00404 00405 // Central jets 00406 if (l1CenJets.isValid()) { 00407 for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin();cj != l1CenJets->end(); cj++) { 00408 l1GctCenJetsRank_->Fill(cj->rank()); 00409 // only plot eta and phi maps for non-zero candidates 00410 if (cj->rank()) { 00411 l1GctCenJetsEtEtaPhi_->Fill(cj->regionId().iphi(),cj->regionId().ieta(),cj->rank()); 00412 l1GctCenJetsOccEtaPhi_->Fill(cj->regionId().iphi(),cj->regionId().ieta()); 00413 l1GctCenJetsOccEta_->Fill(cj->regionId().ieta()); 00414 l1GctCenJetsOccPhi_->Fill(cj->regionId().iphi()); 00415 } 00416 } 00417 if ( l1CenJets->size()==4){ 00418 // Rank for each candidate 00419 l1GctCenJetsRankCand0_->Fill((*l1CenJets).at(0).rank()); 00420 l1GctCenJetsRankCand1_->Fill((*l1CenJets).at(1).rank()); 00421 l1GctCenJetsRankCand2_->Fill((*l1CenJets).at(2).rank()); 00422 l1GctCenJetsRankCand3_->Fill((*l1CenJets).at(3).rank()); 00423 // Differences between candidate ranks 00424 l1GctCenJetsRankDiff01_->Fill((*l1CenJets).at(0).rank()-(*l1CenJets).at(1).rank()); 00425 l1GctCenJetsRankDiff12_->Fill((*l1CenJets).at(1).rank()-(*l1CenJets).at(2).rank()); 00426 l1GctCenJetsRankDiff23_->Fill((*l1CenJets).at(2).rank()-(*l1CenJets).at(3).rank()); 00427 } 00428 } else { 00429 edm::LogWarning("DataNotFound") << " Could not find l1CenJets label was " << gctCenJetsSource_ ; 00430 } 00431 00432 00433 // Forward jets 00434 if (l1ForJets.isValid()) { 00435 for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin(); fj != l1ForJets->end(); fj++) { 00436 l1GctForJetsRank_->Fill(fj->rank()); 00437 // only plot eta and phi maps for non-zero candidates 00438 if (fj->rank()) { 00439 l1GctForJetsEtEtaPhi_->Fill(fj->regionId().iphi(),fj->regionId().ieta(),fj->rank()); 00440 l1GctForJetsOccEtaPhi_->Fill(fj->regionId().iphi(),fj->regionId().ieta()); 00441 l1GctForJetsOccEta_->Fill(fj->regionId().ieta()); 00442 l1GctForJetsOccPhi_->Fill(fj->regionId().iphi()); 00443 } 00444 } 00445 if ( l1ForJets->size()==4){ 00446 // Rank for each candidate 00447 l1GctForJetsRankCand0_->Fill((*l1ForJets).at(0).rank()); 00448 l1GctForJetsRankCand1_->Fill((*l1ForJets).at(1).rank()); 00449 l1GctForJetsRankCand2_->Fill((*l1ForJets).at(2).rank()); 00450 l1GctForJetsRankCand3_->Fill((*l1ForJets).at(3).rank()); 00451 // Differences between candidate ranks 00452 l1GctForJetsRankDiff01_->Fill((*l1ForJets).at(0).rank()-(*l1ForJets).at(1).rank()); 00453 l1GctForJetsRankDiff12_->Fill((*l1ForJets).at(1).rank()-(*l1ForJets).at(2).rank()); 00454 l1GctForJetsRankDiff23_->Fill((*l1ForJets).at(2).rank()-(*l1ForJets).at(3).rank()); 00455 } 00456 } else { 00457 edm::LogWarning("DataNotFound") << " Could not find l1ForJets label was " << gctForJetsSource_ ; 00458 } 00459 00460 // Tau jets 00461 if (l1TauJets.isValid()) { 00462 for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin(); tj != l1TauJets->end(); tj++) { 00463 l1GctTauJetsRank_->Fill(tj->rank()); 00464 // only plot eta and phi maps for non-zero candidates 00465 if (tj->rank()) { 00466 l1GctTauJetsEtEtaPhi_->Fill(tj->regionId().iphi(),tj->regionId().ieta(),tj->rank()); 00467 l1GctTauJetsOccEtaPhi_->Fill(tj->regionId().iphi(),tj->regionId().ieta()); 00468 l1GctTauJetsOccEta_->Fill(tj->regionId().ieta()); 00469 l1GctTauJetsOccPhi_->Fill(tj->regionId().iphi()); 00470 } 00471 } 00472 if (l1TauJets->size()==4){ 00473 // Rank for each candidate 00474 l1GctTauJetsRankCand0_->Fill((*l1TauJets).at(0).rank()); 00475 l1GctTauJetsRankCand1_->Fill((*l1TauJets).at(1).rank()); 00476 l1GctTauJetsRankCand2_->Fill((*l1TauJets).at(2).rank()); 00477 l1GctTauJetsRankCand3_->Fill((*l1TauJets).at(3).rank()); 00478 // Differences between candidate ranks 00479 l1GctTauJetsRankDiff01_->Fill((*l1TauJets).at(0).rank()-(*l1TauJets).at(1).rank()); 00480 l1GctTauJetsRankDiff12_->Fill((*l1TauJets).at(1).rank()-(*l1TauJets).at(2).rank()); 00481 l1GctTauJetsRankDiff23_->Fill((*l1TauJets).at(2).rank()-(*l1TauJets).at(3).rank()); 00482 } 00483 } else { 00484 edm::LogWarning("DataNotFound") << " Could not find l1TauJets label was " << gctTauJetsSource_ ; 00485 } 00486 00487 00488 // Missing ET 00489 if (l1EtMiss.isValid()) { 00490 if (l1EtMiss->size()) { 00491 l1GctEtMiss_->Fill(l1EtMiss->at(0).et()); 00492 l1GctEtMissPhi_->Fill(l1EtMiss->at(0).phi()); 00493 } 00494 } else { 00495 edm::LogWarning("DataNotFound") << " Could not find l1EtMiss label was " << gctEnergySumsSource_ ; 00496 } 00497 00498 // HT 00499 if (l1EtHad.isValid()) { 00500 if (l1EtHad->size()) { 00501 l1GctEtHad_->Fill(l1EtHad->at(0).et()); 00502 } 00503 } else { 00504 edm::LogWarning("DataNotFound") << " Could not find l1EtHad label was " << gctEnergySumsSource_ ; 00505 } 00506 00507 // Total ET 00508 if (l1EtTotal.isValid()) { 00509 if (l1EtTotal->size()) { 00510 l1GctEtTotal_->Fill(l1EtTotal->at(0).et()); 00511 } 00512 } else { 00513 edm::LogWarning("DataNotFound") << " Could not find l1EtTotal label was " << gctEnergySumsSource_ ; 00514 } 00515 00516 //HF Ring Et Sums 00517 if (l1HFSums.isValid()) { 00518 for (L1GctHFRingEtSumsCollection::const_iterator hfs=l1HFSums->begin(); hfs!=l1HFSums->end(); hfs++){ 00519 // Individual ring Et sums 00520 l1GctHFRing0ETSumPosEta_->Fill(hfs->etSum(0)); 00521 l1GctHFRing0ETSumNegEta_->Fill(hfs->etSum(1)); 00522 l1GctHFRing1ETSumPosEta_->Fill(hfs->etSum(2)); 00523 l1GctHFRing1ETSumNegEta_->Fill(hfs->etSum(3)); 00524 // Ratio of ring Et sums 00525 if (hfs->etSum(2)!=0) l1GctHFRingRatioPosEta_->Fill((hfs->etSum(0))/(hfs->etSum(2))); 00526 if (hfs->etSum(3)!=0) l1GctHFRingRatioNegEta_->Fill((hfs->etSum(1))/(hfs->etSum(3))); 00527 // Correlate positive and neagative eta 00528 l1GctHFRing0PosEtaNegEta_->Fill(hfs->etSum(0),hfs->etSum(1)); 00529 l1GctHFRing1PosEtaNegEta_->Fill(hfs->etSum(2),hfs->etSum(3)); 00530 } 00531 } else { 00532 edm::LogWarning("DataNotFound") << " Could not find l1HFSums label was " << gctEnergySumsSource_ ; 00533 } 00534 00535 // HF Ring Counts 00536 if (l1HFCounts.isValid()) { 00537 for (L1GctHFBitCountsCollection::const_iterator hfc=l1HFCounts->begin(); hfc!=l1HFCounts->end(); hfc++){ 00538 // Individual ring counts 00539 l1GctHFRing0TowerCountPosEta_->Fill(hfc->bitCount(0)); 00540 l1GctHFRing0TowerCountNegEta_->Fill(hfc->bitCount(1)); 00541 l1GctHFRing1TowerCountPosEta_->Fill(hfc->bitCount(2)); 00542 l1GctHFRing1TowerCountNegEta_->Fill(hfc->bitCount(3)); 00543 // Correlate positive and negative eta 00544 l1GctHFRing0TowerCountPosEtaNegEta_->Fill(hfc->bitCount(0),hfc->bitCount(1)); 00545 l1GctHFRing1TowerCountPosEtaNegEta_->Fill(hfc->bitCount(2),hfc->bitCount(3)); 00546 } 00547 } else { 00548 edm::LogWarning("DataNotFound") << " Could not find l1HFCounts label was " << gctEnergySumsSource_ ; 00549 } 00550 00551 // Isolated EM 00552 if (l1IsoEm.isValid()) { 00553 for (L1GctEmCandCollection::const_iterator ie=l1IsoEm->begin(); ie!=l1IsoEm->end(); ie++) { 00554 l1GctIsoEmRank_->Fill(ie->rank()); 00555 // only plot eta and phi maps for non-zero candidates 00556 if (ie->rank()){ 00557 l1GctIsoEmRankEtaPhi_->Fill(ie->regionId().iphi(),ie->regionId().ieta(),ie->rank()); 00558 l1GctIsoEmOccEtaPhi_->Fill(ie->regionId().iphi(),ie->regionId().ieta()); 00559 l1GctIsoEmOccEta_->Fill(ie->regionId().ieta()); 00560 l1GctIsoEmOccPhi_->Fill(ie->regionId().iphi()); 00561 } 00562 } 00563 if (l1IsoEm->size()==4){ 00564 // Rank for each candidate 00565 l1GctIsoEmRankCand0_->Fill((*l1IsoEm).at(0).rank()); 00566 l1GctIsoEmRankCand1_->Fill((*l1IsoEm).at(1).rank()); 00567 l1GctIsoEmRankCand2_->Fill((*l1IsoEm).at(2).rank()); 00568 l1GctIsoEmRankCand3_->Fill((*l1IsoEm).at(3).rank()); 00569 // Differences between candidate ranks 00570 l1GctIsoEmRankDiff01_->Fill((*l1IsoEm).at(0).rank()-(*l1IsoEm).at(1).rank()); 00571 l1GctIsoEmRankDiff12_->Fill((*l1IsoEm).at(1).rank()-(*l1IsoEm).at(2).rank()); 00572 l1GctIsoEmRankDiff23_->Fill((*l1IsoEm).at(2).rank()-(*l1IsoEm).at(3).rank()); 00573 } 00574 } else { 00575 edm::LogWarning("DataNotFound") << " Could not find l1IsoEm label was " << gctIsoEmSource_ ; 00576 } 00577 00578 // Non-isolated EM 00579 if (l1NonIsoEm.isValid()) { 00580 for (L1GctEmCandCollection::const_iterator ne=l1NonIsoEm->begin(); ne!=l1NonIsoEm->end(); ne++) { 00581 l1GctNonIsoEmRank_->Fill(ne->rank()); 00582 // only plot eta and phi maps for non-zero candidates 00583 if (ne->rank()){ 00584 l1GctNonIsoEmRankEtaPhi_->Fill(ne->regionId().iphi(),ne->regionId().ieta(),ne->rank()); 00585 l1GctNonIsoEmOccEtaPhi_->Fill(ne->regionId().iphi(),ne->regionId().ieta()); 00586 l1GctNonIsoEmOccEta_->Fill(ne->regionId().ieta()); 00587 l1GctNonIsoEmOccPhi_->Fill(ne->regionId().iphi()); 00588 } 00589 } 00590 if (l1NonIsoEm->size()==4){ 00591 // Rank for each candidate 00592 l1GctNonIsoEmRankCand0_->Fill((*l1NonIsoEm).at(0).rank()); 00593 l1GctNonIsoEmRankCand1_->Fill((*l1NonIsoEm).at(1).rank()); 00594 l1GctNonIsoEmRankCand2_->Fill((*l1NonIsoEm).at(2).rank()); 00595 l1GctNonIsoEmRankCand3_->Fill((*l1NonIsoEm).at(3).rank()); 00596 // Differences between candidate ranks 00597 l1GctNonIsoEmRankDiff01_->Fill((*l1NonIsoEm).at(0).rank()-(*l1NonIsoEm).at(1).rank()); 00598 l1GctNonIsoEmRankDiff12_->Fill((*l1NonIsoEm).at(1).rank()-(*l1NonIsoEm).at(2).rank()); 00599 l1GctNonIsoEmRankDiff23_->Fill((*l1NonIsoEm).at(2).rank()-(*l1NonIsoEm).at(3).rank()); 00600 } 00601 } else { 00602 edm::LogWarning("DataNotFound") << " Could not find l1NonIsoEm label was " << gctNonIsoEmSource_ ; 00603 } 00604 }
void L1TGCT::beginJob | ( | const edm::EventSetup & | c | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 205 of file L1TGCT.cc.
References DQMStore::book1D(), DQMStore::book2D(), dbe, ETABINS, ETAMAX, ETAMIN, l1GctCenJetsEtEtaPhi_, l1GctCenJetsOccEta_, l1GctCenJetsOccEtaPhi_, l1GctCenJetsOccPhi_, l1GctCenJetsRank_, l1GctCenJetsRankCand0_, l1GctCenJetsRankCand1_, l1GctCenJetsRankCand2_, l1GctCenJetsRankCand3_, l1GctCenJetsRankDiff01_, l1GctCenJetsRankDiff12_, l1GctCenJetsRankDiff23_, l1GctEtHad_, l1GctEtMiss_, l1GctEtMissPhi_, l1GctEtTotal_, l1GctForJetsEtEtaPhi_, l1GctForJetsOccEta_, l1GctForJetsOccEtaPhi_, l1GctForJetsOccPhi_, l1GctForJetsRank_, l1GctForJetsRankCand0_, l1GctForJetsRankCand1_, l1GctForJetsRankCand2_, l1GctForJetsRankCand3_, l1GctForJetsRankDiff01_, l1GctForJetsRankDiff12_, l1GctForJetsRankDiff23_, l1GctHFRing0ETSumNegEta_, l1GctHFRing0ETSumPosEta_, l1GctHFRing0PosEtaNegEta_, l1GctHFRing0TowerCountNegEta_, l1GctHFRing0TowerCountPosEta_, l1GctHFRing0TowerCountPosEtaNegEta_, l1GctHFRing1ETSumNegEta_, l1GctHFRing1ETSumPosEta_, l1GctHFRing1PosEtaNegEta_, l1GctHFRing1TowerCountNegEta_, l1GctHFRing1TowerCountPosEta_, l1GctHFRing1TowerCountPosEtaNegEta_, l1GctHFRingRatioNegEta_, l1GctHFRingRatioPosEta_, l1GctIsoEmOccEta_, l1GctIsoEmOccEtaPhi_, l1GctIsoEmOccPhi_, l1GctIsoEmRank_, l1GctIsoEmRankCand0_, l1GctIsoEmRankCand1_, l1GctIsoEmRankCand2_, l1GctIsoEmRankCand3_, l1GctIsoEmRankDiff01_, l1GctIsoEmRankDiff12_, l1GctIsoEmRankDiff23_, l1GctIsoEmRankEtaPhi_, l1GctNonIsoEmOccEta_, l1GctNonIsoEmOccEtaPhi_, l1GctNonIsoEmOccPhi_, l1GctNonIsoEmRank_, l1GctNonIsoEmRankCand0_, l1GctNonIsoEmRankCand1_, l1GctNonIsoEmRankCand2_, l1GctNonIsoEmRankCand3_, l1GctNonIsoEmRankDiff01_, l1GctNonIsoEmRankDiff12_, l1GctNonIsoEmRankDiff23_, l1GctNonIsoEmRankEtaPhi_, l1GctTauJetsEtEtaPhi_, l1GctTauJetsOccEta_, l1GctTauJetsOccEtaPhi_, l1GctTauJetsOccPhi_, l1GctTauJetsRank_, l1GctTauJetsRankCand0_, l1GctTauJetsRankCand1_, l1GctTauJetsRankCand2_, l1GctTauJetsRankCand3_, l1GctTauJetsRankDiff01_, l1GctTauJetsRankDiff12_, l1GctTauJetsRankDiff23_, METPHIBINS, nev_, PHIBINS, PHIMAX, PHIMIN, R12BINS, R12MAX, R12MIN, R3BINS, R3MAX, R3MIN, R5BINS, R5MAX, R5MIN, R6BINS, R6MAX, R6MIN, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
00206 { 00207 00208 nev_ = 0; 00209 00210 // get hold of back-end interface 00211 DQMStore *dbe = 0; 00212 dbe = edm::Service < DQMStore > ().operator->(); 00213 00214 if (dbe) { 00215 dbe->setCurrentFolder("L1T/L1TGCT"); 00216 dbe->rmdir("L1T/L1TGCT"); 00217 } 00218 00219 00220 if (dbe) { 00221 00222 dbe->setCurrentFolder("L1T/L1TGCT"); 00223 00224 // GCT hardware quantities for experts 00225 l1GctCenJetsEtEtaPhi_ = dbe->book2D("CenJetsEtEtaPhi", "CENTRAL JET RANK", 00226 PHIBINS, PHIMIN, PHIMAX, 00227 ETABINS, ETAMIN, ETAMAX); 00228 l1GctForJetsEtEtaPhi_ = dbe->book2D("ForJetsEtEtaPhi", "FORWARD JET RANK", 00229 PHIBINS, PHIMIN, PHIMAX, 00230 ETABINS, ETAMIN, ETAMAX); 00231 l1GctTauJetsEtEtaPhi_ = dbe->book2D("TauJetsEtEtaPhi", "TAU JET RANK", 00232 PHIBINS, PHIMIN, PHIMAX, 00233 ETABINS, ETAMIN, ETAMAX); 00234 l1GctIsoEmRankEtaPhi_ = dbe->book2D("IsoEmRankEtaPhi", "ISO EM RANK", 00235 PHIBINS, PHIMIN, PHIMAX, 00236 ETABINS, ETAMIN, ETAMAX); 00237 l1GctNonIsoEmRankEtaPhi_ = dbe->book2D("NonIsoEmRankEtaPhi", "NON-ISO EM RANK", 00238 PHIBINS, PHIMIN, PHIMAX, 00239 ETABINS, ETAMIN, ETAMAX); 00240 l1GctCenJetsOccEtaPhi_ = dbe->book2D("CenJetsOccEtaPhi", "CENTRAL JET OCCUPANCY", 00241 PHIBINS, PHIMIN, PHIMAX, 00242 ETABINS, ETAMIN, ETAMAX); 00243 l1GctForJetsOccEtaPhi_ = dbe->book2D("ForJetsOccEtaPhi", "FORWARD JET OCCUPANCY", 00244 PHIBINS, PHIMIN, PHIMAX, 00245 ETABINS, ETAMIN, ETAMAX); 00246 l1GctTauJetsOccEtaPhi_ = dbe->book2D("TauJetsOccEtaPhi", "TAU JET OCCUPANCY", 00247 PHIBINS, PHIMIN, PHIMAX, 00248 ETABINS, ETAMIN, ETAMAX); 00249 l1GctIsoEmOccEtaPhi_ = dbe->book2D("IsoEmOccEtaPhi", "ISO EM OCCUPANCY", 00250 PHIBINS, PHIMIN, PHIMAX, 00251 ETABINS, ETAMIN, ETAMAX); 00252 l1GctNonIsoEmOccEtaPhi_ = dbe->book2D("NonIsoEmOccEtaPhi", "NON-ISO EM OCCUPANCY", 00253 PHIBINS, PHIMIN, PHIMAX, 00254 ETABINS, ETAMIN, ETAMAX); 00255 00256 l1GctHFRing0PosEtaNegEta_ = dbe->book2D("HFRing0Corr", "HF RING0 CORRELATION NEG POS ETA", 00257 PHIBINS, PHIMIN, PHIMAX, 00258 ETABINS, ETAMIN, ETAMAX); 00259 l1GctHFRing1PosEtaNegEta_ = dbe->book2D("HFRing1Corr", "HF RING1 CORRELATION NEG POS ETA", 00260 PHIBINS, PHIMIN, PHIMAX, 00261 ETABINS, ETAMIN, ETAMAX); 00262 l1GctHFRing0TowerCountPosEtaNegEta_ = dbe->book2D("HFRing0TowerCountCorr", "HF RING0 TOWER COUNT CORRELATION NEG POS ETA", 00263 PHIBINS, PHIMIN, PHIMAX, 00264 ETABINS, ETAMIN, ETAMAX); 00265 00266 l1GctHFRing1TowerCountPosEtaNegEta_ = dbe->book2D("HFRing1TowerCountCorr", "HF RING1 TOWER COUNT CORRELATION NEG POS ETA", 00267 PHIBINS, PHIMIN, PHIMAX, 00268 ETABINS, ETAMIN, ETAMAX); 00269 00270 00271 // For Qtests need 1D eta and phi histograms (would be better if Qtests ran on 2D histograms too!) 00272 l1GctCenJetsOccEta_ = dbe->book1D("CenJetsOccEta", "CENTRAL JET ETA OCCUPANCY", ETABINS, ETAMIN, ETAMAX); 00273 l1GctCenJetsOccPhi_ = dbe->book1D("CenJetsOccPhi", "CENTRAL JET PHI OCCUPANCY", PHIBINS, PHIMIN, PHIMAX); 00274 l1GctForJetsOccEta_ = dbe->book1D("ForJetsOccEta", "FORWARD JET ETA OCCUPANCY", ETABINS, ETAMIN, ETAMAX); 00275 l1GctForJetsOccPhi_ = dbe->book1D("ForJetsOccPhi", "FORWARD JET PHI OCCUPANCY", PHIBINS, PHIMIN, PHIMAX); 00276 l1GctTauJetsOccEta_ = dbe->book1D("TauJetsOccEta", "TAU JET ETA OCCUPANCY", ETABINS, ETAMIN, ETAMAX); 00277 l1GctTauJetsOccPhi_ = dbe->book1D("TauJetsOccPhi", "TAU JET PHI OCCUPANCY", PHIBINS, PHIMIN, PHIMAX); 00278 l1GctIsoEmOccEta_ = dbe->book1D("IsoEmOccEta", "ISO EM ETA OCCUPANCY", ETABINS, ETAMIN, ETAMAX); 00279 l1GctIsoEmOccPhi_ = dbe->book1D("IsoEmOccPhi", "ISO EM PHI OCCUPANCY", PHIBINS, PHIMIN, PHIMAX); 00280 l1GctNonIsoEmOccEta_ = dbe->book1D("NonIsoEmOccEta", "NON-ISO EM ETA OCCUPANCY", ETABINS, ETAMIN, ETAMAX); 00281 l1GctNonIsoEmOccPhi_ = dbe->book1D("NonIsoEmOccPhi", "NON-ISO EM PHI OCCUPANCY", PHIBINS, PHIMIN, PHIMAX); 00282 00283 //HF Ring stuff 00284 l1GctHFRing0TowerCountPosEta_ = dbe->book1D("HFRing0TowerCountPosEta", "POS ETA RING0 HFRING BIT", R3BINS, R3MIN, R3MAX); 00285 l1GctHFRing0TowerCountNegEta_ = dbe->book1D("HFRing0TowerCountNegEta", "NEG ETA RING0 HFRING BIT", R3BINS, R3MIN, R3MAX); 00286 l1GctHFRing1TowerCountPosEta_ = dbe->book1D("HFRing1TowerCountPosEta", "POS ETA RING1 HFRING BIT", R3BINS, R3MIN, R3MAX); 00287 l1GctHFRing1TowerCountNegEta_ = dbe->book1D("HFRing1TowerCountNegEta", "NEG ETA RING1 HFRING BIT", R3BINS, R3MIN, R3MAX); 00288 00289 l1GctHFRing0ETSumPosEta_ = dbe->book1D("HFRing0ETSumPosEta", "POS ETA RING0 ET SUM", R3BINS, R3MIN, R3MAX); 00290 l1GctHFRing0ETSumNegEta_ = dbe->book1D("HFRing0ETSumNegEta", "NEG ETA RING0 ET SUM", R3BINS, R3MIN, R3MAX); 00291 l1GctHFRing1ETSumPosEta_ = dbe->book1D("HFRing1ETSumPosEta", "POS ETA RING1 ET SUM", R3BINS, R3MIN, R3MAX); 00292 l1GctHFRing1ETSumNegEta_ = dbe->book1D("HFRing1ETSumNegEta", "NEG ETA RING1 ET SUM", R3BINS, R3MIN, R3MAX); 00293 l1GctHFRingRatioPosEta_ = dbe->book1D("HFRingRatioPosEta", "RING RATIO POS ETA", R5BINS, R5MIN, R5MAX); 00294 l1GctHFRingRatioNegEta_ = dbe->book1D("HFRingRatioNegEta", "RING RATIO NEG ETA", R5BINS, R5MIN, R5MAX); 00295 00296 // Rank histograms 00297 l1GctCenJetsRank_ = dbe->book1D("CenJetsRank", "CENTRAL JET RANK", R6BINS, R6MIN, R6MAX); 00298 l1GctForJetsRank_ = dbe->book1D("ForJetsRank", "FORWARD JET RANK", R6BINS, R6MIN, R6MAX); 00299 l1GctTauJetsRank_ = dbe->book1D("TauJetsRank", "TAU JET RANK", R6BINS, R6MIN, R6MAX); 00300 l1GctIsoEmRank_ = dbe->book1D("IsoEmRank", "ISO EM RANK", R6BINS, R6MIN, R6MAX); 00301 l1GctNonIsoEmRank_ = dbe->book1D("NonIsoEmRank", "NON-ISO EM RANK", R6BINS, R6MIN, R6MAX); 00302 00303 // Energy sums 00304 l1GctEtMiss_ = dbe->book1D("EtMiss", "MISSING E_{T}", R12BINS, R12MIN, R12MAX); 00305 l1GctEtMissPhi_ = dbe->book1D("EtMissPhi", "MISSING E_{T} #phi", METPHIBINS, PHIMIN, PHIMAX); 00306 l1GctEtTotal_ = dbe->book1D("EtTotal", "TOTAL E_{T}", R12BINS, R12MIN, R12MAX); 00307 l1GctEtHad_ = dbe->book1D("EtHad", "TOTAL HAD E_{T}", R12BINS, R12MIN, R12MAX); 00308 00309 // More detailed EM quantities 00310 l1GctIsoEmRankCand0_ = dbe->book1D("GctIsoEmRankCand0","ISO EM RANK CAND 0", R6BINS, R6MIN, R6MAX); 00311 l1GctIsoEmRankCand1_ = dbe->book1D("GctIsoEmRankCand1","ISO EM RANK CAND 1", R6BINS, R6MIN, R6MAX); 00312 l1GctIsoEmRankCand2_ = dbe->book1D("GctIsoEmRankCand2","ISO EM RANK CAND 2", R6BINS, R6MIN, R6MAX); 00313 l1GctIsoEmRankCand3_ = dbe->book1D("GctIsoEmRankCand3","ISO EM RANK CAND 3", R6BINS, R6MIN, R6MAX); 00314 00315 l1GctNonIsoEmRankCand0_ = dbe->book1D("GctNonIsoEmRankCand0","NON-ISO EM RANK CAND 0", R6BINS, R6MIN, R6MAX); 00316 l1GctNonIsoEmRankCand1_ = dbe->book1D("GctNonIsoEmRankCand1","NON-ISO EM RANK CAND 1", R6BINS, R6MIN, R6MAX); 00317 l1GctNonIsoEmRankCand2_ = dbe->book1D("GctNonIsoEmRankCand2","NON-ISO EM RANK CAND 2", R6BINS, R6MIN, R6MAX); 00318 l1GctNonIsoEmRankCand3_ = dbe->book1D("GctNonIsoEmRankCand3","NON-ISO EM RANK CAND 3", R6BINS, R6MIN, R6MAX); 00319 00320 l1GctCenJetsRankCand0_ = dbe->book1D("GctCenJetsRankCand0","CEN JET RANK CAND 0", R6BINS, R6MIN, R6MAX); 00321 l1GctCenJetsRankCand1_ = dbe->book1D("GctCenJetsRankCand1","CEN JET RANK CAND 1", R6BINS, R6MIN, R6MAX); 00322 l1GctCenJetsRankCand2_ = dbe->book1D("GctCenJetsRankCand2","CEN JET RANK CAND 2", R6BINS, R6MIN, R6MAX); 00323 l1GctCenJetsRankCand3_ = dbe->book1D("GctCenJetsRankCand3","CEN JET RANK CAND 3", R6BINS, R6MIN, R6MAX); 00324 00325 l1GctForJetsRankCand0_ = dbe->book1D("GctForJetsRankCand0","FOR JET RANK CAND 0", R6BINS, R6MIN, R6MAX); 00326 l1GctForJetsRankCand1_ = dbe->book1D("GctForJetsRankCand1","FOR JET RANK CAND 1", R6BINS, R6MIN, R6MAX); 00327 l1GctForJetsRankCand2_ = dbe->book1D("GctForJetsRankCand2","FOR JET RANK CAND 2", R6BINS, R6MIN, R6MAX); 00328 l1GctForJetsRankCand3_ = dbe->book1D("GctForJetsRankCand3","FOR JET RANK CAND 3", R6BINS, R6MIN, R6MAX); 00329 00330 l1GctTauJetsRankCand0_ = dbe->book1D("GctTauJetsRankCand0","TAU JET RANK CAND 0", R6BINS, R6MIN, R6MAX); 00331 l1GctTauJetsRankCand1_ = dbe->book1D("GctTauJetsRankCand1","TAU JET RANK CAND 1", R6BINS, R6MIN, R6MAX); 00332 l1GctTauJetsRankCand2_ = dbe->book1D("GctTauJetsRankCand2","TAU JET RANK CAND 2", R6BINS, R6MIN, R6MAX); 00333 l1GctTauJetsRankCand3_ = dbe->book1D("GctTauJetsRankCand3","TAU JET RANK CAND 3", R6BINS, R6MIN, R6MAX); 00334 00335 l1GctIsoEmRankDiff01_ = dbe->book1D("GctIsoEmRankDiffCand01","ISO EM RANK CAND 0 - CAND 1", 2*R6BINS, -R6MAX, R6MAX); 00336 l1GctIsoEmRankDiff12_ = dbe->book1D("GctIsoEmRankDiffCand12","ISO EM RANK CAND 1 - CAND 2", 2*R6BINS, -R6MAX, R6MAX); 00337 l1GctIsoEmRankDiff23_ = dbe->book1D("GctIsoEmRankDiffCand23","ISO EM RANK CAND 2 - CAND 3", 2*R6BINS, -R6MAX, R6MAX); 00338 00339 l1GctNonIsoEmRankDiff01_ = dbe->book1D("GctNonIsoEmRankDiffCand01","NON-ISO EM RANK CAND 0 - CAND 1", 2*R6BINS, -R6MAX, R6MAX); 00340 l1GctNonIsoEmRankDiff12_ = dbe->book1D("GctNonIsoEmRankDiffCand12","NON-ISO EM RANK CAND 1 - CAND 2", 2*R6BINS, -R6MAX, R6MAX); 00341 l1GctNonIsoEmRankDiff23_ = dbe->book1D("GctNonIsoEmRankDiffCand23","NON-ISO EM RANK CAND 2 - CAND 3", 2*R6BINS, -R6MAX, R6MAX); 00342 00343 l1GctCenJetsRankDiff01_ = dbe->book1D("GctCenJetsRankDiffCand01","CEN JET RANK CAND 0 - CAND 1", 2*R6BINS, -R6MAX, R6MAX); 00344 l1GctCenJetsRankDiff12_ = dbe->book1D("GctCenJetsRankDiffCand12","CEN JET RANK CAND 1 - CAND 2", 2*R6BINS, -R6MAX, R6MAX); 00345 l1GctCenJetsRankDiff23_ = dbe->book1D("GctCenJetsRankDiffCand23","CEN JET RANK CAND 2 - CAND 3", 2*R6BINS, -R6MAX, R6MAX); 00346 00347 l1GctForJetsRankDiff01_ = dbe->book1D("GctForJetsRankDiffCand01","FOR JET RANK CAND 0 - CAND 1", 2*R6BINS, -R6MAX, R6MAX); 00348 l1GctForJetsRankDiff12_ = dbe->book1D("GctForJetsRankDiffCand12","FOR JET RANK CAND 1 - CAND 2", 2*R6BINS, -R6MAX, R6MAX); 00349 l1GctForJetsRankDiff23_ = dbe->book1D("GctForJetsRankDiffCand23","FOR JET RANK CAND 2 - CAND 3", 2*R6BINS, -R6MAX, R6MAX); 00350 00351 l1GctTauJetsRankDiff01_ = dbe->book1D("GctTauJetsRankDiffCand01","TAU JET RANK CAND 0 - CAND 1", 2*R6BINS, -R6MAX, R6MAX); 00352 l1GctTauJetsRankDiff12_ = dbe->book1D("GctTauJetsRankDiffCand12","TAU JET RANK CAND 1 - CAND 2", 2*R6BINS, -R6MAX, R6MAX); 00353 l1GctTauJetsRankDiff23_ = dbe->book1D("GctTauJetsRankDiffCand23","TAU JET RANK CAND 2 - CAND 3", 2*R6BINS, -R6MAX, R6MAX); 00354 00355 } 00356 00357 }
Reimplemented from edm::EDAnalyzer.
Definition at line 360 of file L1TGCT.cc.
References dbe, lat::endl(), nev_, outputFile_, DQMStore::save(), and verbose_.
00361 { 00362 if (verbose_) 00363 edm::LogInfo("L1TGCT") << "L1TGCT: end job...." << std::endl; 00364 edm::LogInfo("EndJob") << "analyzed " << nev_ << " events"; 00365 00366 if (outputFile_.size() != 0 && dbe) { 00367 dbe->save(outputFile_); 00368 } 00369 00370 return; 00371 }
DQMStore* L1TGCT::dbe [private] |
edm::InputTag L1TGCT::gctCenJetsSource_ [private] |
edm::InputTag L1TGCT::gctEnergySumsSource_ [private] |
edm::InputTag L1TGCT::gctForJetsSource_ [private] |
edm::InputTag L1TGCT::gctIsoEmSource_ [private] |
edm::InputTag L1TGCT::gctNonIsoEmSource_ [private] |
edm::InputTag L1TGCT::gctTauJetsSource_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsEtEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsOccEta_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsOccEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsOccPhi_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRank_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankCand0_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankCand1_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankCand2_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankCand3_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankDiff01_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankDiff12_ [private] |
MonitorElement* L1TGCT::l1GctCenJetsRankDiff23_ [private] |
MonitorElement* L1TGCT::l1GctEtHad_ [private] |
MonitorElement* L1TGCT::l1GctEtMiss_ [private] |
MonitorElement* L1TGCT::l1GctEtMissPhi_ [private] |
MonitorElement* L1TGCT::l1GctEtTotal_ [private] |
MonitorElement* L1TGCT::l1GctForJetsEtEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctForJetsOccEta_ [private] |
MonitorElement* L1TGCT::l1GctForJetsOccEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctForJetsOccPhi_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRank_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankCand0_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankCand1_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankCand2_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankCand3_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankDiff01_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankDiff12_ [private] |
MonitorElement* L1TGCT::l1GctForJetsRankDiff23_ [private] |
MonitorElement* L1TGCT::l1GctHFRing0ETSumNegEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRing0ETSumPosEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRing0PosEtaNegEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRing1ETSumNegEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRing1ETSumPosEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRing1PosEtaNegEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRingRatioNegEta_ [private] |
MonitorElement* L1TGCT::l1GctHFRingRatioPosEta_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmOccEta_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmOccEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmOccPhi_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRank_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankCand0_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankCand1_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankCand2_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankCand3_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankDiff01_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankDiff12_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankDiff23_ [private] |
MonitorElement* L1TGCT::l1GctIsoEmRankEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmOccEta_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmOccEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmOccPhi_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRank_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankCand0_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankCand1_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankCand2_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankCand3_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankDiff01_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankDiff12_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankDiff23_ [private] |
MonitorElement* L1TGCT::l1GctNonIsoEmRankEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsEtEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsOccEta_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsOccEtaPhi_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsOccPhi_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRank_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankCand0_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankCand1_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankCand2_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankCand3_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankDiff01_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankDiff12_ [private] |
MonitorElement* L1TGCT::l1GctTauJetsRankDiff23_ [private] |
ofstream L1TGCT::logFile_ [private] |
bool L1TGCT::monitorDaemon_ [private] |
int L1TGCT::nev_ [private] |
std::string L1TGCT::outputFile_ [private] |
bool L1TGCT::verbose_ [private] |