CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Public Attributes | Private Attributes
PFMETAnalyzer Class Reference

#include <PFMETAnalyzer.h>

Inheritance diagram for PFMETAnalyzer:
PFMETAnalyzerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &)
 Get the analysis. More...
 
void beginJob (DQMStore *dbe)
 Inizialize parameters for histo binning. More...
 
void beginRun (const edm::Run &, const edm::EventSetup &)
 Initialize run-based parameters. More...
 
void bookMESet (std::string)
 
void bookMonitorElement (std::string, bool)
 
void endJob ()
 Finish up a job. More...
 
void endRun (const edm::Run &iRun, const edm::EventSetup &iSetup, DQMStore *dbe)
 Finish up a run. More...
 
void fillMESet (const edm::Event &, std::string, const reco::PFMET &)
 
void fillMonitorElement (const edm::Event &, std::string, std::string, const reco::PFMET &, bool)
 
void makeRatePlot (std::string, double)
 
 PFMETAnalyzer (const edm::ParameterSet &)
 Constructor. More...
 
bool selectHighPtJetEvent (const edm::Event &)
 
bool selectLowPtJetEvent (const edm::Event &)
 
bool selectWElectronEvent (const edm::Event &)
 
bool selectWMuonEvent (const edm::Event &)
 
void setSource (std::string source)
 
void validateMET (const reco::PFMET &, edm::Handle< edm::View< reco::PFCandidate > >)
 
virtual ~PFMETAnalyzer ()
 Destructor. More...
 
- Public Member Functions inherited from PFMETAnalyzerBase
void analyze (const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &, reco::PFMET &pfMET)
 Get the analysis of the muon properties. More...
 
 PFMETAnalyzerBase ()
 Constructor. More...
 
virtual ~PFMETAnalyzerBase ()
 Destructor. More...
 

Public Attributes

int evtCounter
 

Private Attributes

bool _allhist
 
bool _allSelection
 
bool _cleanupSelection
 
DQMStore * _dbe
 
bool _doHLTPhysicsOn
 
bool _doPVCheck
 
GenericTriggerEventFlag * _EleEventFlag
 
double _etThreshold
 
std::vector< std::string > _FolderNames
 
GenericTriggerEventFlag * _HighMETEventFlag
 
double _highPFMETThreshold
 
GenericTriggerEventFlag * _HighPtJetEventFlag
 
double _highPtPFJetThreshold
 
std::string _hlt_PhysDec
 
double _lowPFMETThreshold
 
GenericTriggerEventFlag * _LowPtJetEventFlag
 
double _lowPtPFJetThreshold
 
GenericTriggerEventFlag * _MinBiasEventFlag
 
GenericTriggerEventFlag * _MuonEventFlag
 
int _numPV
 
int _nvtx_min
 
int _nvtxtrks_min
 
std::string _source
 
std::vector< unsigned > _techTrigsAND
 
std::vector< unsigned > _techTrigsNOT
 
std::vector< unsigned > _techTrigsOR
 
bool _tightBHFiltering
 
int _tightJetIDFiltering
 
int _trig_Ele
 
int _trig_HighMET
 
int _trig_HighPtJet
 
int _trig_JetMB
 
int _trig_LowPtJet
 
int _trig_MinBias
 
int _trig_Muon
 
int _trig_PhysDec
 
int _verbose
 
double _vtxchi2_max
 
int _vtxndof_min
 
double _vtxz_max
 
edm::InputTag BeamHaloSummaryTag
 
JetMETDQMDCSFilter * DCSFilter
 
std::vector< std::string > elecExpr_
 
edm::InputTag gtTag
 
edm::InputTag HBHENoiseFilterResultTag
 
edm::InputTag HcalNoiseRBXCollectionTag
 
double hfCalibFactor_
 
std::vector< std::string > highMETExpr_
 
std::vector< std::string > highPtJetExpr_
 
std::vector< std::string > HLTPathsJetMBByName_
 
reco::helper::JetIDHelper * jetID
 
std::vector< std::string > lowPtJetExpr_
 
MonitorElement * meChargedHadronEt
 
MonitorElement * meChargedHadronEt_profile
 
MonitorElement * meChargedHadronEtFraction
 
MonitorElement * meChargedHadronEtFraction_profile
 
MonitorElement * meElectronEt
 
MonitorElement * meElectronEt_profile
 
MonitorElement * meElectronEtFraction
 
MonitorElement * meElectronEtFraction_profile
 
MonitorElement * meHFEMEt
 
MonitorElement * meHFEMEt_profile
 
MonitorElement * meHFEMEtFraction
 
MonitorElement * meHFEMEtFraction_profile
 
MonitorElement * meHFHadronEt
 
MonitorElement * meHFHadronEt_profile
 
MonitorElement * meHFHadronEtFraction
 
MonitorElement * meHFHadronEtFraction_profile
 
MonitorElement * meMuonEt
 
MonitorElement * meMuonEt_profile
 
MonitorElement * meMuonEtFraction
 
MonitorElement * meMuonEtFraction_profile
 
MonitorElement * meNeutralHadronEt
 
MonitorElement * meNeutralHadronEt_profile
 
MonitorElement * meNeutralHadronEtFraction
 
MonitorElement * meNeutralHadronEtFraction_profile
 
MonitorElement * mePfMET
 
MonitorElement * mePfMET_logx
 
MonitorElement * mePfMET_profile
 
MonitorElement * mePfMETPhi
 
MonitorElement * mePfMETRate
 
MonitorElement * mePfMETSig
 
MonitorElement * mePfMEx
 
MonitorElement * mePfMEx_profile
 
MonitorElement * mePfMExLS
 
MonitorElement * mePfMEy
 
MonitorElement * mePfMEy_profile
 
MonitorElement * mePfMEyLS
 
MonitorElement * mePfSumET
 
MonitorElement * mePfSumET_logx
 
MonitorElement * mePfSumET_profile
 
MonitorElement * mePhotonEt
 
MonitorElement * mePhotonEt_profile
 
MonitorElement * mePhotonEtFraction
 
MonitorElement * mePhotonEtFraction_profile
 
MonitorElement * metME
 
std::string metname
 
MonitorElement * meTriggerName_Ele
 
MonitorElement * meTriggerName_HighMET
 
MonitorElement * meTriggerName_HighPtJet
 
MonitorElement * meTriggerName_LowPtJet
 
MonitorElement * meTriggerName_MinBias
 
MonitorElement * meTriggerName_Muon
 
MonitorElement * meTriggerName_PhysDec
 
std::vector< std::string > minbiasExpr_
 
std::vector< std::string > muonExpr_
 
edm::ParameterSet parameters
 
edm::InputTag PFCandidatesTag
 
edm::ParameterSet theCleaningParameters
 
edm::InputTag theJetCollectionLabel
 
edm::InputTag thePfJetCollectionLabel
 
edm::InputTag thePfMETCollectionLabel
 
edm::InputTag vertexTag
 

Detailed Description

DQM monitoring source for PFMET

Date:
2012/05/20 17:57:01
Revision:
1.30
Author
K. Hatakeyama - Rockefeller University A.Apresyan - Caltech

Definition at line 51 of file PFMETAnalyzer.h.

Constructor & Destructor Documentation

PFMETAnalyzer::PFMETAnalyzer ( const edm::ParameterSet &  pSet)

Constructor.

Definition at line 37 of file PFMETAnalyzer.cc.

References edm::ParameterSet::getParameter(), and Parameters::parameters.

37  {
38 
39  parameters = pSet;
40 
41  edm::ParameterSet highptjetparms = parameters.getParameter<edm::ParameterSet>("highPtJetTrigger");
42  edm::ParameterSet lowptjetparms = parameters.getParameter<edm::ParameterSet>("lowPtJetTrigger" );
43  edm::ParameterSet minbiasparms = parameters.getParameter<edm::ParameterSet>("minBiasTrigger" );
44  edm::ParameterSet highmetparms = parameters.getParameter<edm::ParameterSet>("highMETTrigger" );
45  // edm::ParameterSet lowmetparms = parameters.getParameter<edm::ParameterSet>("lowMETTrigger" );
46  edm::ParameterSet eleparms = parameters.getParameter<edm::ParameterSet>("eleTrigger" );
47  edm::ParameterSet muonparms = parameters.getParameter<edm::ParameterSet>("muonTrigger" );
48 
49  //genericTriggerEventFlag_( new GenericTriggerEventFlag( conf_ ) );
50  _HighPtJetEventFlag = new GenericTriggerEventFlag( highptjetparms );
51  _LowPtJetEventFlag = new GenericTriggerEventFlag( lowptjetparms );
52  _MinBiasEventFlag = new GenericTriggerEventFlag( minbiasparms );
53  _HighMETEventFlag = new GenericTriggerEventFlag( highmetparms );
54  // _LowMETEventFlag = new GenericTriggerEventFlag( lowmetparms );
55  _EleEventFlag = new GenericTriggerEventFlag( eleparms );
56  _MuonEventFlag = new GenericTriggerEventFlag( muonparms );
57 
58  highPtJetExpr_ = highptjetparms.getParameter<std::vector<std::string> >("hltPaths");
59  lowPtJetExpr_ = lowptjetparms .getParameter<std::vector<std::string> >("hltPaths");
60  highMETExpr_ = highmetparms .getParameter<std::vector<std::string> >("hltPaths");
61  // lowMETExpr_ = lowmetparms .getParameter<std::vector<std::string> >("hltPaths");
62  muonExpr_ = muonparms .getParameter<std::vector<std::string> >("hltPaths");
63  elecExpr_ = eleparms .getParameter<std::vector<std::string> >("hltPaths");
64  minbiasExpr_ = minbiasparms .getParameter<std::vector<std::string> >("hltPaths");
65 
66 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PFMETAnalyzer::_MuonEventFlag
GenericTriggerEventFlag * _MuonEventFlag
Definition: PFMETAnalyzer.h:126
PFMETAnalyzer::highPtJetExpr_
std::vector< std::string > highPtJetExpr_
Definition: PFMETAnalyzer.h:128
PFMETAnalyzer::parameters
edm::ParameterSet parameters
Definition: PFMETAnalyzer.h:101
PFMETAnalyzer::_HighPtJetEventFlag
GenericTriggerEventFlag * _HighPtJetEventFlag
Definition: PFMETAnalyzer.h:120
GenericTriggerEventFlag
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Definition: GenericTriggerEventFlag.h:36
PFMETAnalyzer::muonExpr_
std::vector< std::string > muonExpr_
Definition: PFMETAnalyzer.h:132
PFMETAnalyzer::_LowPtJetEventFlag
GenericTriggerEventFlag * _LowPtJetEventFlag
Definition: PFMETAnalyzer.h:121
PFMETAnalyzer::lowPtJetExpr_
std::vector< std::string > lowPtJetExpr_
Definition: PFMETAnalyzer.h:129
PFMETAnalyzer::highMETExpr_
std::vector< std::string > highMETExpr_
Definition: PFMETAnalyzer.h:130
PFMETAnalyzer::_EleEventFlag
GenericTriggerEventFlag * _EleEventFlag
Definition: PFMETAnalyzer.h:125
PFMETAnalyzer::_HighMETEventFlag
GenericTriggerEventFlag * _HighMETEventFlag
Definition: PFMETAnalyzer.h:123
PFMETAnalyzer::elecExpr_
std::vector< std::string > elecExpr_
Definition: PFMETAnalyzer.h:133
PFMETAnalyzer::_MinBiasEventFlag
GenericTriggerEventFlag * _MinBiasEventFlag
Definition: PFMETAnalyzer.h:122
PFMETAnalyzer::minbiasExpr_
std::vector< std::string > minbiasExpr_
Definition: PFMETAnalyzer.h:134
PFMETAnalyzer::~PFMETAnalyzer ( )
virtual

Destructor.

Definition at line 69 of file PFMETAnalyzer.cc.

69  {
70 
71  delete _HighPtJetEventFlag;
72  delete _LowPtJetEventFlag;
73  delete _MinBiasEventFlag;
74  delete _HighMETEventFlag;
75  // delete _LowMETEventFlag;
76  delete _EleEventFlag;
77  delete _MuonEventFlag;
78 
79 }
PFMETAnalyzer::_MuonEventFlag
GenericTriggerEventFlag * _MuonEventFlag
Definition: PFMETAnalyzer.h:126
PFMETAnalyzer::_HighPtJetEventFlag
GenericTriggerEventFlag * _HighPtJetEventFlag
Definition: PFMETAnalyzer.h:120
PFMETAnalyzer::_LowPtJetEventFlag
GenericTriggerEventFlag * _LowPtJetEventFlag
Definition: PFMETAnalyzer.h:121
PFMETAnalyzer::_EleEventFlag
GenericTriggerEventFlag * _EleEventFlag
Definition: PFMETAnalyzer.h:125
PFMETAnalyzer::_HighMETEventFlag
GenericTriggerEventFlag * _HighMETEventFlag
Definition: PFMETAnalyzer.h:123
PFMETAnalyzer::_MinBiasEventFlag
GenericTriggerEventFlag * _MinBiasEventFlag
Definition: PFMETAnalyzer.h:122

Member Function Documentation

void PFMETAnalyzer::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
const edm::TriggerResults &  triggerResults 
)

Get the analysis.

Definition at line 458 of file PFMETAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), FamosSequences_cff::caloJets, gather_cfg::cout, reco::BeamHaloSummary::CSCLooseHaloId(), reco::BeamHaloSummary::CSCTightHaloId(), cmsPerfPublish::DirName, reco::BeamHaloSummary::EcalLooseHaloId(), reco::BeamHaloSummary::EcalTightHaloId(), edm::Event::getByLabel(), reco::BeamHaloSummary::GlobalLooseHaloId(), reco::BeamHaloSummary::GlobalTightHaloId(), reco::BeamHaloSummary::HcalLooseHaloId(), reco::BeamHaloSummary::HcalTightHaloId(), i, edm::HandleBase::isValid(), bTagSequences_cff::jetID, LogDebug, LogTrace, metname, reco::tau::pfCandidates(), pfJets_cff::pfJets, edm::Handle< T >::product(), edm::HLTGlobalStatus::size(), edm::TriggerNames::size(), AlCaHLTBitMon_QueryRunRegistry::string, edm::TriggerNames::triggerIndex(), edm::TriggerNames::triggerName(), edm::Event::triggerNames(), findQualityFiles::v, and GoodVertex_cfg::vertexCollection.

459  {
460 
461  if (_verbose) std::cout << "PfMETAnalyzer analyze" << std::endl;
462 
463  LogTrace(metname)<<"[PFMETAnalyzer] Analyze PFMET";
464 
465  metME->Fill(3);
466 
467  // ==========================================================
468  // Trigger information
469  //
470  _trig_JetMB=0;
471  _trig_HighPtJet=0;
472  _trig_LowPtJet=0;
473  _trig_MinBias=0;
474  _trig_HighMET=0;
475  _trig_Ele=0;
476  _trig_Muon=0;
477  _trig_PhysDec=0;
478  if(&triggerResults) {
479 
481 
482  //
483  //
484  // Check how many HLT triggers are in triggerResults
485  int ntrigs = triggerResults.size();
486  if (_verbose) std::cout << "ntrigs=" << ntrigs << std::endl;
487 
488  //
489  //
490  // If index=ntrigs, this HLT trigger doesn't exist in the HLT table for this data.
491  const edm::TriggerNames & triggerNames = iEvent.triggerNames(triggerResults);
492 
493  //
494  //
495  const unsigned int nTrig(triggerNames.size());
496  for (unsigned int i=0;i<nTrig;++i)
497  {
498  if (triggerNames.triggerName(i).find(highPtJetExpr_[0].substr(0,highPtJetExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
499  _trig_HighPtJet=true;
500  else if (triggerNames.triggerName(i).find(lowPtJetExpr_[0].substr(0,lowPtJetExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
501  _trig_LowPtJet=true;
502  else if (triggerNames.triggerName(i).find(highMETExpr_[0].substr(0,highMETExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
503  _trig_HighMET=true;
504  // else if (triggerNames.triggerName(i).find(lowMETExpr_[0].substr(0,lowMETExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
505  // _trig_LowMET=true;
506  else if (triggerNames.triggerName(i).find(muonExpr_[0].substr(0,muonExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
507  _trig_Muon=true;
508  else if (triggerNames.triggerName(i).find(elecExpr_[0].substr(0,elecExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
509  _trig_Ele=true;
510  else if (triggerNames.triggerName(i).find(minbiasExpr_[0].substr(0,minbiasExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
511  _trig_MinBias=true;
512  }
513 
514  // count number of requested Jet or MB HLT paths which have fired
515  for (unsigned int i=0; i!=HLTPathsJetMBByName_.size(); i++) {
516  unsigned int triggerIndex = triggerNames.triggerIndex(HLTPathsJetMBByName_[i]);
517  if (triggerIndex<triggerResults.size()) {
518  if (triggerResults.accept(triggerIndex)) {
519  _trig_JetMB++;
520  }
521  }
522  }
523  // for empty input vectors (n==0), take all HLT triggers!
524  if (HLTPathsJetMBByName_.size()==0) _trig_JetMB=triggerResults.size()-1;
525 
526 
527  if (triggerNames.triggerIndex(_hlt_PhysDec) != triggerNames.size() &&
528  triggerResults.accept(triggerNames.triggerIndex(_hlt_PhysDec))) _trig_PhysDec=1;
529  } else {
530 
531  edm::LogInfo("PFMetAnalyzer") << "TriggerResults::HLT not found, "
532  "automatically select events";
533 
534  // TriggerResults object not found. Look at all events.
535  _trig_JetMB=1;
536  }
537 
538  // ==========================================================
539  // PfMET information
540 
541  // **** Get the MET container
542  edm::Handle<reco::PFMETCollection> pfmetcoll;
543  iEvent.getByLabel(thePfMETCollectionLabel, pfmetcoll);
544 
545  if(!pfmetcoll.isValid()) return;
546 
547  const PFMETCollection *pfmetcol = pfmetcoll.product();
548  const PFMET *pfmet;
549  pfmet = &(pfmetcol->front());
550 
551  LogTrace(metname)<<"[PfMETAnalyzer] Call to the PfMET analyzer";
552 
553  // ==========================================================
554  //
555  edm::Handle<HcalNoiseRBXCollection> HRBXCollection;
556  iEvent.getByLabel(HcalNoiseRBXCollectionTag,HRBXCollection);
557  if (!HRBXCollection.isValid()) {
558  LogDebug("") << "PfMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
559  if (_verbose) std::cout << "PfMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
560  }
561 
562 
563  edm::Handle<bool> HBHENoiseFilterResultHandle;
564  iEvent.getByLabel(HBHENoiseFilterResultTag, HBHENoiseFilterResultHandle);
565  bool HBHENoiseFilterResult = *HBHENoiseFilterResultHandle;
566  if (!HBHENoiseFilterResultHandle.isValid()) {
567  LogDebug("") << "PFMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
568  if (_verbose) std::cout << "PFMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
569  }
570 
571 
572  edm::Handle<reco::CaloJetCollection> caloJets;
573  iEvent.getByLabel(theJetCollectionLabel, caloJets);
574  if (!caloJets.isValid()) {
575  LogDebug("") << "PFMETAnalyzer: Could not find jet product" << std::endl;
576  if (_verbose) std::cout << "PFMETAnalyzer: Could not find jet product" << std::endl;
577  }
578 
579  edm::Handle<edm::View<PFCandidate> > pfCandidates;
580  iEvent.getByLabel(PFCandidatesTag, pfCandidates);
581  if (!pfCandidates.isValid()) {
582  LogDebug("") << "PfMETAnalyzer: Could not find pfcandidates product" << std::endl;
583  if (_verbose) std::cout << "PfMETAnalyzer: Could not find pfcandidates product" << std::endl;
584  }
585 
586  edm::Handle<reco::PFJetCollection> pfJets;
587  iEvent.getByLabel(thePfJetCollectionLabel, pfJets);
588  if (!pfJets.isValid()) {
589  LogDebug("") << "PFMETAnalyzer: Could not find pfjet product" << std::endl;
590  if (_verbose) std::cout << "PFMETAnalyzer: Could not find pfjet product" << std::endl;
591  }
592  // ==========================================================
593  // PfMET sanity check
594 
595  if (_source=="PfMET") validateMET(*pfmet, pfCandidates);
596 
597  // ==========================================================
598  // JetID
599 
600  if (_verbose) std::cout << "JetID starts" << std::endl;
601 
602  //
603  // --- Minimal cuts
604  //
605  bool bJetIDMinimal=true;
606  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
607  cal!=caloJets->end(); ++cal){
608  jetID->calculate(iEvent, *cal);
609  if (cal->pt()>10.){
610  if (fabs(cal->eta())<=2.6 &&
611  cal->emEnergyFraction()<=0.01) bJetIDMinimal=false;
612  }
613  }
614 
615  //
616  // --- Loose cuts, not PF specific for now!
617  //
618  bool bJetIDLoose=true;
619  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
620  cal!=caloJets->end(); ++cal){
621  jetID->calculate(iEvent, *cal);
622  if (_verbose) std::cout << jetID->n90Hits() << " "
623  << jetID->restrictedEMF() << " "
624  << cal->pt() << std::endl;
625  if (cal->pt()>10.){
626  //
627  // for all regions
628  if (jetID->n90Hits()<2) bJetIDLoose=false;
629  if (jetID->fHPD()>=0.98) bJetIDLoose=false;
630  //if (jetID->restrictedEMF()<0.01) bJetIDLoose=false;
631  //
632  // for non-forward
633  if (fabs(cal->eta())<2.55){
634  if (cal->emEnergyFraction()<=0.01) bJetIDLoose=false;
635  }
636  // for forward
637  else {
638  if (cal->emEnergyFraction()<=-0.9) bJetIDLoose=false;
639  if (cal->pt()>80.){
640  if (cal->emEnergyFraction()>= 1.0) bJetIDLoose=false;
641  }
642  } // forward vs non-forward
643  } // pt>10 GeV/c
644  } // calor-jets loop
645 
646  //
647  // --- Tight cuts
648  //
649  bool bJetIDTight=true;
650  bJetIDTight=bJetIDLoose;
651  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
652  cal!=caloJets->end(); ++cal){
653  jetID->calculate(iEvent, *cal);
654  if (cal->pt()>25.){
655  //
656  // for all regions
657  if (jetID->fHPD()>=0.95) bJetIDTight=false;
658  //
659  // for 1.0<|eta|<1.75
660  if (fabs(cal->eta())>=1.00 && fabs(cal->eta())<1.75){
661  if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
662  }
663  //
664  // for 1.75<|eta|<2.55
665  else if (fabs(cal->eta())>=1.75 && fabs(cal->eta())<2.55){
666  if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
667  }
668  //
669  // for 2.55<|eta|<3.25
670  else if (fabs(cal->eta())>=2.55 && fabs(cal->eta())<3.25){
671  if (cal->pt()< 50. && cal->emEnergyFraction()<=-0.3) bJetIDTight=false;
672  if (cal->pt()>=50. && cal->pt()< 80. && cal->emEnergyFraction()<=-0.2) bJetIDTight=false;
673  if (cal->pt()>=80. && cal->pt()<340. && cal->emEnergyFraction()<=-0.1) bJetIDTight=false;
674  if (cal->pt()>=340. && cal->emEnergyFraction()<=-0.1
675  && cal->emEnergyFraction()>=0.95) bJetIDTight=false;
676  }
677  //
678  // for 3.25<|eta|
679  else if (fabs(cal->eta())>=3.25){
680  if (cal->pt()< 50. && cal->emEnergyFraction()<=-0.3
681  && cal->emEnergyFraction()>=0.90) bJetIDTight=false;
682  if (cal->pt()>=50. && cal->pt()<130. && cal->emEnergyFraction()<=-0.2
683  && cal->emEnergyFraction()>=0.80) bJetIDTight=false;
684  if (cal->pt()>=130. && cal->emEnergyFraction()<=-0.1
685  && cal->emEnergyFraction()>=0.70) bJetIDTight=false;
686  }
687  } // pt>10 GeV/c
688  } // calor-jets loop
689 
690  if (_verbose) std::cout << "JetID ends" << std::endl;
691 
692 
693  // ==========================================================
694  // HCAL Noise filter
695 
696  bool bHcalNoiseFilter = HBHENoiseFilterResult;
697 
698  // ==========================================================
699  // Get BeamHaloSummary
700  edm::Handle<BeamHaloSummary> TheBeamHaloSummary ;
701  iEvent.getByLabel(BeamHaloSummaryTag, TheBeamHaloSummary) ;
702 
703  bool bBeamHaloIDTightPass = true;
704  bool bBeamHaloIDLoosePass = true;
705 
706  if(!TheBeamHaloSummary.isValid()) {
707 
708  const BeamHaloSummary TheSummary = (*TheBeamHaloSummary.product() );
709 
710  if( !TheSummary.EcalLooseHaloId() && !TheSummary.HcalLooseHaloId() &&
711  !TheSummary.CSCLooseHaloId() && !TheSummary.GlobalLooseHaloId() )
712  bBeamHaloIDLoosePass = false;
713 
714  if( !TheSummary.EcalTightHaloId() && !TheSummary.HcalTightHaloId() &&
715  !TheSummary.CSCTightHaloId() && !TheSummary.GlobalTightHaloId() )
716  bBeamHaloIDTightPass = false;
717 
718  }
719 
720  // ==========================================================
721  //Vertex information
722 
723  _numPV = 0;
724  bool bPrimaryVertex = true;
725  if(_doPVCheck){
726  bPrimaryVertex = false;
727  Handle<VertexCollection> vertexHandle;
728 
729  iEvent.getByLabel(vertexTag, vertexHandle);
730 
731  if (!vertexHandle.isValid()) {
732  LogDebug("") << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
733  if (_verbose) std::cout << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
734  }
735 
736  if ( vertexHandle.isValid() ){
737  VertexCollection vertexCollection = *(vertexHandle.product());
738  int vertex_number = vertexCollection.size();
739  VertexCollection::const_iterator v = vertexCollection.begin();
740  for ( ; v != vertexCollection.end(); ++v) {
741  double vertex_chi2 = v->normalizedChi2();
742  double vertex_ndof = v->ndof();
743  bool fakeVtx = v->isFake();
744  double vertex_Z = v->z();
745 
746  if ( !fakeVtx
747  && vertex_number>=_nvtx_min
748  && vertex_ndof >_vtxndof_min
749  && vertex_chi2 <_vtxchi2_max
750  && fabs(vertex_Z)<_vtxz_max ) {
751  bPrimaryVertex = true;
752  ++_numPV;
753  }
754  }
755  }
756  }
757  // ==========================================================
758 
759  edm::Handle< L1GlobalTriggerReadoutRecord > gtReadoutRecord;
760  iEvent.getByLabel( gtTag, gtReadoutRecord);
761 
762  if (!gtReadoutRecord.isValid()) {
763  LogDebug("") << "CaloMETAnalyzer: Could not find GT readout record" << std::endl;
764  if (_verbose) std::cout << "CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
765  }
766 
767  bool bTechTriggers = true;
768  bool bTechTriggersAND = true;
769  bool bTechTriggersOR = false;
770  bool bTechTriggersNOT = false;
771 
772  if (gtReadoutRecord.isValid()) {
773  const TechnicalTriggerWord& technicalTriggerWordBeforeMask = gtReadoutRecord->technicalTriggerWord();
774 
775  if (_techTrigsAND.size() == 0)
776  bTechTriggersAND = true;
777  else
778  for (unsigned ttr = 0; ttr != _techTrigsAND.size(); ttr++) {
779  bTechTriggersAND = bTechTriggersAND && technicalTriggerWordBeforeMask.at(_techTrigsAND.at(ttr));
780  }
781 
782  if (_techTrigsAND.size() == 0)
783  bTechTriggersOR = true;
784  else
785  for (unsigned ttr = 0; ttr != _techTrigsOR.size(); ttr++) {
786  bTechTriggersOR = bTechTriggersOR || technicalTriggerWordBeforeMask.at(_techTrigsOR.at(ttr));
787  }
788  if (_techTrigsNOT.size() == 0)
789  bTechTriggersNOT = false;
790  else
791  for (unsigned ttr = 0; ttr != _techTrigsNOT.size(); ttr++) {
792  bTechTriggersNOT = bTechTriggersNOT || technicalTriggerWordBeforeMask.at(_techTrigsNOT.at(ttr));
793  }
794  }
795  else
796  {
797  bTechTriggersAND = true;
798  bTechTriggersOR = true;
799  bTechTriggersNOT = false;
800  }
801 
802  if (_techTrigsAND.size()==0)
803  bTechTriggersAND = true;
804  if (_techTrigsOR.size()==0)
805  bTechTriggersOR = true;
806  if (_techTrigsNOT.size()==0)
807  bTechTriggersNOT = false;
808 
809  bTechTriggers = bTechTriggersAND && bTechTriggersOR && !bTechTriggersNOT;
810 
811  // ==========================================================
812  // Reconstructed MET Information - fill MonitorElements
813 
814  bool bHcalNoise = bHcalNoiseFilter;
815  bool bBeamHaloID = bBeamHaloIDLoosePass;
816  bool bJetID = bJetIDMinimal;
817 
818  bool bPhysicsDeclared = true;
819  if(_doHLTPhysicsOn) bPhysicsDeclared =_trig_PhysDec;
820 
821  if (_tightBHFiltering) bBeamHaloID = bBeamHaloIDTightPass;
822 
823  if (_tightJetIDFiltering==1) bJetID = bJetIDMinimal;
824  else if (_tightJetIDFiltering==2) bJetID = bJetIDLoose;
825  else if (_tightJetIDFiltering==3) bJetID = bJetIDTight;
826  else if (_tightJetIDFiltering==-1) bJetID = true;
827 
828  bool bBasicCleanup = bTechTriggers && bPrimaryVertex && bPhysicsDeclared;
829  bool bExtraCleanup = bBasicCleanup && bHcalNoise && bJetID && bBeamHaloID;
830 
831  std::string DirName = "JetMET/MET/"+_source;
832 
833  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
834  ic != _FolderNames.end(); ic++){
835  if (*ic=="All") fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
836  if (DCSFilter->filter(iEvent, iSetup)) {
837  if (_cleanupSelection){
838  if (*ic=="BasicCleanup" && bBasicCleanup) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
839  if (*ic=="ExtraCleanup" && bExtraCleanup) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
840  }
841  if (_allSelection) {
842  if (*ic=="HcalNoiseFilter" && bHcalNoiseFilter ) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
843  if (*ic=="JetIDMinimal" && bJetIDMinimal) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
844  if (*ic=="JetIDLoose" && bJetIDLoose) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
845  if (*ic=="JetIDTight" && bJetIDTight) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
846  if (*ic=="BeamHaloIDTightPass" && bBeamHaloIDTightPass) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
847  if (*ic=="BeamHaloIDLoosePass" && bBeamHaloIDLoosePass) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
848  if (*ic=="Triggers" && bTechTriggers) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
849  if (*ic=="PV" && bPrimaryVertex) fillMESet(iEvent, DirName+"/"+*ic, *pfmet);
850  }
851  } // DCS
852  }
853 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
reco::BeamHaloSummary::EcalTightHaloId
const bool EcalTightHaloId() const
Definition: BeamHaloSummary.h:38
i
int i
Definition: DBlmapReader.cc:9
PFMETAnalyzer::metname
std::string metname
Definition: PFMETAnalyzer.h:104
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:207
PFMETAnalyzer::_tightJetIDFiltering
int _tightJetIDFiltering
Definition: PFMETAnalyzer.h:147
PFMETAnalyzer::highPtJetExpr_
std::vector< std::string > highPtJetExpr_
Definition: PFMETAnalyzer.h:128
PFMETAnalyzer::DCSFilter
JetMETDQMDCSFilter * DCSFilter
Definition: PFMETAnalyzer.h:182
reco::BeamHaloSummary::HcalTightHaloId
const bool HcalTightHaloId() const
Definition: BeamHaloSummary.h:35
PFMETAnalyzer::_source
std::string _source
Definition: PFMETAnalyzer.h:105
reco::BeamHaloSummary::GlobalLooseHaloId
const bool GlobalLooseHaloId() const
Definition: BeamHaloSummary.h:43
PFMETAnalyzer::validateMET
void validateMET(const reco::PFMET &, edm::Handle< edm::View< reco::PFCandidate > >)
Definition: PFMETAnalyzer.cc:857
PFMETAnalyzer::gtTag
edm::InputTag gtTag
Definition: PFMETAnalyzer.h:115
edm::HLTGlobalStatus::accept
bool accept() const
Has at least one path accepted the event?
Definition: HLTGlobalStatus.h:54
reco::helper::JetIDHelper::fHPD
double fHPD() const
Definition: JetIDHelper.h:33
PFMETAnalyzer::fillMESet
void fillMESet(const edm::Event &, std::string, const reco::PFMET &)
Definition: PFMETAnalyzer.cc:905
PFMETAnalyzer::_techTrigsAND
std::vector< unsigned > _techTrigsAND
Definition: PFMETAnalyzer.h:139
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
PFMETAnalyzer::muonExpr_
std::vector< std::string > muonExpr_
Definition: PFMETAnalyzer.h:132
edm::TriggerNames::size
Strings::size_type size() const
Definition: TriggerNames.cc:39
reco::BeamHaloSummary::GlobalTightHaloId
const bool GlobalTightHaloId() const
Definition: BeamHaloSummary.h:44
reco::helper::JetIDHelper::restrictedEMF
double restrictedEMF() const
Definition: JetIDHelper.h:50
GoodVertex_cfg.vertexCollection
tuple vertexCollection
Definition: GoodVertex_cfg.py:61
PFMETAnalyzer::lowPtJetExpr_
std::vector< std::string > lowPtJetExpr_
Definition: PFMETAnalyzer.h:129
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PFMETAnalyzer::highMETExpr_
std::vector< std::string > highMETExpr_
Definition: PFMETAnalyzer.h:130
PFMETAnalyzer::_vtxchi2_max
double _vtxchi2_max
Definition: PFMETAnalyzer.h:152
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
reco::BeamHaloSummary::HcalLooseHaloId
const bool HcalLooseHaloId() const
Definition: BeamHaloSummary.h:34
reco::tau::pfCandidates
std::vector< PFCandidatePtr > pfCandidates(const PFJet &jet, int particleId, bool sort=true)
Definition: RecoTauCommonUtilities.cc:29
reco::BeamHaloSummary::CSCTightHaloId
const bool CSCTightHaloId() const
Definition: BeamHaloSummary.h:41
PFMETAnalyzer::BeamHaloSummaryTag
edm::InputTag BeamHaloSummaryTag
Definition: PFMETAnalyzer.h:112
PFMETAnalyzer::_hlt_PhysDec
std::string _hlt_PhysDec
Definition: PFMETAnalyzer.h:137
edm::TriggerNames::triggerIndex
unsigned int triggerIndex(std::string const &name) const
Definition: TriggerNames.cc:32
PFMETAnalyzer::jetID
reco::helper::JetIDHelper * jetID
Definition: PFMETAnalyzer.h:179
PFMETAnalyzer::_techTrigsOR
std::vector< unsigned > _techTrigsOR
Definition: PFMETAnalyzer.h:140
edm::HLTGlobalStatus::size
unsigned int size() const
Get number of paths stored.
Definition: HLTGlobalStatus.h:41
PFMETAnalyzer::thePfMETCollectionLabel
edm::InputTag thePfMETCollectionLabel
Definition: PFMETAnalyzer.h:107
PFMETAnalyzer::thePfJetCollectionLabel
edm::InputTag thePfJetCollectionLabel
Definition: PFMETAnalyzer.h:110
TechnicalTriggerWord
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
Definition: L1GlobalTriggerReadoutSetupFwd.h:43
PFMETAnalyzer::PFCandidatesTag
edm::InputTag PFCandidatesTag
Definition: PFMETAnalyzer.h:111
PFMETAnalyzer::elecExpr_
std::vector< std::string > elecExpr_
Definition: PFMETAnalyzer.h:133
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
JetMETDQMDCSFilter::filter
bool filter(const edm::Event &evt, const edm::EventSetup &es)
Definition: JetMETDQMDCSFilter.cc:34
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
PFMETAnalyzer::_tightBHFiltering
bool _tightBHFiltering
Definition: PFMETAnalyzer.h:146
PFMETAnalyzer::_FolderNames
std::vector< std::string > _FolderNames
Definition: PFMETAnalyzer.h:189
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
PFMETAnalyzer::theJetCollectionLabel
edm::InputTag theJetCollectionLabel
Definition: PFMETAnalyzer.h:109
edm::Handle::product
T const * product() const
Definition: Handle.h:74
reco::BeamHaloSummary::CSCLooseHaloId
const bool CSCLooseHaloId() const
Definition: BeamHaloSummary.h:40
PFMETAnalyzer::metME
MonitorElement * metME
Definition: PFMETAnalyzer.h:196
PFMETAnalyzer::HBHENoiseFilterResultTag
edm::InputTag HBHENoiseFilterResultTag
Definition: PFMETAnalyzer.h:113
PFMETAnalyzer::HcalNoiseRBXCollectionTag
edm::InputTag HcalNoiseRBXCollectionTag
Definition: PFMETAnalyzer.h:108
PFMETAnalyzer::_cleanupSelection
bool _cleanupSelection
Definition: PFMETAnalyzer.h:187
FamosSequences_cff.caloJets
tuple caloJets
Definition: FamosSequences_cff.py:68
PFMETAnalyzer::_techTrigsNOT
std::vector< unsigned > _techTrigsNOT
Definition: PFMETAnalyzer.h:141
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
PFMETAnalyzer::minbiasExpr_
std::vector< std::string > minbiasExpr_
Definition: PFMETAnalyzer.h:134
PFMETAnalyzer::HLTPathsJetMBByName_
std::vector< std::string > HLTPathsJetMBByName_
Definition: PFMETAnalyzer.h:118
PFMETAnalyzer::vertexTag
edm::InputTag vertexTag
Definition: PFMETAnalyzer.h:114
pfJets_cff.pfJets
tuple pfJets
Definition: pfJets_cff.py:8
reco::BeamHaloSummary::EcalLooseHaloId
const bool EcalLooseHaloId() const
Definition: BeamHaloSummary.h:37
reco::helper::JetIDHelper::calculate
void calculate(const edm::Event &event, const reco::CaloJet &jet, const int iDbg=0)
Definition: JetIDHelper.cc:86
PFMETCollection
Collection of PF MET.
void PFMETAnalyzer::beginJob ( DQMStore *  dbe)
virtual

Inizialize parameters for histo binning.

Implements PFMETAnalyzerBase.

Definition at line 81 of file PFMETAnalyzer.cc.

References DQMStore::book1D(), cmsPerfPublish::DirName, edm::ParameterSet::getParameter(), bTagSequences_cff::jetID, LogTrace, metname, Parameters::parameters, MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

81  {
82 
83  evtCounter = 0;
84  metname = "pfMETAnalyzer";
85 
86  // trigger information
87  HLTPathsJetMBByName_ = parameters.getParameter<std::vector<std::string > >("HLTPathsJetMB");
88 
89  theCleaningParameters = parameters.getParameter<ParameterSet>("CleaningParameters"),
90 
91  //Trigger parameters
92  gtTag = theCleaningParameters.getParameter<edm::InputTag>("gtLabel");
93  _techTrigsAND = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsAND");
94  _techTrigsOR = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsOR");
95  _techTrigsNOT = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsNOT");
96 
97  _doHLTPhysicsOn = theCleaningParameters.getParameter<bool>("doHLTPhysicsOn");
98  _hlt_PhysDec = theCleaningParameters.getParameter<std::string>("HLT_PhysDec");
99 
100  _tightBHFiltering = theCleaningParameters.getParameter<bool>("tightBHFiltering");
101  _tightJetIDFiltering = theCleaningParameters.getParameter<int>("tightJetIDFiltering");
102 
103 
104  // ==========================================================
105  //DCS information
106  // ==========================================================
107  DCSFilter = new JetMETDQMDCSFilter(parameters.getParameter<ParameterSet>("DCSFilter"));
108 
109  //Vertex requirements
110  _doPVCheck = theCleaningParameters.getParameter<bool>("doPrimaryVertexCheck");
111  vertexTag = theCleaningParameters.getParameter<edm::InputTag>("vertexLabel");
112 
113  if (_doPVCheck) {
114  _nvtx_min = theCleaningParameters.getParameter<int>("nvtx_min");
115  _nvtxtrks_min = theCleaningParameters.getParameter<int>("nvtxtrks_min");
116  _vtxndof_min = theCleaningParameters.getParameter<int>("vtxndof_min");
117  _vtxchi2_max = theCleaningParameters.getParameter<double>("vtxchi2_max");
118  _vtxz_max = theCleaningParameters.getParameter<double>("vtxz_max");
119  }
120 
121 
122  // PFMET information
123  thePfMETCollectionLabel = parameters.getParameter<edm::InputTag>("METCollectionLabel");
124  thePfJetCollectionLabel = parameters.getParameter<edm::InputTag>("PfJetCollectionLabel");
125  _source = parameters.getParameter<std::string>("Source");
126 
127  // Other data collections
128  theJetCollectionLabel = parameters.getParameter<edm::InputTag>("JetCollectionLabel");
129  PFCandidatesTag = parameters.getParameter<edm::InputTag>("PFCandidates");
130  HcalNoiseRBXCollectionTag = parameters.getParameter<edm::InputTag>("HcalNoiseRBXCollection");
131  BeamHaloSummaryTag = parameters.getParameter<edm::InputTag>("BeamHaloSummaryLabel");
132  HBHENoiseFilterResultTag = parameters.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");
133 
134  // misc
135  _verbose = parameters.getParameter<int>("verbose");
136  _etThreshold = parameters.getParameter<double>("etThreshold"); // MET threshold
137  _allhist = parameters.getParameter<bool>("allHist"); // Full set of monitoring histograms
138  _allSelection= parameters.getParameter<bool>("allSelection"); // Plot with all sets of event selection
139  _cleanupSelection= parameters.getParameter<bool>("cleanupSelection"); // Plot with all sets of event selection
140 
141  _highPtPFJetThreshold = parameters.getParameter<double>("HighPtJetThreshold");
142  _lowPtPFJetThreshold = parameters.getParameter<double>("LowPtJetThreshold");
143  _highPFMETThreshold = parameters.getParameter<double>("HighMETThreshold");
144 
145  //
146  jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams"));
147 
148  // DQStore stuff
149  LogTrace(metname)<<"[PFMETAnalyzer] Parameters initialization";
150  std::string DirName = "JetMET/MET/"+_source;
151  dbe->setCurrentFolder(DirName);
152 
153  metME = dbe->book1D("metReco", "metReco", 4, 1, 5);
154  metME->setBinLabel(3,"PFMET",1);
155 
156  _dbe = dbe;
157 
158  _FolderNames.push_back("All");
159  _FolderNames.push_back("BasicCleanup");
160  _FolderNames.push_back("ExtraCleanup");
161  _FolderNames.push_back("HcalNoiseFilter");
162  _FolderNames.push_back("JetIDMinimal");
163  _FolderNames.push_back("JetIDLoose");
164  _FolderNames.push_back("JetIDTight");
165  _FolderNames.push_back("BeamHaloIDTightPass");
166  _FolderNames.push_back("BeamHaloIDLoosePass");
167  _FolderNames.push_back("Triggers");
168  _FolderNames.push_back("PV");
169 
170  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
171  ic != _FolderNames.end(); ic++){
172  if (*ic=="All") bookMESet(DirName+"/"+*ic);
173  if (_cleanupSelection){
174  if (*ic=="BasicCleanup") bookMESet(DirName+"/"+*ic);
175  if (*ic=="ExtraCleanup") bookMESet(DirName+"/"+*ic);
176  }
177  if (_allSelection){
178  if (*ic=="HcalNoiseFilter") bookMESet(DirName+"/"+*ic);
179  if (*ic=="JetIDMinimal") bookMESet(DirName+"/"+*ic);
180  if (*ic=="JetIDLoose") bookMESet(DirName+"/"+*ic);
181  if (*ic=="JetIDTight") bookMESet(DirName+"/"+*ic);
182  if (*ic=="BeamHaloIDTightPass") bookMESet(DirName+"/"+*ic);
183  if (*ic=="BeamHaloIDLoosePass") bookMESet(DirName+"/"+*ic);
184  if (*ic=="Triggers") bookMESet(DirName+"/"+*ic);
185  if (*ic=="PV") bookMESet(DirName+"/"+*ic);
186  }
187  }
188 }
PFMETAnalyzer::_highPFMETThreshold
double _highPFMETThreshold
Definition: PFMETAnalyzer.h:168
PFMETAnalyzer::_highPtPFJetThreshold
double _highPtPFJetThreshold
Definition: PFMETAnalyzer.h:166
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PFMETAnalyzer::metname
std::string metname
Definition: PFMETAnalyzer.h:104
PFMETAnalyzer::_tightJetIDFiltering
int _tightJetIDFiltering
Definition: PFMETAnalyzer.h:147
PFMETAnalyzer::DCSFilter
JetMETDQMDCSFilter * DCSFilter
Definition: PFMETAnalyzer.h:182
PFMETAnalyzer::_source
std::string _source
Definition: PFMETAnalyzer.h:105
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
PFMETAnalyzer::parameters
edm::ParameterSet parameters
Definition: PFMETAnalyzer.h:101
PFMETAnalyzer::gtTag
edm::InputTag gtTag
Definition: PFMETAnalyzer.h:115
MonitorElement::setBinLabel
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)
Definition: MonitorElement.cc:889
PFMETAnalyzer::_techTrigsAND
std::vector< unsigned > _techTrigsAND
Definition: PFMETAnalyzer.h:139
PFMETAnalyzer::theCleaningParameters
edm::ParameterSet theCleaningParameters
Definition: PFMETAnalyzer.h:136
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PFMETAnalyzer::_vtxchi2_max
double _vtxchi2_max
Definition: PFMETAnalyzer.h:152
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
PFMETAnalyzer::BeamHaloSummaryTag
edm::InputTag BeamHaloSummaryTag
Definition: PFMETAnalyzer.h:112
PFMETAnalyzer::_hlt_PhysDec
std::string _hlt_PhysDec
Definition: PFMETAnalyzer.h:137
PFMETAnalyzer::jetID
reco::helper::JetIDHelper * jetID
Definition: PFMETAnalyzer.h:179
PFMETAnalyzer::_lowPtPFJetThreshold
double _lowPtPFJetThreshold
Definition: PFMETAnalyzer.h:167
PFMETAnalyzer::_techTrigsOR
std::vector< unsigned > _techTrigsOR
Definition: PFMETAnalyzer.h:140
PFMETAnalyzer::thePfMETCollectionLabel
edm::InputTag thePfMETCollectionLabel
Definition: PFMETAnalyzer.h:107
PFMETAnalyzer::thePfJetCollectionLabel
edm::InputTag thePfJetCollectionLabel
Definition: PFMETAnalyzer.h:110
PFMETAnalyzer::bookMESet
void bookMESet(std::string)
Definition: PFMETAnalyzer.cc:199
PFMETAnalyzer::PFCandidatesTag
edm::InputTag PFCandidatesTag
Definition: PFMETAnalyzer.h:111
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
PFMETAnalyzer::_tightBHFiltering
bool _tightBHFiltering
Definition: PFMETAnalyzer.h:146
PFMETAnalyzer::_etThreshold
double _etThreshold
Definition: PFMETAnalyzer.h:173
PFMETAnalyzer::_FolderNames
std::vector< std::string > _FolderNames
Definition: PFMETAnalyzer.h:189
PFMETAnalyzer::theJetCollectionLabel
edm::InputTag theJetCollectionLabel
Definition: PFMETAnalyzer.h:109
PFMETAnalyzer::metME
MonitorElement * metME
Definition: PFMETAnalyzer.h:196
PFMETAnalyzer::HBHENoiseFilterResultTag
edm::InputTag HBHENoiseFilterResultTag
Definition: PFMETAnalyzer.h:113
PFMETAnalyzer::HcalNoiseRBXCollectionTag
edm::InputTag HcalNoiseRBXCollectionTag
Definition: PFMETAnalyzer.h:108
PFMETAnalyzer::_cleanupSelection
bool _cleanupSelection
Definition: PFMETAnalyzer.h:187
PFMETAnalyzer::_techTrigsNOT
std::vector< unsigned > _techTrigsNOT
Definition: PFMETAnalyzer.h:141
PFMETAnalyzer::HLTPathsJetMBByName_
std::vector< std::string > HLTPathsJetMBByName_
Definition: PFMETAnalyzer.h:118
PFMETAnalyzer::vertexTag
edm::InputTag vertexTag
Definition: PFMETAnalyzer.h:114
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
void PFMETAnalyzer::beginRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup 
)

Initialize run-based parameters.

Definition at line 342 of file PFMETAnalyzer.cc.

344 {
345  if (_HighPtJetEventFlag->on()) _HighPtJetEventFlag->initRun(iRun, iSetup);
346  if (_LowPtJetEventFlag ->on()) _LowPtJetEventFlag ->initRun(iRun, iSetup);
347  if (_MinBiasEventFlag ->on()) _MinBiasEventFlag ->initRun(iRun, iSetup);
348  if (_HighMETEventFlag ->on()) _HighMETEventFlag ->initRun(iRun, iSetup);
349  if (_EleEventFlag ->on()) _EleEventFlag ->initRun(iRun, iSetup);
350  if (_MuonEventFlag ->on()) _MuonEventFlag ->initRun(iRun, iSetup);
351 
352  if (_HighPtJetEventFlag->on() && _HighPtJetEventFlag->expressionsFromDB(_HighPtJetEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
353  highPtJetExpr_ = _HighPtJetEventFlag->expressionsFromDB(_HighPtJetEventFlag->hltDBKey(), iSetup);
354 
355  if (_LowPtJetEventFlag->on() && _LowPtJetEventFlag->expressionsFromDB(_LowPtJetEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
356  lowPtJetExpr_ = _LowPtJetEventFlag->expressionsFromDB(_LowPtJetEventFlag->hltDBKey(), iSetup);
357 
358  if (_HighMETEventFlag->on() && _HighMETEventFlag->expressionsFromDB(_HighMETEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
359  highMETExpr_ = _HighMETEventFlag->expressionsFromDB(_HighMETEventFlag->hltDBKey(), iSetup);
360 
361  if (_MuonEventFlag->on() && _MuonEventFlag->expressionsFromDB(_MuonEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
362  muonExpr_ = _MuonEventFlag->expressionsFromDB(_MuonEventFlag->hltDBKey(), iSetup);
363 
364  if (_EleEventFlag->on() && _EleEventFlag->expressionsFromDB(_EleEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
365  elecExpr_ = _EleEventFlag->expressionsFromDB(_EleEventFlag->hltDBKey(), iSetup);
366 
367  if (_MinBiasEventFlag->on() && _MinBiasEventFlag->expressionsFromDB(_MinBiasEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
368  minbiasExpr_ = _MinBiasEventFlag->expressionsFromDB(_MinBiasEventFlag->hltDBKey(), iSetup);
369 }
PFMETAnalyzer::_MuonEventFlag
GenericTriggerEventFlag * _MuonEventFlag
Definition: PFMETAnalyzer.h:126
PFMETAnalyzer::highPtJetExpr_
std::vector< std::string > highPtJetExpr_
Definition: PFMETAnalyzer.h:128
PFMETAnalyzer::_HighPtJetEventFlag
GenericTriggerEventFlag * _HighPtJetEventFlag
Definition: PFMETAnalyzer.h:120
PFMETAnalyzer::muonExpr_
std::vector< std::string > muonExpr_
Definition: PFMETAnalyzer.h:132
PFMETAnalyzer::_LowPtJetEventFlag
GenericTriggerEventFlag * _LowPtJetEventFlag
Definition: PFMETAnalyzer.h:121
PFMETAnalyzer::lowPtJetExpr_
std::vector< std::string > lowPtJetExpr_
Definition: PFMETAnalyzer.h:129
PFMETAnalyzer::highMETExpr_
std::vector< std::string > highMETExpr_
Definition: PFMETAnalyzer.h:130
GenericTriggerEventFlag::expressionsFromDB
std::vector< std::string > expressionsFromDB(const std::string &key, const edm::EventSetup &setup)
Reads and returns logical expressions from DB.
Definition: GenericTriggerEventFlag.cc:602
PFMETAnalyzer::_EleEventFlag
GenericTriggerEventFlag * _EleEventFlag
Definition: PFMETAnalyzer.h:125
PFMETAnalyzer::_HighMETEventFlag
GenericTriggerEventFlag * _HighMETEventFlag
Definition: PFMETAnalyzer.h:123
PFMETAnalyzer::elecExpr_
std::vector< std::string > elecExpr_
Definition: PFMETAnalyzer.h:133
PFMETAnalyzer::_MinBiasEventFlag
GenericTriggerEventFlag * _MinBiasEventFlag
Definition: PFMETAnalyzer.h:122
PFMETAnalyzer::minbiasExpr_
std::vector< std::string > minbiasExpr_
Definition: PFMETAnalyzer.h:134
GenericTriggerEventFlag::initRun
void initRun(const edm::Run &run, const edm::EventSetup &setup)
To be called from beginRun() methods.
Definition: GenericTriggerEventFlag.cc:125
void PFMETAnalyzer::bookMESet ( std::string  DirName)

Definition at line 199 of file PFMETAnalyzer.cc.

References gather_cfg::cout.

200 {
201 
202  bool bLumiSecPlot=false;
203  if (DirName.find("All")!=std::string::npos) bLumiSecPlot=true;
204 
205  bookMonitorElement(DirName,bLumiSecPlot);
206 
207  if ( _HighPtJetEventFlag->on() ) {
208  bookMonitorElement(DirName+"/"+"HighPtJet",false);
209  meTriggerName_HighPtJet = _dbe->bookString("triggerName_HighPtJet", highPtJetExpr_[0]);
210  }
211 
212  if ( _LowPtJetEventFlag->on() ) {
213  bookMonitorElement(DirName+"/"+"LowPtJet",false);
214  meTriggerName_LowPtJet = _dbe->bookString("triggerName_LowPtJet", lowPtJetExpr_[0]);
215  }
216 
217  if ( _MinBiasEventFlag->on() ) {
218  bookMonitorElement(DirName+"/"+"MinBias",false);
219  meTriggerName_MinBias = _dbe->bookString("triggerName_MinBias", minbiasExpr_[0]);
220  if (_verbose) std::cout << "_MinBiasEventFlag is on, folder created\n";
221  }
222 
223  if ( _HighMETEventFlag->on() ) {
224  bookMonitorElement(DirName+"/"+"HighMET",false);
225  meTriggerName_HighMET = _dbe->bookString("triggerName_HighMET", highMETExpr_[0]);
226  }
227 
228  if ( _EleEventFlag->on() ) {
229  bookMonitorElement(DirName+"/"+"Ele",false);
230  meTriggerName_Ele = _dbe->bookString("triggerName_Ele", elecExpr_[0]);
231  if (_verbose) std::cout << "_EleEventFlag is on, folder created\n";
232  }
233 
234  if ( _MuonEventFlag->on() ) {
235  bookMonitorElement(DirName+"/"+"Muon",false);
236  meTriggerName_Muon = _dbe->bookString("triggerName_Muon", muonExpr_[0]);
237  if (_verbose) std::cout << "_MuonEventFlag is on, folder created\n";
238  }
239 }
PFMETAnalyzer::meTriggerName_HighPtJet
MonitorElement * meTriggerName_HighPtJet
Definition: PFMETAnalyzer.h:199
PFMETAnalyzer::_MuonEventFlag
GenericTriggerEventFlag * _MuonEventFlag
Definition: PFMETAnalyzer.h:126
PFMETAnalyzer::meTriggerName_Muon
MonitorElement * meTriggerName_Muon
Definition: PFMETAnalyzer.h:204
PFMETAnalyzer::highPtJetExpr_
std::vector< std::string > highPtJetExpr_
Definition: PFMETAnalyzer.h:128
PFMETAnalyzer::meTriggerName_MinBias
MonitorElement * meTriggerName_MinBias
Definition: PFMETAnalyzer.h:201
PFMETAnalyzer::_HighPtJetEventFlag
GenericTriggerEventFlag * _HighPtJetEventFlag
Definition: PFMETAnalyzer.h:120
PFMETAnalyzer::meTriggerName_LowPtJet
MonitorElement * meTriggerName_LowPtJet
Definition: PFMETAnalyzer.h:200
PFMETAnalyzer::muonExpr_
std::vector< std::string > muonExpr_
Definition: PFMETAnalyzer.h:132
PFMETAnalyzer::_LowPtJetEventFlag
GenericTriggerEventFlag * _LowPtJetEventFlag
Definition: PFMETAnalyzer.h:121
PFMETAnalyzer::lowPtJetExpr_
std::vector< std::string > lowPtJetExpr_
Definition: PFMETAnalyzer.h:129
PFMETAnalyzer::highMETExpr_
std::vector< std::string > highMETExpr_
Definition: PFMETAnalyzer.h:130
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
PFMETAnalyzer::_EleEventFlag
GenericTriggerEventFlag * _EleEventFlag
Definition: PFMETAnalyzer.h:125
PFMETAnalyzer::meTriggerName_Ele
MonitorElement * meTriggerName_Ele
Definition: PFMETAnalyzer.h:203
PFMETAnalyzer::_HighMETEventFlag
GenericTriggerEventFlag * _HighMETEventFlag
Definition: PFMETAnalyzer.h:123
DQMStore::bookString
MonitorElement * bookString(const char *name, const char *value)
Book string.
Definition: DQMStore.cc:688
PFMETAnalyzer::elecExpr_
std::vector< std::string > elecExpr_
Definition: PFMETAnalyzer.h:133
PFMETAnalyzer::bookMonitorElement
void bookMonitorElement(std::string, bool)
Definition: PFMETAnalyzer.cc:245
PFMETAnalyzer::meTriggerName_HighMET
MonitorElement * meTriggerName_HighMET
Definition: PFMETAnalyzer.h:202
PFMETAnalyzer::_MinBiasEventFlag
GenericTriggerEventFlag * _MinBiasEventFlag
Definition: PFMETAnalyzer.h:122
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
PFMETAnalyzer::minbiasExpr_
std::vector< std::string > minbiasExpr_
Definition: PFMETAnalyzer.h:134
void PFMETAnalyzer::bookMonitorElement ( std::string  DirName,
bool  bLumiSecPlot = false 
)

Definition at line 245 of file PFMETAnalyzer.cc.

References nbinsPV, PVlow, and PVup.

247 {
248  _dbe->setCurrentFolder(DirName);
249 
250 
251  mePfMEx = _dbe->book1D("METTask_PfMEx", "pfmet.px()", 200, -500, 500);
252  mePfMEy = _dbe->book1D("METTask_PfMEy", "pfmet.py()", 200, -500, 500);
253  mePfMET = _dbe->book1D("METTask_PfMET", "pfmet.pt()", 200, 0, 1000);
254  mePfSumET = _dbe->book1D("METTask_PfSumET", "pfmet.sumEt()", 400, 0, 4000);
255  mePfMETSig = _dbe->book1D("METTask_PfMETSig", "pfmet.mEtSig()", 51, 0, 51);
256  mePfMETPhi = _dbe->book1D("METTask_PfMETPhi", "pfmet.phi()", 60, -3.2, 3.2);
257  mePfMET_logx = _dbe->book1D("METTask_PfMET_logx", "log10(pfmet.pt())", 40, -1, 7);
258  mePfSumET_logx = _dbe->book1D("METTask_PfSumET_logx", "log10(pfmet.sumEt())", 40, -1, 7);
259 
260 
261  mePhotonEtFraction = _dbe->book1D("METTask_PfPhotonEtFraction", "pfmet.photonEtFraction()", 50, 0, 1);
262  mePhotonEt = _dbe->book1D("METTask_PfPhotonEt", "pfmet.photonEt()", 100, 0, 1000);
263  meNeutralHadronEtFraction = _dbe->book1D("METTask_PfNeutralHadronEtFraction", "pfmet.neutralHadronEtFraction()", 50, 0, 1);
264  meNeutralHadronEt = _dbe->book1D("METTask_PfNeutralHadronEt", "pfmet.neutralHadronEt()", 100, 0, 1000);
265  meElectronEtFraction = _dbe->book1D("METTask_PfElectronEtFraction", "pfmet.electronEtFraction()", 50, 0, 1);
266  meElectronEt = _dbe->book1D("METTask_PfElectronEt", "pfmet.electronEt()", 100, 0, 1000);
267  meChargedHadronEtFraction = _dbe->book1D("METTask_PfChargedHadronEtFraction", "pfmet.chargedHadronEtFraction()", 50, 0, 1);
268  meChargedHadronEt = _dbe->book1D("METTask_PfChargedHadronEt", "pfmet.chargedHadronEt()", 100, 0, 1000);
269  meMuonEtFraction = _dbe->book1D("METTask_PfMuonEtFraction", "pfmet.muonEtFraction()", 50, 0, 1);
270  meMuonEt = _dbe->book1D("METTask_PfMuonEt", "pfmet.muonEt()", 100, 0, 1000);
271  meHFHadronEtFraction = _dbe->book1D("METTask_PfHFHadronEtFraction", "pfmet.HFHadronEtFraction()", 50, 0, 1);
272  meHFHadronEt = _dbe->book1D("METTask_PfHFHadronEt", "pfmet.HFHadronEt()", 100, 0, 1000);
273  meHFEMEtFraction = _dbe->book1D("METTask_PfHFEMEtFraction", "pfmet.HFEMEtFraction()", 50, 0, 1);
274  meHFEMEt = _dbe->book1D("METTask_PfHFEMEt", "pfmet.HFEMEt()", 100, 0, 1000);
275 
276 
277  if (_allhist) {
278  if (bLumiSecPlot) {
279 
280  mePfMExLS = _dbe->book2D("METTask_PfMEx_LS", "METTask_PfMEx_LS", 200, -200, 200, 50, 0, 500);
281  mePfMEyLS = _dbe->book2D("METTask_PfMEy_LS", "METTask_PfMEy_LS", 200, -200, 200, 50, 0, 500);
282 
283  mePfMExLS->setAxisTitle("pfmet.px()", 1);
284  mePfMEyLS->setAxisTitle("pfmet.px()", 1);
285 
286  mePfMExLS->setAxisTitle("event.luminosityBlock()", 2);
287  mePfMEyLS->setAxisTitle("event.luminosityBlock()", 2);
288  }
289  }
290 
291 
292  // Book NPV profiles
293  //----------------------------------------------------------------------------
294  mePfMEx_profile = _dbe->bookProfile("METTask_PfMEx_profile", "pfmet.px()", nbinsPV, PVlow, PVup, 200, -500, 500);
295  mePfMEy_profile = _dbe->bookProfile("METTask_PfMEy_profile", "pfmet.py()", nbinsPV, PVlow, PVup, 200, -500, 500);
296  mePfMET_profile = _dbe->bookProfile("METTask_PfMET_profile", "pfmet.pt()", nbinsPV, PVlow, PVup, 200, 0, 1000);
297  mePfSumET_profile = _dbe->bookProfile("METTask_PfSumET_profile", "pfmet.sumEt()", nbinsPV, PVlow, PVup, 400, 0, 4000);
298 
299  mePhotonEtFraction_profile = _dbe->bookProfile("METTask_PfPhotonEtFraction_profile", "pfmet.photonEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
300  mePhotonEt_profile = _dbe->bookProfile("METTask_PfPhotonEt_profile", "pfmet.photonEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
301  meNeutralHadronEtFraction_profile = _dbe->bookProfile("METTask_PfNeutralHadronEtFraction_profile", "pfmet.neutralHadronEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
302  meNeutralHadronEt_profile = _dbe->bookProfile("METTask_PfNeutralHadronEt_profile", "pfmet.neutralHadronEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
303  meElectronEtFraction_profile = _dbe->bookProfile("METTask_PfElectronEtFraction_profile", "pfmet.electronEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
304  meElectronEt_profile = _dbe->bookProfile("METTask_PfElectronEt_profile", "pfmet.electronEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
305  meChargedHadronEtFraction_profile = _dbe->bookProfile("METTask_PfChargedHadronEtFraction_profile", "pfmet.chargedHadronEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
306  meChargedHadronEt_profile = _dbe->bookProfile("METTask_PfChargedHadronEt_profile", "pfmet.chargedHadronEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
307  meMuonEtFraction_profile = _dbe->bookProfile("METTask_PfMuonEtFraction_profile", "pfmet.muonEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
308  meMuonEt_profile = _dbe->bookProfile("METTask_PfMuonEt_profile", "pfmet.muonEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
309  meHFHadronEtFraction_profile = _dbe->bookProfile("METTask_PfHFHadronEtFraction_profile", "pfmet.HFHadronEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
310  meHFHadronEt_profile = _dbe->bookProfile("METTask_PfHFHadronEt_profile", "pfmet.HFHadronEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
311  meHFEMEtFraction_profile = _dbe->bookProfile("METTask_PfHFEMEtFraction_profile", "pfmet.HFEMEtFraction()", nbinsPV, PVlow, PVup, 50, 0, 1);
312  meHFEMEt_profile = _dbe->bookProfile("METTask_PfHFEMEt_profile", "pfmet.HFEMEt()", nbinsPV, PVlow, PVup, 100, 0, 1000);
313 
314 
315  // Set NPV profiles x-axis title
316  //----------------------------------------------------------------------------
317  mePfMEx_profile ->setAxisTitle("nvtx", 1);
318  mePfMEy_profile ->setAxisTitle("nvtx", 1);
319  mePfMET_profile ->setAxisTitle("nvtx", 1);
320  mePfSumET_profile->setAxisTitle("nvtx", 1);
321 
322  mePhotonEtFraction_profile ->setAxisTitle("nvtx", 1);
323  mePhotonEt_profile ->setAxisTitle("nvtx", 1);
324  meNeutralHadronEtFraction_profile->setAxisTitle("nvtx", 1);
325  meNeutralHadronEt_profile ->setAxisTitle("nvtx", 1);
326  meElectronEtFraction_profile ->setAxisTitle("nvtx", 1);
327  meElectronEt_profile ->setAxisTitle("nvtx", 1);
328  meChargedHadronEtFraction_profile->setAxisTitle("nvtx", 1);
329  meChargedHadronEt_profile ->setAxisTitle("nvtx", 1);
330  meMuonEtFraction_profile ->setAxisTitle("nvtx", 1);
331  meMuonEt_profile ->setAxisTitle("nvtx", 1);
332  meHFHadronEtFraction_profile ->setAxisTitle("nvtx", 1);
333  meHFHadronEt_profile ->setAxisTitle("nvtx", 1);
334  meHFEMEtFraction_profile ->setAxisTitle("nvtx", 1);
335  meHFEMEt_profile ->setAxisTitle("nvtx", 1);
336 }
PFMETAnalyzer::mePhotonEtFraction
MonitorElement * mePhotonEtFraction
Definition: PFMETAnalyzer.h:218
PFMETAnalyzer::mePfMET_profile
MonitorElement * mePfMET_profile
Definition: PFMETAnalyzer.h:238
PFMETAnalyzer::meHFEMEt_profile
MonitorElement * meHFEMEt_profile
Definition: PFMETAnalyzer.h:254
PFMETAnalyzer::meHFEMEt
MonitorElement * meHFEMEt
Definition: PFMETAnalyzer.h:231
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
PFMETAnalyzer::meChargedHadronEtFraction_profile
MonitorElement * meChargedHadronEtFraction_profile
Definition: PFMETAnalyzer.h:247
PFMETAnalyzer::mePfMExLS
MonitorElement * mePfMExLS
Definition: PFMETAnalyzer.h:213
PFMETAnalyzer::meHFEMEtFraction_profile
MonitorElement * meHFEMEtFraction_profile
Definition: PFMETAnalyzer.h:253
PFMETAnalyzer::mePfSumET
MonitorElement * mePfSumET
Definition: PFMETAnalyzer.h:212
PFMETAnalyzer::mePfMEy
MonitorElement * mePfMEy
Definition: PFMETAnalyzer.h:208
nbinsPV
const int nbinsPV
Definition: GlobalVariables.h:5
PFMETAnalyzer::meMuonEt_profile
MonitorElement * meMuonEt_profile
Definition: PFMETAnalyzer.h:250
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
PFMETAnalyzer::mePfMEyLS
MonitorElement * mePfMEyLS
Definition: PFMETAnalyzer.h:214
PVup
const double PVup
Definition: GlobalVariables.h:7
PFMETAnalyzer::meHFHadronEtFraction_profile
MonitorElement * meHFHadronEtFraction_profile
Definition: PFMETAnalyzer.h:251
PFMETAnalyzer::meNeutralHadronEt
MonitorElement * meNeutralHadronEt
Definition: PFMETAnalyzer.h:221
PFMETAnalyzer::meChargedHadronEt
MonitorElement * meChargedHadronEt
Definition: PFMETAnalyzer.h:225
PFMETAnalyzer::meHFEMEtFraction
MonitorElement * meHFEMEtFraction
Definition: PFMETAnalyzer.h:230
PFMETAnalyzer::mePfMETSig
MonitorElement * mePfMETSig
Definition: PFMETAnalyzer.h:209
PFMETAnalyzer::mePfMEx
MonitorElement * mePfMEx
Definition: PFMETAnalyzer.h:207
PFMETAnalyzer::mePhotonEt
MonitorElement * mePhotonEt
Definition: PFMETAnalyzer.h:219
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1036
PFMETAnalyzer::meMuonEtFraction_profile
MonitorElement * meMuonEtFraction_profile
Definition: PFMETAnalyzer.h:249
PFMETAnalyzer::meChargedHadronEt_profile
MonitorElement * meChargedHadronEt_profile
Definition: PFMETAnalyzer.h:248
PFMETAnalyzer::meElectronEtFraction_profile
MonitorElement * meElectronEtFraction_profile
Definition: PFMETAnalyzer.h:245
PVlow
const double PVlow
Definition: GlobalVariables.h:6
PFMETAnalyzer::meElectronEtFraction
MonitorElement * meElectronEtFraction
Definition: PFMETAnalyzer.h:222
PFMETAnalyzer::meNeutralHadronEtFraction_profile
MonitorElement * meNeutralHadronEtFraction_profile
Definition: PFMETAnalyzer.h:243
PFMETAnalyzer::mePfSumET_logx
MonitorElement * mePfSumET_logx
Definition: PFMETAnalyzer.h:216
PFMETAnalyzer::mePhotonEtFraction_profile
MonitorElement * mePhotonEtFraction_profile
Definition: PFMETAnalyzer.h:241
PFMETAnalyzer::meChargedHadronEtFraction
MonitorElement * meChargedHadronEtFraction
Definition: PFMETAnalyzer.h:224
PFMETAnalyzer::meHFHadronEtFraction
MonitorElement * meHFHadronEtFraction
Definition: PFMETAnalyzer.h:228
PFMETAnalyzer::mePfMET_logx
MonitorElement * mePfMET_logx
Definition: PFMETAnalyzer.h:215
PFMETAnalyzer::mePhotonEt_profile
MonitorElement * mePhotonEt_profile
Definition: PFMETAnalyzer.h:242
PFMETAnalyzer::mePfMETPhi
MonitorElement * mePfMETPhi
Definition: PFMETAnalyzer.h:211
PFMETAnalyzer::meHFHadronEt
MonitorElement * meHFHadronEt
Definition: PFMETAnalyzer.h:229
PFMETAnalyzer::meElectronEt_profile
MonitorElement * meElectronEt_profile
Definition: PFMETAnalyzer.h:246
PFMETAnalyzer::meMuonEt
MonitorElement * meMuonEt
Definition: PFMETAnalyzer.h:227
PFMETAnalyzer::meHFHadronEt_profile
MonitorElement * meHFHadronEt_profile
Definition: PFMETAnalyzer.h:252
PFMETAnalyzer::mePfMET
MonitorElement * mePfMET
Definition: PFMETAnalyzer.h:210
PFMETAnalyzer::meNeutralHadronEt_profile
MonitorElement * meNeutralHadronEt_profile
Definition: PFMETAnalyzer.h:244
PFMETAnalyzer::mePfSumET_profile
MonitorElement * mePfSumET_profile
Definition: PFMETAnalyzer.h:239
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:850
PFMETAnalyzer::meMuonEtFraction
MonitorElement * meMuonEtFraction
Definition: PFMETAnalyzer.h:226
PFMETAnalyzer::mePfMEy_profile
MonitorElement * mePfMEy_profile
Definition: PFMETAnalyzer.h:237
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:916
PFMETAnalyzer::mePfMEx_profile
MonitorElement * mePfMEx_profile
Definition: PFMETAnalyzer.h:236
PFMETAnalyzer::meNeutralHadronEtFraction
MonitorElement * meNeutralHadronEtFraction
Definition: PFMETAnalyzer.h:220
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
PFMETAnalyzer::meElectronEt
MonitorElement * meElectronEt
Definition: PFMETAnalyzer.h:223
void PFMETAnalyzer::endJob ( void  )

Finish up a job.

Definition at line 191 of file PFMETAnalyzer.cc.

References bTagSequences_cff::jetID.

191  {
192 
193  delete jetID;
194  delete DCSFilter;
195 
196 }
PFMETAnalyzer::DCSFilter
JetMETDQMDCSFilter * DCSFilter
Definition: PFMETAnalyzer.h:182
PFMETAnalyzer::jetID
reco::helper::JetIDHelper * jetID
Definition: PFMETAnalyzer.h:179
void PFMETAnalyzer::endRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup,
DQMStore *  dbe 
)

Finish up a run.

Definition at line 375 of file PFMETAnalyzer.cc.

References cmsPerfPublish::DirName, TrackerOfflineValidation_Dqm_cff::dirName, MonitorElement::getRootObject(), MonitorElement::getTH1F(), i, and AlCaHLTBitMon_QueryRunRegistry::string.

376 {
377 
378  //
379  //--- Check the time length of the Run from the lumi section plots
380 
381  std::string dirName = "JetMET/MET/"+_source+"/";
382  _dbe->setCurrentFolder(dirName);
383 
384  TH1F* tlumisec;
385 
386  MonitorElement *meLumiSec = _dbe->get("aaa");
387  meLumiSec = _dbe->get("JetMET/lumisec");
388 
389  int totlsec=0;
390  double totltime=0.;
391  if ( meLumiSec->getRootObject() ) {
392  tlumisec = meLumiSec->getTH1F();
393  for (int i=0; i<500; i++){
394  if (tlumisec->GetBinContent(i+1)) totlsec++;
395  }
396  totltime = double(totlsec*90); // one lumi sec ~ 90 (sec)
397  }
398 
399  if (totltime==0.) totltime=1.;
400 
401  //
402  //--- Make the integrated plots with rate (Hz)
403 
404  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin(); ic != _FolderNames.end(); ic++)
405  {
406 
407  std::string DirName;
408  DirName = dirName+*ic;
409 
410  makeRatePlot(DirName,totltime);
411  if ( _HighPtJetEventFlag->on() )
412  makeRatePlot(DirName+"/"+"triggerName_HighJetPt",totltime);
413  if ( _LowPtJetEventFlag->on() )
414  makeRatePlot(DirName+"/"+"triggerName_LowJetPt",totltime);
415  if ( _MinBiasEventFlag->on() )
416  makeRatePlot(DirName+"/"+"triggerName_MinBias",totltime);
417  if ( _HighMETEventFlag->on() )
418  makeRatePlot(DirName+"/"+"triggerName_HighMET",totltime);
419  if ( _EleEventFlag->on() )
420  makeRatePlot(DirName+"/"+"triggerName_Ele",totltime);
421  if ( _MuonEventFlag->on() )
422  makeRatePlot(DirName+"/"+"triggerName_Muon",totltime);
423  }
424 }
i
int i
Definition: DBlmapReader.cc:9
PFMETAnalyzer::_MuonEventFlag
GenericTriggerEventFlag * _MuonEventFlag
Definition: PFMETAnalyzer.h:126
PFMETAnalyzer::_source
std::string _source
Definition: PFMETAnalyzer.h:105
PFMETAnalyzer::_HighPtJetEventFlag
GenericTriggerEventFlag * _HighPtJetEventFlag
Definition: PFMETAnalyzer.h:120
PFMETAnalyzer::_LowPtJetEventFlag
GenericTriggerEventFlag * _LowPtJetEventFlag
Definition: PFMETAnalyzer.h:121
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
PFMETAnalyzer::_EleEventFlag
GenericTriggerEventFlag * _EleEventFlag
Definition: PFMETAnalyzer.h:125
PFMETAnalyzer::_HighMETEventFlag
GenericTriggerEventFlag * _HighMETEventFlag
Definition: PFMETAnalyzer.h:123
TrackerOfflineValidation_Dqm_cff.dirName
tuple dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:56
PFMETAnalyzer::makeRatePlot
void makeRatePlot(std::string, double)
Definition: PFMETAnalyzer.cc:428
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1473
PFMETAnalyzer::_FolderNames
std::vector< std::string > _FolderNames
Definition: PFMETAnalyzer.h:189
MonitorElement::getRootObject
TObject * getRootObject(void) const
Definition: MonitorElement.cc:1472
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
PFMETAnalyzer::_MinBiasEventFlag
GenericTriggerEventFlag * _MinBiasEventFlag
Definition: PFMETAnalyzer.h:122
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
void PFMETAnalyzer::fillMESet ( const edm::Event &  iEvent,
std::string  DirName,
const reco::PFMET &  pfmet 
)

Definition at line 905 of file PFMETAnalyzer.cc.

907 {
908 
909  _dbe->setCurrentFolder(DirName);
910 
911  bool bLumiSecPlot=false;
912  if (DirName.find("All")) bLumiSecPlot=true;
913 
914  if (_trig_JetMB)
915  fillMonitorElement(iEvent,DirName,"",pfmet, bLumiSecPlot);
916  if (_trig_HighPtJet)
917  fillMonitorElement(iEvent,DirName,"HighPtJet",pfmet,false);
918  if (_trig_LowPtJet)
919  fillMonitorElement(iEvent,DirName,"LowPtJet",pfmet,false);
920  if (_trig_MinBias)
921  fillMonitorElement(iEvent,DirName,"MinBias",pfmet,false);
922  if (_trig_HighMET)
923  fillMonitorElement(iEvent,DirName,"HighMET",pfmet,false);
924  if (_trig_Ele)
925  fillMonitorElement(iEvent,DirName,"Ele",pfmet,false);
926  if (_trig_Muon)
927  fillMonitorElement(iEvent,DirName,"Muon",pfmet,false);
928 }
PFMETAnalyzer::fillMonitorElement
void fillMonitorElement(const edm::Event &, std::string, std::string, const reco::PFMET &, bool)
Definition: PFMETAnalyzer.cc:934
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
void PFMETAnalyzer::fillMonitorElement ( const edm::Event &  iEvent,
std::string  DirName,
std::string  TriggerTypeName,
const reco::PFMET &  pfmet,
bool  bLumiSecPlot 
)

Definition at line 934 of file PFMETAnalyzer.cc.

References reco::PFMET::chargedHadronEt(), reco::PFMET::chargedHadronEtFraction(), reco::PFMET::electronEt(), reco::PFMET::electronEtFraction(), reco::PFMET::HFEMEt(), reco::PFMET::HFEMEtFraction(), reco::PFMET::HFHadronEt(), reco::PFMET::HFHadronEtFraction(), edm::EventBase::luminosityBlock(), reco::MET::mEtSig(), reco::PFMET::muonEt(), reco::PFMET::muonEtFraction(), reco::PFMET::neutralHadronEt(), reco::PFMET::neutralHadronEtFraction(), pfMET_cfi::pfMET, reco::LeafCandidate::phi(), reco::PFMET::photonEt(), reco::PFMET::photonEtFraction(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::MET::sumEt().

939 {
940  if (TriggerTypeName == "HighPtJet") {
941  if (!selectHighPtJetEvent(iEvent)) return;
942  }
943  else if (TriggerTypeName == "LowPtJet") {
944  if (!selectLowPtJetEvent(iEvent)) return;
945  }
946  else if (TriggerTypeName == "HighMET") {
947  if (pfmet.pt()<_highPFMETThreshold) return;
948  }
949  else if (TriggerTypeName == "Ele") {
950  if (!selectWElectronEvent(iEvent)) return;
951  }
952  else if (TriggerTypeName == "Muon") {
953  if (!selectWMuonEvent(iEvent)) return;
954  }
955 
956  if (TriggerTypeName != "") DirName = DirName + "/" + TriggerTypeName;
957 
958 
959  // Reconstructed PFMET information
960  //----------------------------------------------------------------------------
961  double pfSumET = pfmet.sumEt();
962  double pfMETSig = pfmet.mEtSig();
963  double pfMET = pfmet.pt();
964  double pfMEx = pfmet.px();
965  double pfMEy = pfmet.py();
966  double pfMETPhi = pfmet.phi();
967 
968 
969  // PFMET getters
970  //----------------------------------------------------------------------------
971  double pfPhotonEtFraction = pfmet.photonEtFraction();
972  double pfPhotonEt = pfmet.photonEt();
973  double pfNeutralHadronEtFraction = pfmet.neutralHadronEtFraction();
974  double pfNeutralHadronEt = pfmet.neutralHadronEt();
975  double pfElectronEtFraction = pfmet.electronEtFraction();
976  double pfElectronEt = pfmet.electronEt();
977  double pfChargedHadronEtFraction = pfmet.chargedHadronEtFraction();
978  double pfChargedHadronEt = pfmet.chargedHadronEt();
979  double pfMuonEtFraction = pfmet.muonEtFraction();
980  double pfMuonEt = pfmet.muonEt();
981  double pfHFHadronEtFraction = pfmet.HFHadronEtFraction();
982  double pfHFHadronEt = pfmet.HFHadronEt();
983  double pfHFEMEtFraction = pfmet.HFEMEtFraction();
984  double pfHFEMEt = pfmet.HFEMEt();
985 
986 
987  if (pfSumET > _etThreshold) {
988 
989  mePfMEx = _dbe->get(DirName + "/METTask_PfMEx");
990  mePfMEy = _dbe->get(DirName + "/METTask_PfMEy");
991  mePfMET = _dbe->get(DirName + "/METTask_PfMET");
992  mePfMETPhi = _dbe->get(DirName + "/METTask_PfMETPhi");
993  mePfSumET = _dbe->get(DirName + "/METTask_PfSumET");
994  mePfMETSig = _dbe->get(DirName + "/METTask_PfMETSig");
995  mePfMET_logx = _dbe->get(DirName + "/METTask_PfMET_logx");
996  mePfSumET_logx = _dbe->get(DirName + "/METTask_PfSumET_logx");
997 
998  if (mePfMEx && mePfMEx ->getRootObject()) mePfMEx ->Fill(pfMEx);
999  if (mePfMEy && mePfMEy ->getRootObject()) mePfMEy ->Fill(pfMEy);
1000  if (mePfMET && mePfMET ->getRootObject()) mePfMET ->Fill(pfMET);
1001  if (mePfMETPhi && mePfMETPhi ->getRootObject()) mePfMETPhi ->Fill(pfMETPhi);
1002  if (mePfSumET && mePfSumET ->getRootObject()) mePfSumET ->Fill(pfSumET);
1003  if (mePfMETSig && mePfMETSig ->getRootObject()) mePfMETSig ->Fill(pfMETSig);
1004  if (mePfMET_logx && mePfMET_logx ->getRootObject()) mePfMET_logx ->Fill(log10(pfMET));
1005  if (mePfSumET_logx && mePfSumET_logx->getRootObject()) mePfSumET_logx->Fill(log10(pfSumET));
1006 
1007 
1008  mePhotonEtFraction = _dbe->get(DirName + "/METTask_PfPhotonEtFraction");
1009  mePhotonEt = _dbe->get(DirName + "/METTask_PfPhotonEt");
1010  meNeutralHadronEtFraction = _dbe->get(DirName + "/METTask_PfNeutralHadronEtFraction");
1011  meNeutralHadronEt = _dbe->get(DirName + "/METTask_PfNeutralHadronEt");
1012  meElectronEtFraction = _dbe->get(DirName + "/METTask_PfElectronEtFraction");
1013  meElectronEt = _dbe->get(DirName + "/METTask_PfElectronEt");
1014  meChargedHadronEtFraction = _dbe->get(DirName + "/METTask_PfChargedHadronEtFraction");
1015  meChargedHadronEt = _dbe->get(DirName + "/METTask_PfChargedHadronEt");
1016  meMuonEtFraction = _dbe->get(DirName + "/METTask_PfMuonEtFraction");
1017  meMuonEt = _dbe->get(DirName + "/METTask_PfMuonEt");
1018  meHFHadronEtFraction = _dbe->get(DirName + "/METTask_PfHFHadronEtFraction");
1019  meHFHadronEt = _dbe->get(DirName + "/METTask_PfHFHadronEt");
1020  meHFEMEtFraction = _dbe->get(DirName + "/METTask_PfHFEMEtFraction");
1021  meHFEMEt = _dbe->get(DirName + "/METTask_PfHFEMEt");
1022 
1023  if (mePhotonEtFraction && mePhotonEtFraction ->getRootObject()) mePhotonEtFraction ->Fill(pfPhotonEtFraction);
1024  if (mePhotonEt && mePhotonEt ->getRootObject()) mePhotonEt ->Fill(pfPhotonEt);
1025  if (meNeutralHadronEtFraction && meNeutralHadronEtFraction->getRootObject()) meNeutralHadronEtFraction->Fill(pfNeutralHadronEtFraction);
1026  if (meNeutralHadronEt && meNeutralHadronEt ->getRootObject()) meNeutralHadronEt ->Fill(pfNeutralHadronEt);
1027  if (meElectronEtFraction && meElectronEtFraction ->getRootObject()) meElectronEtFraction ->Fill(pfElectronEtFraction);
1028  if (meElectronEt && meElectronEt ->getRootObject()) meElectronEt ->Fill(pfElectronEt);
1029  if (meChargedHadronEtFraction && meChargedHadronEtFraction->getRootObject()) meChargedHadronEtFraction->Fill(pfChargedHadronEtFraction);
1030  if (meChargedHadronEt && meChargedHadronEt ->getRootObject()) meChargedHadronEt ->Fill(pfChargedHadronEt);
1031  if (meMuonEtFraction && meMuonEtFraction ->getRootObject()) meMuonEtFraction ->Fill(pfMuonEtFraction);
1032  if (meMuonEt && meMuonEt ->getRootObject()) meMuonEt ->Fill(pfMuonEt);
1033  if (meHFHadronEtFraction && meHFHadronEtFraction ->getRootObject()) meHFHadronEtFraction ->Fill(pfHFHadronEtFraction);
1034  if (meHFHadronEt && meHFHadronEt ->getRootObject()) meHFHadronEt ->Fill(pfHFHadronEt);
1035  if (meHFEMEtFraction && meHFEMEtFraction ->getRootObject()) meHFEMEtFraction ->Fill(pfHFEMEtFraction);
1036  if (meHFEMEt && meHFEMEt ->getRootObject()) meHFEMEt ->Fill(pfHFEMEt);
1037 
1038 
1039  if (_allhist) {
1040  if (bLumiSecPlot) {
1041 
1042  mePfMExLS = _dbe->get(DirName + "/METTask_PfMExLS");
1043  mePfMEyLS = _dbe->get(DirName + "/METTask_PfMEyLS");
1044 
1045  if (mePfMExLS && mePfMExLS->getRootObject()) mePfMExLS->Fill(pfMEx, iEvent.luminosityBlock());
1046  if (mePfMEyLS && mePfMEyLS->getRootObject()) mePfMEyLS->Fill(pfMEy, iEvent.luminosityBlock());
1047  }
1048  }
1049 
1050 
1051  // Fill NPV profiles
1052  //--------------------------------------------------------------------------
1053  mePfMEx_profile = _dbe->get(DirName + "/METTask_PfMEx_profile");
1054  mePfMEy_profile = _dbe->get(DirName + "/METTask_PfMEy_profile");
1055  mePfMET_profile = _dbe->get(DirName + "/METTask_PfMET_profile");
1056  mePfSumET_profile = _dbe->get(DirName + "/METTask_PfSumET_profile");
1057 
1058  if (mePfMEx_profile && mePfMEx_profile ->getRootObject()) mePfMEx_profile ->Fill(_numPV, pfMEx);
1059  if (mePfMEy_profile && mePfMEy_profile ->getRootObject()) mePfMEy_profile ->Fill(_numPV, pfMEy);
1060  if (mePfMET_profile && mePfMET_profile ->getRootObject()) mePfMET_profile ->Fill(_numPV, pfMET);
1061  if (mePfSumET_profile && mePfSumET_profile->getRootObject()) mePfSumET_profile->Fill(_numPV, pfSumET);
1062 
1063 
1064  mePhotonEtFraction_profile = _dbe->get(DirName + "/METTask_PfPhotonEtFraction_profile");
1065  mePhotonEt_profile = _dbe->get(DirName + "/METTask_PfPhotonEt_profile");
1066  meNeutralHadronEtFraction_profile = _dbe->get(DirName + "/METTask_PfNeutralHadronEtFraction_profile");
1067  meNeutralHadronEt_profile = _dbe->get(DirName + "/METTask_PfNeutralHadronEt_profile");
1068  meElectronEtFraction_profile = _dbe->get(DirName + "/METTask_PfElectronEtFraction_profile");
1069  meElectronEt_profile = _dbe->get(DirName + "/METTask_PfElectronEt_profile");
1070  meChargedHadronEtFraction_profile = _dbe->get(DirName + "/METTask_PfChargedHadronEtFraction_profile");
1071  meChargedHadronEt_profile = _dbe->get(DirName + "/METTask_PfChargedHadronEt_profile");
1072  meMuonEtFraction_profile = _dbe->get(DirName + "/METTask_PfMuonEtFraction_profile");
1073  meMuonEt_profile = _dbe->get(DirName + "/METTask_PfMuonEt_profile");
1074  meHFHadronEtFraction_profile = _dbe->get(DirName + "/METTask_PfHFHadronEtFraction_profile");
1075  meHFHadronEt_profile = _dbe->get(DirName + "/METTask_PfHFHadronEt_profile");
1076  meHFEMEtFraction_profile = _dbe->get(DirName + "/METTask_PfHFEMEtFraction_profile");
1077  meHFEMEt_profile = _dbe->get(DirName + "/METTask_PfHFEMEt_profile");
1078 
1079  if (mePhotonEtFraction_profile && mePhotonEtFraction_profile ->getRootObject()) mePhotonEtFraction_profile ->Fill(_numPV, pfPhotonEtFraction);
1080  if (mePhotonEt_profile && mePhotonEt_profile ->getRootObject()) mePhotonEt_profile ->Fill(_numPV, pfPhotonEt);
1081  if (meNeutralHadronEtFraction_profile && meNeutralHadronEtFraction_profile->getRootObject()) meNeutralHadronEtFraction_profile->Fill(_numPV, pfNeutralHadronEtFraction);
1082  if (meNeutralHadronEt_profile && meNeutralHadronEt_profile ->getRootObject()) meNeutralHadronEt_profile ->Fill(_numPV, pfNeutralHadronEt);
1083  if (meElectronEtFraction_profile && meElectronEtFraction_profile ->getRootObject()) meElectronEtFraction_profile ->Fill(_numPV, pfElectronEtFraction);
1084  if (meElectronEt_profile && meElectronEt_profile ->getRootObject()) meElectronEt_profile ->Fill(_numPV, pfElectronEt);
1085  if (meChargedHadronEtFraction_profile && meChargedHadronEtFraction_profile->getRootObject()) meChargedHadronEtFraction_profile->Fill(_numPV, pfChargedHadronEtFraction);
1086  if (meChargedHadronEt_profile && meChargedHadronEt_profile ->getRootObject()) meChargedHadronEt_profile ->Fill(_numPV, pfChargedHadronEt);
1087  if (meMuonEtFraction_profile && meMuonEtFraction_profile ->getRootObject()) meMuonEtFraction_profile ->Fill(_numPV, pfMuonEtFraction);
1088  if (meMuonEt_profile && meMuonEt_profile ->getRootObject()) meMuonEt_profile ->Fill(_numPV, pfMuonEt);
1089  if (meHFHadronEtFraction_profile && meHFHadronEtFraction_profile ->getRootObject()) meHFHadronEtFraction_profile ->Fill(_numPV, pfHFHadronEtFraction);
1090  if (meHFHadronEt_profile && meHFHadronEt_profile ->getRootObject()) meHFHadronEt_profile ->Fill(_numPV, pfHFHadronEt);
1091  if (meHFEMEtFraction_profile && meHFEMEtFraction_profile ->getRootObject()) meHFEMEtFraction_profile ->Fill(_numPV, pfHFEMEtFraction);
1092  if (meHFEMEt_profile && meHFEMEt_profile ->getRootObject()) meHFEMEt_profile ->Fill(_numPV, pfHFEMEt);
1093  }
1094 }
PFMETAnalyzer::mePhotonEtFraction
MonitorElement * mePhotonEtFraction
Definition: PFMETAnalyzer.h:218
PFMETAnalyzer::mePfMET_profile
MonitorElement * mePfMET_profile
Definition: PFMETAnalyzer.h:238
PFMETAnalyzer::_highPFMETThreshold
double _highPFMETThreshold
Definition: PFMETAnalyzer.h:168
PFMETAnalyzer::meHFEMEt_profile
MonitorElement * meHFEMEt_profile
Definition: PFMETAnalyzer.h:254
PFMETAnalyzer::meHFEMEt
MonitorElement * meHFEMEt
Definition: PFMETAnalyzer.h:231
reco::PFMET::HFEMEtFraction
double HFEMEtFraction() const
Definition: PFMET.h:47
PFMETAnalyzer::meChargedHadronEtFraction_profile
MonitorElement * meChargedHadronEtFraction_profile
Definition: PFMETAnalyzer.h:247
reco::PFMET::neutralHadronEtFraction
double neutralHadronEtFraction() const
Definition: PFMET.h:32
reco::PFMET::muonEt
double muonEt() const
Definition: PFMET.h:42
PFMETAnalyzer::mePfMExLS
MonitorElement * mePfMExLS
Definition: PFMETAnalyzer.h:213
PFMETAnalyzer::selectWMuonEvent
bool selectWMuonEvent(const edm::Event &)
Definition: PFMETAnalyzer.cc:1158
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:59
PFMETAnalyzer::meHFEMEtFraction_profile
MonitorElement * meHFEMEtFraction_profile
Definition: PFMETAnalyzer.h:253
PFMETAnalyzer::mePfSumET
MonitorElement * mePfSumET
Definition: PFMETAnalyzer.h:212
PFMETAnalyzer::mePfMEy
MonitorElement * mePfMEy
Definition: PFMETAnalyzer.h:208
pfMET_cfi.pfMET
tuple pfMET
Definition: pfMET_cfi.py:7
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:149
PFMETAnalyzer::meMuonEt_profile
MonitorElement * meMuonEt_profile
Definition: PFMETAnalyzer.h:250
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
PFMETAnalyzer::mePfMEyLS
MonitorElement * mePfMEyLS
Definition: PFMETAnalyzer.h:214
PFMETAnalyzer::meHFHadronEtFraction_profile
MonitorElement * meHFHadronEtFraction_profile
Definition: PFMETAnalyzer.h:251
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:155
reco::MET::mEtSig
double mEtSig() const
Definition: MET.h:50
PFMETAnalyzer::meNeutralHadronEt
MonitorElement * meNeutralHadronEt
Definition: PFMETAnalyzer.h:221
reco::PFMET::HFHadronEtFraction
double HFHadronEtFraction() const
Definition: PFMET.h:44
reco::MET::sumEt
double sumEt() const
Definition: MET.h:48
PFMETAnalyzer::meChargedHadronEt
MonitorElement * meChargedHadronEt
Definition: PFMETAnalyzer.h:225
reco::PFMET::muonEtFraction
double muonEtFraction() const
Definition: PFMET.h:41
PFMETAnalyzer::meHFEMEtFraction
MonitorElement * meHFEMEtFraction
Definition: PFMETAnalyzer.h:230
reco::PFMET::HFHadronEt
double HFHadronEt() const
Definition: PFMET.h:45
reco::PFMET::photonEtFraction
double photonEtFraction() const
Definition: PFMET.h:29
PFMETAnalyzer::mePfMETSig
MonitorElement * mePfMETSig
Definition: PFMETAnalyzer.h:209
PFMETAnalyzer::mePfMEx
MonitorElement * mePfMEx
Definition: PFMETAnalyzer.h:207
reco::LeafCandidate::px
virtual double px() const GCC11_FINAL
x coordinate of momentum vector
Definition: LeafCandidate.h:147
PFMETAnalyzer::selectLowPtJetEvent
bool selectLowPtJetEvent(const edm::Event &)
Definition: PFMETAnalyzer.cc:1122
PFMETAnalyzer::mePhotonEt
MonitorElement * mePhotonEt
Definition: PFMETAnalyzer.h:219
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1473
PFMETAnalyzer::meMuonEtFraction_profile
MonitorElement * meMuonEtFraction_profile
Definition: PFMETAnalyzer.h:249
PFMETAnalyzer::meChargedHadronEt_profile
MonitorElement * meChargedHadronEt_profile
Definition: PFMETAnalyzer.h:248
reco::PFMET::HFEMEt
double HFEMEt() const
Definition: PFMET.h:48
PFMETAnalyzer::meElectronEtFraction_profile
MonitorElement * meElectronEtFraction_profile
Definition: PFMETAnalyzer.h:245
PFMETAnalyzer::meElectronEtFraction
MonitorElement * meElectronEtFraction
Definition: PFMETAnalyzer.h:222
PFMETAnalyzer::meNeutralHadronEtFraction_profile
MonitorElement * meNeutralHadronEtFraction_profile
Definition: PFMETAnalyzer.h:243
PFMETAnalyzer::mePfSumET_logx
MonitorElement * mePfSumET_logx
Definition: PFMETAnalyzer.h:216
PFMETAnalyzer::mePhotonEtFraction_profile
MonitorElement * mePhotonEtFraction_profile
Definition: PFMETAnalyzer.h:241
reco::PFMET::electronEt
double electronEt() const
Definition: PFMET.h:36
PFMETAnalyzer::_etThreshold
double _etThreshold
Definition: PFMETAnalyzer.h:173
MonitorElement::getRootObject
TObject * getRootObject(void) const
Definition: MonitorElement.cc:1472
PFMETAnalyzer::meChargedHadronEtFraction
MonitorElement * meChargedHadronEtFraction
Definition: PFMETAnalyzer.h:224
PFMETAnalyzer::selectHighPtJetEvent
bool selectHighPtJetEvent(const edm::Event &)
Definition: PFMETAnalyzer.cc:1100
reco::PFMET::photonEt
double photonEt() const
Definition: PFMET.h:30
PFMETAnalyzer::meHFHadronEtFraction
MonitorElement * meHFHadronEtFraction
Definition: PFMETAnalyzer.h:228
PFMETAnalyzer::mePfMET_logx
MonitorElement * mePfMET_logx
Definition: PFMETAnalyzer.h:215
PFMETAnalyzer::mePhotonEt_profile
MonitorElement * mePhotonEt_profile
Definition: PFMETAnalyzer.h:242
PFMETAnalyzer::mePfMETPhi
MonitorElement * mePfMETPhi
Definition: PFMETAnalyzer.h:211
PFMETAnalyzer::selectWElectronEvent
bool selectWElectronEvent(const edm::Event &)
Definition: PFMETAnalyzer.cc:1145
PFMETAnalyzer::meHFHadronEt
MonitorElement * meHFHadronEt
Definition: PFMETAnalyzer.h:229
PFMETAnalyzer::meElectronEt_profile
MonitorElement * meElectronEt_profile
Definition: PFMETAnalyzer.h:246
reco::PFMET::electronEtFraction
double electronEtFraction() const
Definition: PFMET.h:35
PFMETAnalyzer::meMuonEt
MonitorElement * meMuonEt
Definition: PFMETAnalyzer.h:227
reco::PFMET::chargedHadronEtFraction
double chargedHadronEtFraction() const
Definition: PFMET.h:38
PFMETAnalyzer::meHFHadronEt_profile
MonitorElement * meHFHadronEt_profile
Definition: PFMETAnalyzer.h:252
PFMETAnalyzer::mePfMET
MonitorElement * mePfMET
Definition: PFMETAnalyzer.h:210
PFMETAnalyzer::meNeutralHadronEt_profile
MonitorElement * meNeutralHadronEt_profile
Definition: PFMETAnalyzer.h:244
PFMETAnalyzer::mePfSumET_profile
MonitorElement * mePfSumET_profile
Definition: PFMETAnalyzer.h:239
reco::PFMET::neutralHadronEt
double neutralHadronEt() const
Definition: PFMET.h:33
PFMETAnalyzer::meMuonEtFraction
MonitorElement * meMuonEtFraction
Definition: PFMETAnalyzer.h:226
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:153
PFMETAnalyzer::mePfMEy_profile
MonitorElement * mePfMEy_profile
Definition: PFMETAnalyzer.h:237
PFMETAnalyzer::mePfMEx_profile
MonitorElement * mePfMEx_profile
Definition: PFMETAnalyzer.h:236
PFMETAnalyzer::meNeutralHadronEtFraction
MonitorElement * meNeutralHadronEtFraction
Definition: PFMETAnalyzer.h:220
PFMETAnalyzer::meElectronEt
MonitorElement * meElectronEt
Definition: PFMETAnalyzer.h:223
reco::PFMET::chargedHadronEt
double chargedHadronEt() const
Definition: PFMET.h:39
void PFMETAnalyzer::makeRatePlot ( std::string  DirName,
double  totltime 
)

Definition at line 428 of file PFMETAnalyzer.cc.

References MonitorElement::getRootObject(), MonitorElement::getTH1F(), and i.

429 {
430 
431  _dbe->setCurrentFolder(DirName);
432  MonitorElement *mePfMET = _dbe->get(DirName+"/"+"METTask_PfMET");
433 
434  TH1F* tPfMET;
435  TH1F* tPfMETRate;
436 
437  if ( mePfMET )
438  if ( mePfMET->getRootObject() ) {
439  tPfMET = mePfMET->getTH1F();
440 
441  // Integral plot & convert number of events to rate (hz)
442  tPfMETRate = (TH1F*) tPfMET->Clone("METTask_PfMETRate");
443  for (int i = tPfMETRate->GetNbinsX()-1; i>=0; i--){
444  tPfMETRate->SetBinContent(i+1,tPfMETRate->GetBinContent(i+2)+tPfMET->GetBinContent(i+1));
445  }
446  for (int i = 0; i<tPfMETRate->GetNbinsX(); i++){
447  tPfMETRate->SetBinContent(i+1,tPfMETRate->GetBinContent(i+1)/double(totltime));
448  }
449 
450  tPfMETRate->SetName("METTask_PfMETRate");
451  tPfMETRate->SetTitle("METTask_PfMETRate");
452  mePfMETRate = _dbe->book1D("METTask_PfMETRate",tPfMETRate);
453 
454  }
455 }
i
int i
Definition: DBlmapReader.cc:9
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
PFMETAnalyzer::_dbe
DQMStore * _dbe
Definition: PFMETAnalyzer.h:191
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1473
PFMETAnalyzer::mePfMETRate
MonitorElement * mePfMETRate
Definition: PFMETAnalyzer.h:197
MonitorElement::getRootObject
TObject * getRootObject(void) const
Definition: MonitorElement.cc:1472
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1486
PFMETAnalyzer::mePfMET
MonitorElement * mePfMET
Definition: PFMETAnalyzer.h:210
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
bool PFMETAnalyzer::selectHighPtJetEvent ( const edm::Event &  iEvent)

Definition at line 1100 of file PFMETAnalyzer.cc.

References gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), LogDebug, and pfJets_cff::pfJets.

1100  {
1101 
1102  bool return_value=false;
1103 
1104  edm::Handle<reco::PFJetCollection> pfJets;
1105  iEvent.getByLabel(thePfJetCollectionLabel, pfJets);
1106  if (!pfJets.isValid()) {
1107  LogDebug("") << "PFMETAnalyzer: Could not find pfjet product" << std::endl;
1108  if (_verbose) std::cout << "PFMETAnalyzer: Could not find pfjet product" << std::endl;
1109  }
1110 
1111  for (reco::PFJetCollection::const_iterator pf = pfJets->begin();
1112  pf!=pfJets->end(); ++pf){
1113  if (pf->pt()>_highPtPFJetThreshold){
1114  return_value=true;
1115  }
1116  }
1117 
1118  return return_value;
1119 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
PFMETAnalyzer::_highPtPFJetThreshold
double _highPtPFJetThreshold
Definition: PFMETAnalyzer.h:166
PFMETAnalyzer::thePfJetCollectionLabel
edm::InputTag thePfJetCollectionLabel
Definition: PFMETAnalyzer.h:110
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
pfJets_cff.pfJets
tuple pfJets
Definition: pfJets_cff.py:8
bool PFMETAnalyzer::selectLowPtJetEvent ( const edm::Event &  iEvent)

Definition at line 1122 of file PFMETAnalyzer.cc.

References gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), LogDebug, and pfJets_cff::pfJets.

1122  {
1123 
1124  bool return_value=false;
1125 
1126  edm::Handle<reco::PFJetCollection> pfJets;
1127  iEvent.getByLabel(thePfJetCollectionLabel, pfJets);
1128  if (!pfJets.isValid()) {
1129  LogDebug("") << "PFMETAnalyzer: Could not find jet product" << std::endl;
1130  if (_verbose) std::cout << "PFMETAnalyzer: Could not find jet product" << std::endl;
1131  }
1132 
1133  for (reco::PFJetCollection::const_iterator cal = pfJets->begin();
1134  cal!=pfJets->end(); ++cal){
1135  if (cal->pt()>_lowPtPFJetThreshold){
1136  return_value=true;
1137  }
1138  }
1139 
1140  return return_value;
1141 
1142 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
PFMETAnalyzer::_lowPtPFJetThreshold
double _lowPtPFJetThreshold
Definition: PFMETAnalyzer.h:167
PFMETAnalyzer::thePfJetCollectionLabel
edm::InputTag thePfJetCollectionLabel
Definition: PFMETAnalyzer.h:110
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
pfJets_cff.pfJets
tuple pfJets
Definition: pfJets_cff.py:8
bool PFMETAnalyzer::selectWElectronEvent ( const edm::Event &  iEvent)

Definition at line 1145 of file PFMETAnalyzer.cc.

1145  {
1146 
1147  bool return_value=true;
1148 
1149  /*
1150  W-electron event selection comes here
1151  */
1152 
1153  return return_value;
1154 
1155 }
bool PFMETAnalyzer::selectWMuonEvent ( const edm::Event &  iEvent)

Definition at line 1158 of file PFMETAnalyzer.cc.

1158  {
1159 
1160  bool return_value=true;
1161 
1162  /*
1163  W-muon event selection comes here
1164  */
1165 
1166  return return_value;
1167 
1168 }
void PFMETAnalyzer::setSource ( std::string  source)
inline

Definition at line 92 of file PFMETAnalyzer.h.

References _source, and LaserTracksInput_cfi::source.

92  {
93  _source = source;
94  }
PFMETAnalyzer::_source
std::string _source
Definition: PFMETAnalyzer.h:105
LaserTracksInput_cfi.source
tuple source
Definition: LaserTracksInput_cfi.py:5
void PFMETAnalyzer::validateMET ( const reco::PFMET &  pfmet,
edm::Handle< edm::View< reco::PFCandidate > >  pfCandidates 
)

HF calibration factor (in 31X applied by PFProducer)

Definition at line 857 of file PFMETAnalyzer.cc.

References funct::cos(), gather_cfg::cout, reco::LeafCandidate::energy(), i, reco::tau::pfCandidates(), phi, reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), funct::sin(), mathSSE::sqrt(), reco::MET::sumEt(), reco::LeafCandidate::theta(), and theta().

859 {
860  double sumEx = 0;
861  double sumEy = 0;
862  double sumEt = 0;
863 
864  for( unsigned i=0; i<pfCandidates->size(); i++ ) {
865 
866  const reco::PFCandidate& cand = (*pfCandidates)[i];
867 
868  double E = cand.energy();
869 
871  // if( cand.particleId()==PFCandidate::h_HF ||
872  // cand.particleId()==PFCandidate::egamma_HF )
873  // E *= hfCalibFactor_;
874 
875  double phi = cand.phi();
876  double cosphi = cos(phi);
877  double sinphi = sin(phi);
878 
879  double theta = cand.theta();
880  double sintheta = sin(theta);
881 
882  double et = E*sintheta;
883  double ex = et*cosphi;
884  double ey = et*sinphi;
885 
886  sumEx += ex;
887  sumEy += ey;
888  sumEt += et;
889  }
890 
891  double Et = sqrt( sumEx*sumEx + sumEy*sumEy);
892  XYZTLorentzVector missingEt( -sumEx, -sumEy, 0, Et);
893 
894  if(_verbose)
895  if (sumEt!=pfmet.sumEt() || sumEx!=pfmet.px() || sumEy!=pfmet.py() || missingEt.T()!=pfmet.pt() )
896  {
897  std::cout<<"PFSumEt: " << sumEt <<", "<<"PFMETBlock: "<<pfmet.pt()<<std::endl;
898  std::cout<<"PFMET: " << missingEt.T() <<", "<<"PFMETBlock: "<<pfmet.pt()<<std::endl;
899  std::cout<<"PFMETx: " << missingEt.X() <<", "<<"PFMETBlockx: "<<pfmet.pt()<<std::endl;
900  std::cout<<"PFMETy: " << missingEt.Y() <<", "<<"PFMETBlocky: "<<pfmet.pt()<<std::endl;
901  }
902 }
reco::LeafCandidate::energy
virtual double energy() const GCC11_FINAL
energy
Definition: LeafCandidate.h:134
i
int i
Definition: DBlmapReader.cc:9
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:172
reco::LeafCandidate::theta
virtual double theta() const GCC11_FINAL
momentum polar angle
Definition: LeafCandidate.h:157
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:149
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:30
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:155
reco::MET::sumEt
double sumEt() const
Definition: MET.h:48
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
reco::LeafCandidate::px
virtual double px() const GCC11_FINAL
x coordinate of momentum vector
Definition: LeafCandidate.h:147
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:35
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:153
phi
Definition: DDAxes.h:10

Member Data Documentation

bool PFMETAnalyzer::_allhist
private

Definition at line 185 of file PFMETAnalyzer.h.

bool PFMETAnalyzer::_allSelection
private

Definition at line 186 of file PFMETAnalyzer.h.

bool PFMETAnalyzer::_cleanupSelection
private

Definition at line 187 of file PFMETAnalyzer.h.

DQMStore* PFMETAnalyzer::_dbe
private

Definition at line 191 of file PFMETAnalyzer.h.

bool PFMETAnalyzer::_doHLTPhysicsOn
private

Definition at line 144 of file PFMETAnalyzer.h.

bool PFMETAnalyzer::_doPVCheck
private

Definition at line 143 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_EleEventFlag
private

Definition at line 125 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_etThreshold
private

Definition at line 173 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::_FolderNames
private

Definition at line 189 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_HighMETEventFlag
private

Definition at line 123 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_highPFMETThreshold
private

Definition at line 168 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_HighPtJetEventFlag
private

Definition at line 120 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_highPtPFJetThreshold
private

Definition at line 166 of file PFMETAnalyzer.h.

std::string PFMETAnalyzer::_hlt_PhysDec
private

Definition at line 137 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_lowPFMETThreshold
private

Definition at line 169 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_LowPtJetEventFlag
private

Definition at line 121 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_lowPtPFJetThreshold
private

Definition at line 167 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_MinBiasEventFlag
private

Definition at line 122 of file PFMETAnalyzer.h.

GenericTriggerEventFlag* PFMETAnalyzer::_MuonEventFlag
private

Definition at line 126 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_numPV
private

Definition at line 171 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_nvtx_min
private

Definition at line 149 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_nvtxtrks_min
private

Definition at line 150 of file PFMETAnalyzer.h.

std::string PFMETAnalyzer::_source
private

Definition at line 105 of file PFMETAnalyzer.h.

Referenced by setSource().

std::vector<unsigned > PFMETAnalyzer::_techTrigsAND
private

Definition at line 139 of file PFMETAnalyzer.h.

std::vector<unsigned > PFMETAnalyzer::_techTrigsNOT
private

Definition at line 141 of file PFMETAnalyzer.h.

std::vector<unsigned > PFMETAnalyzer::_techTrigsOR
private

Definition at line 140 of file PFMETAnalyzer.h.

bool PFMETAnalyzer::_tightBHFiltering
private

Definition at line 146 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_tightJetIDFiltering
private

Definition at line 147 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_Ele
private

Definition at line 161 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_HighMET
private

Definition at line 159 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_HighPtJet
private

Definition at line 156 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_JetMB
private

Definition at line 155 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_LowPtJet
private

Definition at line 157 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_MinBias
private

Definition at line 158 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_Muon
private

Definition at line 162 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_trig_PhysDec
private

Definition at line 163 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_verbose
private

Definition at line 102 of file PFMETAnalyzer.h.

Referenced by HiHelperTools.MassSearchReplaceAnyInputTagVisitor::doIt(), helpers.MassSearchReplaceAnyInputTagVisitor::doIt(), HiHelperTools.MassSearchReplaceParamVisitor::enter(), helpers.MassSearchReplaceParamVisitor::enter(), cmsPerfSuite.PerfSuite::mkCandleDir(), cmsPerfSuite.PerfSuite::optionParse(), cmsPerfSuite.PerfSuite::printFlush(), cmsPerfSuite.PerfSuite::runcmd(), and cmsPerfSuite.PerfSuite::runCmdSet().

double PFMETAnalyzer::_vtxchi2_max
private

Definition at line 152 of file PFMETAnalyzer.h.

int PFMETAnalyzer::_vtxndof_min
private

Definition at line 151 of file PFMETAnalyzer.h.

double PFMETAnalyzer::_vtxz_max
private

Definition at line 153 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::BeamHaloSummaryTag
private

Definition at line 112 of file PFMETAnalyzer.h.

JetMETDQMDCSFilter* PFMETAnalyzer::DCSFilter
private

Definition at line 182 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::elecExpr_
private

Definition at line 133 of file PFMETAnalyzer.h.

int PFMETAnalyzer::evtCounter

Definition at line 96 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::gtTag
private

Definition at line 115 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::HBHENoiseFilterResultTag
private

Definition at line 113 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::HcalNoiseRBXCollectionTag
private

Definition at line 108 of file PFMETAnalyzer.h.

double PFMETAnalyzer::hfCalibFactor_
private

Definition at line 176 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::highMETExpr_
private

Definition at line 130 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::highPtJetExpr_
private

Definition at line 128 of file PFMETAnalyzer.h.

std::vector<std::string > PFMETAnalyzer::HLTPathsJetMBByName_
private

Definition at line 118 of file PFMETAnalyzer.h.

reco::helper::JetIDHelper* PFMETAnalyzer::jetID
private

Definition at line 179 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::lowPtJetExpr_
private

Definition at line 129 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meChargedHadronEt
private

Definition at line 225 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meChargedHadronEt_profile
private

Definition at line 248 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meChargedHadronEtFraction
private

Definition at line 224 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meChargedHadronEtFraction_profile
private

Definition at line 247 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meElectronEt
private

Definition at line 223 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meElectronEt_profile
private

Definition at line 246 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meElectronEtFraction
private

Definition at line 222 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meElectronEtFraction_profile
private

Definition at line 245 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFEMEt
private

Definition at line 231 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFEMEt_profile
private

Definition at line 254 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFEMEtFraction
private

Definition at line 230 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFEMEtFraction_profile
private

Definition at line 253 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFHadronEt
private

Definition at line 229 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFHadronEt_profile
private

Definition at line 252 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFHadronEtFraction
private

Definition at line 228 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meHFHadronEtFraction_profile
private

Definition at line 251 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meMuonEt
private

Definition at line 227 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meMuonEt_profile
private

Definition at line 250 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meMuonEtFraction
private

Definition at line 226 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meMuonEtFraction_profile
private

Definition at line 249 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meNeutralHadronEt
private

Definition at line 221 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meNeutralHadronEt_profile
private

Definition at line 244 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meNeutralHadronEtFraction
private

Definition at line 220 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meNeutralHadronEtFraction_profile
private

Definition at line 243 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMET
private

Definition at line 210 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMET_logx
private

Definition at line 215 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMET_profile
private

Definition at line 238 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMETPhi
private

Definition at line 211 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMETRate
private

Definition at line 197 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMETSig
private

Definition at line 209 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMEx
private

Definition at line 207 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMEx_profile
private

Definition at line 236 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMExLS
private

Definition at line 213 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMEy
private

Definition at line 208 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMEy_profile
private

Definition at line 237 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfMEyLS
private

Definition at line 214 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfSumET
private

Definition at line 212 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfSumET_logx
private

Definition at line 216 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePfSumET_profile
private

Definition at line 239 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePhotonEt
private

Definition at line 219 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePhotonEt_profile
private

Definition at line 242 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePhotonEtFraction
private

Definition at line 218 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::mePhotonEtFraction_profile
private

Definition at line 241 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::metME
private

Definition at line 196 of file PFMETAnalyzer.h.

std::string PFMETAnalyzer::metname
private

Definition at line 104 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_Ele
private

Definition at line 203 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_HighMET
private

Definition at line 202 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_HighPtJet
private

Definition at line 199 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_LowPtJet
private

Definition at line 200 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_MinBias
private

Definition at line 201 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_Muon
private

Definition at line 204 of file PFMETAnalyzer.h.

MonitorElement* PFMETAnalyzer::meTriggerName_PhysDec
private

Definition at line 205 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::minbiasExpr_
private

Definition at line 134 of file PFMETAnalyzer.h.

std::vector<std::string> PFMETAnalyzer::muonExpr_
private

Definition at line 132 of file PFMETAnalyzer.h.

edm::ParameterSet PFMETAnalyzer::parameters
private

Definition at line 101 of file PFMETAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::InputTag PFMETAnalyzer::PFCandidatesTag
private

Definition at line 111 of file PFMETAnalyzer.h.

edm::ParameterSet PFMETAnalyzer::theCleaningParameters
private

Definition at line 136 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::theJetCollectionLabel
private

Definition at line 109 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::thePfJetCollectionLabel
private

Definition at line 110 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::thePfMETCollectionLabel
private

Definition at line 107 of file PFMETAnalyzer.h.

edm::InputTag PFMETAnalyzer::vertexTag
private

Definition at line 114 of file PFMETAnalyzer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5