CaloMETAnalyzer Class Reference

#include <CaloMETAnalyzer.h>

Inheritance diagram for CaloMETAnalyzer:

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	bookMonitorElementTriggered (std::string, bool)
	CaloMETAnalyzer (const edm::ParameterSet &, edm::ConsumesCollector &&)
	Constructor. More...
void	computeEmHaMET (edm::Handle< edm::View< reco::Candidate > >)
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::CaloMET &, const reco::CaloMET &)
void	fillMESubSet (const edm::Event &, std::string, const reco::CaloMET &, const reco::CaloMET &)
void	fillMonitorElement (const edm::Event &, std::string, std::string, const reco::CaloMET &, const reco::CaloMET &, bool)
void	fillMonitorElementTriggered (const edm::Event &, std::string, std::string, const reco::CaloMET &, const reco::CaloMET &, bool)
void	makeRatePlot (std::string, double)
bool	selectHighPtJetEvent (const edm::Event &)
bool	selectLowPtJetEvent (const edm::Event &)
bool	selectWElectronEvent (const edm::Event &)
bool	selectWMuonEvent (const edm::Event &)
void	validateMET (const reco::CaloMET &, edm::Handle< edm::View< reco::Candidate > >)
virtual	~CaloMETAnalyzer ()
	Destructor. More...
Public Member Functions inherited from CaloMETAnalyzerBase
void	analyze (const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &, const reco::CaloMET &caloMET)
	Get the analysis of the muon properties. More...
	CaloMETAnalyzerBase ()
	Constructor. More...
virtual	~CaloMETAnalyzerBase ()
	Destructor. More...

Public Attributes
int	evtCounter

Private Attributes
bool	_allhist
bool	_allSelection
bool	_cleanupSelection
DQMStore *	_dbe
bool	_doHLTPhysicsOn
bool	_doPVCheck
GenericTriggerEventFlag *	_EleEventFlag
double	_EmCaloEz
double	_EmMET
double	_EmMetPhi
double	_EmMEx
double	_EmMEy
double	_EmSumEt
double	_etThreshold
std::vector< std::string >	_FolderNames
double	_HaCaloEz
double	_HaMET
double	_HaMetPhi
double	_HaMEx
double	_HaMEy
double	_HaSumEt
GenericTriggerEventFlag *	_HighMETEventFlag
double	_highMETThreshold
GenericTriggerEventFlag *	_HighPtJetEventFlag
double	_highPtJetThreshold
std::string	_hlt_PhysDec
GenericTriggerEventFlag *	_LowPtJetEventFlag
double	_lowPtJetThreshold
GenericTriggerEventFlag *	_MinBiasEventFlag
GenericTriggerEventFlag *	_MuonEventFlag
int	_numPV
int	_nvtx_min
int	_nvtxtrks_min
int	_print
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
edm::Handle< reco::BeamSpot >	beamSpot_h
math::XYZPoint	bspot
edm::Handle< edm::ValueMap < reco::MuonMETCorrectionData > >	corMetGlobalMuons_ValueMap_Handle
JetMETDQMDCSFilter *	DCSFilter
std::vector< std::string >	elecExpr_
edm::InputTag	gtTag
edm::InputTag	HBHENoiseFilterResultTag
edm::InputTag	HcalNoiseRBXCollectionTag
MonitorElement *	hCaloEmEtFraction
MonitorElement *	hCaloEmEtFraction020
MonitorElement *	hCaloEmEtInEB
MonitorElement *	hCaloEmEtInEE
MonitorElement *	hCaloEmEtInHF
MonitorElement *	hCaloEmMET
MonitorElement *	hCaloEmMETPhi
MonitorElement *	hCaloEmMEx
MonitorElement *	hCaloEmMEy
MonitorElement *	hCaloEtFractionHadronic
MonitorElement *	hCaloHadEtInHB
MonitorElement *	hCaloHadEtInHE
MonitorElement *	hCaloHadEtInHF
MonitorElement *	hCaloHadEtInHO
MonitorElement *	hCaloHaMET
MonitorElement *	hCaloHaMETPhi
MonitorElement *	hCaloHaMEx
MonitorElement *	hCaloHaMEy
MonitorElement *	hCaloMaxEtInEmTowers
MonitorElement *	hCaloMaxEtInHadTowers
MonitorElement *	hCaloMET
MonitorElement *	hCaloMET1
MonitorElement *	hCaloMET_logx
MonitorElement *	hCaloMET_profile
MonitorElement *	hCaloMETNoHF
MonitorElement *	hCaloMETNoHF_profile
MonitorElement *	hCaloMETPhi
MonitorElement *	hCaloMETPhi020
MonitorElement *	hCaloMETRate
MonitorElement *	hCaloMETSig
MonitorElement *	hCaloMEx
MonitorElement *	hCaloMEx_profile
MonitorElement *	hCaloMExCorrection
MonitorElement *	hCaloMExLS
MonitorElement *	hCaloMEy
MonitorElement *	hCaloMEy_profile
MonitorElement *	hCaloMEyCorrection
MonitorElement *	hCaloMEyLS
MonitorElement *	hCalomuChi2
MonitorElement *	hCalomuD0
MonitorElement *	hCalomuEta
MonitorElement *	hCalomuNhits
MonitorElement *	hCaloMuonCorrectionFlag
MonitorElement *	hCalomuPt
MonitorElement *	hCaloSumET
MonitorElement *	hCaloSumET_logx
MonitorElement *	hCaloSumET_profile
std::vector< std::string >	highMETExpr_
std::vector< std::string >	highPtJetExpr_
std::vector< std::string >	HLTPathsJetMBByName_
MonitorElement *	hmetME
MonitorElement *	hTriggerName_Ele
MonitorElement *	hTriggerName_HighMET
MonitorElement *	hTriggerName_HighPtJet
MonitorElement *	hTriggerName_LowPtJet
MonitorElement *	hTriggerName_MinBias
MonitorElement *	hTriggerName_Muon
edm::InputTag	inputBeamSpotLabel
edm::InputTag	inputMuonLabel
reco::helper::JetIDHelper *	jetID
std::vector< std::string >	lowPtJetExpr_
std::string	metname
std::vector< std::string >	minbiasExpr_
edm::Handle< reco::MuonCollection >	muon_h
std::vector< std::string >	muonExpr_
edm::ParameterSet	parameters
edm::InputTag	theCaloMETCollectionLabel
edm::InputTag	theCaloTowersLabel
edm::ParameterSet	theCleaningParameters
edm::InputTag	theJetCollectionLabel
edm::InputTag	vertexTag

Detailed Description

DQM monitoring source for CaloMET

Author

F. Chlebana - Fermilab K. Hatakeyama - Rockefeller University

Definition at line 52 of file CaloMETAnalyzer.h.

Constructor & Destructor Documentation

CaloMETAnalyzer::CaloMETAnalyzer	(	const edm::ParameterSet &	pSet,
		edm::ConsumesCollector &&	iC
	)

Constructor.

Definition at line 31 of file CaloMETAnalyzer.cc.

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

                                                                                    {

  parameters = pSet;

  edm::ParameterSet highptjetparms = parameters.getParameter<edm::ParameterSet>("highPtJetTrigger");
  edm::ParameterSet lowptjetparms  = parameters.getParameter<edm::ParameterSet>("lowPtJetTrigger" );
  edm::ParameterSet minbiasparms   = parameters.getParameter<edm::ParameterSet>("minBiasTrigger"  );
  edm::ParameterSet highmetparms   = parameters.getParameter<edm::ParameterSet>("highMETTrigger"  );
  //  edm::ParameterSet lowmetparms    = parameters.getParameter<edm::ParameterSet>("lowMETTrigger"   );
  edm::ParameterSet eleparms       = parameters.getParameter<edm::ParameterSet>("eleTrigger"      );
  edm::ParameterSet muonparms      = parameters.getParameter<edm::ParameterSet>("muonTrigger"     );

  //genericTriggerEventFlag_( new GenericTriggerEventFlag( conf_, iC ) );
  _HighPtJetEventFlag = new GenericTriggerEventFlag( highptjetparms, iC );
  _LowPtJetEventFlag  = new GenericTriggerEventFlag( lowptjetparms, iC );
  _MinBiasEventFlag   = new GenericTriggerEventFlag( minbiasparms,  iC );
  _HighMETEventFlag   = new GenericTriggerEventFlag( highmetparms,  iC );
  //  _LowMETEventFlag    = new GenericTriggerEventFlag( lowmetparms,   iC );
  _EleEventFlag       = new GenericTriggerEventFlag( eleparms,      iC );
  _MuonEventFlag      = new GenericTriggerEventFlag( muonparms,     iC );

  highPtJetExpr_ = highptjetparms.getParameter<std::vector<std::string> >("hltPaths");
  lowPtJetExpr_  = lowptjetparms .getParameter<std::vector<std::string> >("hltPaths");
  highMETExpr_   = highmetparms  .getParameter<std::vector<std::string> >("hltPaths");
  //  lowMETExpr_    = lowmetparms   .getParameter<std::vector<std::string> >("hltPaths");
  muonExpr_      = muonparms     .getParameter<std::vector<std::string> >("hltPaths");
  elecExpr_      = eleparms      .getParameter<std::vector<std::string> >("hltPaths");
  minbiasExpr_   = minbiasparms  .getParameter<std::vector<std::string> >("hltPaths");

}

CaloMETAnalyzer::~CaloMETAnalyzer ( )

virtual

Destructor.

Definition at line 63 of file CaloMETAnalyzer.cc.

                                  {

  delete _HighPtJetEventFlag;
  delete _LowPtJetEventFlag;
  delete _MinBiasEventFlag;
  delete _HighMETEventFlag;
  //  delete _LowMETEventFlag;
  delete _EleEventFlag;
  delete _MuonEventFlag;

}

Member Function Documentation

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

Get the analysis.

Definition at line 506 of file CaloMETAnalyzer.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::EventID::event(), edm::Event::getByLabel(), reco::BeamHaloSummary::GlobalLooseHaloId(), reco::BeamHaloSummary::GlobalTightHaloId(), reco::BeamHaloSummary::HcalLooseHaloId(), reco::BeamHaloSummary::HcalTightHaloId(), i, edm::EventBase::id(), edm::HandleBase::isValid(), bTagSequences_cff::jetID, LogDebug, LogTrace, metname, 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.

                                                           {

  if (_verbose) std::cout << "CaloMETAnalyzer analyze" << std::endl;

  std::string DirName = "JetMET/MET/"+_source;
  if (_print){
  std::cout << " " << std::endl;
  std::cout << "Event = " << iEvent.id().event() << std::endl;
  }

  LogTrace(metname)<<"[CaloMETAnalyzer] Analyze CaloMET";

  hmetME->Fill(1);

  // ==========================================================
  // Trigger information
  //
  _trig_JetMB=0;
  _trig_HighPtJet=0;
  _trig_LowPtJet=0;
  _trig_MinBias=0;
  _trig_HighMET=0;
  //  _trig_LowMET=0;
  _trig_Ele=0;
  _trig_Muon=0;
  _trig_PhysDec=0;
  if(&triggerResults) {


    //
    //
    // Check how many HLT triggers are in triggerResults
    int ntrigs = triggerResults.size();
    if (_verbose) std::cout << "ntrigs=" << ntrigs << std::endl;

    //
    //
    // If index=ntrigs, this HLT trigger doesn't exist in the HLT table for this data.
    const edm::TriggerNames & triggerNames = iEvent.triggerNames(triggerResults);

    //
    //
    const unsigned int nTrig(triggerNames.size());
    for (unsigned int i=0;i<nTrig;++i)
      {
        if (triggerNames.triggerName(i).find(highPtJetExpr_[0].substr(0,highPtJetExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_HighPtJet=true;
        else if (triggerNames.triggerName(i).find(lowPtJetExpr_[0].substr(0,lowPtJetExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_LowPtJet=true;
        else if (triggerNames.triggerName(i).find(highMETExpr_[0].substr(0,highMETExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_HighMET=true;
    //        else if (triggerNames.triggerName(i).find(lowMETExpr_[0].substr(0,lowMETExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
    //    _trig_LowMET=true;
        else if (triggerNames.triggerName(i).find(muonExpr_[0].substr(0,muonExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_Muon=true;
        else if (triggerNames.triggerName(i).find(elecExpr_[0].substr(0,elecExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_Ele=true;
        else if (triggerNames.triggerName(i).find(minbiasExpr_[0].substr(0,minbiasExpr_[0].rfind("_v")+2))!=std::string::npos && triggerResults.accept(i))
      _trig_MinBias=true;
      }

    // count number of requested Jet or MB HLT paths which have fired
    for (unsigned int i=0; i!=HLTPathsJetMBByName_.size(); i++) {
      unsigned int triggerIndex = triggerNames.triggerIndex(HLTPathsJetMBByName_[i]);
      if (triggerIndex<triggerResults.size()) {
    if (triggerResults.accept(triggerIndex)) {
      _trig_JetMB++;
    }
      }
    }
    // for empty input vectors (n==0), take all HLT triggers!
    if (HLTPathsJetMBByName_.size()==0) _trig_JetMB=triggerResults.size()-1;

    /*
      if ( _HighPtJetEventFlag->on() && _HighPtJetEventFlag->accept( iEvent, iSetup) )
      _trig_HighPtJet=1;

      if ( _LowPtJetEventFlag->on() && _LowPtJetEventFlag->accept( iEvent, iSetup) )
      _trig_LowPtJet=1;

      if ( _MinBiasEventFlag->on() && _MinBiasEventFlag->accept( iEvent, iSetup) )
      _trig_MinBias=1;

      if ( _HighMETEventFlag->on() && _HighMETEventFlag->accept( iEvent, iSetup) )
      _trig_HighMET=1;

      if ( _LowMETEventFlag->on() && _LowMETEventFlag->accept( iEvent, iSetup) )
      _trig_LowMET=1;

      if ( _EleEventFlag->on() && _EleEventFlag->accept( iEvent, iSetup) )
      _trig_Ele=1;

      if ( _MuonEventFlag->on() && _MuonEventFlag->accept( iEvent, iSetup) )
      _trig_Muon=1;
    */

    if (triggerNames.triggerIndex(_hlt_PhysDec)   != triggerNames.size() &&
    triggerResults.accept(triggerNames.triggerIndex(_hlt_PhysDec)))   _trig_PhysDec=1;
  } else {

    edm::LogInfo("CaloMetAnalyzer") << "TriggerResults::HLT not found, "
    "automatically select events";
    //
    // TriggerResults object not found. Look at all events.
    _trig_JetMB=1;

  }

  // ==========================================================
  // CaloMET information

  // **** Get the MET container
  edm::Handle<reco::CaloMETCollection> calometcoll;
  iEvent.getByLabel(theCaloMETCollectionLabel, calometcoll);

  if(!calometcoll.isValid()) {
    std::cout<<"Unable to find MET results for CaloMET collection "<<theCaloMETCollectionLabel<<std::endl;
    return;
  }

  const reco::CaloMETCollection *calometcol = calometcoll.product();
  const reco::CaloMET *calomet;
  calomet = &(calometcol->front());

  LogTrace(metname)<<"[CaloMETAnalyzer] Call to the CaloMET analyzer";

  // **** Get the MET no HF container
  edm::Handle<reco::CaloMETCollection> calometnohfcoll;
  iEvent.getByLabel("metNoHF", calometnohfcoll);

  if(!calometnohfcoll.isValid()) {
    std::cout<<"Unable to find MET results for CaloMETNoHF collection metNoHF"<<std::endl;
    return;
  }

  const reco::CaloMETCollection *calometnohfcol = calometnohfcoll.product();
  const reco::CaloMET *calometnohf;
  calometnohf = &(calometnohfcol->front());

  //LogTrace(metname)<<"[CaloMETAnalyzer] Call to the CaloMET analyzer";

  //Only for corMetGlobalMuons
  if (theCaloMETCollectionLabel.label() == "corMetGlobalMuons" ) {

    iEvent.getByLabel("muonMETValueMapProducer" , "muCorrData", corMetGlobalMuons_ValueMap_Handle);
    iEvent.getByLabel("muons", muon_h);
    iEvent.getByLabel(inputBeamSpotLabel, beamSpot_h);

    if(!beamSpot_h.isValid()) edm::LogInfo("OutputInfo") << "falied to retrieve beam spot data require by MET Task";

    bspot = ( beamSpot_h.isValid() ) ? beamSpot_h->position() : math::XYZPoint(0, 0, 0);

  }


  // ==========================================================
  //
  edm::Handle<reco::HcalNoiseRBXCollection> HRBXCollection;
  iEvent.getByLabel(HcalNoiseRBXCollectionTag,HRBXCollection);
  if (!HRBXCollection.isValid()) {
      LogDebug("") << "CaloMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
      if (_verbose) std::cout << "CaloMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
  }


  edm::Handle<bool> HBHENoiseFilterResultHandle;
  iEvent.getByLabel(HBHENoiseFilterResultTag, HBHENoiseFilterResultHandle);
  bool HBHENoiseFilterResult = *HBHENoiseFilterResultHandle;
  if (!HBHENoiseFilterResultHandle.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
  }


  edm::Handle<reco::CaloJetCollection> caloJets;
  iEvent.getByLabel(theJetCollectionLabel, caloJets);
  if (!caloJets.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find jet product" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find jet product" << std::endl;
  }

  edm::Handle<edm::View<reco::Candidate> > towers;
  iEvent.getByLabel(theCaloTowersLabel, towers);
  if (!towers.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find caltower product" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find caltower product" << std::endl;
  }

  // ==========================================================
  // CaloMET sanity check

  if (_source=="CaloMET") validateMET(*calomet,towers);

  // ==========================================================

  if (_allhist) computeEmHaMET(towers);

  // ==========================================================
  // JetID

  if (_verbose) std::cout << "JetID starts" << std::endl;

  //
  // --- Minimal cuts
  //
  bool bJetIDMinimal=true;
  int nj=0;
  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
       cal!=caloJets->end(); ++cal){
    jetID->calculate(iEvent, *cal);
    if (_print && nj<=1) std::cout << "Jet pT = " << cal->pt() << " (GeV) "
                   << " eta = " << cal->eta() << " "
                   << " phi = " << cal->phi() << " "
                   << " emf = " << cal->emEnergyFraction() << std::endl;
    nj++;
    if (cal->pt()>10.){
      if (fabs(cal->eta())<=2.6 &&
      cal->emEnergyFraction()<=0.01) bJetIDMinimal=false;
    }
  }

  //
  // --- Loose cuts
  //
  bool bJetIDLoose=true;
  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
       cal!=caloJets->end(); ++cal){
    jetID->calculate(iEvent, *cal);
    if (_verbose) std::cout << jetID->n90Hits() << " "
                << jetID->restrictedEMF() << " "
                << cal->pt() << std::endl;
    if (cal->pt()>10.){
      //
      // for all regions
      if (jetID->n90Hits()<2)  bJetIDLoose=false;
      if (jetID->fHPD()>=0.98) bJetIDLoose=false;
      //
      // for non-forward
      if (fabs(cal->eta())<2.55){
    if (cal->emEnergyFraction()<=0.01) bJetIDLoose=false;
      }
      // for forward
      else {
    if (cal->emEnergyFraction()<=-0.9) bJetIDLoose=false;
    if (cal->pt()>80.){
    if (cal->emEnergyFraction()>= 1.0) bJetIDLoose=false;
    }
      } // forward vs non-forward
    }   // pt>10 GeV/c
  }     // calor-jets loop

  //
  // --- Tight cuts
  //
  bool bJetIDTight=true;
  bJetIDTight=bJetIDLoose;
  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
       cal!=caloJets->end(); ++cal){
    jetID->calculate(iEvent, *cal);
    if (cal->pt()>25.){
      //
      // for all regions
      if (jetID->fHPD()>=0.95) bJetIDTight=false;
      //
      // for 1.0<|eta|<1.75
      if (fabs(cal->eta())>=1.00 && fabs(cal->eta())<1.75){
    if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
      }
      //
      // for 1.75<|eta|<2.55
      else if (fabs(cal->eta())>=1.75 && fabs(cal->eta())<2.55){
    if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
      }
      //
      // for 2.55<|eta|<3.25
      else if (fabs(cal->eta())>=2.55 && fabs(cal->eta())<3.25){
    if (cal->pt()< 50.                   && cal->emEnergyFraction()<=-0.3) bJetIDTight=false;
    if (cal->pt()>=50. && cal->pt()< 80. && cal->emEnergyFraction()<=-0.2) bJetIDTight=false;
    if (cal->pt()>=80. && cal->pt()<340. && cal->emEnergyFraction()<=-0.1) bJetIDTight=false;
    if (cal->pt()>=340.                  && cal->emEnergyFraction()<=-0.1
                                             && cal->emEnergyFraction()>=0.95) bJetIDTight=false;
      }
      //
      // for 3.25<|eta|
      else if (fabs(cal->eta())>=3.25){
    if (cal->pt()< 50.                   && cal->emEnergyFraction()<=-0.3
                                             && cal->emEnergyFraction()>=0.90) bJetIDTight=false;
    if (cal->pt()>=50. && cal->pt()<130. && cal->emEnergyFraction()<=-0.2
                                             && cal->emEnergyFraction()>=0.80) bJetIDTight=false;
    if (cal->pt()>=130.                  && cal->emEnergyFraction()<=-0.1
                                             && cal->emEnergyFraction()>=0.70) bJetIDTight=false;
      }
    }   // pt>10 GeV/c
  }     // calor-jets loop

  if (_verbose) std::cout << "JetID ends" << std::endl;

  // ==========================================================
  // HCAL Noise filter

  bool bHcalNoiseFilter = HBHENoiseFilterResult;

  // ==========================================================
  // Get BeamHaloSummary
  edm::Handle<reco::BeamHaloSummary> TheBeamHaloSummary ;
  iEvent.getByLabel(BeamHaloSummaryTag, TheBeamHaloSummary) ;

  bool bBeamHaloIDTightPass = true;
  bool bBeamHaloIDLoosePass = true;

  if(TheBeamHaloSummary.isValid()) {

    const reco::BeamHaloSummary TheSummary = (*TheBeamHaloSummary.product() );

    //   std::cout << TheSummary.EcalLooseHaloId() << " "
    //      << TheSummary.HcalLooseHaloId() << " "
    //      << TheSummary.CSCLooseHaloId()  << " "
    //      << TheSummary.GlobalLooseHaloId() << std::endl;

    if( TheSummary.EcalLooseHaloId()  || TheSummary.HcalLooseHaloId() ||
    TheSummary.CSCLooseHaloId()   || TheSummary.GlobalLooseHaloId() )
      bBeamHaloIDLoosePass = false;

    if( TheSummary.EcalTightHaloId()  || TheSummary.HcalTightHaloId() ||
    TheSummary.CSCTightHaloId()   || TheSummary.GlobalTightHaloId() )
      bBeamHaloIDTightPass = false;

  }

  if (_verbose) std::cout << "BeamHaloSummary ends" << std::endl;

  // ==========================================================
  //Vertex information

  _numPV = 0;
  bool bPrimaryVertex = true;
  if(_doPVCheck){
    bPrimaryVertex = false;
    Handle<reco::VertexCollection> vertexHandle;

    iEvent.getByLabel(vertexTag, vertexHandle);

    if (!vertexHandle.isValid()) {
      LogDebug("") << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
      if (_verbose) std::cout << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
    }

    if ( vertexHandle.isValid() ){
      reco::VertexCollection vertexCollection = *(vertexHandle.product());
      int vertex_number     = vertexCollection.size();
      reco::VertexCollection::const_iterator v = vertexCollection.begin();
      for ( ; v != vertexCollection.end(); ++v) {
    double vertex_chi2    = v->normalizedChi2();
    double vertex_ndof    = v->ndof();
    bool   fakeVtx        = v->isFake();
    double vertex_Z       = v->z();

    if (  !fakeVtx
          && vertex_number>=_nvtx_min
          && vertex_ndof   >_vtxndof_min
          && vertex_chi2   <_vtxchi2_max
          && fabs(vertex_Z)<_vtxz_max ) {
      bPrimaryVertex = true;
      ++_numPV;
    }
      }
    }
  }
  // ==========================================================

  edm::Handle< L1GlobalTriggerReadoutRecord > gtReadoutRecord;
  iEvent.getByLabel( gtTag, gtReadoutRecord);

  if (!gtReadoutRecord.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find GT readout record" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
  }

  bool bTechTriggers    = true;
  bool bTechTriggersAND = true;
  bool bTechTriggersOR  = false;
  bool bTechTriggersNOT = false;

  if (gtReadoutRecord.isValid()) {
    const TechnicalTriggerWord&  technicalTriggerWordBeforeMask = gtReadoutRecord->technicalTriggerWord();

    if (_techTrigsAND.size() == 0)
      bTechTriggersAND = true;
    else
      for (unsigned ttr = 0; ttr != _techTrigsAND.size(); ttr++) {
    bTechTriggersAND = bTechTriggersAND && technicalTriggerWordBeforeMask.at(_techTrigsAND.at(ttr));
      }

    if (_techTrigsAND.size() == 0)
      bTechTriggersOR = true;
    else
      for (unsigned ttr = 0; ttr != _techTrigsOR.size(); ttr++) {
    bTechTriggersOR = bTechTriggersOR || technicalTriggerWordBeforeMask.at(_techTrigsOR.at(ttr));
      }
    if (_techTrigsNOT.size() == 0)
      bTechTriggersNOT = false;
    else
      for (unsigned ttr = 0; ttr != _techTrigsNOT.size(); ttr++) {
    bTechTriggersNOT = bTechTriggersNOT || technicalTriggerWordBeforeMask.at(_techTrigsNOT.at(ttr));
      }
  }
  else
    {
      bTechTriggersAND = true;
      bTechTriggersOR  = true;
      bTechTriggersNOT = false;
    }

  if (_techTrigsAND.size()==0)
    bTechTriggersAND = true;
  if (_techTrigsOR.size()==0)
    bTechTriggersOR = true;
  if (_techTrigsNOT.size()==0)
    bTechTriggersNOT = false;

  bTechTriggers = bTechTriggersAND && bTechTriggersOR && !bTechTriggersNOT;

  // ==========================================================
  // Reconstructed MET Information - fill MonitorElements

  bool bHcalNoise  = bHcalNoiseFilter;
  bool bBeamHaloID = bBeamHaloIDLoosePass;
  bool bJetID      = bJetIDMinimal;

  bool bPhysicsDeclared = true;
  if(_doHLTPhysicsOn) bPhysicsDeclared =_trig_PhysDec;

  if      (_tightBHFiltering)       bBeamHaloID = bBeamHaloIDTightPass;

  if      (_tightJetIDFiltering==1)  bJetID      = bJetIDMinimal;
  else if (_tightJetIDFiltering==2)  bJetID      = bJetIDLoose;
  else if (_tightJetIDFiltering==3)  bJetID      = bJetIDTight;
  else if (_tightJetIDFiltering==-1) bJetID      = true;

  bool bBasicCleanup = bTechTriggers && bPrimaryVertex && bPhysicsDeclared;
  bool bExtraCleanup = bBasicCleanup && bHcalNoise && bJetID && bBeamHaloID;

  //std::string DirName = "JetMET/MET/"+_source;

  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
       ic != _FolderNames.end(); ic++){
    if (*ic=="All")                                             fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (DCSFilter->filter(iEvent, iSetup)) {
      if (_cleanupSelection){
    if (*ic=="BasicCleanup"   && bBasicCleanup)                 fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="ExtraCleanup"   && bExtraCleanup)                 fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
      }
      if (_allSelection) {
    if (*ic=="HcalNoiseFilter"      && bHcalNoiseFilter )       fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="JetIDMinimal"         && bJetIDMinimal)           fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="JetIDLoose"           && bJetIDLoose)             fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="JetIDTight"           && bJetIDTight)             fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="BeamHaloIDTightPass"  && bBeamHaloIDTightPass)    fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="BeamHaloIDLoosePass"  && bBeamHaloIDLoosePass)    fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="Triggers"             && bTechTriggers)           fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
    if (*ic=="PV"                   && bPrimaryVertex)          fillMESet(iEvent, DirName+"/"+*ic, *calomet, *calometnohf);
      }
    } // DCS
  }
} // CaloMETAnalyzer::analyze

void CaloMETAnalyzer::beginJob ( DQMStore *  dbe )

virtual

Inizialize parameters for histo binning.

Implements CaloMETAnalyzerBase.

Definition at line 76 of file CaloMETAnalyzer.cc.

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

                                             {

  evtCounter = 0;
  metname = "caloMETAnalyzer";

  // trigger information
  HLTPathsJetMBByName_ = parameters.getParameter<std::vector<std::string > >("HLTPathsJetMB");

  theCleaningParameters = parameters.getParameter<ParameterSet>("CleaningParameters"),

  //Trigger parameters
  gtTag          = theCleaningParameters.getParameter<edm::InputTag>("gtLabel");
  _techTrigsAND  = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsAND");
  _techTrigsOR   = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsOR");
  _techTrigsNOT  = theCleaningParameters.getParameter<std::vector<unsigned > >("techTrigsNOT");

  _doHLTPhysicsOn = theCleaningParameters.getParameter<bool>("doHLTPhysicsOn");
  _hlt_PhysDec    = theCleaningParameters.getParameter<std::string>("HLT_PhysDec");

  _tightBHFiltering     = theCleaningParameters.getParameter<bool>("tightBHFiltering");
  _tightJetIDFiltering  = theCleaningParameters.getParameter<int>("tightJetIDFiltering");

  // ==========================================================
  //DCS information
  // ==========================================================
  DCSFilter = new JetMETDQMDCSFilter(parameters.getParameter<ParameterSet>("DCSFilter"));

  //Vertex requirements
  _doPVCheck          = theCleaningParameters.getParameter<bool>("doPrimaryVertexCheck");
  vertexTag  = theCleaningParameters.getParameter<edm::InputTag>("vertexLabel");

  if (_doPVCheck) {
    _nvtx_min        = theCleaningParameters.getParameter<int>("nvtx_min");
    _nvtxtrks_min    = theCleaningParameters.getParameter<int>("nvtxtrks_min");
    _vtxndof_min     = theCleaningParameters.getParameter<int>("vtxndof_min");
    _vtxchi2_max     = theCleaningParameters.getParameter<double>("vtxchi2_max");
    _vtxz_max        = theCleaningParameters.getParameter<double>("vtxz_max");
  }


  // CaloMET information
  theCaloMETCollectionLabel       = parameters.getParameter<edm::InputTag>("METCollectionLabel");
  _source                         = parameters.getParameter<std::string>("Source");

  if (theCaloMETCollectionLabel.label() == "corMetGlobalMuons" ) {
    inputBeamSpotLabel      = parameters.getParameter<edm::InputTag>("InputBeamSpotLabel");
  }

  // Other data collections
  theCaloTowersLabel          = parameters.getParameter<edm::InputTag>("CaloTowersLabel");
  theJetCollectionLabel       = parameters.getParameter<edm::InputTag>("JetCollectionLabel");
  HcalNoiseRBXCollectionTag   = parameters.getParameter<edm::InputTag>("HcalNoiseRBXCollection");
  BeamHaloSummaryTag          = parameters.getParameter<edm::InputTag>("BeamHaloSummaryLabel");
  HBHENoiseFilterResultTag    = parameters.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");

  // misc
  _verbose     = parameters.getParameter<int>("verbose");
  _print       = parameters.getParameter<int>("printOut");
  _etThreshold = parameters.getParameter<double>("etThreshold"); // MET threshold
  _allhist     = parameters.getParameter<bool>("allHist");       // Full set of monitoring histograms
  _allSelection= parameters.getParameter<bool>("allSelection");  // Plot with all sets of event selection
  _cleanupSelection= parameters.getParameter<bool>("cleanupSelection");  // Plot with all sets of event selection

  _highPtJetThreshold = parameters.getParameter<double>("HighPtJetThreshold"); // High Pt Jet threshold
  _lowPtJetThreshold = parameters.getParameter<double>("LowPtJetThreshold"); // Low Pt Jet threshold
  _highMETThreshold = parameters.getParameter<double>("HighMETThreshold"); // High MET threshold
  //  _lowMETThreshold = parameters.getParameter<double>("LowMETThreshold"); // Low MET threshold

  //
  jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams"));

  // DQStore stuff
  LogTrace(metname)<<"[CaloMETAnalyzer] Parameters initialization";
  std::string DirName = "JetMET/MET/"+_source;
  dbe->setCurrentFolder(DirName);

  hmetME = dbe->book1D("metReco", "metReco", 4, 1, 5);
  hmetME->setBinLabel(1,"CaloMET",1);

  _dbe = dbe;

  _FolderNames.push_back("All");
  _FolderNames.push_back("BasicCleanup");
  _FolderNames.push_back("ExtraCleanup");
  _FolderNames.push_back("HcalNoiseFilter");
  _FolderNames.push_back("JetIDMinimal");
  _FolderNames.push_back("JetIDLoose");
  _FolderNames.push_back("JetIDTight");
  _FolderNames.push_back("BeamHaloIDTightPass");
  _FolderNames.push_back("BeamHaloIDLoosePass");
  _FolderNames.push_back("Triggers");
  _FolderNames.push_back("PV");

  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
       ic != _FolderNames.end(); ic++){
    if (*ic=="All")             bookMESet(DirName+"/"+*ic);
    if (_cleanupSelection){
    if (*ic=="BasicCleanup")    bookMESet(DirName+"/"+*ic);
    if (*ic=="ExtraCleanup")    bookMESet(DirName+"/"+*ic);
    }
    if (_allSelection){
    if (*ic=="HcalNoiseFilter")      bookMESet(DirName+"/"+*ic);
    if (*ic=="JetIDMinimal")         bookMESet(DirName+"/"+*ic);
    if (*ic=="JetIDLoose")           bookMESet(DirName+"/"+*ic);
    if (*ic=="JetIDTight")           bookMESet(DirName+"/"+*ic);
    if (*ic=="BeamHaloIDTightPass")  bookMESet(DirName+"/"+*ic);
    if (*ic=="BeamHaloIDLoosePass")  bookMESet(DirName+"/"+*ic);
    if (*ic=="Triggers")             bookMESet(DirName+"/"+*ic);
    if (*ic=="PV")                   bookMESet(DirName+"/"+*ic);
    }
  }
}

void CaloMETAnalyzer::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

Initialize run-based parameters.

Definition at line 391 of file CaloMETAnalyzer.cc.

{
  if ( _HighPtJetEventFlag->on() ) _HighPtJetEventFlag->initRun( iRun, iSetup );
  if ( _LowPtJetEventFlag ->on() ) _LowPtJetEventFlag ->initRun( iRun, iSetup );
  if ( _MinBiasEventFlag  ->on() ) _MinBiasEventFlag  ->initRun( iRun, iSetup );
  if ( _HighMETEventFlag  ->on() ) _HighMETEventFlag  ->initRun( iRun, iSetup );
  //  if ( _LowMETEventFlag   ->on() ) _LowMETEventFlag   ->initRun( iRun, iSetup );
  if ( _EleEventFlag      ->on() ) _EleEventFlag      ->initRun( iRun, iSetup );
  if ( _MuonEventFlag     ->on() ) _MuonEventFlag     ->initRun( iRun, iSetup );

  if (_HighPtJetEventFlag->on() && _HighPtJetEventFlag->expressionsFromDB(_HighPtJetEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    highPtJetExpr_ = _HighPtJetEventFlag->expressionsFromDB(_HighPtJetEventFlag->hltDBKey(), iSetup);
  if (_LowPtJetEventFlag->on() && _LowPtJetEventFlag->expressionsFromDB(_LowPtJetEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    lowPtJetExpr_  = _LowPtJetEventFlag->expressionsFromDB(_LowPtJetEventFlag->hltDBKey(),   iSetup);
  if (_HighMETEventFlag->on() && _HighMETEventFlag->expressionsFromDB(_HighMETEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    highMETExpr_   = _HighMETEventFlag->expressionsFromDB(_HighMETEventFlag->hltDBKey(),     iSetup);
  //  if (_LowMETEventFlag->on() && _LowMETEventFlag->expressionsFromDB(_LowMETEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
  //    lowMETExpr_    = _LowMETEventFlag->expressionsFromDB(_LowMETEventFlag->hltDBKey(),       iSetup);
  if (_MuonEventFlag->on() && _MuonEventFlag->expressionsFromDB(_MuonEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    muonExpr_      = _MuonEventFlag->expressionsFromDB(_MuonEventFlag->hltDBKey(),           iSetup);
  if (_EleEventFlag->on() && _EleEventFlag->expressionsFromDB(_EleEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    elecExpr_      = _EleEventFlag->expressionsFromDB(_EleEventFlag->hltDBKey(),             iSetup);
  if (_MinBiasEventFlag->on() && _MinBiasEventFlag->expressionsFromDB(_MinBiasEventFlag->hltDBKey(), iSetup)[0] != "CONFIG_ERROR")
    minbiasExpr_   = _MinBiasEventFlag->expressionsFromDB(_MinBiasEventFlag->hltDBKey(),     iSetup);

}

void CaloMETAnalyzer::bookMESet

(

std::string

DirName

)

Definition at line 200 of file CaloMETAnalyzer.cc.

References gather_cfg::cout.

{
  bool bLumiSecPlot=false;
  if (DirName.find("All")!=std::string::npos) bLumiSecPlot=true;

  bookMonitorElement(DirName,bLumiSecPlot);

  if ( _HighPtJetEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"HighPtJet",false);
    hTriggerName_HighPtJet = _dbe->bookString("triggerName_HighPtJet", highPtJetExpr_[0]);
  }

  if ( _LowPtJetEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"LowPtJet",false);
    hTriggerName_LowPtJet = _dbe->bookString("triggerName_LowPtJet", lowPtJetExpr_[0]);
  }

  if ( _MinBiasEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"MinBias",false);
    hTriggerName_MinBias = _dbe->bookString("triggerName_MinBias", minbiasExpr_[0]);
    if (_verbose) std::cout << "_MinBiasEventFlag is on, folder created\n";
  }

  if ( _HighMETEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"HighMET",false);
    hTriggerName_HighMET = _dbe->bookString("triggerName_HighMET", highMETExpr_[0]);
  }

  //  if ( _LowMETEventFlag->on() ) {
  //    bookMonitorElementTriggered(DirName+"/"+"LowMET",false);
  //    hTriggerName_LowMET = _dbe->bookString("triggerName_LowMET", lowMETExpr_[0]);
  //  }

  if ( _EleEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"Ele",false);
    hTriggerName_Ele = _dbe->bookString("triggerName_Ele", elecExpr_[0]);
    if (_verbose) std::cout << "_EleEventFlag is on, folder created\n";
  }

  if ( _MuonEventFlag->on() ) {
    bookMonitorElementTriggered(DirName+"/"+"Muon",false);
    hTriggerName_Muon = _dbe->bookString("triggerName_Muon", muonExpr_[0]);
    if (_verbose) std::cout << "_MuonEventFlag is on, folder created\n";
  }
}

void CaloMETAnalyzer::bookMonitorElement	(	std::string	DirName,
		bool	bLumiSecPlot = false
	)

Definition at line 247 of file CaloMETAnalyzer.cc.

References gather_cfg::cout.

{
  if (_verbose) std::cout << "bookMonitorElement " << DirName << std::endl;

  _dbe->setCurrentFolder(DirName);

  bookMonitorElementTriggered(DirName, bLumiSecPlot);


  if (_allhist) {

    hCaloMaxEtInEmTowers    = _dbe->book1D("METTask_CaloMaxEtInEmTowers",   "METTask_CaloMaxEtInEmTowers"   ,100,0,2000);
    hCaloMaxEtInEmTowers->setAxisTitle("Et(Max) in EM Tower [GeV]",1);
    hCaloMaxEtInHadTowers   = _dbe->book1D("METTask_CaloMaxEtInHadTowers",  "METTask_CaloMaxEtInHadTowers"  ,100,0,2000);
    hCaloMaxEtInHadTowers->setAxisTitle("Et(Max) in Had Tower [GeV]",1);

    hCaloHadEtInHB          = _dbe->book1D("METTask_CaloHadEtInHB","METTask_CaloHadEtInHB",100,0,2000);
    hCaloHadEtInHB->setAxisTitle("Had Et [GeV]",1);
    hCaloHadEtInHO          = _dbe->book1D("METTask_CaloHadEtInHO","METTask_CaloHadEtInHO",25,0,500);
    hCaloHadEtInHO->setAxisTitle("Had Et [GeV]",1);
    hCaloHadEtInHE          = _dbe->book1D("METTask_CaloHadEtInHE","METTask_CaloHadEtInHE",100,0,2000);
    hCaloHadEtInHE->setAxisTitle("Had Et [GeV]",1);
    hCaloHadEtInHF          = _dbe->book1D("METTask_CaloHadEtInHF","METTask_CaloHadEtInHF",50,0,1000);
    hCaloHadEtInHF->setAxisTitle("Had Et [GeV]",1);
    hCaloEmEtInHF           = _dbe->book1D("METTask_CaloEmEtInHF" ,"METTask_CaloEmEtInHF" ,25,0,500);
    hCaloEmEtInHF->setAxisTitle("EM Et [GeV]",1);
    hCaloEmEtInEE           = _dbe->book1D("METTask_CaloEmEtInEE" ,"METTask_CaloEmEtInEE" ,50,0,1000);
    hCaloEmEtInEE->setAxisTitle("EM Et [GeV]",1);
    hCaloEmEtInEB           = _dbe->book1D("METTask_CaloEmEtInEB" ,"METTask_CaloEmEtInEB" ,100,0,2000);
    hCaloEmEtInEB->setAxisTitle("EM Et [GeV]",1);

    hCaloEmMEx= _dbe->book1D("METTask_CaloEmMEx","METTask_CaloEmMEx",200,-500,500);
    hCaloEmMEx->setAxisTitle("EM MEx [GeV]",1);
    hCaloEmMEy= _dbe->book1D("METTask_CaloEmMEy","METTask_CaloEmMEy",200,-500,500);
    hCaloEmMEy->setAxisTitle("EM MEy [GeV]",1);
    hCaloEmMET= _dbe->book1D("METTask_CaloEmMET","METTask_CaloEmMET",200,0,1000);
    hCaloEmMET->setAxisTitle("EM MET [GeV]",1);
    hCaloEmMETPhi= _dbe->book1D("METTask_CaloEmMETPhi","METTask_CaloEmMETPhi", 60, -3.2, 3.2);
    hCaloEmMETPhi->setAxisTitle("EM METPhi [rad]",1);

    hCaloHaMEx= _dbe->book1D("METTask_CaloHaMEx","METTask_CaloHaMEx",200,-500,500);
    hCaloHaMEx->setAxisTitle("HA MEx [GeV]",1);
    hCaloHaMEy= _dbe->book1D("METTask_CaloHaMEy","METTask_CaloHaMEy",200,-500,500);
    hCaloHaMEy->setAxisTitle("HA MEy [GeV]",1);
    hCaloHaMET= _dbe->book1D("METTask_CaloHaMET","METTask_CaloHaMET",200,0,1000);
    hCaloHaMET->setAxisTitle("HA MET [GeV]",1);
    hCaloHaMETPhi= _dbe->book1D("METTask_CaloHaMETPhi","METTask_CaloHaMETPhi", 60, -3.2, 3.2);
    hCaloHaMETPhi->setAxisTitle("HA METPhi [rad]",1);
  }
}

void CaloMETAnalyzer::bookMonitorElementTriggered	(	std::string	DirName,
		bool	bLumiSecPlot = false
	)

Definition at line 300 of file CaloMETAnalyzer.cc.

References gather_cfg::cout, nbinsPV, PVlow, and PVup.

{
  if (_verbose) std::cout << "bookMonitorElementTriggered " << DirName << std::endl;

  _dbe->setCurrentFolder(DirName);


  hCaloMEx        = _dbe->book1D("METTask_CaloMEx",        "METTask_CaloMEx",        200, -500,  500);
  hCaloMEy        = _dbe->book1D("METTask_CaloMEy",        "METTask_CaloMEy",        200, -500,  500);
  hCaloMET        = _dbe->book1D("METTask_CaloMET",        "METTask_CaloMET",        200,    0, 1000);
  hCaloMET1       = _dbe->book1D("METTask_CaloMET1",       "METTask_CaloMET1",        80,    0,  200);
  hCaloMETNoHF    = _dbe->book1D("METTask_CaloMETNoHF",    "METTask_CaloMETNoHF",    200,    0, 1000);
  hCaloSumET      = _dbe->book1D("METTask_CaloSumET",      "METTask_CaloSumET",      400,    0, 4000);
  hCaloMETSig     = _dbe->book1D("METTask_CaloMETSig",     "METTask_CaloMETSig",      51,    0,   51);
  hCaloMETPhi     = _dbe->book1D("METTask_CaloMETPhi",     "METTask_CaloMETPhi",      60, -3.2,  3.2);
  hCaloMETPhi020  = _dbe->book1D("METTask_CaloMETPhi020",  "METTask_CaloMETPhi020",   60, -3.2,  3.2);
  hCaloMET_logx   = _dbe->book1D("METTask_CaloMET_logx",   "METTask_CaloMET_logx",    40,   -1,    7);
  hCaloSumET_logx = _dbe->book1D("METTask_CaloSumET_logx", "METTask_CaloSumET_logx",  40,   -1,    7);

  hCaloMEx       ->setAxisTitle("MEx [GeV]",                 1);
  hCaloMEy       ->setAxisTitle("MEy [GeV]",                 1);
  hCaloMET       ->setAxisTitle("MET [GeV]",                 1);
  hCaloMET1      ->setAxisTitle("MET [GeV]",                 1);
  hCaloMETNoHF   ->setAxisTitle("MET (No HF) [GeV]",         1);
  hCaloSumET     ->setAxisTitle("SumET [GeV]",               1);
  hCaloMETSig    ->setAxisTitle("METSig",                    1);
  hCaloMETPhi    ->setAxisTitle("METPhi [rad]",              1);
  hCaloMETPhi020 ->setAxisTitle("METPhi [rad] (MET>20 GeV)", 1);
  hCaloMET_logx  ->setAxisTitle("log(MET) [GeV]",            1);
  hCaloSumET_logx->setAxisTitle("log(SumET) [GeV]",          1);


  // Book NPV profiles
  //----------------------------------------------------------------------------
  hCaloMEx_profile     = _dbe->bookProfile("METTask_CaloMEx_profile",     "MEx [GeV]",     nbinsPV, PVlow, PVup, 200, -500,  500);
  hCaloMEy_profile     = _dbe->bookProfile("METTask_CaloMEy_profile",     "MEy [GeV]",     nbinsPV, PVlow, PVup, 200, -500,  500);
  hCaloMET_profile     = _dbe->bookProfile("METTask_CaloMET_profile",     "MET [GeV]",     nbinsPV, PVlow, PVup, 200,    0, 1000);
  hCaloMETNoHF_profile = _dbe->bookProfile("METTask_CaloMETNoHF_profile", "METNoHF [GeV]", nbinsPV, PVlow, PVup, 200,    0, 1000);
  hCaloSumET_profile   = _dbe->bookProfile("METTask_CaloSumET_profile",   "SumET [GeV]",   nbinsPV, PVlow, PVup, 400,    0, 4000);


  // Set NPV profiles x-axis title
  //----------------------------------------------------------------------------
  hCaloMEx_profile    ->setAxisTitle("nvtx", 1);
  hCaloMEy_profile    ->setAxisTitle("nvtx", 1);
  hCaloMET_profile    ->setAxisTitle("nvtx", 1);
  hCaloMETNoHF_profile->setAxisTitle("nvtx", 1);
  hCaloSumET_profile  ->setAxisTitle("nvtx", 1);


  if (_allhist){
    if (bLumiSecPlot){
      hCaloMExLS              = _dbe->book2D("METTask_CaloMEx_LS","METTask_CaloMEx_LS",200,-200,200,50,0.,500.);
      hCaloMExLS->setAxisTitle("MEx [GeV]",1);
      hCaloMExLS->setAxisTitle("Lumi Section",2);
      hCaloMEyLS              = _dbe->book2D("METTask_CaloMEy_LS","METTask_CaloMEy_LS",200,-200,200,50,0.,500.);
      hCaloMEyLS->setAxisTitle("MEy [GeV]",1);
      hCaloMEyLS->setAxisTitle("Lumi Section",2);
    }

    //hCaloMaxEtInEmTowers    = _dbe->book1D("METTask_CaloMaxEtInEmTowers",   "METTask_CaloMaxEtInEmTowers"   ,100,0,2000);
    //hCaloMaxEtInEmTowers->setAxisTitle("Et(Max) in EM Tower [GeV]",1);
    //hCaloMaxEtInHadTowers   = _dbe->book1D("METTask_CaloMaxEtInHadTowers",  "METTask_CaloMaxEtInHadTowers"  ,100,0,2000);
    //hCaloMaxEtInHadTowers->setAxisTitle("Et(Max) in Had Tower [GeV]",1);
    hCaloEtFractionHadronic = _dbe->book1D("METTask_CaloEtFractionHadronic","METTask_CaloEtFractionHadronic",100,0,1);
    hCaloEtFractionHadronic->setAxisTitle("Hadronic Et Fraction",1);
    hCaloEmEtFraction       = _dbe->book1D("METTask_CaloEmEtFraction",      "METTask_CaloEmEtFraction"      ,100,0,1);
    hCaloEmEtFraction->setAxisTitle("EM Et Fraction",1);

    //hCaloEmEtFraction002    = _dbe->book1D("METTask_CaloEmEtFraction002",   "METTask_CaloEmEtFraction002"      ,100,0,1);
    //hCaloEmEtFraction002->setAxisTitle("EM Et Fraction (MET>2 GeV)",1);
    //hCaloEmEtFraction010    = _dbe->book1D("METTask_CaloEmEtFraction010",   "METTask_CaloEmEtFraction010"      ,100,0,1);
    //hCaloEmEtFraction010->setAxisTitle("EM Et Fraction (MET>10 GeV)",1);
    hCaloEmEtFraction020    = _dbe->book1D("METTask_CaloEmEtFraction020",   "METTask_CaloEmEtFraction020"      ,100,0,1);
    hCaloEmEtFraction020->setAxisTitle("EM Et Fraction (MET>20 GeV)",1);
  }

  if (theCaloMETCollectionLabel.label() == "corMetGlobalMuons" ) {
    hCalomuPt    = _dbe->book1D("METTask_CalomuonPt", "METTask_CalomuonPt", 50, 0, 500);
    hCalomuEta   = _dbe->book1D("METTask_CalomuonEta", "METTask_CalomuonEta", 60, -3.0, 3.0);
    hCalomuNhits = _dbe->book1D("METTask_CalomuonNhits", "METTask_CalomuonNhits", 50, 0, 50);
    hCalomuChi2  = _dbe->book1D("METTask_CalomuonNormalizedChi2", "METTask_CalomuonNormalizedChi2", 20, 0, 20);
    hCalomuD0    = _dbe->book1D("METTask_CalomuonD0", "METTask_CalomuonD0", 50, -1, 1);
    hCaloMExCorrection       = _dbe->book1D("METTask_CaloMExCorrection", "METTask_CaloMExCorrection", 100, -500.0,500.0);
    hCaloMEyCorrection       = _dbe->book1D("METTask_CaloMEyCorrection", "METTask_CaloMEyCorrection", 100, -500.0,500.0);
    hCaloMuonCorrectionFlag  = _dbe->book1D("METTask_CaloCorrectionFlag","METTask_CaloCorrectionFlag", 5, -0.5, 4.5);
  }

}

void CaloMETAnalyzer::computeEmHaMET

(

edm::Handle< edm::View< reco::Candidate > >

towers

)

Definition at line 975 of file CaloMETAnalyzer.cc.

References funct::cos(), CaloTower::emEnergy(), edm::View< T >::end(), CaloTower::et(), CaloTower::hadEnergy(), phi, reco::Candidate::phi(), funct::pow(), funct::sin(), reco::Candidate::theta(), and theta().

{

  edm::View<reco::Candidate>::const_iterator towerCand = towers->begin();

  double sum_em_et = 0.0;
  double sum_em_ex = 0.0;
  double sum_em_ey = 0.0;
  double sum_em_ez = 0.0;

  double sum_ha_et = 0.0;
  double sum_ha_ex = 0.0;
  double sum_ha_ey = 0.0;
  double sum_ha_ez = 0.0;

  for ( ; towerCand != towers->end(); towerCand++)
    {
      const reco::Candidate* candidate = &(*towerCand);
      if (candidate)
    {
      const CaloTower* calotower = dynamic_cast<const CaloTower*> (candidate);
      if (calotower){
        double Tower_ET = calotower->et();
        if (Tower_ET>0.3) {

          double phi   = candidate->phi();
          double theta = candidate->theta();
          //double e     = candidate->energy();
          double e_em  = calotower->emEnergy();
          double e_ha  = calotower->hadEnergy();
          double et_em = e_em*sin(theta);
          double et_ha = e_ha*sin(theta);

          sum_em_ez += e_em*cos(theta);
          sum_em_et += et_em;
          sum_em_ex += et_em*cos(phi);
          sum_em_ey += et_em*sin(phi);

          sum_ha_ez += e_ha*cos(theta);
          sum_ha_et += et_ha;
          sum_ha_ex += et_ha*cos(phi);
          sum_ha_ey += et_ha*sin(phi);

        } // Et>0.5
      }   // calotower
    }     // candidate
    }         // loop

  //
  _EmMEx = -sum_em_ex;
  _EmMEy = -sum_em_ey;
  _EmMET = pow(_EmMEx*_EmMEx+_EmMEy*_EmMEy,0.5);
  _EmCaloEz = sum_em_ez;
  _EmSumEt  = sum_em_et;
  _EmMetPhi   = atan2( _EmMEy, _EmMEx );
  //
  _HaMEx = -sum_ha_ex;
  _HaMEy = -sum_ha_ey;
  _HaMET = pow(_HaMEx*_HaMEx+_HaMEy*_HaMEy,0.5);
  _HaCaloEz = sum_ha_ez;
  _HaSumEt  = sum_ha_et;
  _HaMetPhi   = atan2( _HaMEy, _HaMEx );

}

void CaloMETAnalyzer::endJob

(

void

)

Finish up a job.

Definition at line 191 of file CaloMETAnalyzer.cc.

References bTagSequences_cff::jetID.

                             {

  delete jetID;
  delete DCSFilter;

}

void CaloMETAnalyzer::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup,
		DQMStore *	dbe
	)

Finish up a run.

Definition at line 420 of file CaloMETAnalyzer.cc.

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

{

  //
  //--- Check the time length of the Run from the lumi section plots

  std::string dirName = "JetMET/MET/"+_source+"/";
  _dbe->setCurrentFolder(dirName);

  TH1F* tlumisec;

  MonitorElement *meLumiSec = _dbe->get("aaa");
  meLumiSec = _dbe->get("JetMET/lumisec");

  int totlsec=0;
  double totltime=0.;
  if ( meLumiSec->getRootObject() ) {
    tlumisec = meLumiSec->getTH1F();
    for (int i=0; i<500; i++){
      if (tlumisec->GetBinContent(i+1)) totlsec++;
    }
    totltime = double(totlsec*90); // one lumi sec ~ 90 (sec)
  }

  if (totltime==0.) totltime=1.;

  //
  //--- Make the integrated plots with rate (Hz)

  for (std::vector<std::string>::const_iterator ic = _FolderNames.begin(); ic != _FolderNames.end(); ic++)
    {

      std::string DirName;
      DirName = dirName+*ic;

      makeRatePlot(DirName,totltime);
      if ( _HighPtJetEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_HighJetPt",totltime);
      if ( _LowPtJetEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_LowJetPt",totltime);
      if ( _MinBiasEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_MinBias",totltime);
      if ( _HighMETEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_HighMET",totltime);
      //      if ( _LowMETEventFlag->on() )
      //    makeRatePlot(DirName+"/"+"triggerName_LowMET",totltime);
      if ( _EleEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_Ele",totltime);
      if ( _MuonEventFlag->on() )
    makeRatePlot(DirName+"/"+"triggerName_Muon",totltime);
    }
}

void CaloMETAnalyzer::fillMESet	(	const edm::Event &	iEvent,
		std::string	DirName,
		const reco::CaloMET &	calomet,
		const reco::CaloMET &	calometnohf
	)

Definition at line 1102 of file CaloMETAnalyzer.cc.

{

  _dbe->setCurrentFolder(DirName);

  bool bLumiSecPlot=false;
  if (DirName.find("All")) bLumiSecPlot=true;

  if (_trig_JetMB)
    fillMonitorElement(iEvent,DirName,"",calomet,calometnohf, bLumiSecPlot);
  if (_trig_HighPtJet)
    fillMonitorElementTriggered(iEvent,DirName,"HighPtJet",calomet,calometnohf,false);
  if (_trig_LowPtJet)
    fillMonitorElementTriggered(iEvent,DirName,"LowPtJet",calomet,calometnohf,false);
  if (_trig_MinBias)
    fillMonitorElementTriggered(iEvent,DirName,"MinBias",calomet,calometnohf,false);
  if (_trig_HighMET)
    fillMonitorElementTriggered(iEvent,DirName,"HighMET",calomet,calometnohf,false);
  //  if (_trig_LowMET)
  //    fillMonitorElementTriggered(iEvent,DirName,"LowMET",calomet,calometnohf,false);
  if (_trig_Ele)
    fillMonitorElementTriggered(iEvent,DirName,"Ele",calomet,calometnohf,false);
  if (_trig_Muon) {
    fillMonitorElementTriggered(iEvent,DirName,"Muon",calomet,calometnohf,false);
  }
}

void CaloMETAnalyzer::fillMESubSet	(	const edm::Event &	,
		std::string	,
		const reco::CaloMET &	,
		const reco::CaloMET &
	)

void CaloMETAnalyzer::fillMonitorElement	(	const edm::Event &	iEvent,
		std::string	DirName,
		std::string	TriggerTypeName,
		const reco::CaloMET &	calomet,
		const reco::CaloMET &	calometnohf,
		bool	bLumiSecPlot
	)

Definition at line 1132 of file CaloMETAnalyzer.cc.

References gather_cfg::cout, reco::CaloMET::emEtInEB(), reco::CaloMET::emEtInEE(), reco::CaloMET::emEtInHF(), reco::CaloMET::hadEtInHB(), reco::CaloMET::hadEtInHE(), reco::CaloMET::hadEtInHF(), reco::CaloMET::hadEtInHO(), reco::CaloMET::maxEtInEmTowers(), reco::CaloMET::maxEtInHadTowers(), and reco::MET::sumEt().

{
  fillMonitorElementTriggered(iEvent, DirName, TriggerTypeName, calomet, calometnohf, bLumiSecPlot);

  double caloSumET = calomet.sumEt();

  /*
  double caloMETSig = calomet.mEtSig();
  //double caloEz     = calomet.e_longitudinal();
  double caloMET    = calomet.pt();
  double caloMEx    = calomet.px();
  double caloMEy    = calomet.py();
  double caloMETPhi = calomet.phi();

  double caloMETNoHF    = calometnohf.pt();
  //double caloMETNoHFPhi = calometnohf.phi();

  if (_verbose) std::cout << _source << " " << caloMET << std::endl;

  double caloEtFractionHadronic = calomet.etFractionHadronic();
  double caloEmEtFraction       = calomet.emEtFraction();
  */
  double caloMaxEtInEMTowers    = calomet.maxEtInEmTowers();
  double caloMaxEtInHadTowers   = calomet.maxEtInHadTowers();

  double caloHadEtInHB = calomet.hadEtInHB();
  double caloHadEtInHO = calomet.hadEtInHO();
  double caloHadEtInHE = calomet.hadEtInHE();
  double caloHadEtInHF = calomet.hadEtInHF();
  double caloEmEtInEB  = calomet.emEtInEB();
  double caloEmEtInEE  = calomet.emEtInEE();
  double caloEmEtInHF  = calomet.emEtInHF();
  //
  //  int myLuminosityBlock;
  //  myLuminosityBlock = (evtCounter++)/1000;
  //  myLuminosityBlock = iEvent.luminosityBlock();
  //

  if (TriggerTypeName!="") DirName = DirName +"/"+TriggerTypeName;

  if (_verbose) std::cout << "_etThreshold = " << _etThreshold << std::endl;
  if (caloSumET>_etThreshold){
  /*
    hCaloMEx    = _dbe->get(DirName+"/"+"METTask_CaloMEx");    if (hCaloMEx     && hCaloMEx->getRootObject() )    hCaloMEx->Fill(caloMEx);
    hCaloMEy    = _dbe->get(DirName+"/"+"METTask_CaloMEy");    if (hCaloMEy     && hCaloMEy->getRootObject() )    hCaloMEy->Fill(caloMEy);
    hCaloMET    = _dbe->get(DirName+"/"+"METTask_CaloMET");    if (hCaloMET     && hCaloMET->getRootObject() )    hCaloMET->Fill(caloMET);
    hCaloMET1   = _dbe->get(DirName+"/"+"METTask_CaloMET1");   if (hCaloMET1    && hCaloMET1->getRootObject() )   hCaloMET1->Fill(caloMET);
    hCaloMETPhi = _dbe->get(DirName+"/"+"METTask_CaloMETPhi"); if (hCaloMETPhi  && hCaloMETPhi->getRootObject() ) hCaloMETPhi->Fill(caloMETPhi);
    hCaloSumET  = _dbe->get(DirName+"/"+"METTask_CaloSumET");  if (hCaloSumET   && hCaloSumET->getRootObject() )  hCaloSumET->Fill(caloSumET);
    hCaloMETSig = _dbe->get(DirName+"/"+"METTask_CaloMETSig"); if (hCaloMETSig  && hCaloMETSig->getRootObject() ) hCaloMETSig->Fill(caloMETSig);
    //hCaloEz     = _dbe->get(DirName+"/"+"METTask_CaloEz");     if (hCaloEz      && hCaloEz->getRootObject() )     hCaloEz->Fill(caloEz);

    hCaloMETNoHF    = _dbe->get(DirName+"/"+"METTask_CaloMETNoHF");    if (hCaloMETNoHF     && hCaloMETNoHF->getRootObject() )    hCaloMETNoHF->Fill(caloMETNoHF);

    hCaloMET_logx   = _dbe->get(DirName+"/"+"METTask_CaloMET_logx");      if (hCaloMET_logx    && hCaloMET_logx->getRootObject() )   hCaloMET_logx->Fill(log10(caloMET));
    hCaloSumET_logx = _dbe->get(DirName+"/"+"METTask_CaloSumET_logx");    if (hCaloSumET_logx  && hCaloSumET_logx->getRootObject() ) hCaloSumET_logx->Fill(log10(caloSumET));

    //  hCaloMETIonFeedbck = _dbe->get(DirName+"/"+"METTask_CaloMETIonFeedbck"); if (hCaloMETIonFeedbck  && hCaloMETIonFeedbck->getRootObject() ) hCaloMETIonFeedbck->Fill(caloMET);
    //  hCaloMETHPDNoise   = _dbe->get(DirName+"/"+"METTask_CaloMETHPDNoise");   if (hCaloMETHPDNoise    && hCaloMETHPDNoise->getRootObject() )   hCaloMETHPDNoise->Fill(caloMET);

    //hCaloMETPhi002 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi002");    if (caloMET>  2. && hCaloMETPhi002  &&  hCaloMETPhi002->getRootObject()) { hCaloMETPhi002->Fill(caloMETPhi);}
    //hCaloMETPhi010 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi010");    if (caloMET> 10. && hCaloMETPhi010  &&  hCaloMETPhi010->getRootObject()) { hCaloMETPhi010->Fill(caloMETPhi);}
    hCaloMETPhi020 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi020");    if (caloMET> 20. && hCaloMETPhi020  &&  hCaloMETPhi020->getRootObject()) { hCaloMETPhi020->Fill(caloMETPhi);}
  */

    if (_allhist){
      /*
      if (bLumiSecPlot){
    hCaloMExLS = _dbe->get(DirName+"/"+"METTask_CaloMEx_LS");   if (hCaloMExLS  &&  hCaloMExLS->getRootObject())    hCaloMExLS->Fill(caloMEx,myLuminosityBlock);
    hCaloMEyLS = _dbe->get(DirName+"/"+"METTask_CaloMEy_LS");   if (hCaloMEyLS  &&  hCaloMEyLS->getRootObject())    hCaloMEyLS->Fill(caloMEy,myLuminosityBlock);
      }
      hCaloEtFractionHadronic = _dbe->get(DirName+"/"+"METTask_CaloEtFractionHadronic"); if (hCaloEtFractionHadronic && hCaloEtFractionHadronic->getRootObject())  hCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
      hCaloEmEtFraction       = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction");       if (hCaloEmEtFraction       && hCaloEmEtFraction->getRootObject())        hCaloEmEtFraction->Fill(caloEmEtFraction);

      //hCaloEmEtFraction002 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction002");       if (caloMET>  2.  &&  hCaloEmEtFraction002    && hCaloEmEtFraction002->getRootObject()) hCaloEmEtFraction002->Fill(caloEmEtFraction);
      //hCaloEmEtFraction010 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction010");       if (caloMET> 10.  &&  hCaloEmEtFraction010    && hCaloEmEtFraction010->getRootObject()) hCaloEmEtFraction010->Fill(caloEmEtFraction);
      hCaloEmEtFraction020 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction020");       if (caloMET> 20.  &&  hCaloEmEtFraction020    && hCaloEmEtFraction020->getRootObject()) hCaloEmEtFraction020->Fill(caloEmEtFraction);
      */

      hCaloMaxEtInEmTowers  = _dbe->get(DirName+"/"+"METTask_CaloMaxEtInEmTowers");   if (hCaloMaxEtInEmTowers  && hCaloMaxEtInEmTowers->getRootObject())   hCaloMaxEtInEmTowers->Fill(caloMaxEtInEMTowers);
      hCaloMaxEtInHadTowers = _dbe->get(DirName+"/"+"METTask_CaloMaxEtInHadTowers");  if (hCaloMaxEtInHadTowers && hCaloMaxEtInHadTowers->getRootObject())  hCaloMaxEtInHadTowers->Fill(caloMaxEtInHadTowers);

      hCaloHadEtInHB = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHB");  if (hCaloHadEtInHB  &&  hCaloHadEtInHB->getRootObject())  hCaloHadEtInHB->Fill(caloHadEtInHB);
      hCaloHadEtInHO = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHO");  if (hCaloHadEtInHO  &&  hCaloHadEtInHO->getRootObject())  hCaloHadEtInHO->Fill(caloHadEtInHO);
      hCaloHadEtInHE = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHE");  if (hCaloHadEtInHE  &&  hCaloHadEtInHE->getRootObject())  hCaloHadEtInHE->Fill(caloHadEtInHE);
      hCaloHadEtInHF = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHF");  if (hCaloHadEtInHF  &&  hCaloHadEtInHF->getRootObject())  hCaloHadEtInHF->Fill(caloHadEtInHF);
      hCaloEmEtInEB  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInEB");   if (hCaloEmEtInEB   &&  hCaloEmEtInEB->getRootObject())   hCaloEmEtInEB->Fill(caloEmEtInEB);
      hCaloEmEtInEE  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInEE");   if (hCaloEmEtInEE   &&  hCaloEmEtInEE->getRootObject())   hCaloEmEtInEE->Fill(caloEmEtInEE);
      hCaloEmEtInHF  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInHF");   if (hCaloEmEtInHF   &&  hCaloEmEtInHF->getRootObject())   hCaloEmEtInHF->Fill(caloEmEtInHF);

      hCaloEmMEx    = _dbe->get(DirName+"/"+"METTask_CaloEmMEx");     if (hCaloEmMEx    && hCaloEmMEx->getRootObject())     hCaloEmMEx->Fill(_EmMEx);
      hCaloEmMEy    = _dbe->get(DirName+"/"+"METTask_CaloEmMEy");     if (hCaloEmMEy    && hCaloEmMEy->getRootObject())     hCaloEmMEy->Fill(_EmMEy);
      //hCaloEmEz     = _dbe->get(DirName+"/"+"METTask_CaloEmEz");      if (hCaloEmEz     && hCaloEmEz->getRootObject())      hCaloEmEz->Fill(_EmCaloEz);
      hCaloEmMET    = _dbe->get(DirName+"/"+"METTask_CaloEmMET");     if (hCaloEmMET    && hCaloEmMET->getRootObject())     hCaloEmMET->Fill(_EmMET);
      hCaloEmMETPhi = _dbe->get(DirName+"/"+"METTask_CaloEmMETPhi");  if (hCaloEmMETPhi && hCaloEmMETPhi->getRootObject())  hCaloEmMETPhi->Fill(_EmMetPhi);
      //hCaloEmSumET  = _dbe->get(DirName+"/"+"METTask_CaloEmSumET");   if (hCaloEmSumET  && hCaloEmSumET->getRootObject())   hCaloEmSumET->Fill(_EmSumEt);

      hCaloHaMEx    = _dbe->get(DirName+"/"+"METTask_CaloHaMEx");     if (hCaloHaMEx    && hCaloHaMEx->getRootObject())     hCaloHaMEx->Fill(_HaMEx);
      hCaloHaMEy    = _dbe->get(DirName+"/"+"METTask_CaloHaMEy");     if (hCaloHaMEy    && hCaloHaMEy->getRootObject())     hCaloHaMEy->Fill(_HaMEy);
      //hCaloHaEz     = _dbe->get(DirName+"/"+"METTask_CaloHaEz");      if (hCaloHaEz     && hCaloHaEz->getRootObject())      hCaloHaEz->Fill(_HaCaloEz);
      hCaloHaMET    = _dbe->get(DirName+"/"+"METTask_CaloHaMET");     if (hCaloHaMET    && hCaloHaMET->getRootObject())     hCaloHaMET->Fill(_HaMET);
      hCaloHaMETPhi = _dbe->get(DirName+"/"+"METTask_CaloHaMETPhi");  if (hCaloHaMETPhi && hCaloHaMETPhi->getRootObject())  hCaloHaMETPhi->Fill(_HaMetPhi);
      //hCaloHaSumET  = _dbe->get(DirName+"/"+"METTask_CaloHaSumET");   if (hCaloHaSumET  && hCaloHaSumET->getRootObject())   hCaloHaSumET->Fill(_HaSumEt);

    } // _allhist
    /*
    if (theCaloMETCollectionLabel.label() == "corMetGlobalMuons" ) {

      for( reco::MuonCollection::const_iterator muonit = muon_h->begin(); muonit != muon_h->end(); muonit++ ) {
    const reco::TrackRef siTrack = muonit->innerTrack();
    hCalomuPt    = _dbe->get(DirName+"/"+"METTask_CalomuPt");     if (hCalomuPt    && hCalomuPt->getRootObject())     hCalomuPt->Fill( muonit->p4().pt() );
    hCalomuEta   = _dbe->get(DirName+"/"+"METTask_CalomuEta");    if (hCalomuEta   && hCalomuEta->getRootObject())    hCalomuEta->Fill( muonit->p4().eta() );
    hCalomuNhits = _dbe->get(DirName+"/"+"METTask_CalomuNhits");  if (hCalomuNhits && hCalomuNhits->getRootObject())  hCalomuNhits->Fill( siTrack.isNonnull() ? siTrack->numberOfValidHits() : -999 );
    hCalomuChi2  = _dbe->get(DirName+"/"+"METTask_CalomuChi2");   if (hCalomuChi2  && hCalomuChi2->getRootObject())   hCalomuChi2->Fill( siTrack.isNonnull() ? siTrack->chi2()/siTrack->ndof() : -999 );
    double d0 = siTrack.isNonnull() ? -1 * siTrack->dxy( bspot) : -999;
    hCalomuD0    = _dbe->get(DirName+"/"+"METTask_CalomuD0");     if (hCalomuD0    && hCalomuD0->getRootObject())  hCalomuD0->Fill( d0 );
      }

      const unsigned int nMuons = muon_h->size();
      for( unsigned int mus = 0; mus < nMuons; mus++ ) {
    reco::MuonRef muref( muon_h, mus);
    reco::MuonMETCorrectionData muCorrData = (*corMetGlobalMuons_ValueMap_Handle)[muref];
    hCaloMExCorrection      = _dbe->get(DirName+"/"+"METTask_CaloMExCorrection");       if (hCaloMExCorrection      && hCaloMExCorrection->getRootObject())       hCaloMExCorrection-> Fill(muCorrData.corrY());
    hCaloMEyCorrection      = _dbe->get(DirName+"/"+"METTask_CaloMEyCorrection");       if (hCaloMEyCorrection      && hCaloMEyCorrection->getRootObject())       hCaloMEyCorrection-> Fill(muCorrData.corrX());
    hCaloMuonCorrectionFlag = _dbe->get(DirName+"/"+"METTask_CaloMuonCorrectionFlag");  if (hCaloMuonCorrectionFlag && hCaloMuonCorrectionFlag->getRootObject())  hCaloMuonCorrectionFlag-> Fill(muCorrData.type());
      }
    }
    */
  } // et threshold cut

}

void CaloMETAnalyzer::fillMonitorElementTriggered	(	const edm::Event &	iEvent,
		std::string	DirName,
		std::string	TriggerTypeName,
		const reco::CaloMET &	calomet,
		const reco::CaloMET &	calometnohf,
		bool	bLumiSecPlot
	)

Definition at line 1269 of file CaloMETAnalyzer.cc.

References reco::MuonMETCorrectionData::corrX(), reco::MuonMETCorrectionData::corrY(), gather_cfg::cout, reco::CaloMET::emEtFraction(), reco::CaloMET::etFractionHadronic(), HcalObjRepresent::Fill(), edm::Ref< C, T, F >::get(), edm::Ref< C, T, F >::isNonnull(), edm::EventBase::luminosityBlock(), reco::MET::mEtSig(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::MET::sumEt(), and reco::MuonMETCorrectionData::type().

{

  if (TriggerTypeName=="HighPtJet") {
    if (!selectHighPtJetEvent(iEvent)) return;
  }
  else if (TriggerTypeName=="LowPtJet") {
    if (!selectLowPtJetEvent(iEvent)) return;
  }
  else if (TriggerTypeName=="HighMET") {
    if (calomet.pt()<_highMETThreshold) return;
  }
  //  else if (TriggerTypeName=="LowMET") {
  //    if (calomet.pt()<_lowMETThreshold) return;
  //  }
  else if (TriggerTypeName=="Ele") {
    if (!selectWElectronEvent(iEvent)) return;
  }
  else if (TriggerTypeName=="Muon") {
    if (!selectWMuonEvent(iEvent)) return;
  }

  double caloSumET  = calomet.sumEt();
  double caloMETSig = calomet.mEtSig();
  //double caloEz     = calomet.e_longitudinal();
  double caloMET    = calomet.pt();
  double caloMEx    = calomet.px();
  double caloMEy    = calomet.py();
  double caloMETPhi = calomet.phi();

  double caloMETNoHF    = calometnohf.pt();
  //double caloMETNoHFPhi = calometnohf.phi();

  if (_verbose) std::cout << _source << " " << caloMET << std::endl;

  double caloEtFractionHadronic = calomet.etFractionHadronic();
  double caloEmEtFraction       = calomet.emEtFraction();

  //double caloMaxEtInEMTowers    = calomet.maxEtInEmTowers();
  //double caloMaxEtInHadTowers   = calomet.maxEtInHadTowers();

  //double caloHadEtInHB = calomet.hadEtInHB();
  //double caloHadEtInHO = calomet.hadEtInHO();
  //double caloHadEtInHE = calomet.hadEtInHE();
  //double caloHadEtInHF = calomet.hadEtInHF();
  //double caloEmEtInEB  = calomet.emEtInEB();
  //double caloEmEtInEE  = calomet.emEtInEE();
  //double caloEmEtInHF  = calomet.emEtInHF();

  //
  int myLuminosityBlock;
  //  myLuminosityBlock = (evtCounter++)/1000;
  myLuminosityBlock = iEvent.luminosityBlock();
  //

  if (TriggerTypeName!="") DirName = DirName +"/"+TriggerTypeName;

  if (_verbose) std::cout << "_etThreshold = " << _etThreshold << std::endl;
  if (caloSumET>_etThreshold){
    hCaloMEx    = _dbe->get(DirName+"/"+"METTask_CaloMEx");    if (hCaloMEx     && hCaloMEx->getRootObject() )    hCaloMEx->Fill(caloMEx);
    hCaloMEy    = _dbe->get(DirName+"/"+"METTask_CaloMEy");    if (hCaloMEy     && hCaloMEy->getRootObject() )    hCaloMEy->Fill(caloMEy);
    hCaloMET    = _dbe->get(DirName+"/"+"METTask_CaloMET");    if (hCaloMET     && hCaloMET->getRootObject() )    hCaloMET->Fill(caloMET);
    hCaloMET1   = _dbe->get(DirName+"/"+"METTask_CaloMET1");   if (hCaloMET1    && hCaloMET1->getRootObject() )   hCaloMET1->Fill(caloMET);
    hCaloMETPhi = _dbe->get(DirName+"/"+"METTask_CaloMETPhi"); if (hCaloMETPhi  && hCaloMETPhi->getRootObject() ) hCaloMETPhi->Fill(caloMETPhi);
    hCaloSumET  = _dbe->get(DirName+"/"+"METTask_CaloSumET");  if (hCaloSumET   && hCaloSumET->getRootObject() )  hCaloSumET->Fill(caloSumET);
    hCaloMETSig = _dbe->get(DirName+"/"+"METTask_CaloMETSig"); if (hCaloMETSig  && hCaloMETSig->getRootObject() ) hCaloMETSig->Fill(caloMETSig);
    //hCaloEz     = _dbe->get(DirName+"/"+"METTask_CaloEz");     if (hCaloEz      && hCaloEz->getRootObject() )     hCaloEz->Fill(caloEz);

    hCaloMETNoHF    = _dbe->get(DirName+"/"+"METTask_CaloMETNoHF");    if (hCaloMETNoHF     && hCaloMETNoHF->getRootObject() )    hCaloMETNoHF->Fill(caloMETNoHF);

    hCaloMET_logx   = _dbe->get(DirName+"/"+"METTask_CaloMET_logx");      if (hCaloMET_logx    && hCaloMET_logx->getRootObject() )   hCaloMET_logx->Fill(log10(caloMET));
    hCaloSumET_logx = _dbe->get(DirName+"/"+"METTask_CaloSumET_logx");    if (hCaloSumET_logx  && hCaloSumET_logx->getRootObject() ) hCaloSumET_logx->Fill(log10(caloSumET));


    // Fill NPV profiles
    //--------------------------------------------------------------------------
    hCaloMEx_profile     = _dbe->get(DirName + "/METTask_CaloMEx_profile");
    hCaloMEy_profile     = _dbe->get(DirName + "/METTask_CaloMEy_profile");
    hCaloMET_profile     = _dbe->get(DirName + "/METTask_CaloMET_profile");
    hCaloMETNoHF_profile = _dbe->get(DirName + "/METTask_CaloMETNoHF_profile");
    hCaloSumET_profile   = _dbe->get(DirName + "/METTask_CaloSumET_profile");

    if (hCaloMEx_profile     && hCaloMEx_profile    ->getRootObject()) hCaloMEx_profile    ->Fill(_numPV, caloMEx);
    if (hCaloMEy_profile     && hCaloMEy_profile    ->getRootObject()) hCaloMEy_profile    ->Fill(_numPV, caloMEy);
    if (hCaloMET_profile     && hCaloMET_profile    ->getRootObject()) hCaloMET_profile    ->Fill(_numPV, caloMET);
    if (hCaloMETNoHF_profile && hCaloMETNoHF_profile->getRootObject()) hCaloMETNoHF_profile->Fill(_numPV, caloMETNoHF);
    if (hCaloSumET_profile   && hCaloSumET_profile  ->getRootObject()) hCaloSumET_profile  ->Fill(_numPV, caloSumET);


    //hCaloEz     = _dbe->get(DirName+"/"+"METTask_CaloEz_profile");     if (hCaloEz      && hCaloEz_profile->getRootObject() )     hCaloEz_profile->Fill(caloEz);
    //  hCaloMETIonFeedbck = _dbe->get(DirName+"/"+"METTask_CaloMETIonFeedbck"); if (hCaloMETIonFeedbck  && hCaloMETIonFeedbck->getRootObject() ) hCaloMETIonFeedbck->Fill(caloMET);
    //  haloMETHPDNoise   = _dbe->get(DirName+"/"+"METTask_CaloMETHPDNoise");   if (hCaloMETHPDNoise    && hCaloMETHPDNoise->getRootObject() )   hCaloMETHPDNoise->Fill(caloMET);

    //hCaloMETPhi002 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi002");    if (caloMET>  2. && hCaloMETPhi002  &&  hCaloMETPhi002->getRootObject()) { hCaloMETPhi002->Fill(caloMETPhi);}
    //hCaloMETPhi010 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi010");    if (caloMET> 10. && hCaloMETPhi010  &&  hCaloMETPhi010->getRootObject()) { hCaloMETPhi010->Fill(caloMETPhi);}
    hCaloMETPhi020 = _dbe->get(DirName+"/"+"METTask_CaloMETPhi020");    if (caloMET> 20. && hCaloMETPhi020  &&  hCaloMETPhi020->getRootObject()) { hCaloMETPhi020->Fill(caloMETPhi);}

    if (_allhist){
      if (bLumiSecPlot){
    hCaloMExLS = _dbe->get(DirName+"/"+"METTask_CaloMEx_LS");   if (hCaloMExLS  &&  hCaloMExLS->getRootObject())    hCaloMExLS->Fill(caloMEx,myLuminosityBlock);
    hCaloMEyLS = _dbe->get(DirName+"/"+"METTask_CaloMEy_LS");   if (hCaloMEyLS  &&  hCaloMEyLS->getRootObject())    hCaloMEyLS->Fill(caloMEy,myLuminosityBlock);
      }

      hCaloEtFractionHadronic = _dbe->get(DirName+"/"+"METTask_CaloEtFractionHadronic"); if (hCaloEtFractionHadronic && hCaloEtFractionHadronic->getRootObject())  hCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
      hCaloEmEtFraction       = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction");       if (hCaloEmEtFraction       && hCaloEmEtFraction->getRootObject())        hCaloEmEtFraction->Fill(caloEmEtFraction);

      //hCaloEmEtFraction002 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction002");       if (caloMET>  2.  &&  hCaloEmEtFraction002    && hCaloEmEtFraction002->getRootObject()) hCaloEmEtFraction002->Fill(caloEmEtFraction);
      //hCaloEmEtFraction010 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction010");       if (caloMET> 10.  &&  hCaloEmEtFraction010    && hCaloEmEtFraction010->getRootObject()) hCaloEmEtFraction010->Fill(caloEmEtFraction);
      hCaloEmEtFraction020 = _dbe->get(DirName+"/"+"METTask_CaloEmEtFraction020");       if (caloMET> 20.  &&  hCaloEmEtFraction020    && hCaloEmEtFraction020->getRootObject()) hCaloEmEtFraction020->Fill(caloEmEtFraction);

      //hCaloMaxEtInEmTowers  = _dbe->get(DirName+"/"+"METTask_CaloMaxEtInEmTowers");   if (hCaloMaxEtInEmTowers  && hCaloMaxEtInEmTowers->getRootObject())   hCaloMaxEtInEmTowers->Fill(caloMaxEtInEMTowers);
      //hCaloMaxEtInHadTowers = _dbe->get(DirName+"/"+"METTask_CaloMaxEtInHadTowers");  if (hCaloMaxEtInHadTowers && hCaloMaxEtInHadTowers->getRootObject())  hCaloMaxEtInHadTowers->Fill(caloMaxEtInHadTowers);
      //
      //hCaloHadEtInHB = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHB");  if (hCaloHadEtInHB  &&  hCaloHadEtInHB->getRootObject())  hCaloHadEtInHB->Fill(caloHadEtInHB);
      //hCaloHadEtInHO = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHO");  if (hCaloHadEtInHO  &&  hCaloHadEtInHO->getRootObject())  hCaloHadEtInHO->Fill(caloHadEtInHO);
      //hCaloHadEtInHE = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHE");  if (hCaloHadEtInHE  &&  hCaloHadEtInHE->getRootObject())  hCaloHadEtInHE->Fill(caloHadEtInHE);
      //hCaloHadEtInHF = _dbe->get(DirName+"/"+"METTask_CaloHadEtInHF");  if (hCaloHadEtInHF  &&  hCaloHadEtInHF->getRootObject())  hCaloHadEtInHF->Fill(caloHadEtInHF);
      //hCaloEmEtInEB  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInEB");   if (hCaloEmEtInEB   &&  hCaloEmEtInEB->getRootObject())   hCaloEmEtInEB->Fill(caloEmEtInEB);
      //hCaloEmEtInEE  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInEE");   if (hCaloEmEtInEE   &&  hCaloEmEtInEE->getRootObject())   hCaloEmEtInEE->Fill(caloEmEtInEE);
      //hCaloEmEtInHF  = _dbe->get(DirName+"/"+"METTask_CaloEmEtInHF");   if (hCaloEmEtInHF   &&  hCaloEmEtInHF->getRootObject())   hCaloEmEtInHF->Fill(caloEmEtInHF);
      //
      //hCaloEmMEx    = _dbe->get(DirName+"/"+"METTask_CaloEmMEx");     if (hCaloEmMEx    && hCaloEmMEx->getRootObject())     hCaloEmMEx->Fill(_EmMEx);
      //hCaloEmMEy    = _dbe->get(DirName+"/"+"METTask_CaloEmMEy");     if (hCaloEmMEy    && hCaloEmMEy->getRootObject())     hCaloEmMEy->Fill(_EmMEy);
      //hCaloEmMET    = _dbe->get(DirName+"/"+"METTask_CaloEmMET");     if (hCaloEmMET    && hCaloEmMET->getRootObject())     hCaloEmMET->Fill(_EmMET);
      //hCaloEmMETPhi = _dbe->get(DirName+"/"+"METTask_CaloEmMETPhi");  if (hCaloEmMETPhi && hCaloEmMETPhi->getRootObject())  hCaloEmMETPhi->Fill(_EmMetPhi);
      //
      //hCaloHaMEx    = _dbe->get(DirName+"/"+"METTask_CaloHaMEx");     if (hCaloHaMEx    && hCaloHaMEx->getRootObject())     hCaloHaMEx->Fill(_HaMEx);
      //hCaloHaMEy    = _dbe->get(DirName+"/"+"METTask_CaloHaMEy");     if (hCaloHaMEy    && hCaloHaMEy->getRootObject())     hCaloHaMEy->Fill(_HaMEy);
      //hCaloHaMET    = _dbe->get(DirName+"/"+"METTask_CaloHaMET");     if (hCaloHaMET    && hCaloHaMET->getRootObject())     hCaloHaMET->Fill(_HaMET);
      //hCaloHaMETPhi = _dbe->get(DirName+"/"+"METTask_CaloHaMETPhi");  if (hCaloHaMETPhi && hCaloHaMETPhi->getRootObject())  hCaloHaMETPhi->Fill(_HaMetPhi);

    } // _allhist
    if (theCaloMETCollectionLabel.label() == "corMetGlobalMuons" ) {

      for( reco::MuonCollection::const_iterator muonit = muon_h->begin(); muonit != muon_h->end(); muonit++ ) {
    const reco::TrackRef siTrack = muonit->innerTrack();
    hCalomuPt    = _dbe->get(DirName+"/"+"METTask_CalomuPt");     if (hCalomuPt    && hCalomuPt->getRootObject())     hCalomuPt->Fill( muonit->p4().pt() );
    hCalomuEta   = _dbe->get(DirName+"/"+"METTask_CalomuEta");    if (hCalomuEta   && hCalomuEta->getRootObject())    hCalomuEta->Fill( muonit->p4().eta() );
    hCalomuNhits = _dbe->get(DirName+"/"+"METTask_CalomuNhits");  if (hCalomuNhits && hCalomuNhits->getRootObject())  hCalomuNhits->Fill( siTrack.isNonnull() ? siTrack->numberOfValidHits() : -999 );
    hCalomuChi2  = _dbe->get(DirName+"/"+"METTask_CalomuChi2");   if (hCalomuChi2  && hCalomuChi2->getRootObject())   hCalomuChi2->Fill( siTrack.isNonnull() ? siTrack->chi2()/siTrack->ndof() : -999 );
    double d0 = siTrack.isNonnull() ? -1 * siTrack->dxy( bspot) : -999;
    hCalomuD0    = _dbe->get(DirName+"/"+"METTask_CalomuD0");     if (hCalomuD0    && hCalomuD0->getRootObject())  hCalomuD0->Fill( d0 );
      }

      const unsigned int nMuons = muon_h->size();
      for( unsigned int mus = 0; mus < nMuons; mus++ ) {
    reco::MuonRef muref( muon_h, mus);
    reco::MuonMETCorrectionData muCorrData = (*corMetGlobalMuons_ValueMap_Handle)[muref];
    hCaloMExCorrection      = _dbe->get(DirName+"/"+"METTask_CaloMExCorrection");       if (hCaloMExCorrection      && hCaloMExCorrection->getRootObject())       hCaloMExCorrection-> Fill(muCorrData.corrY());
    hCaloMEyCorrection      = _dbe->get(DirName+"/"+"METTask_CaloMEyCorrection");       if (hCaloMEyCorrection      && hCaloMEyCorrection->getRootObject())       hCaloMEyCorrection-> Fill(muCorrData.corrX());
    hCaloMuonCorrectionFlag = _dbe->get(DirName+"/"+"METTask_CaloMuonCorrectionFlag");  if (hCaloMuonCorrectionFlag && hCaloMuonCorrectionFlag->getRootObject())  hCaloMuonCorrectionFlag-> Fill(muCorrData.type());
      }
    }
  } // et threshold cut

}

void CaloMETAnalyzer::makeRatePlot	(	std::string	DirName,
		double	totltime
	)

Definition at line 475 of file CaloMETAnalyzer.cc.

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

{

  _dbe->setCurrentFolder(DirName);
  MonitorElement *meCaloMET = _dbe->get(DirName+"/"+"METTask_CaloMET");

  TH1F* tCaloMET;
  TH1F* tCaloMETRate;

  if ( meCaloMET )
    if ( meCaloMET->getRootObject() ) {
      tCaloMET     = meCaloMET->getTH1F();

      // Integral plot & convert number of events to rate (hz)
      tCaloMETRate = (TH1F*) tCaloMET->Clone("METTask_CaloMETRate");
      for (int i = tCaloMETRate->GetNbinsX()-1; i>=0; i--){
    tCaloMETRate->SetBinContent(i+1,tCaloMETRate->GetBinContent(i+2)+tCaloMET->GetBinContent(i+1));
      }
      for (int i = 0; i<tCaloMETRate->GetNbinsX(); i++){
    tCaloMETRate->SetBinContent(i+1,tCaloMETRate->GetBinContent(i+1)/double(totltime));
      }

      tCaloMETRate->SetName("METTask_CaloMETRate");
      tCaloMETRate->SetTitle("METTask_CaloMETRate");
      hCaloMETRate = _dbe->book1D("METTask_CaloMETRate",tCaloMETRate);
      hCaloMETRate->setAxisTitle("MET Threshold [GeV]",1);
    }
}

bool CaloMETAnalyzer::selectHighPtJetEvent

(

const edm::Event &

iEvent

)

Definition at line 1435 of file CaloMETAnalyzer.cc.

References FamosSequences_cff::caloJets, gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), and LogDebug.

                                                                {

  bool return_value=false;

  edm::Handle<reco::CaloJetCollection> caloJets;
  iEvent.getByLabel(theJetCollectionLabel, caloJets);
  if (!caloJets.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find jet product" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find jet product" << std::endl;
  }

  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
       cal!=caloJets->end(); ++cal){
    if (cal->pt()>_highPtJetThreshold){
      return_value=true;
    }
  }

  return return_value;

}

bool CaloMETAnalyzer::selectLowPtJetEvent

(

const edm::Event &

iEvent

)

Definition at line 1458 of file CaloMETAnalyzer.cc.

References FamosSequences_cff::caloJets, gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), and LogDebug.

                                                               {

  bool return_value=false;

  edm::Handle<reco::CaloJetCollection> caloJets;
  iEvent.getByLabel(theJetCollectionLabel, caloJets);
  if (!caloJets.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find jet product" << std::endl;
    if (_verbose) std::cout << "CaloMETAnalyzer: Could not find jet product" << std::endl;
  }

  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
       cal!=caloJets->end(); ++cal){
    if (cal->pt()>_lowPtJetThreshold){
      return_value=true;
    }
  }

  return return_value;

}

bool CaloMETAnalyzer::selectWElectronEvent

(

const edm::Event &

iEvent

)

Definition at line 1481 of file CaloMETAnalyzer.cc.

                                                                {

  bool return_value=true;

  /*
    W-electron event selection comes here
   */

  return return_value;

}

bool CaloMETAnalyzer::selectWMuonEvent

(

const edm::Event &

iEvent

)

Definition at line 1494 of file CaloMETAnalyzer.cc.

                                                            {

  bool return_value=true;

  /*
    W-muon event selection comes here
   */

  return return_value;

}

void CaloMETAnalyzer::validateMET	(	const reco::CaloMET &	calomet,
		edm::Handle< edm::View< reco::Candidate > >	towers
	)

Definition at line 1040 of file CaloMETAnalyzer.cc.

References funct::cos(), gather_cfg::cout, alignCSCRings::e, edm::View< T >::end(), reco::Candidate::energy(), CaloTower::et(), phi, reco::Candidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), funct::sin(), mathSSE::sqrt(), reco::MET::sumEt(), reco::Candidate::theta(), and theta().

{

  edm::View<reco::Candidate>::const_iterator towerCand = towers->begin();

  double sum_et = 0.0;
  double sum_ex = 0.0;
  double sum_ey = 0.0;
  double sum_ez = 0.0;

  for ( ; towerCand != towers->end(); towerCand++)
    {
      const reco::Candidate* candidate = &(*towerCand);
      if (candidate)
    {
      const CaloTower* calotower = dynamic_cast<const CaloTower*> (candidate);
      if (calotower){
        double Tower_ET = calotower->et();
        if (Tower_ET>0.3) {

          double phi   = candidate->phi();
          double theta = candidate->theta();
          double e     = candidate->energy();
          double et    = e*sin(theta);
          sum_ez += e*cos(theta);
          sum_et += et;
          sum_ex += et*cos(phi);
          sum_ey += et*sin(phi);

        } // Et>0.5
      }   // calotower
    }     // candidate
    }         // loop

  double Mex   = -sum_ex;
  double Mey   = -sum_ey;
  //double Mez   = -sum_ez;
  double Met   = sqrt( sum_ex*sum_ex + sum_ey*sum_ey );
  double Sumet = sum_et;
  //double MetPhi   = atan2( -sum_ey, -sum_ex ); // since MET is now a candidate,

  if (_verbose){
    if (Sumet!=calomet.sumEt() || Mex!=calomet.px() || Mey!=calomet.py() || Met!=calomet.pt() ){
      std::cout << _source << std::endl;
      std::cout << "SUMET" << Sumet << " METBlock" << calomet.sumEt() << std::endl;
      std::cout << "MEX"   << Mex   << " METBlock" << calomet.px()    << std::endl;
      std::cout << "MEY"   << Mey   << " METBlock" << calomet.py()    << std::endl;
      std::cout << "MET"   << Met   << " METBlock" << calomet.pt()    << std::endl;
    }
  }

  if (_print){
    std::cout << "SUMET = " << calomet.sumEt() << " (GeV) "
          << "MEX"   << calomet.px() << " (GeV) "
          << "MEY"   << calomet.py() << " (GeV) "
          << "MET"   << calomet.pt() << " (GeV) " << std::endl;
  }

}

Member Data Documentation

bool CaloMETAnalyzer::_allhist

private

Definition at line 186 of file CaloMETAnalyzer.h.

bool CaloMETAnalyzer::_allSelection

private

Definition at line 187 of file CaloMETAnalyzer.h.

bool CaloMETAnalyzer::_cleanupSelection

private

Definition at line 188 of file CaloMETAnalyzer.h.

DQMStore* CaloMETAnalyzer::_dbe

private

Definition at line 217 of file CaloMETAnalyzer.h.

bool CaloMETAnalyzer::_doHLTPhysicsOn

private

Definition at line 150 of file CaloMETAnalyzer.h.

bool CaloMETAnalyzer::_doPVCheck

private

Definition at line 149 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_EleEventFlag

private

Definition at line 131 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmCaloEz

private

Definition at line 198 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmMET

private

Definition at line 197 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmMetPhi

private

Definition at line 200 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmMEx

private

Definition at line 195 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmMEy

private

Definition at line 196 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_EmSumEt

private

Definition at line 199 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_etThreshold

private

Definition at line 177 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::_FolderNames

private

Definition at line 191 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaCaloEz

private

Definition at line 205 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaMET

private

Definition at line 204 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaMetPhi

private

Definition at line 207 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaMEx

private

Definition at line 202 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaMEy

private

Definition at line 203 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_HaSumEt

private

Definition at line 206 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_HighMETEventFlag

private

Definition at line 129 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_highMETThreshold

private

Definition at line 173 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_HighPtJetEventFlag

private

Definition at line 126 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_highPtJetThreshold

private

Definition at line 171 of file CaloMETAnalyzer.h.

std::string CaloMETAnalyzer::_hlt_PhysDec

private

Definition at line 143 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_LowPtJetEventFlag

private

Definition at line 127 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_lowPtJetThreshold

private

Definition at line 172 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_MinBiasEventFlag

private

Definition at line 128 of file CaloMETAnalyzer.h.

GenericTriggerEventFlag* CaloMETAnalyzer::_MuonEventFlag

private

Definition at line 132 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_numPV

private

Definition at line 193 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_nvtx_min

private

Definition at line 155 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_nvtxtrks_min

private

Definition at line 156 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_print

private

Definition at line 105 of file CaloMETAnalyzer.h.

std::string CaloMETAnalyzer::_source

private

Definition at line 108 of file CaloMETAnalyzer.h.

std::vector<unsigned > CaloMETAnalyzer::_techTrigsAND

private

Definition at line 145 of file CaloMETAnalyzer.h.

std::vector<unsigned > CaloMETAnalyzer::_techTrigsNOT

private

Definition at line 147 of file CaloMETAnalyzer.h.

std::vector<unsigned > CaloMETAnalyzer::_techTrigsOR

private

Definition at line 146 of file CaloMETAnalyzer.h.

bool CaloMETAnalyzer::_tightBHFiltering

private

Definition at line 152 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_tightJetIDFiltering

private

Definition at line 153 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_Ele

private

Definition at line 167 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_HighMET

private

Definition at line 165 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_HighPtJet

private

Definition at line 162 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_JetMB

private

Definition at line 161 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_LowPtJet

private

Definition at line 163 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_MinBias

private

Definition at line 164 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_Muon

private

Definition at line 168 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_trig_PhysDec

private

Definition at line 169 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_verbose

private

Definition at line 104 of file CaloMETAnalyzer.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 CaloMETAnalyzer::_vtxchi2_max

private

Definition at line 158 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::_vtxndof_min

private

Definition at line 157 of file CaloMETAnalyzer.h.

double CaloMETAnalyzer::_vtxz_max

private

Definition at line 159 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::BeamHaloSummaryTag

private

Definition at line 115 of file CaloMETAnalyzer.h.

edm::Handle< reco::BeamSpot > CaloMETAnalyzer::beamSpot_h

private

Definition at line 214 of file CaloMETAnalyzer.h.

math::XYZPoint CaloMETAnalyzer::bspot

private

Definition at line 210 of file CaloMETAnalyzer.h.

edm::Handle< edm::ValueMap<reco::MuonMETCorrectionData> > CaloMETAnalyzer::corMetGlobalMuons_ValueMap_Handle

private

Definition at line 212 of file CaloMETAnalyzer.h.

JetMETDQMDCSFilter* CaloMETAnalyzer::DCSFilter

private

Definition at line 183 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::elecExpr_

private

Definition at line 139 of file CaloMETAnalyzer.h.

int CaloMETAnalyzer::evtCounter

Definition at line 97 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::gtTag

private

Definition at line 118 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::HBHENoiseFilterResultTag

private

Definition at line 116 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::HcalNoiseRBXCollectionTag

private

Definition at line 114 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmEtFraction

private

Definition at line 272 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmEtFraction020

private

Definition at line 276 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmEtInEB

private

Definition at line 284 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmEtInEE

private

Definition at line 283 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmEtInHF

private

Definition at line 282 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmMET

private

Definition at line 289 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmMETPhi

private

Definition at line 290 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmMEx

private

Definition at line 286 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEmMEy

private

Definition at line 287 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloEtFractionHadronic

private

Definition at line 271 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHadEtInHB

private

Definition at line 278 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHadEtInHE

private

Definition at line 280 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHadEtInHF

private

Definition at line 281 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHadEtInHO

private

Definition at line 279 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHaMET

private

Definition at line 296 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHaMETPhi

private

Definition at line 297 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHaMEx

private

Definition at line 293 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloHaMEy

private

Definition at line 294 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMaxEtInEmTowers

private

Definition at line 269 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMaxEtInHadTowers

private

Definition at line 270 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMET

private

Definition at line 238 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMET1

private

Definition at line 239 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMET_logx

private

Definition at line 255 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMET_profile

private

Definition at line 250 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETNoHF

private

Definition at line 242 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETNoHF_profile

private

Definition at line 251 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETPhi

private

Definition at line 240 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETPhi020

private

Definition at line 264 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETRate

private

Definition at line 230 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMETSig

private

Definition at line 237 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEx

private

Definition at line 234 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEx_profile

private

Definition at line 248 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMExCorrection

private

Definition at line 305 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMExLS

private

Definition at line 266 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEy

private

Definition at line 235 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEy_profile

private

Definition at line 249 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEyCorrection

private

Definition at line 306 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMEyLS

private

Definition at line 267 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCalomuChi2

private

Definition at line 303 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCalomuD0

private

Definition at line 304 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCalomuEta

private

Definition at line 301 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCalomuNhits

private

Definition at line 302 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloMuonCorrectionFlag

private

Definition at line 307 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCalomuPt

private

Definition at line 300 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloSumET

private

Definition at line 241 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloSumET_logx

private

Definition at line 256 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hCaloSumET_profile

private

Definition at line 252 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::highMETExpr_

private

Definition at line 136 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::highPtJetExpr_

private

Definition at line 134 of file CaloMETAnalyzer.h.

std::vector<std::string > CaloMETAnalyzer::HLTPathsJetMBByName_

private

Definition at line 124 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hmetME

private

Definition at line 232 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_Ele

private

Definition at line 225 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_HighMET

private

Definition at line 223 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_HighPtJet

private

Definition at line 220 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_LowPtJet

private

Definition at line 221 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_MinBias

private

Definition at line 222 of file CaloMETAnalyzer.h.

MonitorElement* CaloMETAnalyzer::hTriggerName_Muon

private

Definition at line 226 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::inputBeamSpotLabel

private

Definition at line 121 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::inputMuonLabel

private

Definition at line 120 of file CaloMETAnalyzer.h.

reco::helper::JetIDHelper* CaloMETAnalyzer::jetID

private

Definition at line 180 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::lowPtJetExpr_

private

Definition at line 135 of file CaloMETAnalyzer.h.

std::string CaloMETAnalyzer::metname

private

Definition at line 107 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::minbiasExpr_

private

Definition at line 140 of file CaloMETAnalyzer.h.

edm::Handle< reco::MuonCollection > CaloMETAnalyzer::muon_h

private

Definition at line 213 of file CaloMETAnalyzer.h.

std::vector<std::string> CaloMETAnalyzer::muonExpr_

private

Definition at line 138 of file CaloMETAnalyzer.h.

edm::ParameterSet CaloMETAnalyzer::parameters

private

Definition at line 102 of file CaloMETAnalyzer.h.

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

edm::InputTag CaloMETAnalyzer::theCaloMETCollectionLabel

private

Definition at line 110 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::theCaloTowersLabel

private

Definition at line 112 of file CaloMETAnalyzer.h.

edm::ParameterSet CaloMETAnalyzer::theCleaningParameters

private

Definition at line 142 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::theJetCollectionLabel

private

Definition at line 113 of file CaloMETAnalyzer.h.

edm::InputTag CaloMETAnalyzer::vertexTag

private

Definition at line 117 of file CaloMETAnalyzer.h.