METTester Class Reference

#include <METTester.h>

Inheritance diagram for METTester:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	FillMETRes ()
	METTester (const edm::ParameterSet &)
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
edm::EDGetTokenT < reco::CaloMETCollection >	caloMETsToken_
edm::EDGetTokenT < reco::GenMETCollection >	genMETsCaloToken_
edm::EDGetTokenT < reco::GenMETCollection >	genMETsToken_
edm::EDGetTokenT < reco::GenMETCollection >	genMETsTrueToken_
edm::InputTag	inputCaloMETLabel_
edm::InputTag	inputMETLabel_
bool	isCaloMET
bool	isGenMET
bool	isMiniAODMET
bool	isPFMET
MonitorElement *	mCaloEmEtFraction
MonitorElement *	mCaloEmEtInEB
MonitorElement *	mCaloEmEtInEE
MonitorElement *	mCaloEmEtInHF
MonitorElement *	mCaloEtFractionHadronic
MonitorElement *	mCaloHadEtInEB
MonitorElement *	mCaloHadEtInEE
MonitorElement *	mCaloHadEtInHB
MonitorElement *	mCaloHadEtInHE
MonitorElement *	mCaloHadEtInHF
MonitorElement *	mCaloHadEtInHO
MonitorElement *	mCaloMaxEtInEmTowers
MonitorElement *	mCaloMaxEtInHadTowers
MonitorElement *	mCaloSETInmHF
MonitorElement *	mCaloSETInpHF
MonitorElement *	mChargedEMEtFraction
MonitorElement *	mChargedHadEtFraction
std::map< std::string, MonitorElement * >	me
std::string	METType_
edm::InputTag	mInputCollection_
MonitorElement *	mInvisibleEtFraction
MonitorElement *	mMET
MonitorElement *	mMET_Nvtx
MonitorElement *	mMETDeltaPhi_GenMETCalo
MonitorElement *	mMETDeltaPhi_GenMETTrue
MonitorElement *	mMETDifference_GenMETCalo
MonitorElement *	mMETDifference_GenMETTrue
MonitorElement *	mMETDifference_GenMETTrue_MET0to20
MonitorElement *	mMETDifference_GenMETTrue_MET100to150
MonitorElement *	mMETDifference_GenMETTrue_MET150to200
MonitorElement *	mMETDifference_GenMETTrue_MET200to300
MonitorElement *	mMETDifference_GenMETTrue_MET20to40
MonitorElement *	mMETDifference_GenMETTrue_MET300to400
MonitorElement *	mMETDifference_GenMETTrue_MET400to500
MonitorElement *	mMETDifference_GenMETTrue_MET40to60
MonitorElement *	mMETDifference_GenMETTrue_MET60to80
MonitorElement *	mMETDifference_GenMETTrue_MET80to100
MonitorElement *	mMETDifference_GenMETTrue_METResolution
MonitorElement *	mMETFine
MonitorElement *	mMETPhi
MonitorElement *	mMETSig
MonitorElement *	mMEx
MonitorElement *	mMEy
MonitorElement *	mMuonEtFraction
MonitorElement *	mNeutralEMEtFraction
MonitorElement *	mNeutralHadEtFraction
MonitorElement *	mNvertex
MonitorElement *	mPFchargedHadronEt
MonitorElement *	mPFchargedHadronEtFraction
MonitorElement *	mPFelectronEt
MonitorElement *	mPFelectronEtFraction
MonitorElement *	mPFHFEMEt
MonitorElement *	mPFHFEMEtFraction
MonitorElement *	mPFHFHadronEt
MonitorElement *	mPFHFHadronEtFraction
MonitorElement *	mPFmuonEt
MonitorElement *	mPFmuonEtFraction
MonitorElement *	mPFneutralHadronEt
MonitorElement *	mPFneutralHadronEtFraction
MonitorElement *	mPFphotonEt
MonitorElement *	mPFphotonEtFraction
MonitorElement *	mSumET
edm::EDGetTokenT < pat::METCollection >	patMETToken_
edm::EDGetTokenT < reco::PFMETCollection >	pfMETsToken_
edm::EDGetTokenT< std::vector < reco::Vertex > >	pvToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 53 of file METTester.h.

Constructor & Destructor Documentation

METTester::METTester ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 32 of file METTester.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

{

  inputMETLabel_           =iConfig.getParameter<edm::InputTag>("InputMETLabel");
  METType_                 =iConfig.getUntrackedParameter<std::string>("METType");

  std::string inputMETCollectionLabel(inputMETLabel_.label());

  isCaloMET = (std::string("calo")==METType_);
  isPFMET   = (std::string("pf")  ==METType_);
  isGenMET  = (std::string("gen") ==METType_);
  isMiniAODMET  = (std::string("miniaod") ==METType_);

  pvToken_ = consumes<std::vector<reco::Vertex> >(iConfig.getParameter<edm::InputTag>("PrimaryVertices"));
  if (isCaloMET)  caloMETsToken_ = consumes<reco::CaloMETCollection> (inputMETLabel_);
  if (isPFMET)    pfMETsToken_ = consumes<reco::PFMETCollection> (inputMETLabel_); 
  if(isMiniAODMET) patMETToken_ = consumes<pat::METCollection> (inputMETLabel_);
  if (isGenMET)   genMETsToken_ = consumes<reco::GenMETCollection> (inputMETLabel_); 
  if(!isMiniAODMET){
    genMETsTrueToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetTrue"));
    genMETsCaloToken_ = consumes<reco::GenMETCollection> (edm::InputTag("genMetCalo"));
  }
  //Events variables
  mNvertex               = 0;

  //Common variables
  mMEx                         = 0;
  mMEy                         = 0;
  mMETSig                      = 0;
  mMET                         = 0;
  mMETFine                     = 0;
  mMET_Nvtx                    = 0;
  mMETPhi                      = 0;
  mSumET                       = 0;
  mMETDifference_GenMETTrue    = 0;
  mMETDeltaPhi_GenMETTrue      = 0;
  mMETDifference_GenMETCalo    = 0;
  mMETDeltaPhi_GenMETCalo      = 0;

  //CaloMET variables
  mCaloMaxEtInEmTowers             = 0;
  mCaloMaxEtInHadTowers            = 0;
  mCaloEtFractionHadronic          = 0;
  mCaloEmEtFraction                = 0;
  mCaloHadEtInHB                   = 0;
  mCaloHadEtInHO                   = 0;
  mCaloHadEtInHE                   = 0;
  mCaloHadEtInHF                   = 0;
  mCaloEmEtInHF                    = 0;
  mCaloSETInpHF                    = 0;
  mCaloSETInmHF                    = 0;
  mCaloEmEtInEE                    = 0;
  mCaloEmEtInEB                    = 0;


  //GenMET variables
  mNeutralEMEtFraction=0;
  mNeutralHadEtFraction=0;
  mChargedEMEtFraction=0;
  mChargedHadEtFraction=0;
  mMuonEtFraction=0; 
  mInvisibleEtFraction=0;
  
  //MET variables
  
  //PFMET variables
  mMETDifference_GenMETTrue_MET0to20=0;
  mMETDifference_GenMETTrue_MET20to40=0;
  mMETDifference_GenMETTrue_MET40to60=0;
  mMETDifference_GenMETTrue_MET60to80=0;
  mMETDifference_GenMETTrue_MET80to100=0;
  mMETDifference_GenMETTrue_MET100to150=0;
  mMETDifference_GenMETTrue_MET150to200=0;
  mMETDifference_GenMETTrue_MET200to300=0;
  mMETDifference_GenMETTrue_MET300to400=0;
  mMETDifference_GenMETTrue_MET400to500=0;
  mMETDifference_GenMETTrue_METResolution=0;
  
 
} 

Member Function Documentation

void METTester::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 197 of file METTester.cc.

References assert(), reco::CaloMET::CaloSETInmHF(), reco::CaloMET::CaloSETInpHF(), reco::GenMET::ChargedEMEtFraction(), pat::MET::ChargedEMEtFraction(), reco::GenMET::ChargedHadEtFraction(), pat::MET::ChargedHadEtFraction(), reco::PFMET::chargedHadronEt(), reco::PFMET::chargedHadronEtFraction(), gather_cfg::cout, reco::PFMET::electronEt(), reco::PFMET::electronEtFraction(), reco::CaloMET::emEtFraction(), reco::CaloMET::emEtInEB(), reco::CaloMET::emEtInEE(), reco::CaloMET::emEtInHF(), reco::CaloMET::etFractionHadronic(), genMetCalo_cfi::genMetCalo, genMetTrue_cff::genMetTrue, edm::Event::getByToken(), reco::CaloMET::hadEtInHB(), reco::CaloMET::hadEtInHE(), reco::CaloMET::hadEtInHF(), reco::CaloMET::hadEtInHO(), reco::PFMET::HFEMEt(), reco::PFMET::HFEMEtFraction(), reco::PFMET::HFHadronEt(), reco::PFMET::HFHadronEtFraction(), reco::GenMET::InvisibleEtFraction(), pat::MET::isPFMET(), edm::HandleBase::isValid(), reco::CaloMET::maxEtInEmTowers(), reco::CaloMET::maxEtInHadTowers(), caloMETBenchmarkGeneric_cfi::met, reco::MET::mEtSig(), reco::PFMET::muonEt(), reco::PFMET::muonEtFraction(), reco::GenMET::MuonEtFraction(), pat::MET::MuonEtFraction(), reco::GenMET::NeutralEMEtFraction(), pat::MET::NeutralEMFraction(), reco::GenMET::NeutralHadEtFraction(), pat::MET::NeutralHadEtFraction(), reco::PFMET::neutralHadronEt(), reco::PFMET::neutralHadronEtFraction(), NULL, reco::LeafCandidate::phi(), reco::PFMET::photonEt(), reco::PFMET::photonEtFraction(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::MET::sumEt(), pat::MET::Type6EtFraction(), and pat::MET::Type7EtFraction().

{ //int counter(0);
  edm::Handle<reco::VertexCollection> pvHandle;
  iEvent.getByToken(pvToken_, pvHandle);
   if (! pvHandle.isValid())
  {
    std::cout << __FUNCTION__ << ":" << __LINE__ << ":pvHandle handle not found!" << std::endl;
    assert(false);
  }
  const int nvtx = pvHandle->size();
  mNvertex->Fill(nvtx);
  //Collections for all MET collections

  edm::Handle<CaloMETCollection> caloMETs;
  edm::Handle<PFMETCollection> pfMETs;
  edm::Handle<GenMETCollection> genMETs;
  edm::Handle<pat::METCollection> patMET;

  if (isCaloMET) iEvent.getByToken(caloMETsToken_, caloMETs);
  if (isPFMET)   iEvent.getByToken(pfMETsToken_,   pfMETs);
  if (isGenMET)  iEvent.getByToken(genMETsToken_,  genMETs);
  if (isMiniAODMET)  iEvent.getByToken(patMETToken_,  patMET);
  if ((isCaloMET) and !caloMETs.isValid()) return; 
  if ((isPFMET)   and !pfMETs.isValid())   return;
  if ((isGenMET)  and !genMETs.isValid())  return;
  if ((isMiniAODMET)  and !patMET.isValid())  return;

  reco::MET met;
  if (isCaloMET) { met = caloMETs->front();}
  if (isPFMET)   { met = pfMETs->front()  ;}
  if (isGenMET)  { met = genMETs->front() ;}
  if (isMiniAODMET)  { met = patMET->front() ;}

  const double SumET = met.sumEt();
  const double METSig = met.mEtSig();
  const double MET = met.pt();
  const double MEx = met.px();
  const double MEy = met.py();
  const double METPhi = met.phi();
  mMEx->Fill(MEx);
  mMEy->Fill(MEy);
  mMET->Fill(MET);
  mMETFine->Fill(MET);
  mMET_Nvtx->Fill((double)nvtx, MET);
  mMETPhi->Fill(METPhi);
  mSumET->Fill(SumET);
  mMETSig->Fill(METSig);

  // Get Generated MET for Resolution plots
  const reco::GenMET *genMetTrue=NULL;
  bool isvalidgenmet=false;

  if(!isMiniAODMET){
    edm::Handle<GenMETCollection> genTrue;
    iEvent.getByToken(genMETsTrueToken_, genTrue);
    if (genTrue.isValid()) {
      isvalidgenmet=true;
      const GenMETCollection *genmetcol = genTrue.product();
      genMetTrue = &(genmetcol->front());
    }
  }else{
    genMetTrue=patMET->front().genMET();
    isvalidgenmet=true;
  }
  if(isvalidgenmet){
    double genMET = genMetTrue->pt();
    double genMETPhi = genMetTrue->phi();

    mMETDifference_GenMETTrue->Fill( MET - genMET );
    mMETDeltaPhi_GenMETTrue->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );

    if(!isGenMET){
      //pfMET resolution in pfMET bins : Sam, Feb, 2012
      if (MET > 0 && MET < 20) mMETDifference_GenMETTrue_MET0to20->Fill( MET - genMET );
      else if (MET > 20 && MET < 40) mMETDifference_GenMETTrue_MET20to40->Fill( MET - genMET );
      else if (MET > 40 && MET < 60) mMETDifference_GenMETTrue_MET40to60->Fill( MET - genMET );
      else if (MET > 60 && MET < 80) mMETDifference_GenMETTrue_MET60to80->Fill( MET - genMET );
      else if (MET > 80 && MET <100) mMETDifference_GenMETTrue_MET80to100->Fill( MET - genMET );
      else if (MET >100 && MET <150) mMETDifference_GenMETTrue_MET100to150->Fill( MET - genMET );
      else if (MET >150 && MET <200) mMETDifference_GenMETTrue_MET150to200->Fill( MET - genMET );
      else if (MET >200 && MET <300) mMETDifference_GenMETTrue_MET200to300->Fill( MET - genMET );
      else if (MET >300 && MET <400) mMETDifference_GenMETTrue_MET300to400->Fill( MET - genMET );
      else if (MET >400 && MET <500) mMETDifference_GenMETTrue_MET400to500->Fill( MET - genMET );
      
    } else {
      edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task:  genMetTrue";
    }    
  }
  if(!isMiniAODMET){
    edm::Handle<GenMETCollection> genCalo;
    iEvent.getByToken(genMETsCaloToken_, genCalo);
    if (genCalo.isValid()) {
      const GenMETCollection *genmetcol = genCalo.product();
      const GenMET  *genMetCalo = &(genmetcol->front());
      const double genMET = genMetCalo->pt();
      const double genMETPhi = genMetCalo->phi();
      
      mMETDifference_GenMETCalo->Fill( MET - genMET );
      mMETDeltaPhi_GenMETCalo->Fill( TMath::ACos( TMath::Cos( METPhi - genMETPhi ) ) );
    } else {
      edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task:  genMetCalo";
    }    
  }
  if ( isCaloMET) { 
    const reco::CaloMET * calomet = &(caloMETs->front());
    // ==========================================================
    // Reconstructed MET Information
    const double caloMaxEtInEMTowers = calomet->maxEtInEmTowers();
    const double caloMaxEtInHadTowers = calomet->maxEtInHadTowers();
    const double caloEtFractionHadronic = calomet->etFractionHadronic();
    const double caloEmEtFraction = calomet->emEtFraction();
    const double caloHadEtInHB = calomet->hadEtInHB();
    const double caloHadEtInHO = calomet->hadEtInHO();
    const double caloHadEtInHE = calomet->hadEtInHE();
    const double caloHadEtInHF = calomet->hadEtInHF();
    const double caloEmEtInEB = calomet->emEtInEB();
    const double caloEmEtInEE = calomet->emEtInEE();
    const double caloEmEtInHF = calomet->emEtInHF();
    const double caloSETInpHF = calomet->CaloSETInpHF();
    const double caloSETInmHF = calomet->CaloSETInmHF();

    mCaloMaxEtInEmTowers->Fill(caloMaxEtInEMTowers);
    mCaloMaxEtInHadTowers->Fill(caloMaxEtInHadTowers);
    mCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
    mCaloEmEtFraction->Fill(caloEmEtFraction);
    mCaloHadEtInHB->Fill(caloHadEtInHB);
    mCaloHadEtInHO->Fill(caloHadEtInHO);
    mCaloHadEtInHE->Fill(caloHadEtInHE);
    mCaloHadEtInHF->Fill(caloHadEtInHF);
    mCaloEmEtInEB->Fill(caloEmEtInEB);
    mCaloEmEtInEE->Fill(caloEmEtInEE);
    mCaloEmEtInHF->Fill(caloEmEtInHF);
    mCaloSETInpHF->Fill(caloSETInpHF);
    mCaloSETInmHF->Fill(caloSETInmHF);
  }  
  if (isGenMET)
  {
    const GenMET *genmet;
    // Get Generated MET
    genmet = &(genMETs->front());

    const double NeutralEMEtFraction = genmet->NeutralEMEtFraction() ;
    const double NeutralHadEtFraction = genmet->NeutralHadEtFraction() ;
    const double ChargedEMEtFraction = genmet->ChargedEMEtFraction () ;
    const double ChargedHadEtFraction = genmet->ChargedHadEtFraction();
    const double MuonEtFraction = genmet->MuonEtFraction() ;
    const double InvisibleEtFraction = genmet->InvisibleEtFraction() ;

    mNeutralEMEtFraction->Fill( NeutralEMEtFraction );
    mNeutralHadEtFraction->Fill( NeutralHadEtFraction );
    mChargedEMEtFraction->Fill( ChargedEMEtFraction );
    mChargedHadEtFraction->Fill( ChargedHadEtFraction );
    mMuonEtFraction->Fill( MuonEtFraction );
    mInvisibleEtFraction->Fill( InvisibleEtFraction );
  } 
  if(isPFMET)
  {
    const reco::PFMET  * pfmet = &(pfMETs->front());
    mPFphotonEtFraction->Fill(pfmet->photonEtFraction());
    mPFphotonEt->Fill(pfmet->photonEt());
    mPFneutralHadronEtFraction->Fill(pfmet->neutralHadronEtFraction());
    mPFneutralHadronEt->Fill(pfmet->neutralHadronEt());
    mPFelectronEtFraction->Fill(pfmet->electronEtFraction());
    mPFelectronEt->Fill(pfmet->electronEt());
    mPFchargedHadronEtFraction->Fill(pfmet->chargedHadronEtFraction());
    mPFchargedHadronEt->Fill(pfmet->chargedHadronEt());
    mPFmuonEtFraction->Fill(pfmet->muonEtFraction());
    mPFmuonEt->Fill(pfmet->muonEt());
    mPFHFHadronEtFraction->Fill(pfmet->HFHadronEtFraction());
    mPFHFHadronEt->Fill(pfmet->HFHadronEt());
    mPFHFEMEtFraction->Fill(pfmet->HFEMEtFraction());
    mPFHFEMEt->Fill(pfmet->HFEMEt());
    // Reconstructed MET Information                                                                                                     

  } 
   if(isMiniAODMET){
    const pat::MET *patmet = &(patMET->front());
    if(patmet->isPFMET()){
      mPFphotonEtFraction->Fill(patmet->NeutralEMFraction());
      mPFneutralHadronEtFraction->Fill(patmet->NeutralHadEtFraction());
      mPFelectronEtFraction->Fill(patmet->ChargedEMEtFraction());
      mPFchargedHadronEtFraction->Fill(patmet->ChargedHadEtFraction());
      mPFmuonEtFraction->Fill(patmet->MuonEtFraction());
      mPFHFHadronEtFraction->Fill(patmet->Type6EtFraction());//HFHadrons
      mPFHFEMEtFraction->Fill(patmet->Type7EtFraction());//HFEMEt  
    }
  }  
  //This is so dirty I could cry. It should be called only ONCE in endJob. But the MonitorElements don't exist then any more.
  FillMETRes();
}

void METTester::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 112 of file METTester.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), and DQMStore::IBooker::setCurrentFolder().

{
    ibooker.setCurrentFolder("JetMET/METValidation/"+inputMETLabel_.label()) ;

    mNvertex                     = ibooker.book1D("Nvertex","Nvertex",80,0,80);
    mMEx                         = ibooker.book1D("MEx","MEx",160,-800,800); 
    mMEy                         = ibooker.book1D("MEy","MEy",160,-800,800);
    mMETSig                      = ibooker.book1D("METSig","METSig",25,0,24.5);
    mMET                         = ibooker.book1D("MET", "MET (20 GeV binning)"           , 100,0,2000);
    mMETFine                     = ibooker.book1D("METFine", "MET (2 GeV binning)"        , 1000,0,2000);
    mMET_Nvtx                    = ibooker.bookProfile("MET_Nvtx", "MET vs. nvtx",    60, 0., 60., 0., 2000., " ");
    mMETPhi                      = ibooker.book1D("METPhi","METPhi",80,-4,4);
    mSumET                       = ibooker.book1D("SumET"            , "SumET"            , 200,0,4000);   //10GeV
    mMETDifference_GenMETTrue    = ibooker.book1D("METDifference_GenMETTrue","METDifference_GenMETTrue", 500,-500,500); 
    mMETDeltaPhi_GenMETTrue      = ibooker.book1D("METDeltaPhi_GenMETTrue","METDeltaPhi_GenMETTrue", 80,0,4); 
    if(!isMiniAODMET){
      mMETDifference_GenMETCalo    = ibooker.book1D("METDifference_GenMETCalo","METDifference_GenMETCalo", 500,-500,500); 
      mMETDeltaPhi_GenMETCalo      = ibooker.book1D("METDeltaPhi_GenMETCalo","METDeltaPhi_GenMETCalo", 80,0,4); 
    }
    if(!isGenMET){
      mMETDifference_GenMETTrue_MET0to20    = ibooker.book1D("METResolution_GenMETTrue_MET0to20"   , "METResolution_GenMETTrue_MET0to20"   , 500,-500,500); 
      mMETDifference_GenMETTrue_MET20to40   = ibooker.book1D("METResolution_GenMETTrue_MET20to40"  , "METResolution_GenMETTrue_MET20to40"  , 500,-500,500); 
      mMETDifference_GenMETTrue_MET40to60   = ibooker.book1D("METResolution_GenMETTrue_MET40to60"  , "METResolution_GenMETTrue_MET40to60"  , 500,-500,500); 
      mMETDifference_GenMETTrue_MET60to80   = ibooker.book1D("METResolution_GenMETTrue_MET60to80"  , "METResolution_GenMETTrue_MET60to80"  , 500,-500,500); 
      mMETDifference_GenMETTrue_MET80to100  = ibooker.book1D("METResolution_GenMETTrue_MET80to100" , "METResolution_GenMETTrue_MET80to100" , 500,-500,500); 
      mMETDifference_GenMETTrue_MET100to150 = ibooker.book1D("METResolution_GenMETTrue_MET100to150", "METResolution_GenMETTrue_MET100to150", 500,-500,500); 
      mMETDifference_GenMETTrue_MET150to200 = ibooker.book1D("METResolution_GenMETTrue_MET150to200", "METResolution_GenMETTrue_MET150to200", 500,-500,500); 
      mMETDifference_GenMETTrue_MET200to300 = ibooker.book1D("METResolution_GenMETTrue_MET200to300", "METResolution_GenMETTrue_MET200to300", 500,-500,500); 
      mMETDifference_GenMETTrue_MET300to400 = ibooker.book1D("METResolution_GenMETTrue_MET300to400", "METResolution_GenMETTrue_MET300to400", 500,-500,500); 
      mMETDifference_GenMETTrue_MET400to500 = ibooker.book1D("METResolution_GenMETTrue_MET400to500", "METResolution_GenMETTrue_MET400to500", 500,-500,500); 
      //this will be filled at the end of the job using info from above hists
      int nBins = 10;
      float bins[] = {0.,20.,40.,60.,80.,100.,150.,200.,300.,400.,500.};
      mMETDifference_GenMETTrue_METResolution     = ibooker.book1D("METResolution_GenMETTrue_InMETBins","METResolution_GenMETTrue_InMETBins",nBins, bins); 
    }
    if ( isCaloMET) { 
      mCaloMaxEtInEmTowers             = ibooker.book1D("CaloMaxEtInEmTowers","CaloMaxEtInEmTowers",300,0,1500);   //5GeV
      mCaloMaxEtInHadTowers            = ibooker.book1D("CaloMaxEtInHadTowers","CaloMaxEtInHadTowers",300,0,1500);  //5GeV
      mCaloEtFractionHadronic          = ibooker.book1D("CaloEtFractionHadronic","CaloEtFractionHadronic",100,0,1);
      mCaloEmEtFraction                = ibooker.book1D("CaloEmEtFraction","CaloEmEtFraction",100,0,1);
      mCaloHadEtInHB                   = ibooker.book1D("CaloHadEtInHB","CaloHadEtInHB", 200, 0, 2000);  //5GeV  
      mCaloHadEtInHE                   = ibooker.book1D("CaloHadEtInHE","CaloHadEtInHE", 100, 0, 500);  //5GeV
      mCaloHadEtInHO                   = ibooker.book1D("CaloHadEtInHO","CaloHadEtInHO", 100, 0, 200);  //5GeV
      mCaloHadEtInHF                   = ibooker.book1D("CaloHadEtInHF","CaloHadEtInHF", 100, 0, 200);  //5GeV
      mCaloSETInpHF                    = ibooker.book1D("CaloSETInpHF","CaloSETInpHF",100, 0, 500);
      mCaloSETInmHF                    = ibooker.book1D("CaloSETInmHF","CaloSETInmHF",100, 0, 500);
      mCaloEmEtInEE                    = ibooker.book1D("CaloEmEtInEE","CaloEmEtInEE",100, 0, 500);    //5GeV
      mCaloEmEtInEB                    = ibooker.book1D("CaloEmEtInEB","CaloEmEtInEB",100, 0, 500);   //5GeV
      mCaloEmEtInHF                    = ibooker.book1D("CaloEmEtInHF","CaloEmEtInHF",100, 0, 500);   //5GeV
    } 

    if(isGenMET){        
      mNeutralEMEtFraction    = ibooker.book1D("GenNeutralEMEtFraction", "GenNeutralEMEtFraction", 120, 0.0, 1.2 );
      mNeutralHadEtFraction   = ibooker.book1D("GenNeutralHadEtFraction", "GenNeutralHadEtFraction", 120, 0.0, 1.2 );
      mChargedEMEtFraction    = ibooker.book1D("GenChargedEMEtFraction", "GenChargedEMEtFraction", 120, 0.0, 1.2);
      mChargedHadEtFraction   = ibooker.book1D("GenChargedHadEtFraction", "GenChargedHadEtFraction", 120, 0.0,1.2);
      mMuonEtFraction         = ibooker.book1D("GenMuonEtFraction", "GenMuonEtFraction", 120, 0.0, 1.2 );
      mInvisibleEtFraction    = ibooker.book1D("GenInvisibleEtFraction", "GenInvisibleEtFraction", 120, 0.0, 1.2 );
    }

    if (isPFMET ||  isMiniAODMET){
      mPFphotonEtFraction = ibooker.book1D("photonEtFraction", "photonEtFraction", 100, 0, 1);
      mPFneutralHadronEtFraction = ibooker.book1D("neutralHadronEtFraction", "neutralHadronEtFraction", 100, 0, 1);
      mPFelectronEtFraction = ibooker.book1D("electronEtFraction", "electronEtFraction", 100, 0, 1);
      mPFchargedHadronEtFraction = ibooker.book1D("chargedHadronEtFraction", "chargedHadronEtFraction", 100, 0, 1);
      mPFHFHadronEtFraction = ibooker.book1D("HFHadronEtFraction", "HFHadronEtFraction", 100, 0, 1);
      mPFmuonEtFraction = ibooker.book1D("muonEtFraction", "muonEtFraction", 100, 0, 1);
      mPFHFEMEtFraction = ibooker.book1D("HFEMEtFraction", "HFEMEtFraction", 100, 0, 1);

      if(!isMiniAODMET){
    mPFphotonEt = ibooker.book1D("photonEt", "photonEt", 100, 0, 1000);
    mPFneutralHadronEt = ibooker.book1D("neutralHadronEt", "neutralHadronEt", 100, 0, 1000);
    mPFelectronEt = ibooker.book1D("electronEt", "electronEt", 100, 0, 1000);
    mPFchargedHadronEt = ibooker.book1D("chargedHadronEt", "chargedHadronEt", 100, 0, 1000);
    mPFmuonEt = ibooker.book1D("muonEt", "muonEt", 100, 0, 1000);
    mPFHFHadronEt = ibooker.book1D("HFHadronEt", "HFHadronEt", 100, 0, 500);
    mPFHFEMEt = ibooker.book1D("HFEMEt", "HFEMEt", 100, 0, 300);
      }

    }
  }

void METTester::FillMETRes

(

)

Definition at line 389 of file METTester.cc.

{
  if(!isGenMET){
    mMETDifference_GenMETTrue_METResolution->setBinContent(1, mMETDifference_GenMETTrue_MET0to20->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(2, mMETDifference_GenMETTrue_MET20to40->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(3, mMETDifference_GenMETTrue_MET40to60->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(4, mMETDifference_GenMETTrue_MET60to80->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(5, mMETDifference_GenMETTrue_MET80to100->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(6, mMETDifference_GenMETTrue_MET100to150->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(7, mMETDifference_GenMETTrue_MET150to200->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(8, mMETDifference_GenMETTrue_MET200to300->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(9, mMETDifference_GenMETTrue_MET300to400->getMean());
    mMETDifference_GenMETTrue_METResolution->setBinContent(10, mMETDifference_GenMETTrue_MET400to500->getMean());
    
    //the error computation should be done in a postProcessor in the harvesting step otherwise the histograms will be just summed
    mMETDifference_GenMETTrue_METResolution->setBinError(1, mMETDifference_GenMETTrue_MET0to20->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(2, mMETDifference_GenMETTrue_MET20to40->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(3, mMETDifference_GenMETTrue_MET40to60->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(4, mMETDifference_GenMETTrue_MET60to80->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(5, mMETDifference_GenMETTrue_MET80to100->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(6, mMETDifference_GenMETTrue_MET100to150->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(7, mMETDifference_GenMETTrue_MET150to200->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(8, mMETDifference_GenMETTrue_MET200to300->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(9, mMETDifference_GenMETTrue_MET300to400->getRMS());
    mMETDifference_GenMETTrue_METResolution->setBinError(10, mMETDifference_GenMETTrue_MET400to500->getRMS());
  }
}

Member Data Documentation

edm::EDGetTokenT<reco::CaloMETCollection> METTester::caloMETsToken_

private

Definition at line 77 of file METTester.h.

edm::EDGetTokenT<reco::GenMETCollection> METTester::genMETsCaloToken_

private

Definition at line 82 of file METTester.h.

edm::EDGetTokenT<reco::GenMETCollection> METTester::genMETsToken_

private

Definition at line 80 of file METTester.h.

edm::EDGetTokenT<reco::GenMETCollection> METTester::genMETsTrueToken_

private

Definition at line 81 of file METTester.h.

edm::InputTag METTester::inputCaloMETLabel_

private

Definition at line 73 of file METTester.h.

edm::InputTag METTester::inputMETLabel_

private

Definition at line 70 of file METTester.h.

bool METTester::isCaloMET

private

Definition at line 160 of file METTester.h.

bool METTester::isGenMET

private

Definition at line 164 of file METTester.h.

bool METTester::isMiniAODMET

private

Definition at line 165 of file METTester.h.

bool METTester::isPFMET

private

Definition at line 163 of file METTester.h.

MonitorElement* METTester::mCaloEmEtFraction

private

Definition at line 107 of file METTester.h.

MonitorElement* METTester::mCaloEmEtInEB

private

Definition at line 118 of file METTester.h.

MonitorElement* METTester::mCaloEmEtInEE

private

Definition at line 117 of file METTester.h.

MonitorElement* METTester::mCaloEmEtInHF

private

Definition at line 114 of file METTester.h.

MonitorElement* METTester::mCaloEtFractionHadronic

private

Definition at line 106 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInEB

private

Definition at line 112 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInEE

private

Definition at line 113 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInHB

private

Definition at line 108 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInHE

private

Definition at line 110 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInHF

private

Definition at line 111 of file METTester.h.

MonitorElement* METTester::mCaloHadEtInHO

private

Definition at line 109 of file METTester.h.

MonitorElement* METTester::mCaloMaxEtInEmTowers

private

Definition at line 104 of file METTester.h.

MonitorElement* METTester::mCaloMaxEtInHadTowers

private

Definition at line 105 of file METTester.h.

MonitorElement* METTester::mCaloSETInmHF

private

Definition at line 116 of file METTester.h.

MonitorElement* METTester::mCaloSETInpHF

private

Definition at line 115 of file METTester.h.

MonitorElement* METTester::mChargedEMEtFraction

private

Definition at line 123 of file METTester.h.

MonitorElement* METTester::mChargedHadEtFraction

private

Definition at line 124 of file METTester.h.

std::map<std::string, MonitorElement*> METTester::me

private

Definition at line 65 of file METTester.h.

std::string METTester::METType_

private

Definition at line 71 of file METTester.h.

edm::InputTag METTester::mInputCollection_

private

Definition at line 69 of file METTester.h.

MonitorElement* METTester::mInvisibleEtFraction

private

Definition at line 126 of file METTester.h.

MonitorElement* METTester::mMET

private

Definition at line 93 of file METTester.h.

MonitorElement* METTester::mMET_Nvtx

private

Definition at line 95 of file METTester.h.

MonitorElement* METTester::mMETDeltaPhi_GenMETCalo

private

Definition at line 101 of file METTester.h.

MonitorElement* METTester::mMETDeltaPhi_GenMETTrue

private

Definition at line 99 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETCalo

private

Definition at line 100 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue

private

Definition at line 98 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET0to20

private

Definition at line 146 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET100to150

private

Definition at line 151 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET150to200

private

Definition at line 152 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET200to300

private

Definition at line 153 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET20to40

private

Definition at line 147 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET300to400

private

Definition at line 154 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET400to500

private

Definition at line 155 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET40to60

private

Definition at line 148 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET60to80

private

Definition at line 149 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_MET80to100

private

Definition at line 150 of file METTester.h.

MonitorElement* METTester::mMETDifference_GenMETTrue_METResolution

private

Definition at line 156 of file METTester.h.

MonitorElement* METTester::mMETFine

private

Definition at line 94 of file METTester.h.

MonitorElement* METTester::mMETPhi

private

Definition at line 96 of file METTester.h.

MonitorElement* METTester::mMETSig

private

Definition at line 92 of file METTester.h.

MonitorElement* METTester::mMEx

private

Definition at line 90 of file METTester.h.

MonitorElement* METTester::mMEy

private

Definition at line 91 of file METTester.h.

MonitorElement* METTester::mMuonEtFraction

private

Definition at line 125 of file METTester.h.

MonitorElement* METTester::mNeutralEMEtFraction

private

Definition at line 121 of file METTester.h.

MonitorElement* METTester::mNeutralHadEtFraction

private

Definition at line 122 of file METTester.h.

MonitorElement* METTester::mNvertex

private

Definition at line 87 of file METTester.h.

MonitorElement* METTester::mPFchargedHadronEt

private

Definition at line 138 of file METTester.h.

MonitorElement* METTester::mPFchargedHadronEtFraction

private

Definition at line 137 of file METTester.h.

MonitorElement* METTester::mPFelectronEt

private

Definition at line 136 of file METTester.h.

MonitorElement* METTester::mPFelectronEtFraction

private

Definition at line 135 of file METTester.h.

MonitorElement* METTester::mPFHFEMEt

private

Definition at line 144 of file METTester.h.

MonitorElement* METTester::mPFHFEMEtFraction

private

Definition at line 143 of file METTester.h.

MonitorElement* METTester::mPFHFHadronEt

private

Definition at line 142 of file METTester.h.

MonitorElement* METTester::mPFHFHadronEtFraction

private

Definition at line 141 of file METTester.h.

MonitorElement* METTester::mPFmuonEt

private

Definition at line 140 of file METTester.h.

MonitorElement* METTester::mPFmuonEtFraction

private

Definition at line 139 of file METTester.h.

MonitorElement* METTester::mPFneutralHadronEt

private

Definition at line 134 of file METTester.h.

MonitorElement* METTester::mPFneutralHadronEtFraction

private

Definition at line 133 of file METTester.h.

MonitorElement* METTester::mPFphotonEt

private

Definition at line 132 of file METTester.h.

MonitorElement* METTester::mPFphotonEtFraction

private

Definition at line 131 of file METTester.h.

MonitorElement* METTester::mSumET

private

Definition at line 97 of file METTester.h.

edm::EDGetTokenT<pat::METCollection> METTester::patMETToken_

private

Definition at line 83 of file METTester.h.

edm::EDGetTokenT<reco::PFMETCollection> METTester::pfMETsToken_

private

Definition at line 78 of file METTester.h.

edm::EDGetTokenT<std::vector<reco::Vertex> > METTester::pvToken_

private

Definition at line 76 of file METTester.h.