METAnalyzer Class Reference

#include <METAnalyzer.h>

Inheritance diagram for METAnalyzer:

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

Public Attributes
int	evtCounter

Private Attributes
bool	_allhist
bool	_allSelection
bool	_cleanupSelection
DQMStore *	_dbe
bool	_doHLTPhysicsOn
bool	_doPVCheck
GenericTriggerEventFlag *	_EleEventFlag
double	_etThreshold
std::string	_FolderName
std::vector< std::string >	_FolderNames
GenericTriggerEventFlag *	_HighMETEventFlag
double	_highMETThreshold
GenericTriggerEventFlag *	_HighPtJetEventFlag
double	_highPtJetThreshold
std::string	_hlt_PhysDec
GenericTriggerEventFlag *	_LowPtJetEventFlag
double	_lowPtJetThreshold
GenericTriggerEventFlag *	_MinBiasEventFlag
GenericTriggerEventFlag *	_MuonEventFlag
int	_nvtx_min
int	_nvtxtrks_min
std::string	_source
std::vector< unsigned >	_techTrigsAND
std::vector< unsigned >	_techTrigsNOT
std::vector< unsigned >	_techTrigsOR
bool	_tightBHFiltering
int	_tightJetIDFiltering
int	_trig_Ele
int	_trig_HighMET
int	_trig_HighPtJet
int	_trig_JetMB
int	_trig_LowPtJet
int	_trig_MinBias
int	_trig_Muon
int	_trig_PhysDec
int	_verbose
double	_vtxchi2_max
int	_vtxndof_min
double	_vtxz_max
edm::InputTag	BeamHaloSummaryTag
edm::Handle< reco::BeamSpot >	beamSpot_h
math::XYZPoint	bspot
JetMETDQMDCSFilter *	DCSFilter
std::vector< std::string >	elecExpr_
edm::Handle< edm::View < reco::GsfElectron > >	electron_h
edm::InputTag	gtTag
edm::InputTag	HBHENoiseFilterResultTag
edm::InputTag	HcalNoiseRBXCollectionTag
MonitorElement *	heleEta
MonitorElement *	heleHoE
MonitorElement *	helePt
double	hfCalibFactor_
std::vector< std::string >	highMETExpr_
std::vector< std::string >	highPtJetExpr_
std::vector< std::string >	HLTPathsJetMBByName_
MonitorElement *	hMET
MonitorElement *	hMET_logx
MonitorElement *	hmetME
MonitorElement *	hMETPhi
MonitorElement *	hMETRate
MonitorElement *	hMETSig
MonitorElement *	hMEx
MonitorElement *	hMExCorrection
MonitorElement *	hMExLS
MonitorElement *	hMEy
MonitorElement *	hMEyCorrection
MonitorElement *	hMEyLS
MonitorElement *	hmuChi2
MonitorElement *	hmuD0
MonitorElement *	hmuEta
MonitorElement *	hmuNhits
MonitorElement *	hMuonCorrectionFlag
MonitorElement *	hmuPt
MonitorElement *	hSumET
MonitorElement *	hSumET_logx
MonitorElement *	hTriggerName_Ele
MonitorElement *	hTriggerName_HighMET
MonitorElement *	hTriggerName_HighPtJet
MonitorElement *	hTriggerName_LowPtJet
MonitorElement *	hTriggerName_MinBias
MonitorElement *	hTriggerName_Muon
MonitorElement *	htrkChi2
MonitorElement *	htrkD0
MonitorElement *	htrkEta
MonitorElement *	htrkNhits
MonitorElement *	htrkPt
edm::InputTag	inputBeamSpotLabel
edm::InputTag	inputElectronLabel
edm::InputTag	inputMuonLabel
edm::InputTag	inputTrackLabel
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	TcCandidatesTag
edm::Handle< edm::ValueMap < reco::MuonMETCorrectionData > >	tcMet_ValueMap_Handle
edm::ParameterSet	theCleaningParameters
edm::InputTag	theJetCollectionLabel
edm::InputTag	theMETCollectionLabel
edm::InputTag	thePfJetCollectionLabel
edm::Handle< edm::View < reco::Track > >	track_h
edm::InputTag	vertexTag

Detailed Description

DQM monitoring source for MET (Mu corrected/TcMET)

Date:

2012/05/20 13:11:46

Revision:

1.26

Author

A.Apresyan - Caltech

Definition at line 54 of file METAnalyzer.h.

Constructor & Destructor Documentation

METAnalyzer::METAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 35 of file METAnalyzer.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_ ) );
  _HighPtJetEventFlag = new GenericTriggerEventFlag( highptjetparms );
  _LowPtJetEventFlag  = new GenericTriggerEventFlag( lowptjetparms  );
  _MinBiasEventFlag   = new GenericTriggerEventFlag( minbiasparms   );
  _HighMETEventFlag   = new GenericTriggerEventFlag( highmetparms   );
  //  _LowMETEventFlag    = new GenericTriggerEventFlag( lowmetparms    );
  _EleEventFlag       = new GenericTriggerEventFlag( eleparms       );
  _MuonEventFlag      = new GenericTriggerEventFlag( muonparms      );

  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");

}

METAnalyzer::~METAnalyzer ( )

virtual

Destructor.

Definition at line 67 of file METAnalyzer.cc.

                          { 

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

}

Member Function Documentation

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

Get the analysis.

Definition at line 423 of file METAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), FamosSequences_cff::caloJets, gather_cfg::cout, reco::BeamHaloSummary::CSCLooseHaloId(), reco::BeamHaloSummary::CSCTightHaloId(), cmsPerfPublish::DirName, reco::BeamHaloSummary::EcalLooseHaloId(), reco::BeamHaloSummary::EcalTightHaloId(), edm::Event::getByLabel(), reco::BeamHaloSummary::GlobalLooseHaloId(), reco::BeamHaloSummary::GlobalTightHaloId(), reco::BeamHaloSummary::HcalLooseHaloId(), reco::BeamHaloSummary::HcalTightHaloId(), i, edm::HandleBase::isValid(), bTagSequences_cff::jetID, LogDebug, LogTrace, CaloMET_cfi::met, 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 << "METAnalyzer analyze" << std::endl;

  std::string DirName = _FolderName+_source;

  LogTrace(metname)<<"[METAnalyzer] Analyze MET";

  hmetME->Fill(2);

  // ==========================================================  
  // 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("MetAnalyzer") << "TriggerResults::HLT not found, "
      "automatically select events"; 

    // TriggerResults object not found. Look at all events.    
    _trig_JetMB=1;
  }

  // ==========================================================
  // MET information
  
  // **** Get the MET container  
  edm::Handle<reco::METCollection> metcoll;
  iEvent.getByLabel(theMETCollectionLabel, metcoll);
  
  if(!metcoll.isValid()) {
    std::cout<<"Unable to find MET results for MET collection "<<theMETCollectionLabel<<std::endl;
    return;
  }

  const METCollection *metcol = metcoll.product();
  const MET *met;
  met = &(metcol->front());
    
  LogTrace(metname)<<"[METAnalyzer] Call to the MET analyzer";

  // ==========================================================
  // TCMET 

  if (theMETCollectionLabel.label() == "tcMet" ) {

    iEvent.getByLabel(inputMuonLabel, muon_h);
    iEvent.getByLabel(inputTrackLabel, track_h);
    iEvent.getByLabel(inputElectronLabel, electron_h);
    iEvent.getByLabel(inputBeamSpotLabel, beamSpot_h);
    iEvent.getByLabel("muonTCMETValueMapProducer" , "muCorrData", tcMet_ValueMap_Handle);

    if(!muon_h.isValid())     edm::LogInfo("OutputInfo") << "falied to retrieve muon data require by MET Task";
    if(!track_h.isValid())    edm::LogInfo("OutputInfo") << "falied to retrieve track data require by MET Task";
    if(!electron_h.isValid()) edm::LogInfo("OutputInfo") << "falied to retrieve electron data require by MET Task";
    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<HcalNoiseRBXCollection> HRBXCollection;
  iEvent.getByLabel(HcalNoiseRBXCollectionTag,HRBXCollection);
  if (!HRBXCollection.isValid()) {
    LogDebug("") << "METAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
    if (_verbose) std::cout << "METAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
  }
  

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


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

  // ==========================================================
  // MET sanity check

  //   if (_source=="MET") validateMET(*met, tcCandidates);
  
  // ==========================================================
  // JetID 

  if (_verbose) std::cout << "JetID starts" << std::endl;
  
  //
  // --- Minimal cuts
  //
  bool bJetIDMinimal=true;
  for (reco::CaloJetCollection::const_iterator cal = caloJets->begin(); 
       cal!=caloJets->end(); ++cal){
    jetID->calculate(iEvent, *cal);
    if (cal->pt()>10.){
      if (fabs(cal->eta())<=2.6 && 
      cal->emEnergyFraction()<=0.01) bJetIDMinimal=false;
    }
  }

  //
  // --- Loose cuts, not  specific for now!
  //
  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; 
      //if (jetID->restrictedEMF()<0.01) 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<BeamHaloSummary> TheBeamHaloSummary ;
  iEvent.getByLabel(BeamHaloSummaryTag, TheBeamHaloSummary) ;

  if (!TheBeamHaloSummary.isValid()) {
    std::cout << "BeamHaloSummary doesn't exist" << std::endl;
  }

  bool bBeamHaloIDTightPass = true;
  bool bBeamHaloIDLoosePass = true;

  if(!TheBeamHaloSummary.isValid()) {

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

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

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

  }

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

  bool bPrimaryVertex = true;
  if(_doPVCheck){
    bPrimaryVertex = false;
    Handle<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() ){
      VertexCollection vertexCollection = *(vertexHandle.product());
      int vertex_number     = vertexCollection.size();
      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;
      }
    }
  }
  // ==========================================================

  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       = true;

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

void METAnalyzer::beginJob ( DQMStore *  dbe )

virtual

Inizialize parameters for histo binning.

Implements METAnalyzerBase.

Definition at line 79 of file METAnalyzer.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 = "METAnalyzer";

  // 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");
  }

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

  if (theMETCollectionLabel.label() == "tcMet" ) {
    inputTrackLabel         = parameters.getParameter<edm::InputTag>("InputTrackLabel");    
    inputMuonLabel          = parameters.getParameter<edm::InputTag>("InputMuonLabel");
    inputElectronLabel      = parameters.getParameter<edm::InputTag>("InputElectronLabel");
    inputBeamSpotLabel      = parameters.getParameter<edm::InputTag>("InputBeamSpotLabel");
  }

  // Other data collections
  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");
  _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

  _FolderName              = parameters.getUntrackedParameter<std::string>("FolderName");

  _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)<<"[METAnalyzer] Parameters initialization";
  std::string DirName = _FolderName+_source;
  dbe->setCurrentFolder(DirName);

  hmetME = dbe->book1D("metReco", "metReco", 4, 1, 5);
  hmetME->setBinLabel(2,"MET",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 METAnalyzer::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

Initialize run-based parameters.

Definition at line 311 of file METAnalyzer.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 METAnalyzer::bookMESet

(

std::string

DirName

)

Definition at line 203 of file METAnalyzer.cc.

References gather_cfg::cout.

{

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

  bookMonitorElement(DirName,bLumiSecPlot);

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

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

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

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

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

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

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

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

Definition at line 251 of file METAnalyzer.cc.

References gather_cfg::cout.

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

  _dbe->setCurrentFolder(DirName);


  hMEx        = _dbe->book1D("METTask_MEx",        "METTask_MEx",        200, -500,  500);
  hMEy        = _dbe->book1D("METTask_MEy",        "METTask_MEy",        200, -500,  500);
  hMET        = _dbe->book1D("METTask_MET",        "METTask_MET",        200,    0, 1000); 
  hSumET      = _dbe->book1D("METTask_SumET",      "METTask_SumET",      400,    0, 4000); 
  hMETSig     = _dbe->book1D("METTask_METSig",     "METTask_METSig",      51,    0,   51);
  hMETPhi     = _dbe->book1D("METTask_METPhi",     "METTask_METPhi",      60, -3.2,  3.2); 
  hMET_logx   = _dbe->book1D("METTask_MET_logx",   "METTask_MET_logx",    40,   -1,    7);
  hSumET_logx = _dbe->book1D("METTask_SumET_logx", "METTask_SumET_logx",  40,   -1,    7);

  hMEx       ->setAxisTitle("MEx [GeV]",        1);
  hMEy       ->setAxisTitle("MEy [GeV]",        1);
  hMET       ->setAxisTitle("MET [GeV]",        1);
  hSumET     ->setAxisTitle("SumET [GeV]",      1);
  hMETSig    ->setAxisTitle("CaloMETSig",       1);
  hMETPhi    ->setAxisTitle("METPhi [rad]",     1);
  hMET_logx  ->setAxisTitle("log(MET) [GeV]",   1);
  hSumET_logx->setAxisTitle("log(SumET) [GeV]", 1);


  if (_allhist){
    if (bLumiSecPlot){
      hMExLS = _dbe->book2D("METTask_MEx_LS","METTask_MEx_LS",200,-200,200,50,0.,500.);
      hMExLS->setAxisTitle("MEx [GeV]",1);
      hMExLS->setAxisTitle("Lumi Section",2);
      hMEyLS = _dbe->book2D("METTask_MEy_LS","METTask_MEy_LS",200,-200,200,50,0.,500.);
      hMEyLS->setAxisTitle("MEy [GeV]",1);
      hMEyLS->setAxisTitle("Lumi Section",2);
    }
  }

  if (theMETCollectionLabel.label() == "tcMet" ) {
    htrkPt    = _dbe->book1D("METTask_trackPt", "METTask_trackPt", 50, 0, 500);
    htrkEta   = _dbe->book1D("METTask_trackEta", "METTask_trackEta", 60, -3.0, 3.0);
    htrkNhits = _dbe->book1D("METTask_trackNhits", "METTask_trackNhits", 50, 0, 50);
    htrkChi2  = _dbe->book1D("METTask_trackNormalizedChi2", "METTask_trackNormalizedChi2", 20, 0, 20);
    htrkD0    = _dbe->book1D("METTask_trackD0", "METTask_trackd0", 50, -1, 1);
    helePt    = _dbe->book1D("METTask_electronPt", "METTask_electronPt", 50, 0, 500);
    heleEta   = _dbe->book1D("METTask_electronEta", "METTask_electronEta", 60, -3.0, 3.0);
    heleHoE   = _dbe->book1D("METTask_electronHoverE", "METTask_electronHoverE", 25, 0, 0.5);
    hmuPt     = _dbe->book1D("METTask_muonPt", "METTask_muonPt", 50, 0, 500);
    hmuEta    = _dbe->book1D("METTask_muonEta", "METTask_muonEta", 60, -3.0, 3.0);
    hmuNhits  = _dbe->book1D("METTask_muonNhits", "METTask_muonNhits", 50, 0, 50);
    hmuChi2   = _dbe->book1D("METTask_muonNormalizedChi2", "METTask_muonNormalizedChi2", 20, 0, 20);
    hmuD0     = _dbe->book1D("METTask_muonD0", "METTask_muonD0", 50, -1, 1);
  }

  hMExCorrection       = _dbe->book1D("METTask_MExCorrection", "METTask_MExCorrection", 100, -500.0,500.0);
  hMEyCorrection       = _dbe->book1D("METTask_MEyCorrection", "METTask_MEyCorrection", 100, -500.0,500.0);
  hMuonCorrectionFlag  = _dbe->book1D("METTask_CorrectionFlag","METTask_CorrectionFlag", 5, -0.5, 4.5);

}

void METAnalyzer::endJob

(

void

)

Finish up a job.

Definition at line 195 of file METAnalyzer.cc.

References bTagSequences_cff::jetID.

                         {

  delete jetID;
  delete DCSFilter;

}

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

Finish up a run.

Definition at line 339 of file METAnalyzer.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 = _FolderName+_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 METAnalyzer::fillMESet	(	const edm::Event &	iEvent,
		std::string	DirName,
		const reco::MET &	met
	)

Definition at line 859 of file METAnalyzer.cc.

{

  _dbe->setCurrentFolder(DirName);

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

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

void METAnalyzer::fillMonitorElement	(	const edm::Event &	iEvent,
		std::string	DirName,
		std::string	TriggerTypeName,
		const reco::MET &	met,
		bool	bLumiSecPlot
	)

Definition at line 887 of file METAnalyzer.cc.

References reco::MuonMETCorrectionData::corrX(), reco::MuonMETCorrectionData::corrY(), gather_cfg::cout, 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 (met.pt()<_highMETThreshold) return;
  }
  //  else if (TriggerTypeName=="LowMET") {
  //    if (met.pt()<_lowMETThreshold) return;
  //  }
  else if (TriggerTypeName=="Ele") {
    if (!selectWElectronEvent(iEvent)) return;
  }
  else if (TriggerTypeName=="Muon") {
    if (!selectWMuonEvent(iEvent)) return;
  }
  
// Reconstructed MET Information
  double SumET  = met.sumEt();
  double METSig = met.mEtSig();
  //double Ez     = met.e_longitudinal();
  double MET    = met.pt();
  double MEx    = met.px();
  double MEy    = met.py();
  double METPhi = met.phi();

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

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

  if (_verbose) std::cout << "_etThreshold = " << _etThreshold << std::endl;
  if (SumET>_etThreshold){
    
    hMEx    = _dbe->get(DirName+"/"+"METTask_MEx");     if (hMEx           && hMEx->getRootObject())     hMEx          ->Fill(MEx);
    hMEy    = _dbe->get(DirName+"/"+"METTask_MEy");     if (hMEy           && hMEy->getRootObject())     hMEy          ->Fill(MEy);
    hMET    = _dbe->get(DirName+"/"+"METTask_MET");     if (hMET           && hMET->getRootObject())     hMET          ->Fill(MET);
    hMETPhi = _dbe->get(DirName+"/"+"METTask_METPhi");  if (hMETPhi        && hMETPhi->getRootObject())  hMETPhi       ->Fill(METPhi);
    hSumET  = _dbe->get(DirName+"/"+"METTask_SumET");   if (hSumET         && hSumET->getRootObject())   hSumET        ->Fill(SumET);
    hMETSig = _dbe->get(DirName+"/"+"METTask_METSig");  if (hMETSig        && hMETSig->getRootObject())  hMETSig       ->Fill(METSig);
    //hEz     = _dbe->get(DirName+"/"+"METTask_Ez");      if (hEz            && hEz->getRootObject())      hEz           ->Fill(Ez);

    hMET_logx   = _dbe->get(DirName+"/"+"METTask_MET_logx");    if (hMET_logx      && hMET_logx->getRootObject())    hMET_logx->Fill(log10(MET));
    hSumET_logx = _dbe->get(DirName+"/"+"METTask_SumET_logx");  if (hSumET_logx    && hSumET_logx->getRootObject())  hSumET_logx->Fill(log10(SumET));

    //hMETIonFeedbck = _dbe->get(DirName+"/"+"METTask_METIonFeedbck");  if (hMETIonFeedbck && hMETIonFeedbck->getRootObject())  hMETIonFeedbck->Fill(MET);
    //hMETHPDNoise   = _dbe->get(DirName+"/"+"METTask_METHPDNoise");    if (hMETHPDNoise   && hMETHPDNoise->getRootObject())    hMETHPDNoise->Fill(MET);
       
    if (_allhist){
      if (bLumiSecPlot){
    hMExLS = _dbe->get(DirName+"/"+"METTask_MExLS"); if (hMExLS  &&  hMExLS->getRootObject())   hMExLS->Fill(MEx,myLuminosityBlock);
    hMEyLS = _dbe->get(DirName+"/"+"METTask_MEyLS"); if (hMEyLS  &&  hMEyLS->getRootObject())   hMEyLS->Fill(MEy,myLuminosityBlock);
      }
    } // _allhist

    if (theMETCollectionLabel.label() == "tcMet" ) {
    
      if(track_h.isValid()) {
    for( edm::View<reco::Track>::const_iterator trkit = track_h->begin(); trkit != track_h->end(); trkit++ ) {
      htrkPt    = _dbe->get(DirName+"/"+"METTask_trackPt");     if (htrkPt    && htrkPt->getRootObject())     htrkPt->Fill( trkit->pt() );
      htrkEta   = _dbe->get(DirName+"/"+"METTask_trackEta");    if (htrkEta   && htrkEta->getRootObject())    htrkEta->Fill( trkit->eta() );
      htrkNhits = _dbe->get(DirName+"/"+"METTask_trackNhits");  if (htrkNhits && htrkNhits->getRootObject())  htrkNhits->Fill( trkit->numberOfValidHits() );
      htrkChi2  = _dbe->get(DirName+"/"+"METTask_trackNormalizedChi2");   if (htrkChi2  && htrkChi2->getRootObject())   htrkChi2->Fill( trkit->chi2() / trkit->ndof() );
      double d0 = -1 * trkit->dxy( bspot );
      htrkD0    = _dbe->get(DirName+"/"+"METTask_trackD0");     if (htrkD0 && htrkD0->getRootObject())        htrkD0->Fill( d0 );
    }
      }
      
      if(electron_h.isValid()) {
    for( edm::View<reco::GsfElectron>::const_iterator eleit = electron_h->begin(); eleit != electron_h->end(); eleit++ ) {
      helePt  = _dbe->get(DirName+"/"+"METTask_electronPt");   if (helePt  && helePt->getRootObject())   helePt->Fill( eleit->p4().pt() );  
      heleEta = _dbe->get(DirName+"/"+"METTask_electronEta");  if (heleEta && heleEta->getRootObject())  heleEta->Fill( eleit->p4().eta() );
      heleHoE = _dbe->get(DirName+"/"+"METTask_electronHoverE");  if (heleHoE && heleHoE->getRootObject())  heleHoE->Fill( eleit->hadronicOverEm() );
    }
      }
      
      if(muon_h.isValid()) {
    for( reco::MuonCollection::const_iterator muonit = muon_h->begin(); muonit != muon_h->end(); muonit++ ) {      
      const reco::TrackRef siTrack = muonit->innerTrack();      
      hmuPt    = _dbe->get(DirName+"/"+"METTask_muonPt");     if (hmuPt    && hmuPt->getRootObject())  hmuPt   ->Fill( muonit->p4().pt() );
      hmuEta   = _dbe->get(DirName+"/"+"METTask_muonEta");    if (hmuEta   && hmuEta->getRootObject())  hmuEta  ->Fill( muonit->p4().eta() );
      hmuNhits = _dbe->get(DirName+"/"+"METTask_muonNhits");  if (hmuNhits && hmuNhits->getRootObject())  hmuNhits->Fill( siTrack.isNonnull() ? siTrack->numberOfValidHits() : -999 );
      hmuChi2  = _dbe->get(DirName+"/"+"METTask_muonNormalizedChi2");   if (hmuChi2  && hmuChi2->getRootObject())  hmuChi2 ->Fill( siTrack.isNonnull() ? siTrack->chi2()/siTrack->ndof() : -999 );
      double d0 = siTrack.isNonnull() ? -1 * siTrack->dxy( bspot) : -999;
      hmuD0    = _dbe->get(DirName+"/"+"METTask_muonD0");     if (hmuD0    && hmuD0->getRootObject())  hmuD0->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 = (*tcMet_ValueMap_Handle)[muref];
      hMExCorrection      = _dbe->get(DirName+"/"+"METTask_MExCorrection");       if (hMExCorrection      && hMExCorrection->getRootObject())       hMExCorrection-> Fill(muCorrData.corrY());
      hMEyCorrection      = _dbe->get(DirName+"/"+"METTask_MEyCorrection");       if (hMEyCorrection      && hMEyCorrection->getRootObject())       hMEyCorrection-> Fill(muCorrData.corrX());
      hMuonCorrectionFlag = _dbe->get(DirName+"/"+"METTask_CorrectionFlag");  if (hMuonCorrectionFlag && hMuonCorrectionFlag->getRootObject())  hMuonCorrectionFlag-> Fill(muCorrData.type());
    }
      }
    }

  } // et threshold cut

}

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

Definition at line 394 of file METAnalyzer.cc.

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

{

  _dbe->setCurrentFolder(DirName);
  MonitorElement *meMET = _dbe->get(DirName+"/"+"METTask_MET");

  TH1F* tMET;
  TH1F* tMETRate;

  if ( meMET )
    if ( meMET->getRootObject() ) {
      tMET     = meMET->getTH1F();
      
      // Integral plot & convert number of events to rate (hz)
      tMETRate = (TH1F*) tMET->Clone("METTask_METRate");
      for (int i = tMETRate->GetNbinsX()-1; i>=0; i--){
    tMETRate->SetBinContent(i+1,tMETRate->GetBinContent(i+2)+tMET->GetBinContent(i+1));
      }
      for (int i = 0; i<tMETRate->GetNbinsX(); i++){
    tMETRate->SetBinContent(i+1,tMETRate->GetBinContent(i+1)/double(totltime));
      }      

      tMETRate->SetName("METTask_METRate");
      tMETRate->SetTitle("METTask_METRate");
      hMETRate      = _dbe->book1D("METTask_METRate",tMETRate);
    }
}

bool METAnalyzer::selectHighPtJetEvent

(

const edm::Event &

iEvent

)

Definition at line 1001 of file METAnalyzer.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("") << "METAnalyzer: Could not find jet product" << std::endl;
    if (_verbose) std::cout << "METAnalyzer: 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 METAnalyzer::selectLowPtJetEvent

(

const edm::Event &

iEvent

)

Definition at line 1023 of file METAnalyzer.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("") << "METAnalyzer: Could not find jet product" << std::endl;
    if (_verbose) std::cout << "METAnalyzer: 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 METAnalyzer::selectWElectronEvent

(

const edm::Event &

iEvent

)

Definition at line 1047 of file METAnalyzer.cc.

                                                            {

  bool return_value=true;

  /*
    W-electron event selection comes here
   */

  return return_value;

}

bool METAnalyzer::selectWMuonEvent

(

const edm::Event &

iEvent

)

Definition at line 1060 of file METAnalyzer.cc.

                                                        {

  bool return_value=true;

  /*
    W-muon event selection comes here
   */

  return return_value;

}

void METAnalyzer::setSource ( std::string  source )

inline

Definition at line 93 of file METAnalyzer.h.

References _source, and LaserTracksInput_cfi::source.

                                   {
    _source = source;
  }

Member Data Documentation

bool METAnalyzer::_allhist

private

Definition at line 195 of file METAnalyzer.h.

bool METAnalyzer::_allSelection

private

Definition at line 196 of file METAnalyzer.h.

bool METAnalyzer::_cleanupSelection

private

Definition at line 197 of file METAnalyzer.h.

DQMStore* METAnalyzer::_dbe

private

Definition at line 212 of file METAnalyzer.h.

bool METAnalyzer::_doHLTPhysicsOn

private

Definition at line 154 of file METAnalyzer.h.

bool METAnalyzer::_doPVCheck

private

Definition at line 153 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_EleEventFlag

private

Definition at line 135 of file METAnalyzer.h.

double METAnalyzer::_etThreshold

private

Definition at line 182 of file METAnalyzer.h.

std::string METAnalyzer::_FolderName

private

Definition at line 109 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::_FolderNames

private

Definition at line 200 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_HighMETEventFlag

private

Definition at line 133 of file METAnalyzer.h.

double METAnalyzer::_highMETThreshold

private

Definition at line 178 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_HighPtJetEventFlag

private

Definition at line 130 of file METAnalyzer.h.

double METAnalyzer::_highPtJetThreshold

private

Definition at line 176 of file METAnalyzer.h.

std::string METAnalyzer::_hlt_PhysDec

private

Definition at line 147 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_LowPtJetEventFlag

private

Definition at line 131 of file METAnalyzer.h.

double METAnalyzer::_lowPtJetThreshold

private

Definition at line 177 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_MinBiasEventFlag

private

Definition at line 132 of file METAnalyzer.h.

GenericTriggerEventFlag* METAnalyzer::_MuonEventFlag

private

Definition at line 136 of file METAnalyzer.h.

int METAnalyzer::_nvtx_min

private

Definition at line 159 of file METAnalyzer.h.

int METAnalyzer::_nvtxtrks_min

private

Definition at line 160 of file METAnalyzer.h.

std::string METAnalyzer::_source

private

Definition at line 107 of file METAnalyzer.h.

Referenced by setSource().

std::vector<unsigned > METAnalyzer::_techTrigsAND

private

Definition at line 149 of file METAnalyzer.h.

std::vector<unsigned > METAnalyzer::_techTrigsNOT

private

Definition at line 151 of file METAnalyzer.h.

std::vector<unsigned > METAnalyzer::_techTrigsOR

private

Definition at line 150 of file METAnalyzer.h.

bool METAnalyzer::_tightBHFiltering

private

Definition at line 156 of file METAnalyzer.h.

int METAnalyzer::_tightJetIDFiltering

private

Definition at line 157 of file METAnalyzer.h.

int METAnalyzer::_trig_Ele

private

Definition at line 171 of file METAnalyzer.h.

int METAnalyzer::_trig_HighMET

private

Definition at line 169 of file METAnalyzer.h.

int METAnalyzer::_trig_HighPtJet

private

Definition at line 166 of file METAnalyzer.h.

int METAnalyzer::_trig_JetMB

private

Definition at line 165 of file METAnalyzer.h.

int METAnalyzer::_trig_LowPtJet

private

Definition at line 167 of file METAnalyzer.h.

int METAnalyzer::_trig_MinBias

private

Definition at line 168 of file METAnalyzer.h.

int METAnalyzer::_trig_Muon

private

Definition at line 172 of file METAnalyzer.h.

int METAnalyzer::_trig_PhysDec

private

Definition at line 173 of file METAnalyzer.h.

int METAnalyzer::_verbose

private

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

private

Definition at line 162 of file METAnalyzer.h.

int METAnalyzer::_vtxndof_min

private

Definition at line 161 of file METAnalyzer.h.

double METAnalyzer::_vtxz_max

private

Definition at line 163 of file METAnalyzer.h.

edm::InputTag METAnalyzer::BeamHaloSummaryTag

private

Definition at line 116 of file METAnalyzer.h.

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

private

Definition at line 209 of file METAnalyzer.h.

math::XYZPoint METAnalyzer::bspot

private

Definition at line 203 of file METAnalyzer.h.

JetMETDQMDCSFilter* METAnalyzer::DCSFilter

private

Definition at line 192 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::elecExpr_

private

Definition at line 143 of file METAnalyzer.h.

edm::Handle< edm::View<reco::GsfElectron > > METAnalyzer::electron_h

private

Definition at line 208 of file METAnalyzer.h.

int METAnalyzer::evtCounter

Definition at line 97 of file METAnalyzer.h.

edm::InputTag METAnalyzer::gtTag

private

Definition at line 119 of file METAnalyzer.h.

edm::InputTag METAnalyzer::HBHENoiseFilterResultTag

private

Definition at line 117 of file METAnalyzer.h.

edm::InputTag METAnalyzer::HcalNoiseRBXCollectionTag

private

Definition at line 112 of file METAnalyzer.h.

MonitorElement* METAnalyzer::heleEta

private

Definition at line 253 of file METAnalyzer.h.

MonitorElement* METAnalyzer::heleHoE

private

Definition at line 254 of file METAnalyzer.h.

MonitorElement* METAnalyzer::helePt

private

Definition at line 252 of file METAnalyzer.h.

double METAnalyzer::hfCalibFactor_

private

Definition at line 185 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::highMETExpr_

private

Definition at line 140 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::highPtJetExpr_

private

Definition at line 138 of file METAnalyzer.h.

std::vector<std::string > METAnalyzer::HLTPathsJetMBByName_

private

Definition at line 128 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMET

private

Definition at line 233 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMET_logx

private

Definition at line 237 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmetME

private

Definition at line 227 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMETPhi

private

Definition at line 234 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMETRate

private

Definition at line 225 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMETSig

private

Definition at line 232 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMEx

private

Definition at line 229 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMExCorrection

private

Definition at line 261 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMExLS

private

Definition at line 244 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMEy

private

Definition at line 230 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMEyCorrection

private

Definition at line 262 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMEyLS

private

Definition at line 245 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmuChi2

private

Definition at line 258 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmuD0

private

Definition at line 259 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmuEta

private

Definition at line 256 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmuNhits

private

Definition at line 257 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hMuonCorrectionFlag

private

Definition at line 263 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hmuPt

private

Definition at line 255 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hSumET

private

Definition at line 235 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hSumET_logx

private

Definition at line 238 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_Ele

private

Definition at line 221 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_HighMET

private

Definition at line 219 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_HighPtJet

private

Definition at line 216 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_LowPtJet

private

Definition at line 217 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_MinBias

private

Definition at line 218 of file METAnalyzer.h.

MonitorElement* METAnalyzer::hTriggerName_Muon

private

Definition at line 222 of file METAnalyzer.h.

MonitorElement* METAnalyzer::htrkChi2

private

Definition at line 250 of file METAnalyzer.h.

MonitorElement* METAnalyzer::htrkD0

private

Definition at line 251 of file METAnalyzer.h.

MonitorElement* METAnalyzer::htrkEta

private

Definition at line 248 of file METAnalyzer.h.

MonitorElement* METAnalyzer::htrkNhits

private

Definition at line 249 of file METAnalyzer.h.

MonitorElement* METAnalyzer::htrkPt

private

Definition at line 247 of file METAnalyzer.h.

edm::InputTag METAnalyzer::inputBeamSpotLabel

private

Definition at line 124 of file METAnalyzer.h.

edm::InputTag METAnalyzer::inputElectronLabel

private

Definition at line 123 of file METAnalyzer.h.

edm::InputTag METAnalyzer::inputMuonLabel

private

Definition at line 122 of file METAnalyzer.h.

edm::InputTag METAnalyzer::inputTrackLabel

private

Definition at line 121 of file METAnalyzer.h.

reco::helper::JetIDHelper* METAnalyzer::jetID

private

Definition at line 188 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::lowPtJetExpr_

private

Definition at line 139 of file METAnalyzer.h.

std::string METAnalyzer::metname

private

Definition at line 106 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::minbiasExpr_

private

Definition at line 144 of file METAnalyzer.h.

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

private

Definition at line 206 of file METAnalyzer.h.

std::vector<std::string> METAnalyzer::muonExpr_

private

Definition at line 142 of file METAnalyzer.h.

edm::ParameterSet METAnalyzer::parameters

private

Definition at line 102 of file METAnalyzer.h.

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

edm::InputTag METAnalyzer::TcCandidatesTag

private

Definition at line 115 of file METAnalyzer.h.

edm::Handle< edm::ValueMap<reco::MuonMETCorrectionData> > METAnalyzer::tcMet_ValueMap_Handle

private

Definition at line 205 of file METAnalyzer.h.

edm::ParameterSet METAnalyzer::theCleaningParameters

private

Definition at line 146 of file METAnalyzer.h.

edm::InputTag METAnalyzer::theJetCollectionLabel

private

Definition at line 113 of file METAnalyzer.h.

edm::InputTag METAnalyzer::theMETCollectionLabel

private

Definition at line 111 of file METAnalyzer.h.

edm::InputTag METAnalyzer::thePfJetCollectionLabel

private

Definition at line 114 of file METAnalyzer.h.

edm::Handle< edm::View<reco::Track> > METAnalyzer::track_h

private

Definition at line 207 of file METAnalyzer.h.

edm::InputTag METAnalyzer::vertexTag

private

Definition at line 118 of file METAnalyzer.h.