|
|
#include <METAnalyzer.h>
|
void | analyze (const edm::Event &, const edm::EventSetup &) override |
| Get the analysis. More...
|
|
void | bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override |
| Inizialize parameters for histo binning. More...
|
|
void | dqmBeginRun (const edm::Run &, const edm::EventSetup &) override |
| Initialize run-based parameters. More...
|
|
void | dqmEndRun (const edm::Run &iRun, const edm::EventSetup &iSetup) override |
| Finish up a run. More...
|
|
void | fillMESet (const edm::Event &, std::string, const reco::MET &, const pat::MET *, const reco::PFMET *, const reco::CaloMET *, const reco::Candidate::PolarLorentzVector &, std::map< std::string, MonitorElement * > &, std::vector< bool >, std::vector< bool >) |
|
void | fillMonitorElement (const edm::Event &, std::string, std::string, const reco::MET &, const pat::MET *, const reco::PFMET *, const reco::CaloMET *, const reco::Candidate::PolarLorentzVector &, std::map< std::string, MonitorElement * > &, bool, bool, std::vector< bool >, std::vector< bool >) |
|
void | makeRatePlot (std::string, double) |
|
| METAnalyzer (const edm::ParameterSet &) |
| Constructor. More...
|
|
| ~METAnalyzer () override |
| Destructor. More...
|
|
void | accumulate (edm::Event const &event, edm::EventSetup const &setup) override |
|
void | beginRun (edm::Run const &run, edm::EventSetup const &setup) final |
|
| DQMOneEDAnalyzer () |
|
void | endRun (edm::Run const &, edm::EventSetup const &) final |
|
void | endRunProduce (edm::Run &run, edm::EventSetup const &setup) final |
|
virtual bool | getCanSaveByLumi () |
|
| EDProducer ()=default |
|
SerialTaskQueue * | globalLuminosityBlocksQueue () final |
|
SerialTaskQueue * | globalRunsQueue () final |
|
bool | hasAbilityToProduceInBeginLumis () const final |
|
bool | hasAbilityToProduceInBeginProcessBlocks () const final |
|
bool | hasAbilityToProduceInBeginRuns () const final |
|
bool | hasAbilityToProduceInEndLumis () const final |
|
bool | hasAbilityToProduceInEndProcessBlocks () const final |
|
bool | hasAbilityToProduceInEndRuns () const final |
|
bool | wantsGlobalLuminosityBlocks () const final |
|
bool | wantsGlobalRuns () const final |
|
bool | wantsInputProcessBlocks () const final |
|
bool | wantsProcessBlocks () const final |
|
| EDProducerBase () |
|
ModuleDescription const & | moduleDescription () const |
|
bool | wantsStreamLuminosityBlocks () const |
|
bool | wantsStreamRuns () const |
|
| ~EDProducerBase () override |
|
void | callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func) |
|
std::vector< edm::ProductResolverIndex > const & | indiciesForPutProducts (BranchType iBranchType) const |
|
| ProducerBase () |
|
std::vector< edm::ProductResolverIndex > const & | putTokenIndexToProductResolverIndex () const |
|
std::vector< bool > const & | recordProvenanceList () const |
|
void | registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &) |
|
std::function< void(BranchDescription const &)> | registrationCallback () const |
| used by the fwk to register list of products More...
|
|
void | resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel) |
|
TypeLabelList const & | typeLabelList () const |
| used by the fwk to register the list of products of this module More...
|
|
| ~ProducerBase () noexcept(false) override |
|
std::vector< ConsumesInfo > | consumesInfo () const |
|
void | convertCurrentProcessAlias (std::string const &processName) |
| Convert "@currentProcess" in InputTag process names to the actual current process name. More...
|
|
| EDConsumerBase () |
|
| EDConsumerBase (EDConsumerBase &&)=default |
|
| EDConsumerBase (EDConsumerBase const &)=delete |
|
ESProxyIndex const * | esGetTokenIndices (edm::Transition iTrans) const |
|
std::vector< ESProxyIndex > const & | esGetTokenIndicesVector (edm::Transition iTrans) const |
|
std::vector< ESRecordIndex > const & | esGetTokenRecordIndicesVector (edm::Transition iTrans) const |
|
ProductResolverIndexAndSkipBit | indexFrom (EDGetToken, BranchType, TypeID const &) const |
|
void | itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const |
|
void | itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const |
|
std::vector< ProductResolverIndexAndSkipBit > const & | itemsToGetFrom (BranchType iType) const |
|
void | labelsForToken (EDGetToken iToken, Labels &oLabels) const |
|
void | modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const |
|
EDConsumerBase & | operator= (EDConsumerBase &&)=default |
|
EDConsumerBase const & | operator= (EDConsumerBase const &)=delete |
|
bool | registeredToConsume (ProductResolverIndex, bool, BranchType) const |
|
bool | registeredToConsumeMany (TypeID const &, BranchType) const |
|
ProductResolverIndexAndSkipBit | uncheckedIndexFrom (EDGetToken) const |
|
void | updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet) |
|
void | updateLookup (eventsetup::ESRecordsToProxyIndices const &) |
|
virtual | ~EDConsumerBase () noexcept(false) |
|
|
typedef dqm::reco::DQMStore | DQMStore |
|
typedef dqm::reco::MonitorElement | MonitorElement |
|
typedef EDProducerBase | ModuleType |
|
using | ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > > |
|
typedef ProductRegistryHelper::TypeLabelList | TypeLabelList |
|
typedef ProductLabels | Labels |
|
static const std::string & | baseType () |
|
static void | fillDescriptions (ConfigurationDescriptions &descriptions) |
|
static void | prevalidate (ConfigurationDescriptions &descriptions) |
|
template<class ProductType > |
BranchAliasSetterT< ProductType > | produces () |
| declare what type of product will make and with which optional label More...
|
|
template<typename ProductType , BranchType B> |
BranchAliasSetterT< ProductType > | produces () |
|
template<typename ProductType , Transition B> |
BranchAliasSetterT< ProductType > | produces () |
|
BranchAliasSetter | produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true) |
|
template<BranchType B> |
BranchAliasSetter | produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true) |
|
template<Transition B> |
BranchAliasSetter | produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true) |
|
template<typename ProductType , Transition B> |
BranchAliasSetterT< ProductType > | produces (std::string instanceName) |
|
template<class ProductType > |
BranchAliasSetterT< ProductType > | produces (std::string instanceName) |
|
template<typename ProductType , BranchType B> |
BranchAliasSetterT< ProductType > | produces (std::string instanceName) |
|
ProducesCollector | producesCollector () |
|
EDGetToken | consumes (const TypeToGet &id, edm::InputTag const &tag) |
|
template<typename ProductType , BranchType B = InEvent> |
EDGetTokenT< ProductType > | consumes (edm::InputTag const &tag) |
|
template<BranchType B = InEvent> |
EDConsumerBaseAdaptor< B > | consumes (edm::InputTag tag) noexcept |
|
template<BranchType B> |
EDGetToken | consumes (TypeToGet const &id, edm::InputTag const &tag) |
|
ConsumesCollector | consumesCollector () |
| Use a ConsumesCollector to gather consumes information from helper functions. More...
|
|
template<typename ProductType , BranchType B = InEvent> |
void | consumesMany () |
|
void | consumesMany (const TypeToGet &id) |
|
template<BranchType B> |
void | consumesMany (const TypeToGet &id) |
|
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event> |
auto | esConsumes () |
|
template<Transition Tr = Transition::Event> |
constexpr auto | esConsumes () noexcept |
|
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event> |
auto | esConsumes (ESInputTag const &tag) |
|
template<Transition Tr = Transition::Event> |
auto | esConsumes (ESInputTag tag) noexcept |
|
template<Transition Tr = Transition::Event> |
ESGetTokenGeneric | esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey) |
| Used with EventSetupRecord::doGet. More...
|
|
EDGetToken | mayConsume (const TypeToGet &id, edm::InputTag const &tag) |
|
template<BranchType B> |
EDGetToken | mayConsume (const TypeToGet &id, edm::InputTag const &tag) |
|
template<typename ProductType , BranchType B = InEvent> |
EDGetTokenT< ProductType > | mayConsume (edm::InputTag const &tag) |
|
edm::EDPutTokenT< DQMToken > | runToken_ |
|
DQM MET analysis monitoring
- Author
- F. Chlebana - Fermilab K. Hatakeyama - Rockefeller University
Jan. '14: modified by
M. Artur Weber R. Schoefbeck V. Sordini
Definition at line 89 of file METAnalyzer.h.
◆ METAnalyzer()
◆ ~METAnalyzer()
METAnalyzer::~METAnalyzer |
( |
| ) |
|
|
override |
◆ analyze()
Get the analysis.
Reimplemented from DQMOneEDAnalyzer<>.
Definition at line 1415 of file METAnalyzer.cc.
1417 int myLuminosityBlock;
1418 myLuminosityBlock =
iEvent.luminosityBlock();
1431 std::cout <<
"METAnalyzer analyze" << std::endl;
1447 std::vector<int> triggerFolderDecisions;
1456 int ntrigs = (*triggerResults).size();
1458 std::cout <<
"ntrigs=" << ntrigs << std::endl;
1464 bool fd = (*it)->accept(
iEvent, iSetup);
1503 const MET*
met =
nullptr;
1505 const PFMET* pfmet =
nullptr;
1506 const CaloMET* calomet =
nullptr;
1511 assert(!pfmetcoll->empty());
1512 met = &(pfmetcoll->front());
1513 pfmet = &(pfmetcoll->front());
1516 assert(!calometcoll->empty());
1517 met = &(calometcoll->front());
1518 calomet = &(calometcoll->front());
1521 assert(!patmetcoll->empty());
1522 met = &(patmetcoll->front());
1523 patmet = &(patmetcoll->front());
1526 LogTrace(
"METAnalyzer") <<
"[METAnalyzer] Call to the MET analyzer";
1548 bool bJetID =
false;
1549 bool bDiJetID =
false;
1563 LogDebug(
"") <<
"METAnalyzer: Could not find calojet product" << std::endl;
1565 std::cout <<
"METAnalyzer: Could not find calojet product" << std::endl;
1589 LogDebug(
"") <<
"METAnalyzer: Could not find patjet product" << std::endl;
1591 std::cout <<
"METAnalyzer: Could not find patjet product" << std::endl;
1599 LogDebug(
"") <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1601 std::cout <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1603 collsize =
pfJets->size();
1606 unsigned int ind1 = -1;
1608 bool pass_jetID1 =
false;
1609 unsigned int ind2 = -1;
1611 bool pass_jetID2 =
false;
1614 bool pass_correction_flag =
false;
1618 pass_correction_flag =
true;
1621 pass_correction_flag =
true;
1624 for (
int ijet = 0; ijet < collsize; ijet++) {
1625 double pt_jet = -10;
1627 bool iscleaned =
false;
1628 if (pass_correction_flag) {
1640 pt_jet =
scale * (*caloJets)[ijet].pt();
1644 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
1667 pt_jet =
scale * (*pfJets)[ijet].pt();
1673 pt_jet = (*patJets)[ijet].pt();
1685 pass_jetID2 = pass_jetID1;
1688 pass_jetID1 = iscleaned;
1689 }
else if (pt_jet >
pt2) {
1692 pass_jetID2 = iscleaned;
1695 if (pass_jetID1 && pass_jetID2) {
1709 if (dphi > acos(-1.)) {
1710 dphi = 2 * acos(-1.) - dphi;
1722 if (!vertexHandle.
isValid()) {
1723 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1725 std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1734 bool bZJets =
false;
1741 double pt_muon0 = -1;
1742 double pt_muon1 = -1;
1747 for (
unsigned int i = 0;
i <
Muons->size();
i++) {
1748 bool pass_muon_id =
false;
1749 bool pass_muon_iso =
false;
1751 double dz = fabs((*
Muons)[
i].muonBestTrack()->
dz());
1757 if ((*
Muons)[
i].isGlobalMuon() && (*
Muons)[
i].isPFMuon() &&
1759 (*
Muons)[
i].numberOfMatchedStations() > 1 &&
dxy < 0.2 && (*
Muons)[
i].numberOfMatchedStations() > 1 &&
1761 (*
Muons)[
i].
innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5) {
1762 pass_muon_id =
true;
1766 ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt +
1771 if (muonIsoPFdb < 0.12) {
1772 pass_muon_iso =
true;
1775 if (pass_muon_id && pass_muon_iso) {
1777 mu_index1 = mu_index0;
1778 pt_muon1 = pt_muon0;
1780 pt_muon0 = (*Muons)[
i].pt();
1781 }
else if ((*
Muons)[
i].
pt() > pt_muon1) {
1783 pt_muon1 = (*Muons)[
i].pt();
1788 if (mu_index0 >= 0 && mu_index1 >= 0) {
1790 zCand = (*Muons)[mu_index0].polarP4() + (*Muons)[mu_index1].polarP4();
1791 if (fabs(zCand.M() - 91.) < 20) {
1803 bool techTriggerResultBxM2 =
false;
1804 bool techTriggerResultBxM1 =
false;
1805 bool techTriggerResultBx0 =
false;
1807 if (!gtReadoutRecord.
isValid()) {
1808 LogDebug(
"") <<
"METAnalyzer: Could not find GT readout record" << std::endl;
1810 std::cout <<
"METAnalyzer: Could not find GT readout record product" << std::endl;
1820 if (techTriggerResultBx0 != 0) {
1821 techTriggerResultBxM2 = technicalTriggerWordBeforeMaskBxM2.at(
m_bitAlgTechTrig_);
1822 techTriggerResultBxM1 = technicalTriggerWordBeforeMaskBxM1.at(
m_bitAlgTechTrig_);
1827 std::vector<bool> trigger_flag(4,
false);
1828 if (techTriggerResultBx0 && techTriggerResultBxM2 &&
1829 techTriggerResultBxM1) {
1830 trigger_flag[0] =
true;
1832 if (techTriggerResultBx0 && techTriggerResultBxM1) {
1833 trigger_flag[1] =
true;
1835 if (techTriggerResultBx0 && !techTriggerResultBxM1) {
1836 trigger_flag[2] =
true;
1838 if (techTriggerResultBx0 && !techTriggerResultBxM2 &&
1839 !techTriggerResultBxM1) {
1840 trigger_flag[3] =
true;
1842 std::vector<bool> filter_decisions(
1848 if (!HBHENoiseFilterResultHandle.
isValid()) {
1849 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1851 std::cout <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1853 filter_decisions[0] = *HBHENoiseFilterResultHandle;
1856 if (!CSCTightHaloFilterResultHandle.
isValid()) {
1857 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHaloFilterResultHandle" << std::endl;
1859 std::cout <<
"METAnalyzer: CSCTightHaloFilterResultHandle" << std::endl;
1861 filter_decisions[1] = *CSCTightHaloFilterResultHandle;
1864 if (!eeBadScFilterResultHandle.
isValid()) {
1865 LogDebug(
"") <<
"METAnalyzer: Could not find eeBadScFilterResultHandle" << std::endl;
1867 std::cout <<
"METAnalyzer: eeBadScFilterResultHandle" << std::endl;
1869 filter_decisions[2] = *eeBadScFilterResultHandle;
1872 if (!HBHENoiseIsoFilterResultHandle.
isValid()) {
1873 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1875 std::cout <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1877 filter_decisions[3] = *HBHENoiseIsoFilterResultHandle;
1880 if (!CSCTightHalo2015FilterResultHandle.
isValid()) {
1881 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1883 std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1885 filter_decisions[4] = *CSCTightHalo2015FilterResultHandle;
1888 if (!EcalDeadCellTriggerFilterResultHandle.
isValid()) {
1889 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellTriggerFilterResultHandle" << std::endl;
1891 std::cout <<
"METAnalyzer: EcalDeadCellTriggerFilterResultHandle" << std::endl;
1893 filter_decisions[5] = *EcalDeadCellTriggerFilterResultHandle;
1896 if (!EcalDeadCellBoundaryHandle.
isValid()) {
1897 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellBoundaryHandle" << std::endl;
1899 std::cout <<
"METAnalyzer: EcalDeadCellBoundaryHandle" << std::endl;
1901 filter_decisions[6] = *EcalDeadCellBoundaryHandle;
1904 if (!HcalStripHaloFilterHandle.
isValid()) {
1905 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1907 std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1909 filter_decisions[7] = *HcalStripHaloFilterHandle;
1914 if (metFilterResults.
isValid()) {
1941 if (metFilterResults.
isValid()) {
1969 bool HBHENoiseFilterResultFlag = filter_decisions[0];
1972 bool bHBHENoiseFilter = HBHENoiseFilterResultFlag;
1980 bool pass_selection =
false;
1981 if ((*ic ==
"Uncleaned") && (
isCaloMet_ || bPrimaryVertex)) {
1983 DirName_old +
"/" + *ic,
1992 pass_selection =
true;
1995 if ((*ic ==
"Cleaned") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bJetID) {
1997 DirName_old +
"/" + *ic,
2006 pass_selection =
true;
2008 if ((*ic ==
"DiJet") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bDiJetID) {
2010 DirName_old +
"/" + *ic,
2019 pass_selection =
true;
2021 if ((*ic ==
"ZJets") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bZJets) {
2023 DirName_old +
"/" + *ic,
2032 pass_selection =
true;
2035 DirName = DirName_old +
"/" + *ic;
References edm::HLTGlobalStatus::accept(), cms::cuda::assert(), isolatedTracks_cfi::caloJets, ALCARECOTkAlJpsiMuMu_cff::charge, reco::JetCorrector::correction(), gather_cfg::cout, fed_dqm_sourceclient-live_cfg::DirName, PVValHelper::dxy, PVValHelper::dz, PVValHelper::eta, ztee::fd, cutBasedMuonId_MuonPOG_V0_cff::globalTrack, mps_fire::i, iEvent, susybsm::HSCParticleType::innerTrack, edm::HandleBase::isValid(), LogDebug, LogTrace, SiStripPI::max, BTaggingMonitor_cfi::met, HLTTauDQMOffline_cfi::Muons, TrackCollections2monitor_cff::numberOfValidPixelHits, hemisphereProducer_cfi::patJets, muons_cfi::pfIsolationR04, pfJetBenchmark_cfi::pfJets, position, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, HLT_FULL_cff::pt1, HLT_FULL_cff::pt2, HLT_FULL_cff::scale, AlCaHLTBitMon_QueryRunRegistry::string, L1GlobalTriggerReadoutRecord::technicalTriggerWord(), edm::triggerResults(), and spclusmultinvestigator_cfi::vertexCollection.
◆ bookHistograms()
◆ bookMESet()
Definition at line 224 of file METAnalyzer.cc.
229 bool fillPFCandidatePlots =
false;
230 bool fillZPlots =
false;
232 if (
DirName.find(
"Cleaned") != std::string::npos) {
233 fillPFCandidatePlots =
true;
241 }
else if (
DirName.find(
"ZJets") != std::string::npos) {
242 fillPFCandidatePlots =
false;
References fed_dqm_sourceclient-live_cfg::DirName.
◆ bookMonitorElement()
void METAnalyzer::bookMonitorElement |
( |
std::string |
DirName, |
|
|
DQMStore::IBooker & |
ibooker, |
|
|
std::map< std::string, MonitorElement * > & |
map_of_MEs, |
|
|
bool |
bLumiSecPlot = false , |
|
|
bool |
fillPFCandPlots = false , |
|
|
bool |
fillZPlots = false |
|
) |
| |
|
private |
Definition at line 251 of file METAnalyzer.cc.
291 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"u_parallel_Z_inc",
meZJets_u_par));
307 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"u_perp_Z_inc",
meZJets_u_perp));
325 hTrigger = ibooker.
book1D(
"triggerResults",
"triggerResults", 500, 0, 500);
335 hMEx = ibooker.
book1D(
"MEx",
"MEx", 200, -500, 500);
336 hMEy = ibooker.
book1D(
"MEy",
"MEy", 200, -500, 500);
337 hMET = ibooker.
book1D(
"MET",
"MET", 200, 0, 1000);
338 hMET_2 = ibooker.
book1D(
"MET_2",
"MET Range 2", 200, 0, 2000);
339 hSumET = ibooker.
book1D(
"SumET",
"SumET", 400, 0, 4000);
360 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"triggerResults",
hTrigger));
361 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MEx",
hMEx));
362 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MEy",
hMEy));
363 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET",
hMET));
364 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_2",
hMET_2));
365 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"SumET",
hSumET));
366 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METSig",
hMETSig));
367 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhi",
hMETPhi));
368 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_logx",
hMET_logx));
369 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"SumET_logx",
hSumET_logx));
376 ibooker.
book1D(
"MET_CSCTightHalo2015Filter",
"MET_CSCTightHalo2015Filtered", 200, 0, 1000);
378 ibooker.
book1D(
"MET_EcalDeadCellTriggerFilter",
"MET_EcalDeadCellTriggerFiltered", 200, 0, 1000);
380 ibooker.
book1D(
"MET_EcalDeadCellBoundaryFilter",
"MET_EcalDeadCellBoundaryFiltered", 200, 0, 1000);
390 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_CSCTightHalo2015Filter",
392 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_EcalDeadCellTriggerFilter",
394 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_EcalDeadCellBoundaryFilter",
412 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MEx_profile",
meMEx_profile));
413 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MEy_profile",
meMEy_profile));
414 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_profile",
meMET_profile));
415 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"SumET_profile",
meSumET_profile));
433 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloHadEtInHO",
hCaloHadEtInHO));
434 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloHadEtInHF",
hCaloHadEtInHF));
435 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloHadEtInHE",
hCaloHadEtInHE));
436 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloHadEtInHB",
hCaloHadEtInHB));
437 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloEmEtInHF",
hCaloEmEtInHF));
438 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloEmEtInEE",
hCaloEmEtInEE));
439 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloEmEtInEB",
hCaloEmEtInEB));
452 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"CaloMETPhi020",
hCaloMETPhi020));
461 if (fillPFCandPlots &&
464 meCHF_Barrel = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel",
"chargedHadronEtFractionBarrel", 50, 0, 1);
466 ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus",
"chargedHadronEtFractionEndcapPlus", 50, 0, 1);
468 ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus",
"chargedHadronEtFractionEndcapMinus", 50, 0, 1);
470 "PfChargedHadronEtFractionBarrel_BXm1Empty",
"chargedHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
472 "chargedHadronEtFractionEndcapPlus prev empty bunch",
477 "chargedHadronEtFractionEndcapMinus prev empty bunch",
482 "chargedHadronEtFractionBarrel prev filled 2 bunches",
487 "chargedHadronEtFractionEndcapPlus prev filled bunch",
492 "chargedHadronEtFractionEndcapMinus prev filled bunch",
498 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEtFractionBarrel",
meCHF_Barrel));
499 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus",
501 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
503 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
505 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
507 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
512 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
514 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
516 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
522 meNHF_Barrel = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel",
"neutralHadronEtFractionBarrel", 50, 0, 1);
524 ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus",
"neutralHadronEtFractionEndcapPlus", 50, 0, 1);
526 ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus",
"neutralHadronEtFractionEndcapMinus", 50, 0, 1);
528 "PfNeutralHadronEtFractionBarrel_BXm1Empty",
"neutralHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
530 "neutralHadronEtFractionEndcapPlus prev empty bunch",
535 "neutralHadronEtFractionEndcapMinus prev empty bunch",
543 "neutralHadronEtFractionBarrel prev filled 2 bunches",
548 "neutralHadronEtFractionEndcapPlus prev filled bunch",
553 "neutralHadronEtFractionEndcapMinus prev filled bunch",
562 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEtFractionBarrel",
meNHF_Barrel));
563 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus",
565 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
567 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
569 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
571 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
576 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
578 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
580 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
586 mePhF_Barrel = ibooker.
book1D(
"PfPhotonEtFractionBarrel",
"photonEtFractionBarrel", 50, 0, 1);
587 mePhF_EndcapPlus = ibooker.
book1D(
"PfPhotonEtFractionEndcapPlus",
"photonEtFractionEndcapPlus", 50, 0, 1);
590 ibooker.
book1D(
"PfPhotonEtFractionBarrel_BXm1Empty",
"photonEtFractionBarrel prev empty bunch", 50, 0, 1);
592 "PfPhotonEtFractionEndcapPlus_BXm1Empty",
"photonEtFractionEndcapPlus prev empty bunch", 50, 0, 1);
594 "PfPhotonEtFractionEndcapMinus_BXm1Empty",
"photonEtFractionEndcapMinus prev empty bunch", 50, 0, 1);
599 "PfPhotonEtFractionBarrel_BXm1Filled",
"photonEtFractionBarrel prev filled 2 bunches", 50, 0, 1);
601 "PfPhotonEtFractionEndcapPlus_BXm1Filled",
"photonEtFractionEndcapPlus prev filled bunch", 50, 0, 1);
603 "PfPhotonEtFractionEndcapMinus_BXm1Filled",
"photonEtFractionEndcapMinus prev filled bunch", 50, 0, 1);
609 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFractionBarrel",
mePhF_Barrel));
612 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFractionEndcapMinus",
614 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Empty",
616 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
618 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
623 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Filled",
625 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
627 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
633 meHFHadF_Plus = ibooker.
book1D(
"PfHFHadronEtFractionPlus",
"HFHadronEtFractionPlus", 50, 0, 1);
634 meHFHadF_Minus = ibooker.
book1D(
"PfHFHadronEtFractionMinus",
"HFHadronEtFractionMinus", 50, 0, 1);
636 ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Empty",
"HFHadronEtFractionPlus prev empty bunch", 50, 0, 1);
638 ibooker.
book1D(
"PfHFHadronEtFractionMinus_BXm1Empty",
"HFHadronEtFractionMinus prev empty bunch", 50, 0, 1);
642 ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Filled",
"HFHadronEtFractionPlus prev filled bunch", 50, 0, 1);
644 "PfHFHadronEtFractionMinus_BXm1Filled",
"HFHadronEtFractionMinus prev filled bunch", 50, 0, 1);
649 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFractionPlus",
meHFHadF_Plus));
651 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFractionMinus",
meHFHadF_Minus));
652 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Empty",
654 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFractionMinus_BXm1Empty",
658 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Filled",
660 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
665 meHFEMF_Plus = ibooker.
book1D(
"PfHFEMEtFractionPlus",
"HFEMEtFractionPlus", 50, 0, 1);
668 ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Empty",
"HFEMEtFractionPlus prev empty bunch", 50, 0, 1);
670 ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Empty",
"HFEMEtFractionMinus prev empty bunch", 50, 0, 1);
674 ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Filled",
"HFEMEtFractionPlus prev filled bunch", 50, 0, 1);
676 ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Filled",
"HFEMEtFractionMinus prev filled bunch", 50, 0, 1);
681 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionPlus",
meHFEMF_Plus));
683 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionMinus",
meHFEMF_Minus));
684 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Empty",
686 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Empty",
690 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Filled",
692 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Filled",
755 ibooker.
book1D(
"PfPhotonEtFraction_BXm1Empty",
"photonEtFraction() prev empty bunch", 50, 0, 1);
757 ibooker.
book1D(
"PfPhotonEtFraction_BXm1Filled",
"photonEtFraction() prev filled bunch", 50, 0, 1);
759 "PfNeutralHadronEtFraction_BXm1Empty",
"neutralHadronEtFraction() prev empty bunch", 50, 0, 1);
761 "PfNeutralHadronEtFraction_BXm1Filled",
"neutralHadronEtFraction() prev filled bunch", 50, 0, 1);
763 "PfChargedHadronEtFraction_BXm1Empty",
"chargedHadronEtFraction() prev empty bunch", 50, 0, 1);
765 "PfChargedHadronEtFraction_BXm1Filled",
"chargedHadronEtFraction() prev filled bunch", 50, 0, 1);
771 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFraction_BXm1Empty",
773 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFraction_BXm1Filled",
775 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEtFraction_BXm1Empty",
777 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralEtFraction_BXm1Filled",
779 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEtFraction_BXm1Empty",
781 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedEtFraction_BXm1Filled",
783 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_BXm1Empty",
meMET_BXm1Empty));
784 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MET_BXm1Filled",
meMET_BXm1Filled));
791 "METPhiChargedHadronsBarrel_BXm1Filled",
"METPhi_PFChargedHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
793 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Filled",
794 "METPhi_PFChargedHadronsEndcapPlus prev bunch filled",
799 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Filled",
800 "METPhi_PFChargedHadronsEndcapMinus prev bunch filled",
805 "METPhiNeutralHadronsBarrel_BXm1Filled",
"METPhi_PFNeutralHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
807 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Filled",
808 "METPhi_PFNeutralHadronsEndcapPlus prev bunch filled",
813 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Filled",
814 "METPhi_PFNeutralHadronsEndcapMinus prev bunch filled",
819 "METPhiPhotonsBarrel_BXm1Filled",
"METPhi_PFPhotonsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
821 "METPhiPhotonsEndcapPlus_BXm1Filled",
"METPhi_PFPhotonsEndcapPlus prev bunch filled", 50, -
M_PI,
M_PI);
823 "METPhiPhotonsEndcapMinus_BXm1Filled",
"METPhi_PFPhotonsEndcapMinus prev bunch filled", 50, -
M_PI,
M_PI);
825 "METPhiHFHadronsPlus_BXm1Filled",
"METPhi_PFHFHadronsPlus prev bunch filled", 50, -
M_PI,
M_PI);
827 "METPhiHFHadronsMinus_BXm1Filled",
"METPhi_PFHFHadronsMinus prev bunch filled", 50, -
M_PI,
M_PI);
829 "METPhiHFEGammasPlus_BXm1Filled",
"METPhi_PFHFEGammasPlus prev bunch filled", 50, -
M_PI,
M_PI);
831 "METPhiHFEGammasMinus_BXm1Filled",
"METPhi_PFHFEGammasMinus prev bunch filled", 50, -
M_PI,
M_PI);
833 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
835 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
838 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Filled",
840 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
842 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
845 std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Filled",
847 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Filled",
849 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Filled",
851 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Filled",
853 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Filled",
855 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Filled",
857 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Filled",
859 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Filled",
863 "METPhiChargedHadronsBarrel_BXm1Empty",
"METPhi_PFChargedHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
865 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Empty",
866 "METPhi_PFChargedHadronsEndcapPlus prev bunch empty",
871 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Empty",
872 "METPhi_PFChargedHadronsEndcapMinus prev bunch empty",
877 "METPhiNeutralHadronsBarrel_BXm1Empty",
"METPhi_PFNeutralHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
879 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Empty",
880 "METPhi_PFNeutralHadronsEndcapPlus prev bunch empty",
885 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Empty",
886 "METPhi_PFNeutralHadronsEndcapMinus prev bunch empty",
891 ibooker.
book1D(
"METPhiPhotonsBarrel_BXm1Empty",
"METPhi_PFPhotonsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
893 "METPhiPhotonsEndcapPlus_BXm1Empty",
"METPhi_PFPhotonsEndcapPlus prev bunch empty", 50, -
M_PI,
M_PI);
895 "METPhiPhotonsEndcapMinus_BXm1Empty",
"METPhi_PFPhotonsEndcapMinus prev bunch empty", 50, -
M_PI,
M_PI);
897 ibooker.
book1D(
"METPhiHFHadronsPlus_BXm1Empty",
"METPhi_PFHFHadronsPlus prev bunch empty", 50, -
M_PI,
M_PI);
899 "METPhiHFHadronsMinus_BXm1Empty",
"METPhi_PFHFHadronsMinus prev bunch empty", 50, -
M_PI,
M_PI);
901 ibooker.
book1D(
"METPhiHFEGammasPlus_BXm1Empty",
"METPhi_PFHFEGammasPlus prev bunch empty", 50, -
M_PI,
M_PI);
903 "METPhiHFEGammasMinus_BXm1Empty",
"METPhi_PFHFEGammasMinus prev bunch empty", 50, -
M_PI,
M_PI);
905 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
907 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
909 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
911 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
913 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
915 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
917 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Empty",
919 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Empty",
921 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Empty",
923 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Empty",
925 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Empty",
927 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Empty",
929 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Empty",
933 ibooker.
book1D(
"METPhiChargedHadronsBarrel",
"METPhi_PFChargedHadronsBarrel", 50, -
M_PI,
M_PI);
935 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus",
"METPhi_PFChargedHadronsEndcapPlus", 50, -
M_PI,
M_PI);
937 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus",
"METPhi_PFChargedHadronsEndcapMinus", 50, -
M_PI,
M_PI);
939 ibooker.
book1D(
"METPhiNeutralHadronsBarrel",
"METPhi_PFNeutralHadronsBarrel", 50, -
M_PI,
M_PI);
941 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus",
"METPhi_PFNeutralHadronsEndcapPlus", 50, -
M_PI,
M_PI);
943 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus",
"METPhi_PFNeutralHadronsEndcapMinus", 50, -
M_PI,
M_PI);
946 ibooker.
book1D(
"METPhiPhotonsEndcapPlus",
"METPhi_PFPhotonsEndcapPlus", 50, -
M_PI,
M_PI);
948 ibooker.
book1D(
"METPhiPhotonsEndcapMinus",
"METPhi_PFPhotonsEndcapMinus", 50, -
M_PI,
M_PI);
954 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiChargedHadronsBarrel",
956 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiChargedHadronsEndcapPlus",
958 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiChargedHadronsEndcapMinus",
960 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiNeutralHadronsBarrel",
962 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus",
964 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus",
968 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapPlus",
970 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METPhiPhotonsEndcapMinus",
998 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
999 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
1000 int nMinPFCand =
v->getParameter<
int>(
"nMin");
1001 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
1002 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
1021 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
1033 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
1040 if (fillPFCandPlots &&
1053 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
1054 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
1075 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEtFraction",
1077 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEtFraction",
1081 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFraction",
meHFEMEtFraction));
1099 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEtFraction_profile",
1101 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEtFraction_profile",
1103 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEtFraction_profile",
1105 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEtFraction_profile",
1107 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEtFraction_profile",
1118 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfPhotonEt",
mePhotonEt));
1121 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfElectronEt",
meElectronEt));
1124 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfMuonEt",
meMuonEt));
1125 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFHadronEt",
meHFHadronEt));
1126 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEt",
meHFEMEt));
1146 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfNeutralHadronEt_profile",
1148 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfChargedHadronEt_profile",
1152 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"PfHFEMEt_profile",
meHFEMEt_profile));
1157 hMExLS = ibooker.
book2D(
"MExLS",
"MEx_LS", 200, -200, 200, 250, 0., 2500.);
1161 hMEyLS = ibooker.
book2D(
"MEyLS",
"MEy_LS", 200, -200, 200, 250, 0., 2500.);
1165 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MExLS",
hMExLS));
1166 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"MEyLS",
hMEyLS));
1171 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"METRate",
hMETRate));
1175 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
"JetMET/lumisec",
lumisecME));
References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), gather_cfg::cout, fed_dqm_sourceclient-live_cfg::DirName, mps_fire::i, M_PI, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and findQualityFiles::v.
◆ dqmBeginRun()
Initialize run-based parameters.
Reimplemented from DQMOneEDAnalyzer<>.
Definition at line 1180 of file METAnalyzer.cc.
1184 for (
CItAlgo techTrig =
menu->gtTechnicalTriggerMap().begin(); techTrig !=
menu->gtTechnicalTriggerMap().end();
1207 std::cout <<
"HLTEventAnalyzerAOD::analyze:"
1223 (*it)->initRun(iRun, iSetup);
1225 if ((*it)->expressionsFromDB((*it)->hltDBKey(), iSetup)[0] !=
"CONFIG_ERROR")
1231 bool changed_filter =
true;
1232 std::vector<int> initializeFilter(8, -1);
1242 if (
found != std::string::npos) {
1247 if (
found != std::string::npos) {
1252 if (
found != std::string::npos) {
1257 if (
found != std::string::npos) {
1262 if (
found != std::string::npos) {
1267 if (
found != std::string::npos) {
1272 if (
found != std::string::npos) {
1277 if (
found != std::string::npos) {
1289 if (
found != std::string::npos) {
1294 if (
found != std::string::npos) {
1299 if (
found != std::string::npos) {
1304 if (
found != std::string::npos) {
1309 if (
found != std::string::npos) {
1314 if (
found != std::string::npos) {
1319 if (
found != std::string::npos) {
1324 if (
found != std::string::npos) {
1330 <<
"nothing found with both RECO and reRECO label" << std::endl;
References gather_cfg::cout, newFWLiteAna::found, edm::EventSetup::get(), edm::get(), mps_fire::i, optionsL1T::menu, edm::ESHandle< T >::product(), cond::persistency::search(), electronDataDiscovery::search2(), and AlCaHLTBitMon_QueryRunRegistry::string.
◆ dqmEndRun()
◆ fillMESet()
void METAnalyzer::fillMESet |
( |
const edm::Event & |
, |
|
|
std::string |
, |
|
|
const reco::MET & |
, |
|
|
const pat::MET * |
, |
|
|
const reco::PFMET * |
, |
|
|
const reco::CaloMET * |
, |
|
|
const reco::Candidate::PolarLorentzVector & |
, |
|
|
std::map< std::string, MonitorElement * > & |
, |
|
|
std::vector< bool > |
, |
|
|
std::vector< bool > |
|
|
) |
| |
Definition at line 2041 of file METAnalyzer.cc.
2052 bool fillPFCandidatePlots =
false;
2053 if (
DirName.find(
"Cleaned") != std::string::npos) {
2054 fillPFCandidatePlots =
true;
2065 fillPFCandidatePlots,
2069 fillPFCandidatePlots =
false;
2081 fillPFCandidatePlots,
2086 }
else if (
DirName.find(
"DiJet") != std::string::npos) {
2097 fillPFCandidatePlots,
2112 fillPFCandidatePlots,
2116 }
else if (
DirName.find(
"ZJets") != std::string::npos) {
2127 fillPFCandidatePlots,
2141 fillPFCandidatePlots,
References fed_dqm_sourceclient-live_cfg::DirName, mps_fire::i, iEvent, BTaggingMonitor_cfi::met, and AlCaHLTBitMon_QueryRunRegistry::string.
◆ fillMonitorElement()
void METAnalyzer::fillMonitorElement |
( |
const edm::Event & |
, |
|
|
std::string |
, |
|
|
std::string |
, |
|
|
const reco::MET & |
, |
|
|
const pat::MET * |
, |
|
|
const reco::PFMET * |
, |
|
|
const reco::CaloMET * |
, |
|
|
const reco::Candidate::PolarLorentzVector & |
, |
|
|
std::map< std::string, MonitorElement * > & |
, |
|
|
bool |
, |
|
|
bool |
, |
|
|
std::vector< bool > |
, |
|
|
std::vector< bool > |
|
|
) |
| |
Definition at line 2148 of file METAnalyzer.cc.
2161 bool do_only_Z_histograms =
false;
2162 if (
DirName.find(
"ZJets") != std::string::npos) {
2163 do_only_Z_histograms =
true;
2166 double u_x = -
met.px() - zCand.Px();
2167 double u_y = -
met.py() - zCand.Py();
2171 double u_perp =
sqrt(u_x * u_x + u_y * u_y);
2174 if (zCand.Pt() != 0) {
2175 e_Z_x = zCand.Px() / zCand.Pt();
2176 e_Z_y = zCand.Py() / zCand.Pt();
2178 u_par = u_x * e_Z_x + u_y * e_Z_y;
2179 u_perp = -e_Z_y * u_x + e_Z_x * u_y;
2184 if (zCand.Pt() < 15) {
2188 }
else if (zCand.Pt() < 30) {
2192 }
else if (zCand.Pt() < 55) {
2196 }
else if (zCand.Pt() < 75) {
2200 }
else if (zCand.Pt() < 150) {
2204 }
else if (zCand.Pt() < 290) {
2217 if (zCand.Pt() < 15) {
2221 }
else if (zCand.Pt() < 30) {
2225 }
else if (zCand.Pt() < 55) {
2229 }
else if (zCand.Pt() < 75) {
2233 }
else if (zCand.Pt() < 150) {
2237 }
else if (zCand.Pt() < 290) {
2247 if (!do_only_Z_histograms) {
2249 double SumET =
met.sumEt();
2250 double METSig =
met.mEtSig();
2253 double MEx =
met.px();
2254 double MEy =
met.py();
2255 double METPhi =
met.phi();
2257 int myLuminosityBlock;
2258 myLuminosityBlock =
iEvent.luminosityBlock();
2261 if (!subFolderName.empty()) {
2296 bool HBHENoiseFilterResult =
false;
2297 bool CSCTightHaloFilterResult =
false;
2298 bool eeBadScFilterResult =
false;
2299 bool HBHEIsoNoiseFilterResult =
false;
2300 bool CSCTightHalo2015FilterResult =
false;
2301 bool EcalDeadCellTriggerFilterResult =
false;
2302 bool EcalDeadCellBoundaryFilterResult =
false;
2303 bool HcalStripHaloFilterResult =
false;
2304 HBHENoiseFilterResult = METFilterDecision[0];
2305 if (HBHENoiseFilterResult) {
2310 CSCTightHaloFilterResult = METFilterDecision[1];
2311 if (CSCTightHaloFilterResult) {
2316 eeBadScFilterResult = METFilterDecision[2];
2317 if (eeBadScFilterResult) {
2322 HBHEIsoNoiseFilterResult = METFilterDecision[3];
2323 if (HBHEIsoNoiseFilterResult) {
2328 CSCTightHalo2015FilterResult = METFilterDecision[4];
2329 if (CSCTightHalo2015FilterResult) {
2334 EcalDeadCellTriggerFilterResult = METFilterDecision[5];
2335 if (EcalDeadCellTriggerFilterResult) {
2340 EcalDeadCellBoundaryFilterResult = METFilterDecision[6];
2341 if (EcalDeadCellBoundaryFilterResult) {
2346 HcalStripHaloFilterResult = METFilterDecision[7];
2347 if (HcalStripHaloFilterResult) {
2389 double caloEtFractionHadronic = calomet->etFractionHadronic();
2390 double caloEmEtFraction = calomet->emEtFraction();
2392 double caloHadEtInHB = calomet->hadEtInHB();
2393 double caloHadEtInHO = calomet->hadEtInHO();
2394 double caloHadEtInHE = calomet->hadEtInHE();
2395 double caloHadEtInHF = calomet->hadEtInHF();
2396 double caloEmEtInEB = calomet->emEtInEB();
2397 double caloEmEtInEE = calomet->emEtInEE();
2398 double caloEmEtInHF = calomet->emEtInHF();
2449 float pt_sum_CHF_Barrel = 0;
2450 float pt_sum_CHF_Endcap_plus = 0;
2451 float pt_sum_CHF_Endcap_minus = 0;
2452 float pt_sum_NHF_Barrel = 0;
2453 float pt_sum_NHF_Endcap_plus = 0;
2454 float pt_sum_NHF_Endcap_minus = 0;
2455 float pt_sum_PhF_Barrel = 0;
2456 float pt_sum_PhF_Endcap_plus = 0;
2457 float pt_sum_PhF_Endcap_minus = 0;
2458 float pt_sum_HFH_plus = 0;
2459 float pt_sum_HFH_minus = 0;
2460 float pt_sum_HFE_plus = 0;
2461 float pt_sum_HFE_minus = 0;
2463 float px_chargedHadronsBarrel = 0;
2464 float py_chargedHadronsBarrel = 0;
2465 float px_chargedHadronsEndcapPlus = 0;
2466 float py_chargedHadronsEndcapPlus = 0;
2467 float px_chargedHadronsEndcapMinus = 0;
2468 float py_chargedHadronsEndcapMinus = 0;
2469 float px_neutralHadronsBarrel = 0;
2470 float py_neutralHadronsBarrel = 0;
2471 float px_neutralHadronsEndcapPlus = 0;
2472 float py_neutralHadronsEndcapPlus = 0;
2473 float px_neutralHadronsEndcapMinus = 0;
2474 float py_neutralHadronsEndcapMinus = 0;
2475 float px_PhotonsBarrel = 0;
2476 float py_PhotonsBarrel = 0;
2477 float px_PhotonsEndcapPlus = 0;
2478 float py_PhotonsEndcapPlus = 0;
2479 float px_PhotonsEndcapMinus = 0;
2480 float py_PhotonsEndcapMinus = 0;
2481 float px_HFHadronsPlus = 0;
2482 float py_HFHadronsPlus = 0;
2483 float px_HFHadronsMinus = 0;
2484 float py_HFHadronsMinus = 0;
2485 float px_HFEGammasPlus = 0;
2486 float py_HFEGammasPlus = 0;
2487 float px_HFEGammasMinus = 0;
2488 float py_HFEGammasMinus = 0;
2491 if (
c.particleId() == 1) {
2493 if (
c.eta() > (-3.0) &&
c.eta() < (-1.392)) {
2494 px_chargedHadronsEndcapMinus -=
c.px();
2495 py_chargedHadronsEndcapMinus -=
c.py();
2496 pt_sum_CHF_Endcap_minus +=
c.et();
2497 }
else if (
c.eta() >= (-1.392) &&
c.eta() <= 1.392) {
2498 px_chargedHadronsBarrel -=
c.px();
2499 py_chargedHadronsBarrel -=
c.py();
2500 pt_sum_CHF_Barrel +=
c.et();
2501 }
else if (
c.eta() > 1.392 &&
c.eta() < 3.0) {
2502 px_chargedHadronsEndcapPlus -=
c.px();
2503 py_chargedHadronsEndcapPlus -=
c.py();
2504 pt_sum_CHF_Endcap_plus +=
c.et();
2507 if (
c.particleId() == 5) {
2509 if (
c.eta() > (-3.0) &&
c.eta() < (-1.392)) {
2510 px_neutralHadronsEndcapMinus -=
c.px();
2511 py_neutralHadronsEndcapMinus -=
c.py();
2512 pt_sum_NHF_Endcap_minus +=
c.et();
2513 }
else if (
c.eta() >= (-1.392) &&
c.eta() <= 1.392) {
2514 px_neutralHadronsBarrel -=
c.px();
2515 py_neutralHadronsBarrel -=
c.py();
2516 pt_sum_NHF_Barrel +=
c.et();
2517 }
else if (
c.eta() > 1.392 &&
c.eta() < 3.0) {
2518 px_neutralHadronsEndcapPlus -=
c.px();
2519 py_neutralHadronsEndcapPlus -=
c.py();
2520 pt_sum_NHF_Endcap_plus +=
c.et();
2523 if (
c.particleId() == 4) {
2525 if (
c.eta() > (-3.0) &&
c.eta() < (-1.479)) {
2526 px_PhotonsEndcapMinus -=
c.px();
2527 py_PhotonsEndcapMinus -=
c.py();
2528 pt_sum_PhF_Endcap_minus +=
c.et();
2529 }
else if (
c.eta() >= (-1.479) &&
c.eta() <= 1.479) {
2530 px_PhotonsBarrel -=
c.px();
2531 py_PhotonsBarrel -=
c.py();
2532 pt_sum_PhF_Barrel +=
c.et();
2533 }
else if (
c.eta() > 1.479 &&
c.eta() < 3.0) {
2534 px_PhotonsEndcapPlus -=
c.px();
2535 py_PhotonsEndcapPlus -=
c.py();
2536 pt_sum_PhF_Endcap_plus +=
c.et();
2539 if (
c.particleId() == 6) {
2541 if (
c.eta() > (-5.20) &&
c.eta() < -2.901376) {
2542 pt_sum_HFH_minus +=
c.et();
2543 px_HFHadronsMinus -=
c.px();
2544 py_HFHadronsMinus -=
c.py();
2545 }
else if (
c.eta() > 2.901376 &&
c.eta() < 5.20) {
2546 px_HFHadronsPlus -=
c.px();
2547 py_HFHadronsPlus -=
c.py();
2548 pt_sum_HFH_plus +=
c.et();
2551 if (
c.particleId() == 7) {
2553 if (
c.eta() > (-5.20) &&
c.eta() < -2.901376) {
2554 pt_sum_HFE_minus +=
c.et();
2555 px_HFEGammasMinus -=
c.px();
2556 py_HFEGammasMinus -=
c.py();
2557 }
else if (
c.eta() > 2.901376 &&
c.eta() < 5.20) {
2558 px_HFEGammasPlus -=
c.px();
2559 py_HFEGammasPlus -=
c.py();
2560 pt_sum_HFE_plus +=
c.et();
2626 if (pt_sum_CHF_Barrel) {
2631 if (pt_sum_CHF_Endcap_plus) {
2636 if (pt_sum_CHF_Endcap_minus) {
2641 if (pt_sum_NHF_Barrel) {
2646 if (pt_sum_NHF_Endcap_plus) {
2651 if (pt_sum_NHF_Endcap_minus) {
2656 if (pt_sum_PhF_Barrel) {
2661 if (pt_sum_PhF_Endcap_plus) {
2666 if (pt_sum_PhF_Endcap_minus) {
2671 if (pt_sum_HFH_plus) {
2676 if (pt_sum_HFH_minus) {
2681 if (pt_sum_HFE_plus) {
2686 if (pt_sum_HFE_minus) {
2751 if (techTriggerCase[1]) {
2806 if (pt_sum_CHF_Barrel) {
2808 map_of_MEs[
DirName +
"/" +
"METPhiChargedHadronsBarrel_BXm1Filled"];
2812 if (pt_sum_CHF_Endcap_plus) {
2814 map_of_MEs[
DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Filled"];
2818 atan2(py_chargedHadronsEndcapPlus, px_chargedHadronsEndcapPlus));
2820 if (pt_sum_CHF_Endcap_minus) {
2822 map_of_MEs[
DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Filled"];
2826 atan2(py_chargedHadronsEndcapMinus, px_chargedHadronsEndcapMinus));
2828 if (pt_sum_NHF_Barrel) {
2830 map_of_MEs[
DirName +
"/" +
"METPhiNeutralHadronsBarrel_BXm1Filled"];
2834 if (pt_sum_NHF_Endcap_plus) {
2836 map_of_MEs[
DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Filled"];
2840 atan2(py_neutralHadronsEndcapPlus, px_neutralHadronsEndcapPlus));
2842 if (pt_sum_NHF_Endcap_minus) {
2844 map_of_MEs[
DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Filled"];
2848 atan2(py_neutralHadronsEndcapMinus, px_neutralHadronsEndcapMinus));
2850 if (pt_sum_PhF_Barrel) {
2855 if (pt_sum_PhF_Endcap_plus) {
2860 if (pt_sum_PhF_Endcap_minus) {
2865 if (pt_sum_HFH_plus) {
2870 if (pt_sum_HFH_minus) {
2875 if (pt_sum_HFE_plus) {
2880 if (pt_sum_HFE_minus) {
2945 if (techTriggerCase[2]) {
3000 if (pt_sum_CHF_Barrel) {
3005 if (pt_sum_CHF_Endcap_plus) {
3007 map_of_MEs[
DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Empty"];
3011 atan2(py_chargedHadronsEndcapPlus, px_chargedHadronsEndcapPlus));
3013 if (pt_sum_CHF_Endcap_minus) {
3015 map_of_MEs[
DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Empty"];
3019 atan2(py_chargedHadronsEndcapMinus, px_chargedHadronsEndcapMinus));
3021 if (pt_sum_NHF_Barrel) {
3026 if (pt_sum_NHF_Endcap_plus) {
3028 map_of_MEs[
DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Empty"];
3032 atan2(py_neutralHadronsEndcapPlus, px_neutralHadronsEndcapPlus));
3034 if (pt_sum_NHF_Endcap_minus) {
3036 map_of_MEs[
DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Empty"];
3040 atan2(py_neutralHadronsEndcapMinus, px_neutralHadronsEndcapMinus));
3042 if (pt_sum_PhF_Barrel) {
3047 if (pt_sum_PhF_Endcap_plus) {
3052 if (pt_sum_PhF_Endcap_minus) {
3057 if (pt_sum_HFH_plus) {
3062 if (pt_sum_HFH_minus) {
3067 if (pt_sum_HFE_plus) {
3072 if (pt_sum_HFE_minus) {
3083 double pfPhotonEtFraction = pfmet->photonEtFraction();
3084 double pfPhotonEt = pfmet->photonEt();
3085 double pfNeutralHadronEtFraction = pfmet->neutralHadronEtFraction();
3086 double pfNeutralHadronEt = pfmet->neutralHadronEt();
3087 double pfElectronEt = pfmet->electronEt();
3088 double pfChargedHadronEtFraction = pfmet->chargedHadronEtFraction();
3089 double pfChargedHadronEt = pfmet->chargedHadronEt();
3090 double pfMuonEt = pfmet->muonEt();
3091 double pfHFHadronEtFraction = pfmet->HFHadronEtFraction();
3092 double pfHFHadronEt = pfmet->HFHadronEt();
3093 double pfHFEMEtFraction = pfmet->HFEMEtFraction();
3094 double pfHFEMEt = pfmet->HFEMEt();
3221 meHFEMEt->
Fill(patmet->Type7EtFraction() * patmet->sumEt());
3223 meMuonEt->
Fill(patmet->MuonEtFraction() * patmet->sumEt());
References HltBtagPostValidation_cff::c, pat::MET::ChargedEMEtFraction(), pat::MET::ChargedHadEtFraction(), reco::PFMET::chargedHadronEt(), reco::PFMET::chargedHadronEtFraction(), fed_dqm_sourceclient-live_cfg::DirName, reco::PFMET::electronEt(), reco::CaloMET::emEtFraction(), reco::CaloMET::emEtInEB(), reco::CaloMET::emEtInEE(), reco::CaloMET::emEtInHF(), reco::CaloMET::etFractionHadronic(), HcalObjRepresent::Fill(), reco::CaloMET::hadEtInHB(), reco::CaloMET::hadEtInHE(), reco::CaloMET::hadEtInHF(), reco::CaloMET::hadEtInHO(), reco::PFMET::HFEMEt(), reco::PFMET::HFEMEtFraction(), reco::PFMET::HFHadronEt(), reco::PFMET::HFHadronEtFraction(), mps_fire::i, iEvent, dqmiolumiharvest::j, BTaggingMonitor_cfi::met, reco::PFMET::muonEt(), pat::MET::MuonEtFraction(), pat::MET::NeutralEMFraction(), pat::MET::NeutralHadEtFraction(), reco::PFMET::neutralHadronEt(), reco::PFMET::neutralHadronEtFraction(), pfLinker_cff::particleFlow, reco::PFMET::photonEt(), reco::PFMET::photonEtFraction(), reco::LeafCandidate::pt(), mathSSE::sqrt(), reco::MET::sumEt(), pat::MET::Type6EtFraction(), and pat::MET::Type7EtFraction().
◆ makeRatePlot()
void METAnalyzer::makeRatePlot |
( |
std::string |
DirName, |
|
|
double |
totltime |
|
) |
| |
◆ allTriggerDecisions_
std::vector<int> METAnalyzer::allTriggerDecisions_ |
|
private |
◆ allTriggerNames_
std::vector<std::string> METAnalyzer::allTriggerNames_ |
|
private |
◆ beamSpot_
◆ bypassAllDCSChecks_
bool METAnalyzer::bypassAllDCSChecks_ |
|
private |
◆ bypassAllPVChecks_
bool METAnalyzer::bypassAllPVChecks_ |
|
private |
◆ caloJetsToken_
◆ caloMetToken_
◆ cleaningParameters_
◆ countsPFCand_
std::vector<int> METAnalyzer::countsPFCand_ |
|
private |
◆ CSCHalo2015ResultTag_
◆ CSCHalo2015ResultToken_
◆ CSCHaloResultTag_
◆ CSCHaloResultToken_
◆ DCSFilter_
◆ diagnosticsParameters_
◆ EcalDeadCellBoundaryTag_
◆ EcalDeadCellBoundaryToken_
◆ EcalDeadCellTriggerTag_
◆ EcalDeadCellTriggerToken_
◆ eeBadScFilterTag_
◆ eeBadScFilterToken_
◆ etaMaxPFCand_
std::vector<double> METAnalyzer::etaMaxPFCand_ |
|
private |
◆ etaMinPFCand_
std::vector<double> METAnalyzer::etaMinPFCand_ |
|
private |
◆ etaNBinsPFCand_
std::vector<int> METAnalyzer::etaNBinsPFCand_ |
|
private |
◆ fill_met_high_level_histo
bool METAnalyzer::fill_met_high_level_histo |
|
private |
◆ fillCandidateMap_histos
bool METAnalyzer::fillCandidateMap_histos |
|
private |
◆ FilterhltConfig_
◆ FolderName_
std::string METAnalyzer::FolderName_ |
|
private |
◆ folderNames_
std::vector<std::string> METAnalyzer::folderNames_ |
|
private |
◆ gtTag_
◆ gtToken_
◆ hbheIsoNoiseFilterResultToken_
◆ HBHEIsoNoiseStringMiniAOD
std::string METAnalyzer::HBHEIsoNoiseStringMiniAOD |
|
private |
◆ hbheNoiseFilterResultTag_
◆ hbheNoiseFilterResultToken_
◆ hbheNoiseIsoFilterResultTag_
◆ HBHENoiseStringMiniAOD
std::string METAnalyzer::HBHENoiseStringMiniAOD |
|
private |
◆ hcalNoiseRBXCollectionTag_
◆ hCaloEmEtFraction
◆ hCaloEmEtFraction020
◆ hCaloEmEtInEB
◆ hCaloEmEtInEE
◆ hCaloEmEtInHF
◆ hCaloEmMET
◆ hCaloEmMETPhi
◆ hCaloEmMEx
◆ hCaloEmMEy
◆ hCaloEtFractionHadronic
◆ hCaloHadEtInHB
◆ hCaloHadEtInHE
◆ hCaloHadEtInHF
◆ hCaloHadEtInHO
◆ hCaloHaMET
◆ hCaloHaMETPhi
◆ hCaloHaMEx
◆ hCaloHaMEy
◆ hCaloMETPhi020
◆ HcalStripHaloTag_
◆ HcalStripHaloToken_
◆ hfCalibFactor_
double METAnalyzer::hfCalibFactor_ |
|
private |
◆ hltConfig_
◆ hltPhysDec_
std::string METAnalyzer::hltPhysDec_ |
|
private |
◆ hMET
◆ hMET_2
◆ hMET_CSCTightHalo2015Filter
◆ hMET_CSCTightHaloFilter
◆ hMET_EcalDeadCellBoundaryFilter
◆ hMET_EcalDeadCellTriggerFilter
◆ hMET_eeBadScFilter
◆ hMET_HBHEIsoNoiseFilter
◆ hMET_HBHENoiseFilter
◆ hMET_HcalStripHaloFilter
◆ hMET_logx
◆ hMETPhi
◆ hMETRate
◆ hMETSig
◆ hMEx
◆ hMExLS
◆ hMEy
◆ hMEyLS
◆ hSumET
◆ hSumET_logx
◆ hTrigger
◆ hTriggerLabelsIsSet_
bool METAnalyzer::hTriggerLabelsIsSet_ |
|
private |
◆ inputJetIDValueMap
◆ isCaloMet_
bool METAnalyzer::isCaloMet_ |
|
private |
◆ isMiniAODMet_
bool METAnalyzer::isMiniAODMet_ |
|
private |
◆ isPFMet_
bool METAnalyzer::isPFMet_ |
|
private |
◆ jetCollectionLabel_
◆ jetCorrectorToken_
◆ jetID_ValueMapToken_
◆ jetIDFunctorLoose
◆ LSBegin_
int METAnalyzer::LSBegin_ |
|
private |
◆ LSEnd_
◆ lumisecME
◆ m_bitAlgTechTrig_
int METAnalyzer::m_bitAlgTechTrig_ |
|
private |
◆ m_l1algoname_
std::string METAnalyzer::m_l1algoname_ |
|
private |
◆ map_dijet_MEs
◆ meChargedHadronEt
◆ meChargedHadronEt_profile
◆ meChargedHadronEtFraction
◆ meChargedHadronEtFraction_BXm1Empty
◆ meChargedHadronEtFraction_BXm1Filled
◆ meChargedHadronEtFraction_profile
◆ meCHF_Barrel
◆ meCHF_Barrel_BXm1Empty
◆ meCHF_Barrel_BXm1Filled
◆ meCHF_EndcapMinus
◆ meCHF_EndcapMinus_BXm1Empty
◆ meCHF_EndcapMinus_BXm1Filled
◆ meCHF_EndcapPlus
◆ meCHF_EndcapPlus_BXm1Empty
◆ meCHF_EndcapPlus_BXm1Filled
◆ meElectronEt
◆ meHFEMEt
◆ meHFEMEt_profile
◆ meHFEMEtFraction
◆ meHFEMEtFraction_profile
◆ meHFEMF_Minus
◆ meHFEMF_Minus_BXm1Empty
◆ meHFEMF_Minus_BXm1Filled
◆ meHFEMF_Plus
◆ meHFEMF_Plus_BXm1Empty
◆ meHFEMF_Plus_BXm1Filled
◆ meHFHadF_Minus
◆ meHFHadF_Minus_BXm1Empty
◆ meHFHadF_Minus_BXm1Filled
◆ meHFHadF_Plus
◆ meHFHadF_Plus_BXm1Empty
◆ meHFHadF_Plus_BXm1Filled
◆ meHFHadronEt
◆ meHFHadronEt_profile
◆ meHFHadronEtFraction
◆ meHFHadronEtFraction_profile
◆ meMET_BXm1Empty
◆ meMET_BXm1Filled
◆ meMET_profile
◆ meMETPhiChargedHadronsBarrel
◆ meMETPhiChargedHadronsBarrel_BXm1Empty
◆ meMETPhiChargedHadronsBarrel_BXm1Filled
◆ meMETPhiChargedHadronsEndcapMinus
◆ meMETPhiChargedHadronsEndcapMinus_BXm1Empty
MonitorElement* METAnalyzer::meMETPhiChargedHadronsEndcapMinus_BXm1Empty |
|
private |
◆ meMETPhiChargedHadronsEndcapMinus_BXm1Filled
MonitorElement* METAnalyzer::meMETPhiChargedHadronsEndcapMinus_BXm1Filled |
|
private |
◆ meMETPhiChargedHadronsEndcapPlus
◆ meMETPhiChargedHadronsEndcapPlus_BXm1Empty
MonitorElement* METAnalyzer::meMETPhiChargedHadronsEndcapPlus_BXm1Empty |
|
private |
◆ meMETPhiChargedHadronsEndcapPlus_BXm1Filled
MonitorElement* METAnalyzer::meMETPhiChargedHadronsEndcapPlus_BXm1Filled |
|
private |
◆ meMETPhiHFEGammasMinus
◆ meMETPhiHFEGammasMinus_BXm1Empty
◆ meMETPhiHFEGammasMinus_BXm1Filled
◆ meMETPhiHFEGammasPlus
◆ meMETPhiHFEGammasPlus_BXm1Empty
◆ meMETPhiHFEGammasPlus_BXm1Filled
◆ meMETPhiHFHadronsMinus
◆ meMETPhiHFHadronsMinus_BXm1Empty
◆ meMETPhiHFHadronsMinus_BXm1Filled
◆ meMETPhiHFHadronsPlus
◆ meMETPhiHFHadronsPlus_BXm1Empty
◆ meMETPhiHFHadronsPlus_BXm1Filled
◆ meMETPhiNeutralHadronsBarrel
◆ meMETPhiNeutralHadronsBarrel_BXm1Empty
◆ meMETPhiNeutralHadronsBarrel_BXm1Filled
◆ meMETPhiNeutralHadronsEndcapMinus
◆ meMETPhiNeutralHadronsEndcapMinus_BXm1Empty
MonitorElement* METAnalyzer::meMETPhiNeutralHadronsEndcapMinus_BXm1Empty |
|
private |
◆ meMETPhiNeutralHadronsEndcapMinus_BXm1Filled
MonitorElement* METAnalyzer::meMETPhiNeutralHadronsEndcapMinus_BXm1Filled |
|
private |
◆ meMETPhiNeutralHadronsEndcapPlus
◆ meMETPhiNeutralHadronsEndcapPlus_BXm1Empty
MonitorElement* METAnalyzer::meMETPhiNeutralHadronsEndcapPlus_BXm1Empty |
|
private |
◆ meMETPhiNeutralHadronsEndcapPlus_BXm1Filled
MonitorElement* METAnalyzer::meMETPhiNeutralHadronsEndcapPlus_BXm1Filled |
|
private |
◆ meMETPhiPhotonsBarrel
◆ meMETPhiPhotonsBarrel_BXm1Empty
◆ meMETPhiPhotonsBarrel_BXm1Filled
◆ meMETPhiPhotonsEndcapMinus
◆ meMETPhiPhotonsEndcapMinus_BXm1Empty
◆ meMETPhiPhotonsEndcapMinus_BXm1Filled
◆ meMETPhiPhotonsEndcapPlus
◆ meMETPhiPhotonsEndcapPlus_BXm1Empty
◆ meMETPhiPhotonsEndcapPlus_BXm1Filled
◆ meMEx_profile
◆ meMEy_profile
◆ meMuonEt
◆ meNeutralHadronEt
◆ meNeutralHadronEt_profile
◆ meNeutralHadronEtFraction
◆ meNeutralHadronEtFraction_BXm1Empty
◆ meNeutralHadronEtFraction_BXm1Filled
◆ meNeutralHadronEtFraction_profile
◆ meNHF_Barrel
◆ meNHF_Barrel_BXm1Empty
◆ meNHF_Barrel_BXm1Filled
◆ meNHF_EndcapMinus
◆ meNHF_EndcapMinus_BXm1Empty
◆ meNHF_EndcapMinus_BXm1Filled
◆ meNHF_EndcapPlus
◆ meNHF_EndcapPlus_BXm1Empty
◆ meNHF_EndcapPlus_BXm1Filled
◆ mePhF_Barrel
◆ mePhF_Barrel_BXm1Empty
◆ mePhF_Barrel_BXm1Filled
◆ mePhF_EndcapMinus
◆ mePhF_EndcapMinus_BXm1Empty
◆ mePhF_EndcapMinus_BXm1Filled
◆ mePhF_EndcapPlus
◆ mePhF_EndcapPlus_BXm1Empty
◆ mePhF_EndcapPlus_BXm1Filled
◆ mePhotonEt
◆ mePhotonEt_profile
◆ mePhotonEtFraction
◆ mePhotonEtFraction_BXm1Empty
◆ mePhotonEtFraction_BXm1Filled
◆ mePhotonEtFraction_profile
◆ meSumET_BXm1Empty
◆ meSumET_BXm1Filled
◆ meSumET_profile
◆ metCollectionLabel_
◆ METFilterMiniAODLabel2_
◆ METFilterMiniAODLabel_
◆ METFilterMiniAODToken2_
◆ METFilterMiniAODToken_
◆ MetType_
std::string METAnalyzer::MetType_ |
|
private |
◆ MExPFCand_
std::vector<double> METAnalyzer::MExPFCand_ |
|
private |
◆ MEyPFCand_
std::vector<double> METAnalyzer::MEyPFCand_ |
|
private |
◆ meZJets_u_par
◆ meZJets_u_par_ZPt_0_15
◆ meZJets_u_par_ZPt_150_290
◆ meZJets_u_par_ZPt_15_30
◆ meZJets_u_par_ZPt_290
◆ meZJets_u_par_ZPt_30_55
◆ meZJets_u_par_ZPt_55_75
◆ meZJets_u_par_ZPt_75_150
◆ meZJets_u_perp
◆ meZJets_u_perp_ZPt_0_15
◆ meZJets_u_perp_ZPt_150_290
◆ meZJets_u_perp_ZPt_15_30
◆ meZJets_u_perp_ZPt_290
◆ meZJets_u_perp_ZPt_30_55
◆ meZJets_u_perp_ZPt_55_75
◆ meZJets_u_perp_ZPt_75_150
◆ miniaodfilterdec
int METAnalyzer::miniaodfilterdec |
|
private |
◆ miniaodFilterIndex_
std::vector<int> METAnalyzer::miniaodFilterIndex_ |
|
private |
◆ MuonsToken_
◆ nbinsPV_
int METAnalyzer::nbinsPV_ |
|
private |
◆ nPVMax_
double METAnalyzer::nPVMax_ |
|
private |
◆ nPVMin_
double METAnalyzer::nPVMin_ |
|
private |
◆ numPV_
◆ onlyCleaned_
bool METAnalyzer::onlyCleaned_ |
|
private |
◆ parameters
◆ patJetsToken_
◆ patMetToken_
◆ pfjetIDFunctorLoose
◆ pfJetsToken_
◆ pflowPackedToken_
◆ pflowToken_
◆ pfMetToken_
◆ profilePFCand_x_
◆ profilePFCand_x_name_
std::vector<std::string> METAnalyzer::profilePFCand_x_name_ |
|
private |
◆ profilePFCand_y_
◆ profilePFCand_y_name_
std::vector<std::string> METAnalyzer::profilePFCand_y_name_ |
|
private |
◆ ptMinCand_
double METAnalyzer::ptMinCand_ |
|
private |
◆ ptThreshold_
double METAnalyzer::ptThreshold_ |
|
private |
◆ runcosmics_
bool METAnalyzer::runcosmics_ |
|
private |
◆ triggerFolderDecisions_
std::vector<int> METAnalyzer::triggerFolderDecisions_ |
|
private |
◆ triggerFolderEventFlag_
◆ triggerFolderExpr_
std::vector<std::vector<std::string> > METAnalyzer::triggerFolderExpr_ |
|
private |
◆ triggerFolderLabels_
std::vector<std::string> METAnalyzer::triggerFolderLabels_ |
|
private |
◆ triggerResultsLabel_
◆ triggerResultsToken_
◆ triggerSelectedSubFolders_
◆ typePFCand_
std::vector<int> METAnalyzer::typePFCand_ |
|
private |
◆ verbose_
int METAnalyzer::verbose_ |
|
private |
◆ vertexTag_
◆ vertexToken_
T const * product() const
bool filter(const edm::Event &evt, const edm::EventSetup &es)
T const * product() const
virtual MonitorElementData::Scope setScope(MonitorElementData::Scope newscope)
virtual void setOption(const char *option)
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
numberOfValidPixelHits
at most 999 lost hits
std::vector< Vertex > VertexCollection
collection of Vertex objects
static const std::string triggerResults
virtual void setCurrentFolder(std::string const &fullpath)
std::vector< ParameterSet > VParameterSet
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
double correction(const LorentzVector &fJet) const
get correction using Jet information only
T getUntrackedParameter(std::string const &, T const &) const
virtual TH1F * getTH1F() const
Log< level::Warning, false > LogWarning
PF Jet selector for pat::Jets.
TObject * getRootObject() const override
Jet selector for pat::Jets and for CaloJets.
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
std::vector< T >::const_iterator search(const cond::Time_t &val, const std::vector< T > &container)
unsigned int size() const
number of trigger paths in trigger table
pat::strbitset getBitTemplate() const
Get an empty bitset with the proper names.
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
static int position[264][3]
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Analysis-level MET class.
UseScope< MonitorElementData::Scope::LUMI > UseLumiScope
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
virtual const TObject * getRootObject() const
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
virtual int getNbinsX() const
get # of bins in X-axis
T getParameter(std::string const &) const
bool accept() const
Has at least one path accepted the event?
Particle reconstructed by the particle flow algorithm.
const std::string & triggerName(unsigned int triggerIndex) const
const TechnicalTriggerWord & technicalTriggerWord(int bxInEventValue) const
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())