JetTester Class Reference

#include <JetTester.h>

Inheritance diagram for JetTester:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	JetTester (const edm::ParameterSet &)
	~JetTester ()
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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	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 Member Functions
void	fillMatchHists (const double GenEta, const double GenPhi, const double GenPt, const double RecoEta, const double RecoPhi, const double RecoPt)

Private Attributes
edm::EDGetTokenT < reco::CaloJetCollection >	caloJetsToken_
MonitorElement *	chargedEmEnergy
MonitorElement *	chargedEmEnergyFraction
MonitorElement *	chargedHadronEnergy
MonitorElement *	chargedHadronEnergyFraction
MonitorElement *	chargedHadronMultiplicity
MonitorElement *	chargedMuEnergy
MonitorElement *	chargedMuEnergyFraction
MonitorElement *	chargedMultiplicity
MonitorElement *	electronEnergy
MonitorElement *	electronEnergyFraction
MonitorElement *	electronMultiplicity
MonitorElement *	emEnergyFraction
MonitorElement *	emEnergyInEB
MonitorElement *	emEnergyInEE
MonitorElement *	emEnergyInHF
MonitorElement *	energyFractionHadronic
edm::EDGetTokenT < GenEventInfoProduct >	evtToken_
edm::EDGetTokenT < reco::GenJetCollection >	genJetsToken_
MonitorElement *	hadEnergyInHB
MonitorElement *	hadEnergyInHE
MonitorElement *	hadEnergyInHF
MonitorElement *	hadEnergyInHO
MonitorElement *	HFEMEnergy
MonitorElement *	HFEMEnergyFraction
MonitorElement *	HFEMMultiplicity
MonitorElement *	HFHadronEnergy
MonitorElement *	HFHadronEnergyFraction
MonitorElement *	HFHadronMultiplicity
MonitorElement *	HOEnergy
MonitorElement *	HOEnergyFraction
bool	isCaloJet
bool	isMiniAODJet
bool	isPFJet
edm::EDGetTokenT < reco::JetCorrector >	jetCorrectorToken_
std::string	JetType
MonitorElement *	maxEInEmTowers
MonitorElement *	maxEInHadTowers
MonitorElement *	mConstituents
MonitorElement *	mCorrJetEta
MonitorElement *	mCorrJetEta_Pt40
MonitorElement *	mCorrJetPhi
MonitorElement *	mCorrJetPhi_Pt40
MonitorElement *	mCorrJetPt
MonitorElement *	mDeltaEta
MonitorElement *	mDeltaPhi
MonitorElement *	mDeltaPt
MonitorElement *	mEnergy
MonitorElement *	mEta
MonitorElement *	mEtaFirst
MonitorElement *	mGenEta
MonitorElement *	mGenEtaFirst
MonitorElement *	mGenPhi
MonitorElement *	mGenPhiFirst
MonitorElement *	mGenPt
edm::InputTag	mInputCollection
edm::InputTag	mInputGenCollection
MonitorElement *	mJetArea
edm::InputTag	mJetCorrector
MonitorElement *	mMass
double	mMatchGenPtThreshold
MonitorElement *	mMjj
MonitorElement *	mNJets1
MonitorElement *	mNJets2
MonitorElement *	mNJets_40
MonitorElement *	mNJetsEta_B_20_40
MonitorElement *	mNJetsEta_B_40
MonitorElement *	mNJetsEta_E_20_40
MonitorElement *	mNJetsEta_E_40
MonitorElement *	mNvtx
MonitorElement *	mP
MonitorElement *	mPhi
MonitorElement *	mPhiFirst
MonitorElement *	mPt
MonitorElement *	mPtCorrOverGen_GenEta_1500_3500
MonitorElement *	mPtCorrOverGen_GenEta_200_600
MonitorElement *	mPtCorrOverGen_GenEta_20_40
MonitorElement *	mPtCorrOverGen_GenEta_3500
MonitorElement *	mPtCorrOverGen_GenEta_3500_5000
MonitorElement *	mPtCorrOverGen_GenEta_40_200
MonitorElement *	mPtCorrOverGen_GenEta_5000_6500
MonitorElement *	mPtCorrOverGen_GenEta_600_1500
MonitorElement *	mPtCorrOverGen_GenPt_B
MonitorElement *	mPtCorrOverGen_GenPt_E
MonitorElement *	mPtCorrOverGen_GenPt_F
MonitorElement *	mPtCorrOverReco_Eta_1500_3500
MonitorElement *	mPtCorrOverReco_Eta_200_600
MonitorElement *	mPtCorrOverReco_Eta_20_40
MonitorElement *	mPtCorrOverReco_Eta_3500
MonitorElement *	mPtCorrOverReco_Eta_3500_5000
MonitorElement *	mPtCorrOverReco_Eta_40_200
MonitorElement *	mPtCorrOverReco_Eta_5000_6500
MonitorElement *	mPtCorrOverReco_Eta_600_1500
MonitorElement *	mPtCorrOverReco_Pt_B
MonitorElement *	mPtCorrOverReco_Pt_E
MonitorElement *	mPtCorrOverReco_Pt_F
MonitorElement *	mPtFirst
MonitorElement *	mPtHat
MonitorElement *	mPtRecoOverGen_B_1500_3500
MonitorElement *	mPtRecoOverGen_B_200_600
MonitorElement *	mPtRecoOverGen_B_20_40
MonitorElement *	mPtRecoOverGen_B_3500
MonitorElement *	mPtRecoOverGen_B_3500_5000
MonitorElement *	mPtRecoOverGen_B_40_200
MonitorElement *	mPtRecoOverGen_B_5000_6500
MonitorElement *	mPtRecoOverGen_B_600_1500
MonitorElement *	mPtRecoOverGen_E_1500_3500
MonitorElement *	mPtRecoOverGen_E_200_600
MonitorElement *	mPtRecoOverGen_E_20_40
MonitorElement *	mPtRecoOverGen_E_3500
MonitorElement *	mPtRecoOverGen_E_3500_5000
MonitorElement *	mPtRecoOverGen_E_40_200
MonitorElement *	mPtRecoOverGen_E_5000_6500
MonitorElement *	mPtRecoOverGen_E_600_1500
MonitorElement *	mPtRecoOverGen_F_1500_3500
MonitorElement *	mPtRecoOverGen_F_200_600
MonitorElement *	mPtRecoOverGen_F_20_40
MonitorElement *	mPtRecoOverGen_F_3500
MonitorElement *	mPtRecoOverGen_F_40_200
MonitorElement *	mPtRecoOverGen_F_600_1500
MonitorElement *	mPtRecoOverGen_GenEta_1500_3500
MonitorElement *	mPtRecoOverGen_GenEta_200_600
MonitorElement *	mPtRecoOverGen_GenEta_20_40
MonitorElement *	mPtRecoOverGen_GenEta_3500
MonitorElement *	mPtRecoOverGen_GenEta_3500_5000
MonitorElement *	mPtRecoOverGen_GenEta_40_200
MonitorElement *	mPtRecoOverGen_GenEta_5000_6500
MonitorElement *	mPtRecoOverGen_GenEta_600_1500
MonitorElement *	mPtRecoOverGen_GenPhi_B
MonitorElement *	mPtRecoOverGen_GenPhi_E
MonitorElement *	mPtRecoOverGen_GenPhi_F
MonitorElement *	mPtRecoOverGen_GenPt_B
MonitorElement *	mPtRecoOverGen_GenPt_E
MonitorElement *	mPtRecoOverGen_GenPt_F
double	mRecoJetPtThreshold
double	mRThreshold
MonitorElement *	muonEnergy
MonitorElement *	muonEnergyFraction
MonitorElement *	muonMultiplicity
MonitorElement *	n60
MonitorElement *	n90
MonitorElement *	neutralEmEnergy
MonitorElement *	neutralEmEnergyFraction
MonitorElement *	neutralHadronEnergy
MonitorElement *	neutralHadronEnergyFraction
MonitorElement *	neutralHadronMultiplicity
MonitorElement *	neutralMultiplicity
edm::EDGetTokenT < pat::JetCollection >	patJetsToken_
edm::EDGetTokenT < reco::PFJetCollection >	pfJetsToken_
MonitorElement *	photonEnergy
MonitorElement *	photonEnergyFraction
MonitorElement *	photonMultiplicity
edm::EDGetTokenT< std::vector < reco::Vertex > >	pvToken_
MonitorElement *	towersArea

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
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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::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 43 of file JetTester.h.

Constructor & Destructor Documentation

JetTester::JetTester

(

const edm::ParameterSet &

iConfig

)

returns the maximum energy deposited in ECAL towers

returns the maximum energy deposited in HCAL towers

returns the jet hadronic energy fraction

returns the jet electromagnetic energy fraction

returns the jet hadronic energy in HB

returns the jet hadronic energy in HO

returns the jet hadronic energy in HE

returns the jet hadronic energy in HF

returns the jet electromagnetic energy in EB

returns the jet electromagnetic energy in EE

returns the jet electromagnetic energy extracted from HF

returns area of contributing towers

returns the number of constituents carrying a 90% of the total Jet energy*/

returns the number of constituents carrying a 60% of the total Jet energy*/

chargedMultiplicity

muonMultiplicity

chargedMultiplicity

chargedEmEnergy

neutralEmEnergy

chargedHadronEnergy

neutralHadronEnergy

chargedHadronEnergyFraction (relative to uncorrected jet energy)

neutralHadronEnergyFraction (relative to uncorrected jet energy)

chargedEmEnergyFraction (relative to uncorrected jet energy)

neutralEmEnergyFraction (relative to uncorrected jet energy)

photonEnergy

photonEnergyFraction (relative to corrected jet energy)

electronEnergy

electronEnergyFraction (relative to corrected jet energy)

muonEnergy

muonEnergyFraction (relative to corrected jet energy)

HFHadronEnergy

HFHadronEnergyFraction (relative to corrected jet energy)

HFEMEnergy

HFEMEnergyFraction (relative to corrected jet energy)

chargedHadronMultiplicity

neutralHadronMultiplicity

photonMultiplicity

electronMultiplicity

HFHadronMultiplicity

HFEMMultiplicity

chargedMuEnergy

chargedMuEnergyFraction

neutralMultiplicity

HOEnergy

HOEnergyFraction (relative to corrected jet energy)

Definition at line 12 of file JetTester.cc.

References caloJetsToken_, chargedEmEnergy, chargedEmEnergyFraction, chargedHadronEnergy, chargedHadronEnergyFraction, chargedHadronMultiplicity, chargedMuEnergy, chargedMuEnergyFraction, chargedMultiplicity, electronEnergy, electronEnergyFraction, electronMultiplicity, emEnergyFraction, emEnergyInEB, emEnergyInEE, emEnergyInHF, energyFractionHadronic, evtToken_, genJetsToken_, edm::ParameterSet::getParameter(), hadEnergyInHB, hadEnergyInHE, hadEnergyInHF, hadEnergyInHO, HFEMEnergy, HFEMEnergyFraction, HFEMMultiplicity, HFHadronEnergy, HFHadronEnergyFraction, HFHadronMultiplicity, HOEnergy, HOEnergyFraction, HLT_25ns14e33_v1_cff::InputTag, isCaloJet, isMiniAODJet, isPFJet, jetCorrectorToken_, JetType, edm::InputTag::label(), maxEInEmTowers, maxEInHadTowers, mConstituents, mCorrJetEta, mCorrJetEta_Pt40, mCorrJetPhi, mCorrJetPhi_Pt40, mCorrJetPt, mDeltaEta, mDeltaPhi, mDeltaPt, mEnergy, mEta, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mInputCollection, mInputGenCollection, mJetArea, mJetCorrector, mMass, mMjj, mNJets1, mNJets2, mNJetsEta_B_20_40, mNJetsEta_B_40, mNJetsEta_E_20_40, mNJetsEta_E_40, mNvtx, mP, mPhi, mPhiFirst, mPt, mPtCorrOverGen_GenEta_1500_3500, mPtCorrOverGen_GenEta_200_600, mPtCorrOverGen_GenEta_20_40, mPtCorrOverGen_GenEta_3500, mPtCorrOverGen_GenEta_3500_5000, mPtCorrOverGen_GenEta_40_200, mPtCorrOverGen_GenEta_5000_6500, mPtCorrOverGen_GenEta_600_1500, mPtCorrOverGen_GenPt_B, mPtCorrOverGen_GenPt_E, mPtCorrOverGen_GenPt_F, mPtCorrOverReco_Eta_1500_3500, mPtCorrOverReco_Eta_200_600, mPtCorrOverReco_Eta_20_40, mPtCorrOverReco_Eta_3500, mPtCorrOverReco_Eta_3500_5000, mPtCorrOverReco_Eta_40_200, mPtCorrOverReco_Eta_5000_6500, mPtCorrOverReco_Eta_600_1500, mPtCorrOverReco_Pt_B, mPtCorrOverReco_Pt_E, mPtCorrOverReco_Pt_F, mPtFirst, mPtHat, mPtRecoOverGen_B_1500_3500, mPtRecoOverGen_B_200_600, mPtRecoOverGen_B_20_40, mPtRecoOverGen_B_3500, mPtRecoOverGen_B_3500_5000, mPtRecoOverGen_B_40_200, mPtRecoOverGen_B_5000_6500, mPtRecoOverGen_B_600_1500, mPtRecoOverGen_E_1500_3500, mPtRecoOverGen_E_200_600, mPtRecoOverGen_E_20_40, mPtRecoOverGen_E_3500, mPtRecoOverGen_E_3500_5000, mPtRecoOverGen_E_40_200, mPtRecoOverGen_E_5000_6500, mPtRecoOverGen_E_600_1500, mPtRecoOverGen_F_1500_3500, mPtRecoOverGen_F_200_600, mPtRecoOverGen_F_20_40, mPtRecoOverGen_F_3500, mPtRecoOverGen_F_40_200, mPtRecoOverGen_F_600_1500, mPtRecoOverGen_GenEta_1500_3500, mPtRecoOverGen_GenEta_200_600, mPtRecoOverGen_GenEta_20_40, mPtRecoOverGen_GenEta_3500, mPtRecoOverGen_GenEta_3500_5000, mPtRecoOverGen_GenEta_40_200, mPtRecoOverGen_GenEta_5000_6500, mPtRecoOverGen_GenEta_600_1500, mPtRecoOverGen_GenPhi_B, mPtRecoOverGen_GenPhi_E, mPtRecoOverGen_GenPhi_F, mPtRecoOverGen_GenPt_B, mPtRecoOverGen_GenPt_E, mPtRecoOverGen_GenPt_F, muonEnergy, muonEnergyFraction, muonMultiplicity, n60, n90, neutralEmEnergy, neutralEmEnergyFraction, neutralHadronEnergy, neutralHadronEnergyFraction, neutralHadronMultiplicity, neutralMultiplicity, patJetsToken_, pfJetsToken_, photonEnergy, photonEnergyFraction, photonMultiplicity, pvToken_, AlCaHLTBitMon_QueryRunRegistry::string, and towersArea.

                                                   :
  mInputCollection               (iConfig.getParameter<edm::InputTag>       ("src")),
//  rhoTag                         (iConfig.getParameter<edm::InputTag>       ("srcRho")), 
  JetType                        (iConfig.getUntrackedParameter<std::string>("JetType")),
  mRecoJetPtThreshold            (iConfig.getParameter<double>              ("recoJetPtThreshold")),
  mMatchGenPtThreshold           (iConfig.getParameter<double>              ("matchGenPtThreshold")),
  mRThreshold                    (iConfig.getParameter<double>              ("RThreshold"))
{
  std::string inputCollectionLabel(mInputCollection.label());

  isCaloJet = (std::string("calo")==JetType);
  isPFJet   = (std::string("pf")  ==JetType);
  isMiniAODJet = (std::string("miniaod")  ==JetType);
  if(!isMiniAODJet){
    mJetCorrector                  =iConfig.getParameter<edm::InputTag>("JetCorrections"); 
  } 

  //consumes
  pvToken_ = consumes<std::vector<reco::Vertex> >(iConfig.getParameter<edm::InputTag>("primVertex"));
  if (isCaloJet) caloJetsToken_  = consumes<reco::CaloJetCollection>(mInputCollection);
  if (isPFJet)   pfJetsToken_    = consumes<reco::PFJetCollection>(mInputCollection);
  if(isMiniAODJet)patJetsToken_ = consumes<pat::JetCollection>(mInputCollection);
  mInputGenCollection            =iConfig.getParameter<edm::InputTag>("srcGen");
  genJetsToken_ = consumes<reco::GenJetCollection>(edm::InputTag(mInputGenCollection));
  evtToken_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
  if(!isMiniAODJet && !mJetCorrector.label().empty()){
    jetCorrectorToken_ = consumes<reco::JetCorrector>(mJetCorrector);
  }
  
  // Events variables
  mNvtx           = 0;

  // Jet parameters
  mEta          = 0;
  mPhi          = 0;
  mEnergy       = 0;
  mP            = 0;
  mPt           = 0;
  mMass         = 0;
  mConstituents = 0;
  mJetArea      = 0;
//  mRho          = 0;

  // Corrected jets
  mCorrJetPt  = 0;
  mCorrJetEta = 0;
  mCorrJetPhi = 0;
  mCorrJetEta_Pt40 = 0;
  mCorrJetPhi_Pt40 = 0;

  // Corrected jets profiles
  mPtCorrOverReco_Pt_B                = 0;
  mPtCorrOverReco_Pt_E                = 0;
  mPtCorrOverReco_Pt_F                = 0;
  mPtCorrOverReco_Eta_20_40          = 0;
  mPtCorrOverReco_Eta_40_200          = 0;
  mPtCorrOverReco_Eta_200_600         = 0;
  mPtCorrOverReco_Eta_600_1500        = 0;
  mPtCorrOverReco_Eta_1500_3500       = 0;
  mPtCorrOverReco_Eta_3500_5000       = 0;
  mPtCorrOverReco_Eta_5000_6500       = 0;
  mPtCorrOverReco_Eta_3500       = 0;
  mPtCorrOverGen_GenPt_B          = 0;
  mPtCorrOverGen_GenPt_E          = 0;
  mPtCorrOverGen_GenPt_F          = 0;
  mPtCorrOverGen_GenEta_20_40    = 0;
  mPtCorrOverGen_GenEta_40_200    = 0;
  mPtCorrOverGen_GenEta_200_600   = 0;
  mPtCorrOverGen_GenEta_600_1500  = 0;
  mPtCorrOverGen_GenEta_1500_3500 = 0;
  mPtCorrOverGen_GenEta_3500_5000 = 0;
  mPtCorrOverGen_GenEta_5000_6500 = 0;
  mPtCorrOverGen_GenEta_3500 = 0;

  // Generation
  mGenEta      = 0;
  mGenPhi      = 0;
  mGenPt       = 0;
  mGenEtaFirst = 0;
  mGenPhiFirst = 0;
  mPtHat       = 0;
  mDeltaEta    = 0;
  mDeltaPhi    = 0;
  mDeltaPt     = 0;

  mPtRecoOverGen_B_20_40    = 0;
  mPtRecoOverGen_E_20_40    = 0;
  mPtRecoOverGen_F_20_40    = 0;
  mPtRecoOverGen_B_40_200    = 0;
  mPtRecoOverGen_E_40_200    = 0;
  mPtRecoOverGen_F_40_200    = 0;
  mPtRecoOverGen_B_200_600   = 0;
  mPtRecoOverGen_E_200_600   = 0;
  mPtRecoOverGen_F_200_600   = 0;
  mPtRecoOverGen_B_600_1500  = 0;
  mPtRecoOverGen_E_600_1500  = 0;
  mPtRecoOverGen_F_600_1500  = 0;
  mPtRecoOverGen_B_1500_3500 = 0;
  mPtRecoOverGen_E_1500_3500 = 0;
  mPtRecoOverGen_F_1500_3500 = 0;
  mPtRecoOverGen_B_3500_5000 = 0;
  mPtRecoOverGen_E_3500_5000 = 0;
  mPtRecoOverGen_B_5000_6500 = 0;
  mPtRecoOverGen_E_5000_6500 = 0;
  mPtRecoOverGen_B_3500 = 0;
  mPtRecoOverGen_E_3500 = 0;
  mPtRecoOverGen_F_3500 = 0;

  // Generation profiles
  mPtRecoOverGen_GenPt_B          = 0;
  mPtRecoOverGen_GenPt_E          = 0;
  mPtRecoOverGen_GenPt_F          = 0;
  mPtRecoOverGen_GenPhi_B         = 0;
  mPtRecoOverGen_GenPhi_E         = 0;
  mPtRecoOverGen_GenPhi_F         = 0;
  mPtRecoOverGen_GenEta_20_40    = 0;
  mPtRecoOverGen_GenEta_40_200    = 0;
  mPtRecoOverGen_GenEta_200_600   = 0;
  mPtRecoOverGen_GenEta_600_1500  = 0;
  mPtRecoOverGen_GenEta_1500_3500 = 0;
  mPtRecoOverGen_GenEta_3500_5000 = 0;
  mPtRecoOverGen_GenEta_5000_6500 = 0;
  mPtRecoOverGen_GenEta_3500 = 0;

  // Some jet algebra
  mEtaFirst   = 0;
  mPhiFirst   = 0;
  mPtFirst    = 0;
  mMjj        = 0;
  mNJetsEta_B_20_40 = 0;
  mNJetsEta_E_20_40 = 0;
  mNJetsEta_B_40 = 0;
  mNJetsEta_E_40 = 0;
  mNJets1     = 0;
  mNJets2     = 0;

//  // PFJet specific
//  mHadEnergyInHF       = 0;
//  mEmEnergyInHF        = 0;
//  mChargedEmEnergy     = 0;
//  mChargedHadronEnergy = 0;
//  mNeutralEmEnergy     = 0;
//  mNeutralHadronEnergy = 0;

  // ---- Calo Jet specific information ----
  maxEInEmTowers = 0;
  maxEInHadTowers = 0;
  energyFractionHadronic = 0;
  emEnergyFraction = 0;
  hadEnergyInHB = 0;
  hadEnergyInHO = 0;
  hadEnergyInHE = 0;
  hadEnergyInHF = 0;
  emEnergyInEB = 0;
  emEnergyInEE = 0;
  emEnergyInHF = 0;
  towersArea = 0;
  n90 = 0;
  n60 = 0;

  // ---- JPT Jet specific information ----
//  elecMultiplicity = 0;

  // ---- JPT or PF Jet specific information ----
  muonMultiplicity = 0;
  chargedMultiplicity = 0;
  chargedEmEnergy = 0;
  neutralEmEnergy = 0;
  chargedHadronEnergy = 0;
  neutralHadronEnergy = 0;
  chargedHadronEnergyFraction = 0;
  neutralHadronEnergyFraction = 0;
  chargedEmEnergyFraction = 0;
  neutralEmEnergyFraction = 0;

  // ---- PF Jet specific information ----
  photonEnergy = 0;
  photonEnergyFraction = 0;
  electronEnergy = 0;
  electronEnergyFraction = 0;
  muonEnergy = 0;
  muonEnergyFraction = 0;
  HFHadronEnergy = 0;
  HFHadronEnergyFraction = 0;
  HFEMEnergy = 0;
  HFEMEnergyFraction = 0;
  chargedHadronMultiplicity = 0;
  neutralHadronMultiplicity = 0;
  photonMultiplicity = 0;
  electronMultiplicity = 0;
  HFHadronMultiplicity = 0;
  HFEMMultiplicity = 0;
  chargedMuEnergy = 0;
  chargedMuEnergyFraction = 0;
  neutralMultiplicity = 0;

  HOEnergy = 0;
  HOEnergyFraction = 0;
}

JetTester::~JetTester

(

)

Definition at line 450 of file JetTester.cc.

{}

Member Function Documentation

void JetTester::analyze ( const edm::Event &  mEvent, const edm::EventSetup &  mSetup  )

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 456 of file JetTester.cc.

References caloJetsToken_, chargedEmEnergy, chargedEmEnergyFraction, chargedHadronEnergy, chargedHadronEnergyFraction, chargedHadronMultiplicity, chargedMuEnergy, chargedMuEnergyFraction, chargedMultiplicity, deltaR(), electronEnergy, electronEnergyFraction, electronMultiplicity, emEnergyFraction, emEnergyInEB, emEnergyInEE, emEnergyInHF, relval_parameters_module::energy, energyFractionHadronic, eta, reco::LeafCandidate::eta(), evtToken_, MonitorElement::Fill(), HcalObjRepresent::Fill(), fillMatchHists(), genJetsToken_, edm::Event::getByToken(), hadEnergyInHB, hadEnergyInHE, hadEnergyInHF, hadEnergyInHO, HFEMEnergy, HFEMEnergyFraction, HFEMMultiplicity, HFHadronEnergy, HFHadronEnergyFraction, HFHadronMultiplicity, HOEnergy, HOEnergyFraction, i, isCaloJet, isMiniAODJet, isPFJet, edm::EventBase::isRealData(), edm::HandleBase::isValid(), jetCorrectorToken_, edm::InputTag::label(), maxEInEmTowers, maxEInHadTowers, mConstituents, mCorrJetEta, mCorrJetEta_Pt40, mCorrJetPhi, mCorrJetPhi_Pt40, mCorrJetPt, mEnergy, mEta, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mInputGenCollection, mJetArea, mJetCorrector, mMass, mMatchGenPtThreshold, mMjj, mNJets1, mNJets2, mNJets_40, mNJetsEta_B_20_40, mNJetsEta_B_40, mNJetsEta_E_20_40, mNJetsEta_E_40, mNvtx, mP, mPhi, mPhiFirst, mPt, mPtCorrOverGen_GenEta_1500_3500, mPtCorrOverGen_GenEta_200_600, mPtCorrOverGen_GenEta_20_40, mPtCorrOverGen_GenEta_3500, mPtCorrOverGen_GenEta_3500_5000, mPtCorrOverGen_GenEta_40_200, mPtCorrOverGen_GenEta_5000_6500, mPtCorrOverGen_GenEta_600_1500, mPtCorrOverGen_GenPt_B, mPtCorrOverGen_GenPt_E, mPtCorrOverGen_GenPt_F, mPtCorrOverReco_Eta_1500_3500, mPtCorrOverReco_Eta_200_600, mPtCorrOverReco_Eta_20_40, mPtCorrOverReco_Eta_3500, mPtCorrOverReco_Eta_3500_5000, mPtCorrOverReco_Eta_40_200, mPtCorrOverReco_Eta_5000_6500, mPtCorrOverReco_Eta_600_1500, mPtCorrOverReco_Pt_B, mPtCorrOverReco_Pt_E, mPtCorrOverReco_Pt_F, mPtFirst, mPtHat, mRecoJetPtThreshold, mRThreshold, muonEnergy, muonEnergyFraction, muonMultiplicity, n60, n90, neutralEmEnergy, neutralEmEnergyFraction, neutralHadronEnergy, neutralHadronEnergyFraction, neutralHadronMultiplicity, neutralMultiplicity, AlCaHLTBitMon_ParallelJobs::p, jetProducer_cfi::patJets, patJetsToken_, pfJets_cff::pfJets, pfJetsToken_, phi, reco::LeafCandidate::phi(), photonEnergy, photonEnergyFraction, photonMultiplicity, position, EnergyCorrector::pt, reco::LeafCandidate::pt(), pvToken_, RecoJets_cff::recoJets, rho, pileupReCalc_HLTpaths::scale, reco::Jet::scaleEnergy(), towersArea, and z.

{
  // Get the primary vertices
  //----------------------------------------------------------------------------
  edm::Handle<vector<reco::Vertex> > pvHandle;
  mEvent.getByToken(pvToken_, pvHandle);

  int nGoodVertices = 0;

  if (pvHandle.isValid())
    {
      for (unsigned i=0; i<pvHandle->size(); i++)
    {
      if ((*pvHandle)[i].ndof() > 4 &&
          (fabs((*pvHandle)[i].z()) <= 24) &&
          (fabs((*pvHandle)[i].position().rho()) <= 2))
        nGoodVertices++;
    }
    }

  mNvtx->Fill(nGoodVertices);

//  // Get the jet rho
//  //----------------------------------------------------------------------------
//  edm::Handle<double> pRho;
//  mEvent.getByToken(rhoTag, pRho);
//
//  if (pRho.isValid())
//    {
//      double jetRho = *pRho;
//
//      if (mRho) mRho->Fill(jetRho);
//    }


  // Get the Jet collection
  //----------------------------------------------------------------------------
  math::XYZTLorentzVector p4tmp[2];

  std::vector<Jet> recoJets;
  recoJets.clear();

  edm::Handle<CaloJetCollection> caloJets;
  edm::Handle<PFJetCollection>   pfJets;
//  edm::Handle<JPTJetCollection>  jptJets;
  edm::Handle<pat::JetCollection> patJets;

  if (isCaloJet) mEvent.getByToken(caloJetsToken_, caloJets);
  if (isPFJet)   mEvent.getByToken(pfJetsToken_, pfJets);
//  if (isJPTJet)  mEvent.getByToken(jptJetsToken_, jptJets);
  if (isMiniAODJet)   mEvent.getByToken(patJetsToken_, patJets);

  if (isCaloJet && !caloJets.isValid()) return;
  if (isPFJet   && !pfJets.isValid())   return;
//  if (isJPTJet  && !jptJets.isValid())  return;
 if (isMiniAODJet   && !patJets.isValid())   return;


  if (isCaloJet)
    {
      for (unsigned ijet=0; ijet<caloJets->size(); ijet++)
    recoJets.push_back((*caloJets)[ijet]);
    }

/*  if (isJPTJet)
    {
      for (unsigned ijet=0; ijet<jptJets->size(); ijet++)
    recoJets.push_back((*jptJets)[ijet]);
    }*/

  if (isPFJet) {
    for (unsigned ijet=0; ijet<pfJets->size(); ijet++)
      recoJets.push_back((*pfJets)[ijet]);
  }
  if (isMiniAODJet) {
    for (unsigned ijet=0; ijet<patJets->size(); ijet++)
      recoJets.push_back((*patJets)[ijet]);
  }

  int nJet      = 0;
  int nJet_E_20_40 = 0;
  int nJet_B_20_40 = 0;
  int nJet_E_40 = 0;
  int nJet_B_40 = 0;
  int nJet_40 = 0;

  int index_first_jet=-1;
  double pt_first=-1;

 int index_second_jet=-1;
 double pt_second=-1;

  for (unsigned ijet=0; ijet<recoJets.size(); ijet++) {
    bool pass_lowjet=false;
    bool pass_mediumjet = false;
    if(!isMiniAODJet){
      if (  (recoJets[ijet].pt() > 20.) &&  (recoJets[ijet].pt() < mRecoJetPtThreshold)) {
    pass_lowjet=true;
      }
    }
    if(isMiniAODJet){
      if((recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected"))>20. && ((recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected")) < mRecoJetPtThreshold)){
    pass_lowjet=true;
      }
    }
    if (pass_lowjet) {
      if (fabs(recoJets[ijet].eta()) > 1.5)
        nJet_E_20_40++;
      else
        nJet_B_20_40++;   
    }
    if(!isMiniAODJet){
      if (recoJets[ijet].pt() > mRecoJetPtThreshold) {
    pass_mediumjet = true;
      }
    }
    if(isMiniAODJet){
      if((recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected")) > mRecoJetPtThreshold){
    pass_mediumjet=true;
      }
    }
    if (pass_mediumjet) {
      if(isMiniAODJet){
    if( (recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected"))>pt_first){
      pt_second=pt_first;
      pt_first=recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected");
      index_second_jet=index_first_jet;
      index_first_jet=ijet;
    }else if( (recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected"))>pt_second){
      index_second_jet=ijet;
      pt_second=recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected");
    }
      }
      //counting forward and barrel jets
      if (fabs(recoJets[ijet].eta()) > 1.5)
        nJet_E_40++;
      else
        nJet_B_40++;      
      nJet_40++;

      if (mEta) mEta->Fill(recoJets[ijet].eta());

      if (mJetArea)      mJetArea     ->Fill(recoJets[ijet].jetArea());
      if (mPhi)          mPhi         ->Fill(recoJets[ijet].phi());
      if(!isMiniAODJet){
    if (mEnergy)       mEnergy      ->Fill(recoJets[ijet].energy());
    if (mP)            mP           ->Fill(recoJets[ijet].p());
    if (mPt)           mPt          ->Fill(recoJets[ijet].pt());
    if (mMass)         mMass        ->Fill(recoJets[ijet].mass());
      }else{
    if (mEnergy)       mEnergy      ->Fill(recoJets[ijet].energy()*(*patJets)[ijet].jecFactor("Uncorrected"));
    if (mP)            mP           ->Fill(recoJets[ijet].p()*(*patJets)[ijet].jecFactor("Uncorrected"));
    if (mPt)           mPt          ->Fill(recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected"));
    if (mMass)         mMass        ->Fill(recoJets[ijet].mass()*(*patJets)[ijet].jecFactor("Uncorrected"));
      }
      if (mConstituents) mConstituents->Fill(recoJets[ijet].nConstituents());
      if(!isMiniAODJet){
    if (ijet == 0) {
      if (mEtaFirst) mEtaFirst->Fill(recoJets[ijet].eta());
      if (mPhiFirst) mPhiFirst->Fill(recoJets[ijet].phi());
      if (mPtFirst)  mPtFirst ->Fill(recoJets[ijet].pt());
    }

    if (ijet == 0) {nJet++; p4tmp[0] = recoJets[ijet].p4();}
    if (ijet == 1) {nJet++; p4tmp[1] = recoJets[ijet].p4();}
      }
  //    if (isPFJet || isCaloJet) {
  //      if (mHadEnergyInHF)       mHadEnergyInHF      ->Fill((*pfJets)[ijet].HFHadronEnergy());
  //      if (mEmEnergyInHF)        mEmEnergyInHF       ->Fill((*pfJets)[ijet].HFEMEnergy());
  //      if (mChargedEmEnergy)     mChargedEmEnergy    ->Fill((*pfJets)[ijet].chargedEmEnergy());
  //      if (mChargedHadronEnergy) mChargedHadronEnergy->Fill((*pfJets)[ijet].chargedHadronEnergy());
  //      if (mNeutralEmEnergy)     mNeutralEmEnergy    ->Fill((*pfJets)[ijet].neutralEmEnergy());
  //      if (mNeutralHadronEnergy) mNeutralHadronEnergy->Fill((*pfJets)[ijet].neutralHadronEnergy());
  //    }


      // ---- Calo Jet specific information ----
      if (isCaloJet) {
        maxEInEmTowers              ->Fill((*caloJets)[ijet].maxEInEmTowers());
        maxEInHadTowers             ->Fill((*caloJets)[ijet].maxEInHadTowers());
        energyFractionHadronic      ->Fill((*caloJets)[ijet].energyFractionHadronic());
        emEnergyFraction            ->Fill((*caloJets)[ijet].emEnergyFraction());
        hadEnergyInHB               ->Fill((*caloJets)[ijet].hadEnergyInHB());
        hadEnergyInHO               ->Fill((*caloJets)[ijet].hadEnergyInHO());
        hadEnergyInHE               ->Fill((*caloJets)[ijet].hadEnergyInHE());
        hadEnergyInHF               ->Fill((*caloJets)[ijet].hadEnergyInHF());
        emEnergyInEB                ->Fill((*caloJets)[ijet].emEnergyInEB());
        emEnergyInEE                ->Fill((*caloJets)[ijet].emEnergyInEE());
        emEnergyInHF                ->Fill((*caloJets)[ijet].emEnergyInHF());
        towersArea                  ->Fill((*caloJets)[ijet].towersArea());
        n90                         ->Fill((*caloJets)[ijet].n90());
        n60                         ->Fill((*caloJets)[ijet].n60());
      }
      // ---- PF Jet specific information ----
      if (isPFJet) {
        muonMultiplicity ->Fill((*pfJets)[ijet].muonMultiplicity());
        chargedMultiplicity ->Fill((*pfJets)[ijet].chargedMultiplicity());
        chargedEmEnergy ->Fill((*pfJets)[ijet].chargedEmEnergy());
        neutralEmEnergy ->Fill((*pfJets)[ijet].neutralEmEnergy());
        chargedHadronEnergy ->Fill((*pfJets)[ijet].chargedHadronEnergy());
        neutralHadronEnergy ->Fill((*pfJets)[ijet].neutralHadronEnergy());
        chargedHadronEnergyFraction ->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
        neutralHadronEnergyFraction ->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
        chargedEmEnergyFraction ->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
        neutralEmEnergyFraction ->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
        photonEnergy ->Fill((*pfJets)[ijet].photonEnergy());
        photonEnergyFraction ->Fill((*pfJets)[ijet].photonEnergyFraction());
        electronEnergy ->Fill((*pfJets)[ijet].electronEnergy());
        electronEnergyFraction ->Fill((*pfJets)[ijet].electronEnergyFraction());
        muonEnergy ->Fill((*pfJets)[ijet].muonEnergy());
        muonEnergyFraction ->Fill((*pfJets)[ijet].muonEnergyFraction());
        HFHadronEnergy ->Fill((*pfJets)[ijet].HFHadronEnergy());
        HFHadronEnergyFraction ->Fill((*pfJets)[ijet].HFHadronEnergyFraction());
        HFEMEnergy ->Fill((*pfJets)[ijet].HFEMEnergy());
        HFEMEnergyFraction ->Fill((*pfJets)[ijet].HFEMEnergyFraction());
        chargedHadronMultiplicity ->Fill((*pfJets)[ijet].chargedHadronMultiplicity());
        neutralHadronMultiplicity ->Fill((*pfJets)[ijet].neutralHadronMultiplicity());
        photonMultiplicity ->Fill((*pfJets)[ijet].photonMultiplicity());
        electronMultiplicity ->Fill((*pfJets)[ijet].electronMultiplicity());
        HFHadronMultiplicity ->Fill((*pfJets)[ijet].HFHadronMultiplicity());
        HFEMMultiplicity ->Fill((*pfJets)[ijet].HFEMMultiplicity());
        chargedMuEnergy ->Fill((*pfJets)[ijet].chargedMuEnergy());
        chargedMuEnergyFraction ->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
        neutralMultiplicity ->Fill((*pfJets)[ijet].neutralMultiplicity()); 
        HOEnergy ->Fill((*pfJets)[ijet].hoEnergy());
        HOEnergyFraction ->Fill((*pfJets)[ijet].hoEnergyFraction());
     }
   if (isMiniAODJet && (*patJets)[ijet].isPFJet()) {
        muonMultiplicity ->Fill((*patJets)[ijet].muonMultiplicity());
        chargedMultiplicity ->Fill((*patJets)[ijet].chargedMultiplicity());
        chargedEmEnergy ->Fill((*patJets)[ijet].chargedEmEnergy());
        neutralEmEnergy ->Fill((*patJets)[ijet].neutralEmEnergy());
        chargedHadronEnergy ->Fill((*patJets)[ijet].chargedHadronEnergy());
        neutralHadronEnergy ->Fill((*patJets)[ijet].neutralHadronEnergy());
        chargedHadronEnergyFraction ->Fill((*patJets)[ijet].chargedHadronEnergyFraction());
        neutralHadronEnergyFraction ->Fill((*patJets)[ijet].neutralHadronEnergyFraction());
        chargedEmEnergyFraction ->Fill((*patJets)[ijet].chargedEmEnergyFraction());
        neutralEmEnergyFraction ->Fill((*patJets)[ijet].neutralEmEnergyFraction());
        photonEnergy ->Fill((*patJets)[ijet].photonEnergy());
        photonEnergyFraction ->Fill((*patJets)[ijet].photonEnergyFraction());
        electronEnergy ->Fill((*patJets)[ijet].electronEnergy());
        electronEnergyFraction ->Fill((*patJets)[ijet].electronEnergyFraction());
        muonEnergy ->Fill((*patJets)[ijet].muonEnergy());
        muonEnergyFraction ->Fill((*patJets)[ijet].muonEnergyFraction());
        HFHadronEnergy ->Fill((*patJets)[ijet].HFHadronEnergy());
        HFHadronEnergyFraction ->Fill((*patJets)[ijet].HFHadronEnergyFraction());
        HFEMEnergy ->Fill((*patJets)[ijet].HFEMEnergy());
        HFEMEnergyFraction ->Fill((*patJets)[ijet].HFEMEnergyFraction());
        chargedHadronMultiplicity ->Fill((*patJets)[ijet].chargedHadronMultiplicity());
        neutralHadronMultiplicity ->Fill((*patJets)[ijet].neutralHadronMultiplicity());
        photonMultiplicity ->Fill((*patJets)[ijet].photonMultiplicity());
        electronMultiplicity ->Fill((*patJets)[ijet].electronMultiplicity());
        HFHadronMultiplicity ->Fill((*patJets)[ijet].HFHadronMultiplicity());
        HFEMMultiplicity ->Fill((*patJets)[ijet].HFEMMultiplicity());
        chargedMuEnergy ->Fill((*patJets)[ijet].chargedMuEnergy());
        chargedMuEnergyFraction ->Fill((*patJets)[ijet].chargedMuEnergyFraction());
        neutralMultiplicity ->Fill((*patJets)[ijet].neutralMultiplicity()); 
        HOEnergy ->Fill((*patJets)[ijet].hoEnergy());
        HOEnergyFraction ->Fill((*patJets)[ijet].hoEnergyFraction());
      }
    }//fill quantities for medium jets
  }

  if (mNJetsEta_B_20_40) mNJetsEta_B_20_40->Fill(nJet_B_20_40);
  if (mNJetsEta_E_20_40) mNJetsEta_E_20_40->Fill(nJet_E_20_40);
  if (mNJetsEta_B_40) mNJetsEta_B_40->Fill(nJet_B_40);
  if (mNJetsEta_E_40) mNJetsEta_E_40->Fill(nJet_E_40);
  if (mNJets_40) mNJets_40->Fill(nJet_40); 
  if(!isMiniAODJet){
    if (nJet >= 2)
      {
    if (mMjj) mMjj->Fill((p4tmp[0]+p4tmp[1]).mass());
      }
  }else{
    if(index_first_jet>-1){
      if (mEtaFirst) mEtaFirst->Fill(recoJets[index_first_jet].eta());
      if (mPhiFirst) mPhiFirst->Fill(recoJets[index_first_jet].phi());
      if (mPtFirst)  mPtFirst ->Fill(recoJets[index_first_jet].pt()*(*patJets)[index_first_jet].jecFactor("Uncorrected"));
      nJet++; p4tmp[0] = recoJets[index_first_jet].p4()*(*patJets)[index_first_jet].jecFactor("Uncorrected");
    }
    if(index_second_jet>-1){
      nJet++; p4tmp[1] = recoJets[index_second_jet].p4()*(*patJets)[index_second_jet].jecFactor("Uncorrected");
    }
    if (nJet >= 2)
      {
    if (mMjj) mMjj->Fill((p4tmp[0]+p4tmp[1]).mass());
      }
  }

  // Count jets above pt cut
  //----------------------------------------------------------------------------
  for (int istep=0; istep<100; ++istep)
    {
      int njets1 = 0;
      int njets2 = 0;

      float ptStep1 = (istep * ( 200. / 100.));
      float ptStep2 = (istep * (4000. / 100.));

      for (unsigned ijet=0; ijet<recoJets.size(); ijet++) {
    if(!isMiniAODJet){
      if (recoJets[ijet].pt() > ptStep1) njets1++;
      if (recoJets[ijet].pt() > ptStep2) njets2++;
    }else{
      if ((recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected")) > ptStep1) njets1++;
      if ((recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected")) > ptStep2) njets2++;
    }
    mNJets1->Fill(ptStep1, njets1);
    mNJets2->Fill(ptStep2, njets2);
      }
    }


  // Corrected jets
  //----------------------------------------------------------------------------
  double scale = -999;
  edm::Handle<reco::JetCorrector> jetCorr;
  bool pass_correction_flag =false;
  if(!isMiniAODJet && !mJetCorrector.label().empty()){
    mEvent.getByToken(jetCorrectorToken_, jetCorr);
    if (jetCorr.isValid()){
      pass_correction_flag=true;
    }
  }
  if(isMiniAODJet){
    pass_correction_flag =true;
  }
  for (unsigned ijet=0; ijet<recoJets.size(); ijet++) {
    Jet correctedJet = recoJets[ijet];
    if(pass_correction_flag){
      if (isCaloJet) scale = jetCorr->correction((*caloJets)[ijet]); 
      if (isPFJet)   scale = jetCorr->correction((*pfJets)[ijet]); 
      //if (isJPTJet)  scale = jetCorr->correction((*jptJets)[ijet]);
      if(!isMiniAODJet){
    correctedJet.scaleEnergy(scale); 
      }
      
      if (correctedJet.pt() < 20) continue;
      
      mCorrJetEta->Fill(correctedJet.eta());
      mCorrJetPhi->Fill(correctedJet.phi());
      mCorrJetPt ->Fill(correctedJet.pt());
      if (correctedJet.pt() >= 40) {
    mCorrJetEta_Pt40->Fill(correctedJet.eta()); 
    mCorrJetPhi_Pt40->Fill(correctedJet.phi());
      }
      
      double ijetEta = recoJets[ijet].eta();
      double ijetPt  = recoJets[ijet].pt();
      if(isMiniAODJet){
    ijetPt=recoJets[ijet].pt()*(*patJets)[ijet].jecFactor("Uncorrected");
      }
      double ratio   = correctedJet.pt() / ijetPt;
      if(isMiniAODJet){
    ratio =1./(*patJets)[ijet].jecFactor("Uncorrected");
      }
      
      if      (fabs(ijetEta) < 1.5) mPtCorrOverReco_Pt_B->Fill(log10(ijetPt), ratio);
      else if (fabs(ijetEta) < 3.0) mPtCorrOverReco_Pt_E->Fill(log10(ijetPt), ratio);
      else if (fabs(ijetEta) < 6.0) mPtCorrOverReco_Pt_F->Fill(log10(ijetPt), ratio);
      
      if      (ijetPt <  40) mPtCorrOverReco_Eta_20_40   ->Fill(ijetEta, ratio);
      else if (ijetPt <  200) mPtCorrOverReco_Eta_40_200  ->Fill(ijetEta, ratio);
      else if (ijetPt <  600) mPtCorrOverReco_Eta_200_600  ->Fill(ijetEta, ratio);
      else if (ijetPt < 1500) mPtCorrOverReco_Eta_600_1500 ->Fill(ijetEta, ratio);
      else if (ijetPt < 3500) mPtCorrOverReco_Eta_1500_3500->Fill(ijetEta, ratio);
      else if (ijetPt < 5000) mPtCorrOverReco_Eta_3500_5000->Fill(ijetEta, ratio);
      else if (ijetPt < 6500) mPtCorrOverReco_Eta_5000_6500->Fill(ijetEta, ratio);
      if (ijetPt > 3500) mPtCorrOverReco_Eta_3500->Fill(ijetEta, ratio);
    }
  }

  //----------------------------------------------------------------------------
  //
  // Generation
  //
  //----------------------------------------------------------------------------
  if (!mEvent.isRealData())
    {
      // Get ptHat
      //------------------------------------------------------------------------
      edm::Handle<GenEventInfoProduct> myGenEvt;
      mEvent.getByToken(evtToken_, myGenEvt);

      if (myGenEvt.isValid()) {
        if(myGenEvt->hasBinningValues()){
      double ptHat = myGenEvt->binningValues()[0];
      if (mPtHat) mPtHat->Fill(ptHat);
        }
      }
      // Gen jets
      //------------------------------------------------------------------------
      edm::Handle<GenJetCollection> genJets;
      mEvent.getByToken(genJetsToken_, genJets);

      if (!genJets.isValid()) return;
      
      for (GenJetCollection::const_iterator gjet=genJets->begin();  gjet!=genJets->end(); gjet++)   {
    //for MiniAOD we have here intrinsic thresholds, introduce also threshold for RECO
    if(gjet->pt() > mMatchGenPtThreshold){
      if (mGenEta) mGenEta->Fill(gjet->eta());
      if (mGenPhi) mGenPhi->Fill(gjet->phi());
      if (mGenPt)  mGenPt ->Fill(gjet->pt());
      if (gjet == genJets->begin()) {
        if (mGenEtaFirst) mGenEtaFirst->Fill(gjet->eta());
        if (mGenPhiFirst) mGenPhiFirst->Fill(gjet->phi());
      }
    }
      }

      if (!(mInputGenCollection.label().empty())) {
      for (GenJetCollection::const_iterator gjet=genJets->begin(); gjet!=genJets->end(); gjet++) {
        if (fabs(gjet->eta()) > 6.) continue;  // Out of the detector 
        if (gjet->pt() < mMatchGenPtThreshold) continue;
        if (recoJets.size() <= 0) continue;
        // pt response
        //------------------------------------------------------------
    int iMatch    =   -1;
    double CorrdeltaRBest = 999;
    double CorrJetPtBest  =   0;
    for (unsigned ijet=0; ijet<recoJets.size(); ++ijet) {
      Jet correctedJet = recoJets[ijet];
      if(pass_correction_flag && !isMiniAODJet){
        if (isCaloJet) scale = jetCorr->correction((*caloJets)[ijet]); 
        if (isPFJet)   scale = jetCorr->correction((*pfJets)[ijet]); 
        correctedJet.scaleEnergy(scale);
      }
      double CorrJetPt = correctedJet.pt();
      if (CorrJetPt > 10) {
        double CorrdR = deltaR(gjet->eta(), gjet->phi(), correctedJet.eta(), correctedJet.phi());
        if (CorrdR < CorrdeltaRBest) {
          CorrdeltaRBest = CorrdR;
          CorrJetPtBest  = CorrJetPt;
          iMatch = ijet;
        }
      }
    }
    if (iMatch<0) continue;
    if(!isMiniAODJet){
        fillMatchHists(gjet->eta(),  gjet->phi(),  gjet->pt(), recoJets[iMatch].eta(), recoJets[iMatch].phi(),  recoJets[iMatch].pt());
      }else{
      fillMatchHists(gjet->eta(),  gjet->phi(),  gjet->pt(), (*patJets)[iMatch].eta(), (*patJets)[iMatch].phi(),(*patJets)[iMatch].pt()*(*patJets)[iMatch].jecFactor("Uncorrected"));
    }
        if (pass_correction_flag) {//fill only for corrected jets
      if (CorrdeltaRBest < mRThreshold) {
            double response = CorrJetPtBest / gjet->pt();
            
            if      (fabs(gjet->eta()) < 1.5) mPtCorrOverGen_GenPt_B->Fill(log10(gjet->pt()), response);
            else if (fabs(gjet->eta()) < 3.0) mPtCorrOverGen_GenPt_E->Fill(log10(gjet->pt()), response);   
            else if (fabs(gjet->eta()) < 6.0) mPtCorrOverGen_GenPt_F->Fill(log10(gjet->pt()), response);
            
            if (gjet->pt() > 20) {
              if      (gjet->pt() <  40) mPtCorrOverGen_GenEta_20_40   ->Fill(gjet->eta(), response);
              else if (gjet->pt() <  200) mPtCorrOverGen_GenEta_40_200   ->Fill(gjet->eta(), response);
              else if (gjet->pt() <  600) mPtCorrOverGen_GenEta_200_600  ->Fill(gjet->eta(), response);
              else if (gjet->pt() < 1500) mPtCorrOverGen_GenEta_600_1500 ->Fill(gjet->eta(), response);
              else if (gjet->pt() < 3500) mPtCorrOverGen_GenEta_1500_3500->Fill(gjet->eta(), response);
              else if (gjet->pt() < 5000) mPtCorrOverGen_GenEta_3500_5000->Fill(gjet->eta(), response);
              else if (gjet->pt() < 6500) mPtCorrOverGen_GenEta_5000_6500->Fill(gjet->eta(), response);
              if (gjet->pt() > 3500) mPtCorrOverGen_GenEta_3500->Fill(gjet->eta(), response);
            }
          }
        }
      }
    }
    }
}

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 258 of file JetTester.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), chargedEmEnergy, chargedEmEnergyFraction, chargedHadronEnergy, chargedHadronEnergyFraction, chargedHadronMultiplicity, chargedMuEnergy, chargedMuEnergyFraction, chargedMultiplicity, electronEnergy, electronEnergyFraction, electronMultiplicity, emEnergyFraction, emEnergyInEB, emEnergyInEE, emEnergyInHF, energyFractionHadronic, MuonErrorMatrixValues_cff::etaRange, hadEnergyInHB, hadEnergyInHE, hadEnergyInHF, hadEnergyInHO, HFEMEnergy, HFEMEnergyFraction, HFEMMultiplicity, HFHadronEnergy, HFHadronEnergyFraction, HFHadronMultiplicity, HOEnergy, HOEnergyFraction, isCaloJet, isMiniAODJet, isPFJet, edm::InputTag::label(), maxEInEmTowers, maxEInHadTowers, mConstituents, mCorrJetEta, mCorrJetEta_Pt40, mCorrJetPhi, mCorrJetPhi_Pt40, mCorrJetPt, mDeltaEta, mDeltaPhi, mDeltaPt, mEnergy, mEta, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mInputCollection, mJetArea, mJetCorrector, mMass, mMjj, mNJets1, mNJets2, mNJets_40, mNJetsEta_B_20_40, mNJetsEta_B_40, mNJetsEta_E_20_40, mNJetsEta_E_40, mNvtx, mP, mPhi, mPhiFirst, mPt, mPtCorrOverGen_GenEta_1500_3500, mPtCorrOverGen_GenEta_200_600, mPtCorrOverGen_GenEta_20_40, mPtCorrOverGen_GenEta_3500, mPtCorrOverGen_GenEta_3500_5000, mPtCorrOverGen_GenEta_40_200, mPtCorrOverGen_GenEta_5000_6500, mPtCorrOverGen_GenEta_600_1500, mPtCorrOverGen_GenPt_B, mPtCorrOverGen_GenPt_E, mPtCorrOverGen_GenPt_F, mPtCorrOverReco_Eta_1500_3500, mPtCorrOverReco_Eta_200_600, mPtCorrOverReco_Eta_20_40, mPtCorrOverReco_Eta_3500, mPtCorrOverReco_Eta_3500_5000, mPtCorrOverReco_Eta_40_200, mPtCorrOverReco_Eta_5000_6500, mPtCorrOverReco_Eta_600_1500, mPtCorrOverReco_Pt_B, mPtCorrOverReco_Pt_E, mPtCorrOverReco_Pt_F, mPtFirst, mPtHat, mPtRecoOverGen_B_1500_3500, mPtRecoOverGen_B_200_600, mPtRecoOverGen_B_20_40, mPtRecoOverGen_B_3500, mPtRecoOverGen_B_3500_5000, mPtRecoOverGen_B_40_200, mPtRecoOverGen_B_5000_6500, mPtRecoOverGen_B_600_1500, mPtRecoOverGen_E_1500_3500, mPtRecoOverGen_E_200_600, mPtRecoOverGen_E_20_40, mPtRecoOverGen_E_3500, mPtRecoOverGen_E_3500_5000, mPtRecoOverGen_E_40_200, mPtRecoOverGen_E_5000_6500, mPtRecoOverGen_E_600_1500, mPtRecoOverGen_F_1500_3500, mPtRecoOverGen_F_200_600, mPtRecoOverGen_F_20_40, mPtRecoOverGen_F_3500, mPtRecoOverGen_F_40_200, mPtRecoOverGen_F_600_1500, mPtRecoOverGen_GenEta_1500_3500, mPtRecoOverGen_GenEta_200_600, mPtRecoOverGen_GenEta_20_40, mPtRecoOverGen_GenEta_3500, mPtRecoOverGen_GenEta_3500_5000, mPtRecoOverGen_GenEta_40_200, mPtRecoOverGen_GenEta_5000_6500, mPtRecoOverGen_GenEta_600_1500, mPtRecoOverGen_GenPhi_B, mPtRecoOverGen_GenPhi_E, mPtRecoOverGen_GenPhi_F, mPtRecoOverGen_GenPt_B, mPtRecoOverGen_GenPt_E, mPtRecoOverGen_GenPt_F, muonEnergy, muonEnergyFraction, muonMultiplicity, n60, n90, neutralEmEnergy, neutralEmEnergyFraction, neutralHadronEnergy, neutralHadronEnergyFraction, neutralHadronMultiplicity, neutralMultiplicity, photonEnergy, photonEnergyFraction, photonMultiplicity, DQMStore::IBooker::setCurrentFolder(), and towersArea.

                                                          { 

  ibooker.setCurrentFolder("JetMET/JetValidation/"+mInputCollection.label());  

    double log10PtMin  = 0.50;
    double log10PtMax  = 3.75;
    int    log10PtBins = 26; 

    //if eta range changed here need change in JetTesterPostProcessor as well
    double etaRange[91] = {-6.0, -5.8, -5.6, -5.4, -5.2, -5.0, -4.8, -4.6, -4.4, -4.2,
             -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.9, -2.8, -2.7, -2.6,
             -2.5, -2.4, -2.3, -2.2, -2.1, -2.0, -1.9, -1.8, -1.7, -1.6,
             -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.7, -0.6,
             -0.5, -0.4, -0.3, -0.2, -0.1,
             0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
             1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9,
             2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9,
             3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8,
             5.0, 5.2, 5.4, 5.6, 5.8, 6.0};

    // Event variables
    mNvtx           = ibooker.book1D("Nvtx",           "number of vertices", 60, 0, 60);

    // Jet parameters
    mEta          = ibooker.book1D("Eta",          "Eta",          120,   -6,    6); 
    mPhi          = ibooker.book1D("Phi",          "Phi",           70, -3.5,  3.5); 
    mPt           = ibooker.book1D("Pt",           "Pt",           100,    0,  1000); 
    mP            = ibooker.book1D("P",            "P",            100,    0,  1000); 
    mEnergy       = ibooker.book1D("Energy",       "Energy",       100,    0,  1000); 
    mMass         = ibooker.book1D("Mass",         "Mass",         100,    0,  200); 
    mConstituents = ibooker.book1D("Constituents", "Constituents", 100,    0,  100); 
    mJetArea      = ibooker.book1D("JetArea",      "JetArea",       100,   0, 4);

    // Corrected jets
    if (isMiniAODJet || !mJetCorrector.label().empty()) {//if correction label is filled, but fill also for MiniAOD though
      mCorrJetPt  = ibooker.book1D("CorrJetPt",  "CorrJetPt",  150,    0, 1500);
      mCorrJetEta = ibooker.book1D("CorrJetEta", "CorrJetEta Pt>20", 60,   -6,   6);
      mCorrJetPhi = ibooker.book1D("CorrJetPhi", "CorrJetPhi Pt>20",  70, -3.5, 3.5);
      mCorrJetEta_Pt40 = ibooker.book1D("CorrJetEta_Pt40", "CorrJetEta Pt>40", 60,   -6,   6);
      mCorrJetPhi_Pt40 = ibooker.book1D("CorrJetPhi_Pt40", "CorrJetPhi Pt>40",  70, -3.5, 3.5);

      // Corrected jets profiles
      mPtCorrOverReco_Pt_B = ibooker.bookProfile("PtCorrOverReco_Pt_B", "0<|eta|<1.5", log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");
      mPtCorrOverReco_Pt_E = ibooker.bookProfile("PtCorrOverReco_Pt_E", "1.5<|eta|<3", log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");
      mPtCorrOverReco_Pt_F = ibooker.bookProfile("PtCorrOverReco_Pt_F", "3<|eta|<6",   log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");

      mPtCorrOverReco_Eta_20_40    = ibooker.bookProfile("PtCorrOverReco_Eta_20_40",    "20<genPt<40",    90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_40_200    = ibooker.bookProfile("PtCorrOverReco_Eta_40_200",    "40<genPt<200",    90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_200_600   = ibooker.bookProfile("PtCorrOverReco_Eta_200_600",   "200<genPt<600",   90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_600_1500  = ibooker.bookProfile("PtCorrOverReco_Eta_600_1500",  "600<genPt<1500",  90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_1500_3500 = ibooker.bookProfile("PtCorrOverReco_Eta_1500_3500", "1500<genPt<3500", 90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_3500_5000 = ibooker.bookProfile("PtCorrOverReco_Eta_3500_5000", "3500<genPt<5000", 90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_5000_6500 = ibooker.bookProfile("PtCorrOverReco_Eta_5000_6500", "5000<genPt<6500", 90, etaRange, 0, 5, " ");
      mPtCorrOverReco_Eta_3500      = ibooker.bookProfile("PtCorrOverReco_Eta_3500",      "genPt>3500",      90, etaRange, 0, 5, " ");

      mPtCorrOverGen_GenPt_B = ibooker.bookProfile("PtCorrOverGen_GenPt_B", "0<|eta|<1.5", log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
      mPtCorrOverGen_GenPt_E = ibooker.bookProfile("PtCorrOverGen_GenPt_E", "1.5<|eta|<3", log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
      mPtCorrOverGen_GenPt_F = ibooker.bookProfile("PtCorrOverGen_GenPt_F", "3<|eta|<6",   log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
      //if eta range changed here need change in JetTesterPostProcessor as well
      mPtCorrOverGen_GenEta_20_40    = ibooker.bookProfile("PtCorrOverGen_GenEta_20_40",      "20<genPt<40;#eta",    90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_40_200    = ibooker.bookProfile("PtCorrOverGen_GenEta_40_200",    "40<genPt<200;#eta",    90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_200_600   = ibooker.bookProfile("PtCorrOverGen_GenEta_200_600",   "200<genPt<600;#eta",   90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_600_1500  = ibooker.bookProfile("PtCorrOverGen_GenEta_600_1500",  "600<genPt<1500;#eta",  90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_1500_3500 = ibooker.bookProfile("PtCorrOverGen_GenEta_1500_3500", "1500<genPt<3500;#eta", 90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_3500_5000 = ibooker.bookProfile("PtCorrOverGen_GenEta_3500_5000", "3500<genPt<5000;#eta", 90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_5000_6500 = ibooker.bookProfile("PtCorrOverGen_GenEta_5000_6500", "5000<genPt<6500;#eta", 90, etaRange, 0.8, 1.2, " ");
      mPtCorrOverGen_GenEta_3500      = ibooker.bookProfile("PtCorrOverGen_GenEta_3500",      "genPt>3500;#eta",      90, etaRange, 0.8, 1.2, " ");
    }

    mGenEta      = ibooker.book1D("GenEta",      "GenEta",      120,   -6,    6);
    mGenPhi      = ibooker.book1D("GenPhi",      "GenPhi",       70, -3.5,  3.5);
    mGenPt       = ibooker.book1D("GenPt",       "GenPt",       100,    0,  1000);
    mGenEtaFirst = ibooker.book1D("GenEtaFirst", "GenEtaFirst", 120,   -6,    6);
    mGenPhiFirst = ibooker.book1D("GenPhiFirst", "GenPhiFirst",  70, -3.5,  3.5);
    mPtHat       = ibooker.book1D("PtHat",       "PtHat",       100,    0, 1000); 
    mDeltaEta    = ibooker.book1D("DeltaEta",    "DeltaEta",    100, -0.5,  0.5);
    mDeltaPhi    = ibooker.book1D("DeltaPhi",    "DeltaPhi",    100, -0.5,  0.5);
    mDeltaPt     = ibooker.book1D("DeltaPt",     "DeltaPt",     100, -1.0,  1.0);
    
    mPtRecoOverGen_B_20_40    = ibooker.book1D("PtRecoOverGen_B_20_40",    "20<genpt<40",    90, 0, 2);
    mPtRecoOverGen_E_20_40    = ibooker.book1D("PtRecoOverGen_E_20_40",    "20<genpt<40",    90, 0, 2);
    mPtRecoOverGen_F_20_40    = ibooker.book1D("PtRecoOverGen_F_20_40",    "20<genpt<40",    90, 0, 2);
    mPtRecoOverGen_B_40_200    = ibooker.book1D("PtRecoOverGen_B_40_200",    "40<genpt<200",    90, 0, 2);
    mPtRecoOverGen_E_40_200    = ibooker.book1D("PtRecoOverGen_E_40_200",    "40<genpt<200",    90, 0, 2);
    mPtRecoOverGen_F_40_200    = ibooker.book1D("PtRecoOverGen_F_40_200",    "40<genpt<200",    90, 0, 2);
    mPtRecoOverGen_B_200_600   = ibooker.book1D("PtRecoOverGen_B_200_600",   "200<genpt<600",   90, 0, 2);
    mPtRecoOverGen_E_200_600   = ibooker.book1D("PtRecoOverGen_E_200_600",   "200<genpt<600",   90, 0, 2);
    mPtRecoOverGen_F_200_600   = ibooker.book1D("PtRecoOverGen_F_200_600",   "200<genpt<600",   90, 0, 2);
    mPtRecoOverGen_B_600_1500  = ibooker.book1D("PtRecoOverGen_B_600_1500",  "600<genpt<1500",  90, 0, 2);
    mPtRecoOverGen_E_600_1500  = ibooker.book1D("PtRecoOverGen_E_600_1500",  "600<genpt<1500",  90, 0, 2);
    mPtRecoOverGen_F_600_1500  = ibooker.book1D("PtRecoOverGen_F_600_1500",  "600<genpt<1500",  90, 0, 2);
    mPtRecoOverGen_B_1500_3500 = ibooker.book1D("PtRecoOverGen_B_1500_3500", "1500<genpt<3500", 90, 0, 2);
    mPtRecoOverGen_E_1500_3500 = ibooker.book1D("PtRecoOverGen_E_1500_3500", "1500<genpt<3500", 90, 0, 2);
    mPtRecoOverGen_F_1500_3500 = ibooker.book1D("PtRecoOverGen_F_1500_3500", "1500<genpt<3500", 90, 0, 2);
    mPtRecoOverGen_B_3500_5000 = ibooker.book1D("PtRecoOverGen_B_3500_5000", "3500<genpt<5000", 90, 0, 2);
    mPtRecoOverGen_E_3500_5000 = ibooker.book1D("PtRecoOverGen_E_3500_5000", "3500<genpt<5000", 90, 0, 2);
    mPtRecoOverGen_B_5000_6500 = ibooker.book1D("PtRecoOverGen_B_5000_6500", "5000<genpt<6500", 90, 0, 2);
    mPtRecoOverGen_E_5000_6500 = ibooker.book1D("PtRecoOverGen_E_5000_6500", "5000<genpt<6500", 90, 0, 2);
    mPtRecoOverGen_B_3500      = ibooker.book1D("PtRecoOverGen_B_3500",      "genpt>3500",      90, 0, 2);
    mPtRecoOverGen_E_3500      = ibooker.book1D("PtRecoOverGen_E_3500",      "genpt>3500",      90, 0, 2);
    mPtRecoOverGen_F_3500      = ibooker.book1D("PtRecoOverGen_F_3500",      "genpt>3500",      90, 0, 2);

    // Generation profiles
    mPtRecoOverGen_GenPt_B          = ibooker.bookProfile("PtRecoOverGen_GenPt_B",          "0<|eta|<1.5",     log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mPtRecoOverGen_GenPt_E          = ibooker.bookProfile("PtRecoOverGen_GenPt_E",          "1.5<|eta|<3",     log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mPtRecoOverGen_GenPt_F          = ibooker.bookProfile("PtRecoOverGen_GenPt_F",          "3<|eta|<6",       log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mPtRecoOverGen_GenPhi_B         = ibooker.bookProfile("PtRecoOverGen_GenPhi_B",         "0<|eta|<1.5",     70, -3.5, 3.5, 0, 2, " ");
    mPtRecoOverGen_GenPhi_E         = ibooker.bookProfile("PtRecoOverGen_GenPhi_E",         "1.5<|eta|<3",     70, -3.5, 3.5, 0, 2, " ");
    mPtRecoOverGen_GenPhi_F         = ibooker.bookProfile("PtRecoOverGen_GenPhi_F",         "3<|eta|<6",       70, -3.5, 3.5, 0, 2, " ");
    //if eta range changed here need change in JetTesterPostProcessor as well
    mPtRecoOverGen_GenEta_20_40    = ibooker.bookProfile("PtRecoOverGen_GenEta_20_40",    "20<genpt<40",    90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_40_200    = ibooker.bookProfile("PtRecoOverGen_GenEta_40_200",    "40<genpt<200",    90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_200_600   = ibooker.bookProfile("PtRecoOverGen_GenEta_200_600",   "200<genpt<600",   90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_600_1500  = ibooker.bookProfile("PtRecoOverGen_GenEta_600_1500",  "600<genpt<1500",  90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_1500_3500 = ibooker.bookProfile("PtRecoOverGen_GenEta_1500_3500", "1500<genpt<3500", 90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_3500_5000 = ibooker.bookProfile("PtRecoOverGen_GenEta_3500_5000", "3500<genpt<5000", 90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_5000_6500 = ibooker.bookProfile("PtRecoOverGen_GenEta_5000_6500", "5000<genpt<6500", 90, etaRange, 0, 2, " ");
    mPtRecoOverGen_GenEta_3500      = ibooker.bookProfile("PtRecoOverGen_GenEta_3500",      "genpt>3500",      90, etaRange, 0, 2, " "); 
  
    // Some jet algebra
    //------------------------------------------------------------------------
    mEtaFirst   = ibooker.book1D("EtaFirst",   "EtaFirst",   120,   -6,    6); 
    mPhiFirst   = ibooker.book1D("PhiFirst",   "PhiFirst",    70, -3.5,  3.5);      
    mPtFirst    = ibooker.book1D("PtFirst",    "PtFirst",    50,    0,  1000); 
    mMjj        = ibooker.book1D("Mjj",        "Mjj",        100,    0, 2000); 
    mNJetsEta_B_20_40 = ibooker.book1D("NJetsEta_B_20_40", "NJetsEta_B 20<Pt<40",  15,    0,   15);
    mNJetsEta_E_20_40 = ibooker.book1D("NJetsEta_E_20_40", "NJetsEta_E 20<Pt<40",  15,    0,   15);
    mNJetsEta_B_40 = ibooker.book1D("NJetsEta_B", "NJetsEta_B 40<Pt",  15,    0,   15);
    mNJetsEta_E_40 = ibooker.book1D("NJetsEta_E", "NJetsEta_E 40<Pt",  15,    0,   15);
    mNJets_40 = ibooker.book1D("NJets", "NJets 40>Pt",  15,    0,   15);
    mNJets1 = ibooker.bookProfile("NJets1", "Number of jets above Pt threshold", 100, 0,  200, 100, 0, 50, "s");
    mNJets2 = ibooker.bookProfile("NJets2", "Number of jets above Pt threshold", 100, 0, 4000, 100, 0, 50, "s");


    if (isCaloJet) {
      maxEInEmTowers              = ibooker.book1D("maxEInEmTowers", "maxEInEmTowers", 50,0,500);
      maxEInHadTowers             = ibooker.book1D("maxEInHadTowers", "maxEInHadTowers", 50,0,500);
      energyFractionHadronic      = ibooker.book1D("energyFractionHadronic", "energyFractionHadronic", 50,0,1);
      emEnergyFraction            = ibooker.book1D("emEnergyFraction", "emEnergyFraction", 50,0,1);
      hadEnergyInHB               = ibooker.book1D("hadEnergyInHB", "hadEnergyInHB", 50,0,500);
      hadEnergyInHO               = ibooker.book1D("hadEnergyInHO", "hadEnergyInHO", 50,0,500);
      hadEnergyInHE               = ibooker.book1D("hadEnergyInHE", "hadEnergyInHE", 50,0,500);
      hadEnergyInHF               = ibooker.book1D("hadEnergyInHF", "hadEnergyInHF", 50,0,500);
      emEnergyInEB                = ibooker.book1D("emEnergyInEB", "emEnergyInEB", 50,0,500);
      emEnergyInEE                = ibooker.book1D("emEnergyInEE", "emEnergyInEE", 50,0,500);
      emEnergyInHF                = ibooker.book1D("emEnergyInHF", "emEnergyInHF", 50,0,500);
      towersArea                  = ibooker.book1D("towersArea", "towersArea", 50,0,1);
      n90                         = ibooker.book1D("n90", "n90", 30,0,30);
      n60                         = ibooker.book1D("n60", "n60", 30,0,30);
    }

    if (isPFJet || isMiniAODJet) {
      muonMultiplicity = ibooker.book1D("muonMultiplicity", "muonMultiplicity", 10,0,10);
      chargedMultiplicity = ibooker.book1D("chargedMultiplicity", "chargedMultiplicity", 100,0,100);
      chargedEmEnergy = ibooker.book1D("chargedEmEnergy", "chargedEmEnergy", 100,0,500);
      neutralEmEnergy = ibooker.book1D("neutralEmEnergy", "neutralEmEnergy", 100,0,500);
      chargedHadronEnergy = ibooker.book1D("chargedHadronEnergy", "chargedHadronEnergy", 100,0,500);
      neutralHadronEnergy = ibooker.book1D("neutralHadronEnergy", "neutralHadronEnergy", 100,0,500);
      chargedHadronEnergyFraction = ibooker.book1D("chargedHadronEnergyFraction", "chargedHadronEnergyFraction", 50,0,1);
      neutralHadronEnergyFraction = ibooker.book1D("neutralHadronEnergyFraction", "neutralHadronEnergyFraction", 50,0,1);
      chargedEmEnergyFraction = ibooker.book1D("chargedEmEnergyFraction", "chargedEmEnergyFraction", 50,0,1);
      neutralEmEnergyFraction = ibooker.book1D("neutralEmEnergyFraction", "neutralEmEnergyFraction", 50,0,1);
      photonEnergy = ibooker.book1D("photonEnergy", "photonEnergy", 50,0,500);
      photonEnergyFraction = ibooker.book1D("photonEnergyFraction", "photonEnergyFraction", 50,0,1);
      electronEnergy = ibooker.book1D("electronEnergy", "electronEnergy", 50,0,500);
      electronEnergyFraction = ibooker.book1D("electronEnergyFraction", "electronEnergyFraction", 50,0,1);
      muonEnergy = ibooker.book1D("muonEnergy", "muonEnergy", 50,0,500);
      muonEnergyFraction = ibooker.book1D("muonEnergyFraction", "muonEnergyFraction", 50,0,1);
      HFHadronEnergy = ibooker.book1D("HFHadronEnergy", "HFHadronEnergy", 50,0,500);
      HFHadronEnergyFraction = ibooker.book1D("HFHadronEnergyFraction", "HFHadronEnergyFraction", 50,0,1);
      HFEMEnergy = ibooker.book1D("HFEMEnergy", "HFEMEnergy", 50,0,500);
      HFEMEnergyFraction = ibooker.book1D("HFEMEnergyFraction", "HFEMEnergyFraction", 50,0,1);
      chargedHadronMultiplicity = ibooker.book1D("chargedHadronMultiplicity", "chargedHadronMultiplicity", 50,0,50);
      neutralHadronMultiplicity = ibooker.book1D("neutralHadronMultiplicity", "neutralHadronMultiplicity", 50,0,50);
      photonMultiplicity = ibooker.book1D("photonMultiplicity", "photonMultiplicity", 10,0,10);
      electronMultiplicity = ibooker.book1D("electronMultiplicity", "electronMultiplicity", 10,0,10);
      HFHadronMultiplicity = ibooker.book1D("HFHadronMultiplicity", "HFHadronMultiplicity", 50,0,50);
      HFEMMultiplicity = ibooker.book1D("HFEMMultiplicity", "HFEMMultiplicity", 50,0,50);
      chargedMuEnergy = ibooker.book1D("chargedMuEnergy", "chargedMuEnergy", 50,0,500);
      chargedMuEnergyFraction = ibooker.book1D("chargedMuEnergyFraction", "chargedMuEnergyFraction", 50,0,1);
      neutralMultiplicity = ibooker.book1D("neutralMultiplicity", "neutralMultiplicity", 50,0,50);
      HOEnergy = ibooker.book1D("HOEnergy", "HOEnergy", 50,0,500);
      HOEnergyFraction = ibooker.book1D("HOEnergyFraction", "HOEnergyFraction", 50,0,1);
    }
}

void JetTester::fillMatchHists ( const double  GenEta, const double  GenPhi, const double  GenPt, const double  RecoEta, const double  RecoPhi, const double  RecoPt  )

private

Definition at line 928 of file JetTester.cc.

References MonitorElement::Fill(), mDeltaEta, mDeltaPhi, mDeltaPt, mMatchGenPtThreshold, mPtRecoOverGen_B_1500_3500, mPtRecoOverGen_B_200_600, mPtRecoOverGen_B_20_40, mPtRecoOverGen_B_3500, mPtRecoOverGen_B_3500_5000, mPtRecoOverGen_B_40_200, mPtRecoOverGen_B_5000_6500, mPtRecoOverGen_B_600_1500, mPtRecoOverGen_E_1500_3500, mPtRecoOverGen_E_200_600, mPtRecoOverGen_E_20_40, mPtRecoOverGen_E_3500, mPtRecoOverGen_E_3500_5000, mPtRecoOverGen_E_40_200, mPtRecoOverGen_E_5000_6500, mPtRecoOverGen_E_600_1500, mPtRecoOverGen_F_1500_3500, mPtRecoOverGen_F_200_600, mPtRecoOverGen_F_20_40, mPtRecoOverGen_F_3500, mPtRecoOverGen_F_40_200, mPtRecoOverGen_F_600_1500, mPtRecoOverGen_GenEta_1500_3500, mPtRecoOverGen_GenEta_200_600, mPtRecoOverGen_GenEta_20_40, mPtRecoOverGen_GenEta_3500, mPtRecoOverGen_GenEta_3500_5000, mPtRecoOverGen_GenEta_40_200, mPtRecoOverGen_GenEta_5000_6500, mPtRecoOverGen_GenEta_600_1500, mPtRecoOverGen_GenPhi_B, mPtRecoOverGen_GenPhi_E, mPtRecoOverGen_GenPhi_F, mPtRecoOverGen_GenPt_B, mPtRecoOverGen_GenPt_E, and mPtRecoOverGen_GenPt_F.

Referenced by analyze().

{
  if (GenPt > mMatchGenPtThreshold) {
    mDeltaEta->Fill(GenEta - RecoEta);
    mDeltaPhi->Fill(GenPhi - RecoPhi);
    mDeltaPt ->Fill((GenPt - RecoPt) / GenPt);
  }

  if (fabs(GenEta) < 1.5)
    {
      mPtRecoOverGen_GenPt_B ->Fill(log10(GenPt),  RecoPt / GenPt);
      mPtRecoOverGen_GenPhi_B->Fill(GenPhi, RecoPt / GenPt);
    
      if (GenPt > 20 && GenPt < 40) mPtRecoOverGen_B_20_40   ->Fill(RecoPt / GenPt);
      else if (GenPt <  200)         mPtRecoOverGen_B_40_200  ->Fill(RecoPt / GenPt);
      else if (GenPt <  600)         mPtRecoOverGen_B_200_600  ->Fill(RecoPt / GenPt);
      else if (GenPt < 1500)         mPtRecoOverGen_B_600_1500 ->Fill(RecoPt / GenPt);
      else if (GenPt < 3500)         mPtRecoOverGen_B_1500_3500->Fill(RecoPt / GenPt);
      else if (GenPt < 5000)         mPtRecoOverGen_B_3500_5000->Fill(RecoPt / GenPt);
      else if (GenPt < 6500)         mPtRecoOverGen_B_5000_6500->Fill(RecoPt / GenPt);
      if (GenPt>3500)         mPtRecoOverGen_B_3500->Fill(RecoPt / GenPt);
    }
  else if (fabs(GenEta) < 3.0)
    {
      mPtRecoOverGen_GenPt_E ->Fill(log10(GenPt),  RecoPt / GenPt);
      mPtRecoOverGen_GenPhi_E->Fill(GenPhi, RecoPt / GenPt);
    
      if (GenPt > 20 && GenPt < 40) mPtRecoOverGen_E_20_40   ->Fill(RecoPt / GenPt);
      else if (GenPt <  200)         mPtRecoOverGen_E_40_200  ->Fill(RecoPt / GenPt);
      else if (GenPt <  600)         mPtRecoOverGen_E_200_600  ->Fill(RecoPt / GenPt);
      else if (GenPt < 1500)         mPtRecoOverGen_E_600_1500 ->Fill(RecoPt / GenPt);
      else if (GenPt < 3500)         mPtRecoOverGen_E_1500_3500->Fill(RecoPt / GenPt);
      else if (GenPt < 5000)         mPtRecoOverGen_E_3500_5000->Fill(RecoPt / GenPt);
      else if (GenPt < 6500)         mPtRecoOverGen_E_5000_6500->Fill(RecoPt / GenPt);
      if (GenPt>3500)         mPtRecoOverGen_E_3500->Fill(RecoPt / GenPt);
    }
  else if (fabs(GenEta) < 6.0)
    {
      mPtRecoOverGen_GenPt_F ->Fill (log10(GenPt),  RecoPt / GenPt);
      mPtRecoOverGen_GenPhi_F->Fill (GenPhi, RecoPt / GenPt);
    
      if (GenPt > 20 && GenPt < 40) mPtRecoOverGen_F_20_40   ->Fill(RecoPt / GenPt);
      else if (GenPt <  200)         mPtRecoOverGen_F_40_200  ->Fill(RecoPt / GenPt);
      else if (GenPt <  600)         mPtRecoOverGen_F_200_600  ->Fill(RecoPt / GenPt);
      else if (GenPt < 1500)         mPtRecoOverGen_F_600_1500 ->Fill(RecoPt / GenPt);
      else if (GenPt < 3500)         mPtRecoOverGen_F_1500_3500->Fill(RecoPt / GenPt);
      if (GenPt>3500)                mPtRecoOverGen_F_3500->Fill(RecoPt / GenPt);
    }

  if (GenPt > 20 && GenPt < 40) mPtRecoOverGen_GenEta_20_40   ->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt <  200)         mPtRecoOverGen_GenEta_40_200  ->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt <  600)         mPtRecoOverGen_GenEta_200_600  ->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt < 1500)         mPtRecoOverGen_GenEta_600_1500 ->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt < 3500)         mPtRecoOverGen_GenEta_1500_3500->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt < 5000)         mPtRecoOverGen_GenEta_3500_5000->Fill(GenEta, RecoPt / GenPt);
  else if (GenPt < 6500)         mPtRecoOverGen_GenEta_5000_6500->Fill(GenEta, RecoPt / GenPt);
  if (GenPt > 3500)              mPtRecoOverGen_GenEta_3500->Fill(GenEta, RecoPt / GenPt);
}

Member Data Documentation

edm::EDGetTokenT<reco::CaloJetCollection> JetTester::caloJetsToken_

private

Definition at line 64 of file JetTester.h.

Referenced by analyze(), and JetTester().

MonitorElement* JetTester::chargedEmEnergy

private

Definition at line 199 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::chargedEmEnergyFraction

private

Definition at line 205 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::chargedHadronEnergy

private

Definition at line 201 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::chargedHadronEnergyFraction

private

Definition at line 203 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::chargedHadronMultiplicity

private

Definition at line 219 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::chargedMuEnergy

private

Definition at line 225 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::chargedMuEnergyFraction

private

Definition at line 226 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::chargedMultiplicity

private

Definition at line 198 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::electronEnergy

private

Definition at line 211 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::electronEnergyFraction

private

Definition at line 212 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::electronMultiplicity

private

Definition at line 222 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::emEnergyFraction

private

Definition at line 184 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::emEnergyInEB

private

Definition at line 189 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::emEnergyInEE

private

Definition at line 190 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::emEnergyInHF

private

Definition at line 191 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::energyFractionHadronic

private

Definition at line 183 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

edm::EDGetTokenT<GenEventInfoProduct> JetTester::evtToken_

private

Definition at line 67 of file JetTester.h.

Referenced by analyze(), and JetTester().

edm::EDGetTokenT<reco::GenJetCollection> JetTester::genJetsToken_

private

Definition at line 66 of file JetTester.h.

Referenced by analyze(), and JetTester().

MonitorElement* JetTester::hadEnergyInHB

private

Definition at line 185 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::hadEnergyInHE

private

Definition at line 187 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::hadEnergyInHF

private

Definition at line 188 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::hadEnergyInHO

private

Definition at line 186 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFEMEnergy

private

Definition at line 217 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFEMEnergyFraction

private

Definition at line 218 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFEMMultiplicity

private

Definition at line 224 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFHadronEnergy

private

Definition at line 215 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFHadronEnergyFraction

private

Definition at line 216 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HFHadronMultiplicity

private

Definition at line 223 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HOEnergy

private

Definition at line 228 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::HOEnergyFraction

private

Definition at line 229 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

bool JetTester::isCaloJet

private

Definition at line 235 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

bool JetTester::isMiniAODJet

private

Definition at line 237 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

bool JetTester::isPFJet

private

Definition at line 236 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

edm::EDGetTokenT<reco::JetCorrector> JetTester::jetCorrectorToken_

private

Definition at line 69 of file JetTester.h.

Referenced by analyze(), and JetTester().

std::string JetTester::JetType

private

Definition at line 60 of file JetTester.h.

Referenced by JetTester().

MonitorElement* JetTester::maxEInEmTowers

private

Definition at line 181 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::maxEInHadTowers

private

Definition at line 182 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mConstituents

private

Definition at line 81 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mCorrJetEta

private

Definition at line 87 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mCorrJetEta_Pt40

private

Definition at line 89 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mCorrJetPhi

private

Definition at line 88 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mCorrJetPhi_Pt40

private

Definition at line 90 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mCorrJetPt

private

Definition at line 86 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mDeltaEta

private

Definition at line 122 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mDeltaPhi

private

Definition at line 123 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mDeltaPt

private

Definition at line 124 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mEnergy

private

Definition at line 79 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mEta

private

Definition at line 75 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mEtaFirst

private

Definition at line 168 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mGenEta

private

Definition at line 116 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mGenEtaFirst

private

Definition at line 119 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mGenPhi

private

Definition at line 117 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mGenPhiFirst

private

Definition at line 120 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mGenPt

private

Definition at line 118 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

edm::InputTag JetTester::mInputCollection

private

Definition at line 57 of file JetTester.h.

Referenced by bookHistograms(), and JetTester().

edm::InputTag JetTester::mInputGenCollection

private

Definition at line 58 of file JetTester.h.

Referenced by analyze(), and JetTester().

MonitorElement* JetTester::mJetArea

private

Definition at line 82 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

edm::InputTag JetTester::mJetCorrector

private

Definition at line 59 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mMass

private

Definition at line 80 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

double JetTester::mMatchGenPtThreshold

private

Definition at line 233 of file JetTester.h.

Referenced by analyze(), and fillMatchHists().

MonitorElement* JetTester::mMjj

private

Definition at line 171 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJets1

private

Definition at line 177 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJets2

private

Definition at line 178 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJets_40

private

Definition at line 176 of file JetTester.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* JetTester::mNJetsEta_B_20_40

private

Definition at line 172 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJetsEta_B_40

private

Definition at line 174 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJetsEta_E_20_40

private

Definition at line 173 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNJetsEta_E_40

private

Definition at line 175 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mNvtx

private

Definition at line 72 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mP

private

Definition at line 78 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPhi

private

Definition at line 76 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPhiFirst

private

Definition at line 169 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPt

private

Definition at line 77 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_1500_3500

private

Definition at line 110 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_200_600

private

Definition at line 108 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_20_40

private

Definition at line 106 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_3500

private

Definition at line 113 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_3500_5000

private

Definition at line 111 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_40_200

private

Definition at line 107 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_5000_6500

private

Definition at line 112 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenEta_600_1500

private

Definition at line 109 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenPt_B

private

Definition at line 103 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenPt_E

private

Definition at line 104 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverGen_GenPt_F

private

Definition at line 105 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_1500_3500

private

Definition at line 99 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_200_600

private

Definition at line 97 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_20_40

private

Definition at line 95 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_3500

private

Definition at line 102 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_3500_5000

private

Definition at line 100 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_40_200

private

Definition at line 96 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_5000_6500

private

Definition at line 101 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Eta_600_1500

private

Definition at line 98 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Pt_B

private

Definition at line 92 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Pt_E

private

Definition at line 93 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtCorrOverReco_Pt_F

private

Definition at line 94 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtFirst

private

Definition at line 170 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtHat

private

Definition at line 121 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_1500_3500

private

Definition at line 138 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_200_600

private

Definition at line 132 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_20_40

private

Definition at line 126 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_3500

private

Definition at line 146 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_3500_5000

private

Definition at line 142 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_40_200

private

Definition at line 129 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_5000_6500

private

Definition at line 144 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_B_600_1500

private

Definition at line 135 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_1500_3500

private

Definition at line 139 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_200_600

private

Definition at line 133 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_20_40

private

Definition at line 127 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_3500

private

Definition at line 147 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_3500_5000

private

Definition at line 143 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_40_200

private

Definition at line 130 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_5000_6500

private

Definition at line 145 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_E_600_1500

private

Definition at line 136 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_1500_3500

private

Definition at line 140 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_200_600

private

Definition at line 134 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_20_40

private

Definition at line 128 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_3500

private

Definition at line 148 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_40_200

private

Definition at line 131 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_F_600_1500

private

Definition at line 137 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_1500_3500

private

Definition at line 161 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_200_600

private

Definition at line 159 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_20_40

private

Definition at line 157 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_3500

private

Definition at line 164 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_3500_5000

private

Definition at line 162 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_40_200

private

Definition at line 158 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_5000_6500

private

Definition at line 163 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenEta_600_1500

private

Definition at line 160 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPhi_B

private

Definition at line 154 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPhi_E

private

Definition at line 155 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPhi_F

private

Definition at line 156 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPt_B

private

Definition at line 151 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPt_E

private

Definition at line 152 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

MonitorElement* JetTester::mPtRecoOverGen_GenPt_F

private

Definition at line 153 of file JetTester.h.

Referenced by bookHistograms(), fillMatchHists(), and JetTester().

double JetTester::mRecoJetPtThreshold

private

Definition at line 232 of file JetTester.h.

Referenced by analyze().

double JetTester::mRThreshold

private

Definition at line 234 of file JetTester.h.

Referenced by analyze().

MonitorElement* JetTester::muonEnergy

private

Definition at line 213 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::muonEnergyFraction

private

Definition at line 214 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::muonMultiplicity

private

Definition at line 197 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::n60

private

Definition at line 194 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::n90

private

Definition at line 193 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::neutralEmEnergy

private

Definition at line 200 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::neutralEmEnergyFraction

private

Definition at line 206 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::neutralHadronEnergy

private

Definition at line 202 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

MonitorElement* JetTester::neutralHadronEnergyFraction

private

Definition at line 204 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::neutralHadronMultiplicity

private

Definition at line 220 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::neutralMultiplicity

private

Definition at line 227 of file JetTester.h.

Referenced by analyze(), bookHistograms(), Jet.Jet::jetID(), and JetTester().

edm::EDGetTokenT<pat::JetCollection> JetTester::patJetsToken_

private

Definition at line 68 of file JetTester.h.

Referenced by analyze(), and JetTester().

edm::EDGetTokenT<reco::PFJetCollection> JetTester::pfJetsToken_

private

Definition at line 65 of file JetTester.h.

Referenced by analyze(), and JetTester().

MonitorElement* JetTester::photonEnergy

private

Definition at line 209 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::photonEnergyFraction

private

Definition at line 210 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

MonitorElement* JetTester::photonMultiplicity

private

Definition at line 221 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().

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

private

Definition at line 63 of file JetTester.h.

Referenced by analyze(), and JetTester().

MonitorElement* JetTester::towersArea

private

Definition at line 192 of file JetTester.h.

Referenced by analyze(), bookHistograms(), and JetTester().