MCEmbeddingValidationAnalyzer Class Reference

#include <MCEmbeddingValidationAnalyzer.h>

Inheritance diagram for MCEmbeddingValidationAnalyzer:

Classes
struct	electronDistributionExtra
struct	l1ExtraObjectDistributionT
struct	leptonDistributionT
struct	leptonEfficiencyT
struct	leptonL1TriggerEfficiencyT1T2
struct	metDistributionType
struct	metL1TriggerEfficiencyType
struct	plotEntryTypeEvtWeight
struct	plotEntryTypeL1ETM
struct	plotEntryTypeMuonRadCorrUncertainty
struct	tauDistributionExtra

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	beginJob ()
	MCEmbeddingValidationAnalyzer (const edm::ParameterSet &)
	~MCEmbeddingValidationAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	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
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
typedef std::vector < edm::InputTag >	vInputTag

Private Member Functions
template<typename T >
void	bookHistograms (std::vector< T * > collection, DQMStore &dqmStore)
template<typename T >
void	cleanCollection (std::vector< T * > collection)
std::string	dqmDirectory_full (const std::string &dqmSubDirectory)
template<typename T >
void	fillHistograms (std::vector< T * > collection, int numJets, const edm::Event &evt, double evtWeight)
template<typename T >
void	fillHistograms (std::vector< T * > collection, int numJets, const edm::Event &evt, const edm::EventSetup &es, double evtWeight)
void	setupElectronDistributionExtra (int, int, const edm::ParameterSet &, const std::string &, std::vector< electronDistributionExtra * > &)
template<typename T >
void	setupL1ExtraObjectDistribution (int, int, const edm::ParameterSet &, const std::string &, std::vector< l1ExtraObjectDistributionT< T > * > &)
template<typename T >
void	setupLeptonDistribution (int, int, const edm::ParameterSet &, const std::string &, std::vector< leptonDistributionT< T > * > &)
template<typename T >
void	setupLeptonEfficiency (int, int, const edm::ParameterSet &, const std::string &, std::vector< leptonEfficiencyT< T > * > &)
template<typename T1 , typename T2 >
void	setupLeptonL1TriggerEfficiency (int, int, const edm::ParameterSet &, const std::string &, std::vector< leptonL1TriggerEfficiencyT1T2< T1, T2 > * > &)
void	setupMEtDistribution (int, int, const edm::ParameterSet &, const std::string &, std::vector< metDistributionType * > &)
void	setupMEtL1TriggerEfficiency (int, int, const edm::ParameterSet &, const std::string &, std::vector< metL1TriggerEfficiencyType * > &)
void	setupTauDistributionExtra (int, int, const edm::ParameterSet &, const std::string &, std::vector< tauDistributionExtra * > &)

Private Attributes
std::string	dqmDirectory_
double	dRminSeparation_
edm::EDGetTokenT < EcalRecHitCollection >	ebRHToken_
edm::EDGetTokenT < EcalRecHitCollection >	eeRHToken_
std::vector < leptonDistributionT < pat::Electron > * >	electronDistributions_
std::vector < electronDistributionExtra * >	electronDistributionsExtra_
std::vector< leptonEfficiencyT < pat::Electron > * >	electronEfficiencies_
std::vector < leptonL1TriggerEfficiencyT1T2 < pat::Electron, l1extra::L1EmParticle > * >	electronL1TriggerEfficiencies_
std::vector < plotEntryTypeEvtWeight * >	evtWeightPlotEntries_
std::vector< leptonEfficiencyT < reco::GsfElectron > * >	gsfElectronEfficiencies_
MonitorElement *	histogramBeamSpotX_
MonitorElement *	histogramBeamSpotY_
MonitorElement *	histogramCorrJetEtaPtGt20_
MonitorElement *	histogramCorrJetEtaPtGt30_
MonitorElement *	histogramCorrJetPt_
MonitorElement *	histogramCorrJetPtAbsEta2_5to4_5_
MonitorElement *	histogramCorrJetPtAbsEtaLt2_5_
MonitorElement *	histogramEventCounter_
MonitorElement *	histogramGenCaloMEt_
MonitorElement *	histogramGenDeltaPhiLeg1Leg2_
MonitorElement *	histogramGenDiTauDecayAngle_
MonitorElement *	histogramGenDiTauEta_
MonitorElement *	histogramGenDiTauMass_
MonitorElement *	histogramGenDiTauPhi_
MonitorElement *	histogramGenDiTauPt_
MonitorElement *	histogramGenFilterEfficiency_
MonitorElement *	histogramGenLeg1Eta_
MonitorElement *	histogramGenLeg1Mt_
MonitorElement *	histogramGenLeg1Phi_
MonitorElement *	histogramGenLeg1Pt_
MonitorElement *	histogramGenLeg1X_
MonitorElement *	histogramGenLeg1XforGenLeg2X0_00to0_25_
MonitorElement *	histogramGenLeg1XforGenLeg2X0_25to0_50_
MonitorElement *	histogramGenLeg1XforGenLeg2X0_50to0_75_
MonitorElement *	histogramGenLeg1XforGenLeg2X0_75to1_00_
MonitorElement *	histogramGenLeg2Eta_
MonitorElement *	histogramGenLeg2Mt_
MonitorElement *	histogramGenLeg2Phi_
MonitorElement *	histogramGenLeg2Pt_
MonitorElement *	histogramGenLeg2X_
MonitorElement *	histogramGenLeg2XforGenLeg1X0_00to0_25_
MonitorElement *	histogramGenLeg2XforGenLeg1X0_25to0_50_
MonitorElement *	histogramGenLeg2XforGenLeg1X0_50to0_75_
MonitorElement *	histogramGenLeg2XforGenLeg1X0_75to1_00_
MonitorElement *	histogramGenPFMEt_
MonitorElement *	histogramGenTau1Eta_
MonitorElement *	histogramGenTau1Phi_
MonitorElement *	histogramGenTau1Pt_
MonitorElement *	histogramGenTau2Eta_
MonitorElement *	histogramGenTau2Phi_
MonitorElement *	histogramGenTau2Pt_
MonitorElement *	histogramGenVisDeltaPhiLeg1Leg2_
MonitorElement *	histogramGenVisDiTauEta_
MonitorElement *	histogramGenVisDiTauMass_
MonitorElement *	histogramGenVisDiTauPhi_
MonitorElement *	histogramGenVisDiTauPt_
MonitorElement *	histogramNumChargedPFCandsPtGt10_
MonitorElement *	histogramNumChargedPFCandsPtGt20_
MonitorElement *	histogramNumChargedPFCandsPtGt30_
MonitorElement *	histogramNumChargedPFCandsPtGt40_
MonitorElement *	histogramNumChargedPFCandsPtGt5_
MonitorElement *	histogramNumGlobalMuons_
MonitorElement *	histogramNumJetsCorrPtGt20_
MonitorElement *	histogramNumJetsCorrPtGt20AbsEta2_5to4_5_
MonitorElement *	histogramNumJetsCorrPtGt20AbsEtaLt2_5_
MonitorElement *	histogramNumJetsCorrPtGt30_
MonitorElement *	histogramNumJetsCorrPtGt30AbsEta2_5to4_5_
MonitorElement *	histogramNumJetsCorrPtGt30AbsEtaLt2_5_
MonitorElement *	histogramNumJetsRawPtGt20_
MonitorElement *	histogramNumJetsRawPtGt20AbsEta2_5to4_5_
MonitorElement *	histogramNumJetsRawPtGt20AbsEtaLt2_5_
MonitorElement *	histogramNumJetsRawPtGt30_
MonitorElement *	histogramNumJetsRawPtGt30AbsEta2_5to4_5_
MonitorElement *	histogramNumJetsRawPtGt30AbsEtaLt2_5_
MonitorElement *	histogramNumNeutralPFCandsPtGt10_
MonitorElement *	histogramNumNeutralPFCandsPtGt20_
MonitorElement *	histogramNumNeutralPFCandsPtGt30_
MonitorElement *	histogramNumNeutralPFCandsPtGt40_
MonitorElement *	histogramNumNeutralPFCandsPtGt5_
MonitorElement *	histogramNumPFMuons_
MonitorElement *	histogramNumRecVertices_
MonitorElement *	histogramNumRecVerticesWithBS_
MonitorElement *	histogramNumStandAloneMuons_
MonitorElement *	histogramNumTracksPtGt10_
MonitorElement *	histogramNumTracksPtGt20_
MonitorElement *	histogramNumTracksPtGt30_
MonitorElement *	histogramNumTracksPtGt40_
MonitorElement *	histogramNumTracksPtGt5_
MonitorElement *	histogramPhiRotLegMinus_
MonitorElement *	histogramPhiRotLegPlus_
MonitorElement *	histogramRawJetEtaPtGt20_
MonitorElement *	histogramRawJetEtaPtGt30_
MonitorElement *	histogramRawJetPt_
MonitorElement *	histogramRawJetPtAbsEta2_5to4_5_
MonitorElement *	histogramRawJetPtAbsEtaLt2_5_
MonitorElement *	histogramRecCaloMEtECAL_
MonitorElement *	histogramRecCaloMEtHCAL_
MonitorElement *	histogramRecCaloMEtHF_
MonitorElement *	histogramRecCaloMEtHO_
MonitorElement *	histogramRecCaloSumEtECAL_
MonitorElement *	histogramRecCaloSumEtHCAL_
MonitorElement *	histogramRecCaloSumEtHF_
MonitorElement *	histogramRecCaloSumEtHO_
MonitorElement *	histogramRecLeg1PFMt_
MonitorElement *	histogramRecLeg1X_
MonitorElement *	histogramRecLeg2PFMt_
MonitorElement *	histogramRecLeg2X_
MonitorElement *	histogramRecVertexWithBSx_
MonitorElement *	histogramRecVertexWithBSy_
MonitorElement *	histogramRecVertexWithBSz_
MonitorElement *	histogramRecVertexX_
MonitorElement *	histogramRecVertexY_
MonitorElement *	histogramRecVertexZ_
MonitorElement *	histogramRecVisDeltaPhiLeg1Leg2_
MonitorElement *	histogramRecVisDiTauEta_
MonitorElement *	histogramRecVisDiTauMass_
MonitorElement *	histogramRecVisDiTauPhi_
MonitorElement *	histogramRecVisDiTauPt_
MonitorElement *	histogramRotationAngleMatrix_
MonitorElement *	histogramRotationLegMinusDeltaR_
MonitorElement *	histogramRotationLegPlusDeltaR_
MonitorElement *	histogramSumGenParticlePt_
MonitorElement *	histogramSumGenParticlePt_charged_
MonitorElement *	histogramTheRecVertexX_
MonitorElement *	histogramTheRecVertexY_
MonitorElement *	histogramTheRecVertexZ_
MonitorElement *	histogramWarning_recMuonNearReplacedMuon_
MonitorElement *	histogramWarning_recPFCandNearReplacedMuon_
MonitorElement *	histogramWarning_recTrackNearReplacedMuon_
std::vector < l1ExtraObjectDistributionT < l1extra::L1JetParticle > * >	l1CentralJetDistributions_
std::vector < l1ExtraObjectDistributionT < l1extra::L1EmParticle > * >	l1ElectronDistributions_
std::vector< plotEntryTypeL1ETM * >	l1ETMplotEntries_
std::vector < l1ExtraObjectDistributionT < l1extra::L1JetParticle > * >	l1ForwardJetDistributions_
std::vector < l1ExtraObjectDistributionT < l1extra::L1MuonParticle > * >	l1MuonDistributions_
std::vector < l1ExtraObjectDistributionT < l1extra::L1JetParticle > * >	l1TauDistributions_
std::vector < metDistributionType * >	metDistributions_
std::vector < metL1TriggerEfficiencyType * >	metL1TriggerEfficiencies_
std::string	moduleLabel_
std::vector < leptonDistributionT < pat::Muon > * >	muonDistributions_
std::vector< leptonEfficiencyT < pat::Muon > * >	muonEfficiencies_
std::vector < leptonL1TriggerEfficiencyT1T2 < pat::Muon, l1extra::L1MuonParticle > * >	muonL1TriggerEfficiencies_
int	muonRadCorrUncertainty_maxWarnings_
int	muonRadCorrUncertainty_numWarnings_
std::vector < plotEntryTypeMuonRadCorrUncertainty * >	muonRadCorrUncertaintyPlotEntries_afterRad_
std::vector < plotEntryTypeMuonRadCorrUncertainty * >	muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_
std::vector < plotEntryTypeMuonRadCorrUncertainty * >	muonRadCorrUncertaintyPlotEntries_beforeRad_
double	replacedMuonPtThresholdHigh_
double	replacedMuonPtThresholdLow_
edm::InputTag	srcBeamSpot_
edm::InputTag	srcCaloTowers_
edm::InputTag	srcGenCaloMEt_
edm::InputTag	srcGenDiTaus_
edm::InputTag	srcGenFilterInfo_
edm::InputTag	srcGenLeg1_
edm::InputTag	srcGenLeg2_
edm::InputTag	srcGenParticles_
edm::InputTag	srcGenPFMEt_
edm::InputTag	srcL1ETM_
edm::InputTag	srcMuonRadCorrWeight_
edm::InputTag	srcMuonRadCorrWeightDown_
edm::InputTag	srcMuonRadCorrWeightUp_
edm::InputTag	srcMuonsAfterRad_
edm::InputTag	srcMuonsBeforeRad_
vInputTag	srcOtherWeights_
edm::InputTag	srcRecCaloMEt_
edm::InputTag	srcRecJets_
edm::InputTag	srcRecLeg1_
edm::InputTag	srcRecLeg2_
edm::InputTag	srcRecMuons_
edm::InputTag	srcRecPFCandidates_
edm::InputTag	srcRecPFMEt_
edm::InputTag	srcRecTracks_
edm::InputTag	srcRecVertices_
edm::InputTag	srcRecVerticesWithBS_
edm::InputTag	srcReplacedMuons_
edm::InputTag	srcTheRecVertex_
std::vector < leptonDistributionT < pat::Tau > * >	tauDistributions_
std::vector < tauDistributionExtra * >	tauDistributionsExtra_
std::vector< leptonEfficiencyT < pat::Tau > * >	tauEfficiencies_
int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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

Compare Ztautau events produced via MCEmbedding to Ztautau events produced via direct Monte Carlo production in terms of: o kinematic distributions of electrons, muons, taus and MET o electron, muon and tau reconstruction and identification efficiencies o MET resolution o trigger efficiencies

Author

Christian Veelken, LLR

Version

Revision:

1.20

Id:

MCEmbeddingValidationAnalyzer.h,v 1.20 2013/06/04 13:44:44 veelken Exp

Definition at line 197 of file MCEmbeddingValidationAnalyzer.h.

Member Typedef Documentation

typedef std::vector<edm::InputTag> MCEmbeddingValidationAnalyzer::vInputTag

private

Definition at line 244 of file MCEmbeddingValidationAnalyzer.h.

Constructor & Destructor Documentation

MCEmbeddingValidationAnalyzer::MCEmbeddingValidationAnalyzer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 80 of file MCEmbeddingValidationAnalyzer.cc.

References ebRHToken_, eeRHToken_, electronDistributions_, electronDistributionsExtra_, electronEfficiencies_, electronL1TriggerEfficiencies_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), gsfElectronEfficiencies_, l1CentralJetDistributions_, l1ElectronDistributions_, l1ForwardJetDistributions_, l1MuonDistributions_, l1TauDistributions_, metDistributions_, metL1TriggerEfficiencies_, muonDistributions_, muonEfficiencies_, muonL1TriggerEfficiencies_, setupElectronDistributionExtra(), setupL1ExtraObjectDistribution(), setupLeptonDistribution(), setupLeptonEfficiency(), setupLeptonL1TriggerEfficiency(), setupMEtDistribution(), setupMEtL1TriggerEfficiency(), setupTauDistributionExtra(), tauDistributions_, tauDistributionsExtra_, tauEfficiencies_, and verbosity_.

  : moduleLabel_(cfg.getParameter<std::string>("@module_label")),    
    srcReplacedMuons_(cfg.getParameter<edm::InputTag>("srcReplacedMuons")),
    srcRecMuons_(cfg.getParameter<edm::InputTag>("srcRecMuons")),
    srcRecTracks_(cfg.getParameter<edm::InputTag>("srcRecTracks")),
    srcCaloTowers_(cfg.getParameter<edm::InputTag>("srcCaloTowers")),
    srcRecPFCandidates_(cfg.getParameter<edm::InputTag>("srcRecPFCandidates")),
    srcRecJets_(cfg.getParameter<edm::InputTag>("srcRecJets")),
    srcTheRecVertex_(cfg.getParameter<edm::InputTag>("srcTheRecVertex")),
    srcRecVertices_(cfg.getParameter<edm::InputTag>("srcRecVertices")),
    srcRecVerticesWithBS_(cfg.getParameter<edm::InputTag>("srcRecVerticesWithBS")),
    srcBeamSpot_(cfg.getParameter<edm::InputTag>("srcBeamSpot")),
    srcGenDiTaus_(cfg.getParameter<edm::InputTag>("srcGenDiTaus")),
    dRminSeparation_(cfg.getParameter<double>("dRminSeparation")),
    srcGenLeg1_(cfg.getParameter<edm::InputTag>("srcGenLeg1")),
    srcRecLeg1_(cfg.getParameter<edm::InputTag>("srcRecLeg1")),
    srcGenLeg2_(cfg.getParameter<edm::InputTag>("srcGenLeg2")),
    srcRecLeg2_(cfg.getParameter<edm::InputTag>("srcRecLeg2")),
    srcGenParticles_(cfg.getParameter<edm::InputTag>("srcGenParticles")),
    srcL1ETM_(cfg.getParameter<edm::InputTag>("srcL1ETM")),
    srcGenCaloMEt_(cfg.getParameter<edm::InputTag>("srcGenCaloMEt")),
    srcGenPFMEt_(cfg.getParameter<edm::InputTag>("srcGenPFMEt")),
    srcRecCaloMEt_(cfg.getParameter<edm::InputTag>("srcRecCaloMEt")),
    srcRecPFMEt_(cfg.getParameter<edm::InputTag>("srcRecPFMEt")),
    srcMuonsBeforeRad_(cfg.getParameter<edm::InputTag>("srcMuonsBeforeRad")),
    srcMuonsAfterRad_(cfg.getParameter<edm::InputTag>("srcMuonsAfterRad")),
    srcMuonRadCorrWeight_(cfg.getParameter<edm::InputTag>("srcMuonRadCorrWeight")),
    srcMuonRadCorrWeightUp_(cfg.getParameter<edm::InputTag>("srcMuonRadCorrWeightUp")),
    srcMuonRadCorrWeightDown_(cfg.getParameter<edm::InputTag>("srcMuonRadCorrWeightDown")),
    srcOtherWeights_(cfg.getParameter<vInputTag>("srcOtherWeights")),
    srcGenFilterInfo_(cfg.getParameter<edm::InputTag>("srcGenFilterInfo")),
    dqmDirectory_(cfg.getParameter<std::string>("dqmDirectory")),
    replacedMuonPtThresholdHigh_(cfg.getParameter<double>("replacedMuonPtThresholdHigh")),
    replacedMuonPtThresholdLow_(cfg.getParameter<double>("replacedMuonPtThresholdLow"))
{
  typedef std::pair<int, int> pint;
  std::vector<pint> jetBins;
  jetBins.push_back(pint(-1, -1));
  jetBins.push_back(pint(0, 0));
  jetBins.push_back(pint(1, 1));
  jetBins.push_back(pint(2, 2));
  jetBins.push_back(pint(3, 1000));
  for ( std::vector<pint>::const_iterator jetBin = jetBins.begin();
    jetBin != jetBins.end(); ++jetBin ) {
//--- setup electron Pt, eta and phi distributions;
//    electron id & isolation and trigger efficiencies
    setupLeptonDistribution(jetBin->first, jetBin->second, cfg, "electronDistributions", electronDistributions_);
    setupElectronDistributionExtra(jetBin->first, jetBin->second, cfg, "electronDistributions", electronDistributionsExtra_);
    setupLeptonEfficiency(jetBin->first, jetBin->second, cfg, "electronEfficiencies", electronEfficiencies_);
    setupLeptonEfficiency(jetBin->first, jetBin->second, cfg, "gsfElectronEfficiencies", gsfElectronEfficiencies_);
    setupLeptonL1TriggerEfficiency(jetBin->first, jetBin->second, cfg, "electronL1TriggerEfficiencies", electronL1TriggerEfficiencies_);

//--- setup muon Pt, eta and phi distributions;
//    muon id & isolation and trigger efficiencies
    setupLeptonDistribution(jetBin->first, jetBin->second, cfg, "muonDistributions", muonDistributions_);
    setupLeptonEfficiency(jetBin->first, jetBin->second, cfg, "muonEfficiencies", muonEfficiencies_);
    setupLeptonL1TriggerEfficiency(jetBin->first, jetBin->second, cfg, "muonL1TriggerEfficiencies", muonL1TriggerEfficiencies_);
  
//--- setup tau Pt, eta and phi distributions;
//    tau id efficiency
    setupLeptonDistribution(jetBin->first, jetBin->second, cfg, "tauDistributions", tauDistributions_);
    setupTauDistributionExtra(jetBin->first, jetBin->second, cfg, "tauDistributions", tauDistributionsExtra_);
    setupLeptonEfficiency(jetBin->first, jetBin->second, cfg, "tauEfficiencies", tauEfficiencies_);

//--- setup Pt, eta and phi distributions of L1Extra objects
//   (electrons, muons, tau-jet, central and forward jets)
    setupL1ExtraObjectDistribution(jetBin->first, jetBin->second, cfg, "l1ElectronDistributions", l1ElectronDistributions_);
    setupL1ExtraObjectDistribution(jetBin->first, jetBin->second, cfg, "l1MuonDistributions", l1MuonDistributions_);
    setupL1ExtraObjectDistribution(jetBin->first, jetBin->second, cfg, "l1TauDistributions", l1TauDistributions_);
    setupL1ExtraObjectDistribution(jetBin->first, jetBin->second, cfg, "l1CentralJetDistributions", l1CentralJetDistributions_);
    setupL1ExtraObjectDistribution(jetBin->first, jetBin->second, cfg, "l1ForwardJetDistributions", l1ForwardJetDistributions_);

//--- setup MET Pt and phi distributions;
//    efficiency of L1 (Calo)MET trigger requirement
    setupMEtDistribution(jetBin->first, jetBin->second, cfg, "metDistributions", metDistributions_);
    setupMEtL1TriggerEfficiency(jetBin->first, jetBin->second, cfg, "metL1TriggerEfficiencies", metL1TriggerEfficiencies_);
  }

  ebRHToken_ = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
  eeRHToken_ = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));

  verbosity_ = ( cfg.exists("verbosity") ) ?
    cfg.getParameter<int>("verbosity") : 0;
}

MCEmbeddingValidationAnalyzer::~MCEmbeddingValidationAnalyzer

(

)

Definition at line 165 of file MCEmbeddingValidationAnalyzer.cc.

References cleanCollection(), electronDistributions_, electronDistributionsExtra_, electronEfficiencies_, electronL1TriggerEfficiencies_, evtWeightPlotEntries_, gsfElectronEfficiencies_, l1CentralJetDistributions_, l1ElectronDistributions_, l1ETMplotEntries_, l1ForwardJetDistributions_, l1MuonDistributions_, l1TauDistributions_, metDistributions_, metL1TriggerEfficiencies_, muonDistributions_, muonEfficiencies_, muonL1TriggerEfficiencies_, muonRadCorrUncertaintyPlotEntries_afterRad_, muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_, muonRadCorrUncertaintyPlotEntries_beforeRad_, tauDistributions_, tauDistributionsExtra_, and tauEfficiencies_.

{
  for ( std::vector<plotEntryTypeEvtWeight*>::iterator it = evtWeightPlotEntries_.begin();
    it != evtWeightPlotEntries_.end(); ++it ) {
    delete (*it);
  }

  for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator it = muonRadCorrUncertaintyPlotEntries_beforeRad_.begin();
    it != muonRadCorrUncertaintyPlotEntries_beforeRad_.end(); ++it ) {
    delete (*it);
  }
  for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator it = muonRadCorrUncertaintyPlotEntries_afterRad_.begin();
    it != muonRadCorrUncertaintyPlotEntries_afterRad_.end(); ++it ) {
    delete (*it);
  }
  for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator it = muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_.begin();
    it != muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_.end(); ++it ) {
    delete (*it);
  }

  for ( std::vector<plotEntryTypeL1ETM*>::iterator it = l1ETMplotEntries_.begin();
    it != l1ETMplotEntries_.end(); ++it ) {
    delete (*it);
  }

  cleanCollection(electronDistributions_);
  cleanCollection(electronDistributionsExtra_);
  cleanCollection(electronEfficiencies_);
  cleanCollection(gsfElectronEfficiencies_);
  cleanCollection(electronL1TriggerEfficiencies_);
  cleanCollection(muonDistributions_);
  cleanCollection(muonEfficiencies_);
  cleanCollection(muonL1TriggerEfficiencies_);
  cleanCollection(tauDistributions_);
  cleanCollection(tauDistributionsExtra_);
  cleanCollection(tauEfficiencies_);
  cleanCollection(l1ElectronDistributions_);
  cleanCollection(l1MuonDistributions_);
  cleanCollection(l1TauDistributions_);
  cleanCollection(l1CentralJetDistributions_);
  cleanCollection(l1ForwardJetDistributions_);
  cleanCollection(metDistributions_);
  cleanCollection(metL1TriggerEfficiencies_);
}

Member Function Documentation

void MCEmbeddingValidationAnalyzer::analyze ( const edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDAnalyzer.

Definition at line 777 of file MCEmbeddingValidationAnalyzer.cc.

References Abs(), SiPixelRawToDigiRegional_cfi::beamSpot, HiRecoJets_cff::caloTowers, reco::Candidate::charge(), gather_cfg::cout, reco::CompositeCandidate::daughter(), deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRminSeparation_, alignCSCRings::e, electronDistributions_, electronDistributionsExtra_, electronEfficiencies_, electronL1TriggerEfficiencies_, reco::Candidate::eta(), evtWeightPlotEntries_, MonitorElement::Fill(), fillHistograms(), findMuons(), genParticleCandidates2GenParticles_cfi::genParticles, reco::Candidate::get(), edm::Event::getByLabel(), getSelMuons(), gsfElectronEfficiencies_, reco::Candidate::hasMasterClone(), histogramBeamSpotX_, histogramBeamSpotY_, histogramCorrJetEtaPtGt20_, histogramCorrJetPt_, histogramCorrJetPtAbsEta2_5to4_5_, histogramCorrJetPtAbsEtaLt2_5_, histogramEventCounter_, histogramGenCaloMEt_, histogramGenDeltaPhiLeg1Leg2_, histogramGenDiTauDecayAngle_, histogramGenDiTauEta_, histogramGenDiTauMass_, histogramGenDiTauPhi_, histogramGenDiTauPt_, histogramGenFilterEfficiency_, histogramGenLeg1Eta_, histogramGenLeg1Mt_, histogramGenLeg1Phi_, histogramGenLeg1Pt_, histogramGenLeg1X_, histogramGenLeg1XforGenLeg2X0_00to0_25_, histogramGenLeg1XforGenLeg2X0_25to0_50_, histogramGenLeg1XforGenLeg2X0_50to0_75_, histogramGenLeg1XforGenLeg2X0_75to1_00_, histogramGenLeg2Eta_, histogramGenLeg2Mt_, histogramGenLeg2Phi_, histogramGenLeg2Pt_, histogramGenLeg2X_, histogramGenLeg2XforGenLeg1X0_00to0_25_, histogramGenLeg2XforGenLeg1X0_25to0_50_, histogramGenLeg2XforGenLeg1X0_50to0_75_, histogramGenLeg2XforGenLeg1X0_75to1_00_, histogramGenPFMEt_, histogramGenTau1Eta_, histogramGenTau1Phi_, histogramGenTau1Pt_, histogramGenTau2Eta_, histogramGenTau2Phi_, histogramGenTau2Pt_, histogramGenVisDeltaPhiLeg1Leg2_, histogramGenVisDiTauEta_, histogramGenVisDiTauMass_, histogramGenVisDiTauPhi_, histogramGenVisDiTauPt_, histogramNumChargedPFCandsPtGt10_, histogramNumChargedPFCandsPtGt20_, histogramNumChargedPFCandsPtGt30_, histogramNumChargedPFCandsPtGt40_, histogramNumChargedPFCandsPtGt5_, histogramNumGlobalMuons_, histogramNumJetsCorrPtGt20_, histogramNumJetsCorrPtGt20AbsEta2_5to4_5_, histogramNumJetsCorrPtGt20AbsEtaLt2_5_, histogramNumJetsCorrPtGt30_, histogramNumJetsCorrPtGt30AbsEta2_5to4_5_, histogramNumJetsCorrPtGt30AbsEtaLt2_5_, histogramNumJetsRawPtGt20_, histogramNumJetsRawPtGt20AbsEta2_5to4_5_, histogramNumJetsRawPtGt20AbsEtaLt2_5_, histogramNumJetsRawPtGt30_, histogramNumJetsRawPtGt30AbsEta2_5to4_5_, histogramNumJetsRawPtGt30AbsEtaLt2_5_, histogramNumNeutralPFCandsPtGt10_, histogramNumNeutralPFCandsPtGt20_, histogramNumNeutralPFCandsPtGt30_, histogramNumNeutralPFCandsPtGt40_, histogramNumNeutralPFCandsPtGt5_, histogramNumPFMuons_, histogramNumRecVertices_, histogramNumRecVerticesWithBS_, histogramNumStandAloneMuons_, histogramNumTracksPtGt10_, histogramNumTracksPtGt20_, histogramNumTracksPtGt30_, histogramNumTracksPtGt40_, histogramNumTracksPtGt5_, histogramPhiRotLegMinus_, histogramPhiRotLegPlus_, histogramRawJetEtaPtGt20_, histogramRawJetPt_, histogramRawJetPtAbsEta2_5to4_5_, histogramRawJetPtAbsEtaLt2_5_, histogramRecCaloMEtECAL_, histogramRecCaloMEtHCAL_, histogramRecCaloMEtHF_, histogramRecCaloMEtHO_, histogramRecCaloSumEtECAL_, histogramRecCaloSumEtHCAL_, histogramRecCaloSumEtHF_, histogramRecCaloSumEtHO_, histogramRecLeg1PFMt_, histogramRecLeg1X_, histogramRecLeg2PFMt_, histogramRecLeg2X_, histogramRecVertexWithBSx_, histogramRecVertexWithBSy_, histogramRecVertexWithBSz_, histogramRecVertexX_, histogramRecVertexY_, histogramRecVertexZ_, histogramRecVisDeltaPhiLeg1Leg2_, histogramRecVisDiTauEta_, histogramRecVisDiTauMass_, histogramRecVisDiTauPhi_, histogramRecVisDiTauPt_, histogramRotationAngleMatrix_, histogramRotationLegMinusDeltaR_, histogramRotationLegPlusDeltaR_, histogramSumGenParticlePt_, histogramTheRecVertexX_, histogramTheRecVertexY_, histogramTheRecVertexZ_, histogramWarning_recMuonNearReplacedMuon_, histogramWarning_recPFCandNearReplacedMuon_, histogramWarning_recTrackNearReplacedMuon_, edm::InputTag::instance(), edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, l1CentralJetDistributions_, l1ElectronDistributions_, l1ETMplotEntries_, l1ForwardJetDistributions_, l1MuonDistributions_, l1TauDistributions_, edm::InputTag::label(), reco::Candidate::masterClone(), metDistributions_, metL1TriggerEfficiencies_, moduleLabel_, metsig::muon, muonDistributions_, muonEfficiencies_, muonL1TriggerEfficiencies_, muonRadCorrUncertainty_maxWarnings_, muonRadCorrUncertainty_numWarnings_, muonRadCorrUncertaintyPlotEntries_afterRad_, muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_, muonRadCorrUncertaintyPlotEntries_beforeRad_, patZpeak::muons, normalizedPhi(), reco::CompositeCandidate::numberOfDaughters(), reco::Candidate::p4(), reco::tau::pfCandidates(), reco::Candidate::phi(), SVfit_namespace::phiLabFromLabMomenta(), reco::Candidate::pt(), replacedMuonPtThresholdHigh_, replacedMuonPtThresholdLow_, srcBeamSpot_, srcCaloTowers_, srcGenCaloMEt_, srcGenDiTaus_, srcGenFilterInfo_, srcGenLeg1_, srcGenLeg2_, srcGenParticles_, srcGenPFMEt_, srcL1ETM_, srcMuonRadCorrWeight_, srcMuonRadCorrWeightDown_, srcMuonRadCorrWeightUp_, srcMuonsAfterRad_, srcMuonsBeforeRad_, srcOtherWeights_, srcRecCaloMEt_, srcRecJets_, srcRecLeg1_, srcRecLeg2_, srcRecMuons_, srcRecPFCandidates_, srcRecPFMEt_, srcRecTracks_, srcRecVertices_, srcRecVerticesWithBS_, srcReplacedMuons_, srcTheRecVertex_, AlCaHLTBitMon_QueryRunRegistry::string, tauDistributions_, tauDistributionsExtra_, tauEfficiencies_, testEve_cfg::tracks, verbosity_, and histoStyle::weight.

{
  if ( verbosity_ ) {
    std::cout << "<MCEmbeddingValidationAnalyzer::analyze>:" << std::endl;
    std::cout << " moduleLabel = " << moduleLabel_ << std::endl;
  }

  const reco::Candidate* replacedMuonPlus  = 0;
  const reco::Candidate* replacedMuonMinus = 0;
  if ( srcReplacedMuons_.label() != "" ) {
    std::vector<reco::CandidateBaseRef> replacedMuons = getSelMuons(evt, srcReplacedMuons_);
    for ( std::vector<reco::CandidateBaseRef>::const_iterator replacedMuon = replacedMuons.begin();
      replacedMuon != replacedMuons.end(); ++replacedMuon ) {
      if      ( (*replacedMuon)->charge() > +0.5 ) replacedMuonPlus  = replacedMuon->get();
      else if ( (*replacedMuon)->charge() < -0.5 ) replacedMuonMinus = replacedMuon->get();
    }
  }
  if ( verbosity_ ) {
    if ( replacedMuonPlus ) std::cout << "replacedMuonPlus: Pt = " << replacedMuonPlus->pt() << ", eta = " << replacedMuonPlus->eta() << ", phi = " << replacedMuonPlus->phi() << std::endl;
    if ( replacedMuonMinus ) std::cout << "replacedMuonMinus: Pt = " << replacedMuonMinus->pt() << ", eta = " << replacedMuonMinus->eta() << ", phi = " << replacedMuonMinus->phi() << std::endl;
  }
  
  typedef edm::View<reco::Candidate> CandidateView;
  edm::Handle<CandidateView> genDiTaus;
  evt.getByLabel(srcGenDiTaus_, genDiTaus);

  const reco::Candidate* genTauPlus  = 0;
  const reco::Candidate* genTauMinus = 0;
  for ( CandidateView::const_iterator genDiTau = genDiTaus->begin();
    genDiTau != genDiTaus->end(); ++genDiTau ) {
    const reco::CompositeCandidate* genDiTau_composite = dynamic_cast<const reco::CompositeCandidate*>(&(*genDiTau));
    if ( !(genDiTau_composite && genDiTau_composite->numberOfDaughters() == 2) ) continue;
    size_t numDaughters = genDiTau_composite->numberOfDaughters();
    for ( size_t iDaughter = 0; iDaughter < numDaughters; ++iDaughter ) {
      const reco::Candidate* daughter = genDiTau_composite->daughter(iDaughter);
      if      ( daughter->charge() > +0.5 ) genTauPlus  = daughter;
      else if ( daughter->charge() < -0.5 ) genTauMinus = daughter;
    }
  }
  if ( verbosity_ ) {
    if ( genTauPlus ) std::cout << "genTauPlus: Pt = " << genTauPlus->pt() << ", eta = " << genTauPlus->eta() << ", phi = " << genTauPlus->phi() << std::endl;
    if ( genTauMinus ) std::cout << "genTauMinus: Pt = " << genTauMinus->pt() << ", eta = " << genTauMinus->eta() << ", phi = " << genTauMinus->phi() << std::endl;
  }

  double dRmin = 1.e+3;
  if ( replacedMuonPlus && genTauPlus ) {
    double dR = deltaR(genTauPlus->p4(),  replacedMuonPlus->p4());
    if ( verbosity_ ) std::cout << " dR(replacedMuonPlus, genTauPlus) = " << dR << std::endl;
    if ( dR < dRmin ) dRmin = dR;
  }
  if ( replacedMuonPlus && genTauMinus ) {
    double dR = deltaR(genTauMinus->p4(),  replacedMuonPlus->p4());
    if ( verbosity_ ) std::cout << " dR(replacedMuonPlus, genTauMinus) = " << dR << std::endl;
    if ( dR < dRmin ) dRmin = dR;
  }
  if ( replacedMuonMinus && genTauPlus ) {
    double dR = deltaR(genTauPlus->p4(),  replacedMuonMinus->p4());
    if ( verbosity_ ) std::cout << " dR(replacedMuonMinus, genTauPlus) = " << dR << std::endl;
    if ( dR < dRmin ) dRmin = dR;
  }
  if ( replacedMuonMinus && genTauMinus ) {
    double dR = deltaR(genTauMinus->p4(),  replacedMuonMinus->p4());
    if ( verbosity_ ) std::cout << " dR(replacedMuonMinus, genTauMinus) = " << dR << std::endl;
    if ( dR < dRmin ) dRmin = dR;
  }
  if ( verbosity_ ) std::cout << "--> dRmin = " << dRmin << std::endl;
  if ( dRmin < dRminSeparation_ ) return;
  if ( verbosity_ ) std::cout << " cut on dRminSeparation = " << dRminSeparation_ << " passed." << std::endl;

//--- compute event weight
  double evtWeight = 1.0;
  std::map<std::string, double> evtWeightMap;
  for ( vInputTag::const_iterator srcWeight = srcOtherWeights_.begin();
    srcWeight != srcOtherWeights_.end(); ++srcWeight ) {
    edm::Handle<double> weight;
    evt.getByLabel(*srcWeight, weight);
    //std::cout << "weight(" << srcWeight->label().data() << ":" << srcWeight->instance().data() << ") = " << (*weight) << std::endl;
    evtWeight *= (*weight);
    evtWeightMap[Form("%s_%s", srcWeight->label().data(), srcWeight->instance().data())] = (*weight);
  }
  if ( srcGenFilterInfo_.label() != "" ) {
    edm::Handle<GenFilterInfo> genFilterInfo;
    evt.getByLabel(srcGenFilterInfo_, genFilterInfo);
    //std::cout << "genFilterInfo: numEventsTried = " << genFilterInfo->numEventsTried() << ", numEventsPassed = " << genFilterInfo->numEventsPassed() << std::endl;
    if ( genFilterInfo->numEventsTried() > 0 ) {
      double weight = genFilterInfo->filterEfficiency();
      //std::cout << "weight(genFilterInfo) = " << weight << std::endl;
      histogramGenFilterEfficiency_->Fill(weight, evtWeight);
      evtWeight *= weight;
      evtWeightMap[Form("%s_%s", srcGenFilterInfo_.label().data(), srcGenFilterInfo_.instance().data())] = weight;
    }
  }

  if ( evtWeight < 1.e-3 || evtWeight > 1.e+3 || TMath::IsNaN(evtWeight) ) return;

  histogramEventCounter_->Fill(1., evtWeight);

  double muonRadCorrWeight     = 1.;
  double muonRadCorrWeightUp   = 1.;
  double muonRadCorrWeightDown = 1.;
  if ( srcMuonRadCorrWeight_.label() != "" && srcMuonRadCorrWeightUp_.label() != "" && srcMuonRadCorrWeightDown_.label() != "" ) {
    edm::Handle<double> muonRadCorrWeight_handle;
    evt.getByLabel(srcMuonRadCorrWeight_, muonRadCorrWeight_handle);
    muonRadCorrWeight = (*muonRadCorrWeight_handle);
    evtWeightMap["muonRadCorrWeight"] = muonRadCorrWeight;
    edm::Handle<double> muonRadCorrWeightUp_handle;
    evt.getByLabel(srcMuonRadCorrWeightUp_, muonRadCorrWeightUp_handle);
    muonRadCorrWeightUp = (*muonRadCorrWeightUp_handle);
    edm::Handle<double> muonRadCorrWeightDown_handle;
    evt.getByLabel(srcMuonRadCorrWeightDown_, muonRadCorrWeightDown_handle);
    muonRadCorrWeightDown = (*muonRadCorrWeightDown_handle);
  }
  //std::cout << " muonRadCorrWeight = " << muonRadCorrWeight 
  //          << " + " << (muonRadCorrWeightUp - muonRadCorrWeight)
  //          << " - " << (muonRadCorrWeight - muonRadCorrWeightDown) << std::endl;
  
  for ( std::vector<plotEntryTypeEvtWeight*>::iterator evtWeightPlotEntry = evtWeightPlotEntries_.begin();
    evtWeightPlotEntry != evtWeightPlotEntries_.end(); ++evtWeightPlotEntry ) {
    (*evtWeightPlotEntry)->fillHistograms(evt, evtWeightMap);
  }
  
//--- fill all histograms
  edm::Handle<reco::TrackCollection> tracks;
  evt.getByLabel(srcRecTracks_, tracks);
  int numTracksPtGt5  = 0;
  int numTracksPtGt10 = 0;
  int numTracksPtGt20 = 0;
  int numTracksPtGt30 = 0;
  int numTracksPtGt40 = 0;
  for ( reco::TrackCollection::const_iterator track = tracks->begin();
    track != tracks->end(); ++track ) {
    if ( track->pt() >  5. ) ++numTracksPtGt5;
    if ( track->pt() > 10. ) ++numTracksPtGt10;
    if ( track->pt() > 20. ) ++numTracksPtGt20;
    if ( track->pt() > 30. ) ++numTracksPtGt30;
    if ( track->pt() > 40. ) ++numTracksPtGt40;
  }
  histogramNumTracksPtGt5_->Fill(numTracksPtGt5, muonRadCorrWeight*evtWeight);
  histogramNumTracksPtGt10_->Fill(numTracksPtGt10, muonRadCorrWeight*evtWeight);
  histogramNumTracksPtGt20_->Fill(numTracksPtGt20, muonRadCorrWeight*evtWeight);
  histogramNumTracksPtGt30_->Fill(numTracksPtGt30, muonRadCorrWeight*evtWeight);
  histogramNumTracksPtGt40_->Fill(numTracksPtGt40, muonRadCorrWeight*evtWeight);
  
 edm::Handle<reco::PFCandidateCollection> pfCandidates;
  evt.getByLabel(srcRecPFCandidates_, pfCandidates);
  int numChargedPFCandsPtGt5  = 0;
  int numChargedPFCandsPtGt10 = 0;
  int numChargedPFCandsPtGt20 = 0;
  int numChargedPFCandsPtGt30 = 0;
  int numChargedPFCandsPtGt40 = 0;
  int numNeutralPFCandsPtGt5  = 0;
  int numNeutralPFCandsPtGt10 = 0;
  int numNeutralPFCandsPtGt20 = 0;
  int numNeutralPFCandsPtGt30 = 0;
  int numNeutralPFCandsPtGt40 = 0;
  for ( reco::PFCandidateCollection::const_iterator pfCandidate = pfCandidates->begin();
    pfCandidate != pfCandidates->end(); ++pfCandidate ) {
    if ( TMath::Abs(pfCandidate->charge()) > 0.5 ) {
      if ( pfCandidate->pt() >  5. ) ++numChargedPFCandsPtGt5;
      if ( pfCandidate->pt() > 10. ) ++numChargedPFCandsPtGt10;
      if ( pfCandidate->pt() > 20. ) ++numChargedPFCandsPtGt20;
      if ( pfCandidate->pt() > 30. ) ++numChargedPFCandsPtGt30;
      if ( pfCandidate->pt() > 40. ) ++numChargedPFCandsPtGt40;
    } else {
      if ( pfCandidate->pt() >  5. ) ++numNeutralPFCandsPtGt5;
      if ( pfCandidate->pt() > 10. ) ++numNeutralPFCandsPtGt10;
      if ( pfCandidate->pt() > 20. ) ++numNeutralPFCandsPtGt20;
      if ( pfCandidate->pt() > 30. ) ++numNeutralPFCandsPtGt30;
      if ( pfCandidate->pt() > 40. ) ++numNeutralPFCandsPtGt40;
    }
  }
  histogramNumChargedPFCandsPtGt5_->Fill(numChargedPFCandsPtGt5, muonRadCorrWeight*evtWeight);
  histogramNumChargedPFCandsPtGt10_->Fill(numChargedPFCandsPtGt10, muonRadCorrWeight*evtWeight);
  histogramNumChargedPFCandsPtGt20_->Fill(numChargedPFCandsPtGt20, muonRadCorrWeight*evtWeight);
  histogramNumChargedPFCandsPtGt30_->Fill(numChargedPFCandsPtGt30, muonRadCorrWeight*evtWeight);
  histogramNumChargedPFCandsPtGt40_->Fill(numChargedPFCandsPtGt40, muonRadCorrWeight*evtWeight);
  histogramNumNeutralPFCandsPtGt5_->Fill(numNeutralPFCandsPtGt5, muonRadCorrWeight*evtWeight);
  histogramNumNeutralPFCandsPtGt10_->Fill(numNeutralPFCandsPtGt10, muonRadCorrWeight*evtWeight);
  histogramNumNeutralPFCandsPtGt20_->Fill(numNeutralPFCandsPtGt20, muonRadCorrWeight*evtWeight);
  histogramNumNeutralPFCandsPtGt30_->Fill(numNeutralPFCandsPtGt30, muonRadCorrWeight*evtWeight);
  histogramNumNeutralPFCandsPtGt40_->Fill(numNeutralPFCandsPtGt40, muonRadCorrWeight*evtWeight);

  edm::Handle<pat::JetCollection> jets;
  evt.getByLabel(srcRecJets_, jets);
  int numJetsRawPtGt20                = 0;
  int numJetsRawPtGt20AbsEtaLt2_5     = 0;
  int numJetsRawPtGt20AbsEta2_5to4_5  = 0;
  int numJetsCorrPtGt20               = 0;
  int numJetsCorrPtGt20AbsEtaLt2_5    = 0;
  int numJetsCorrPtGt20AbsEta2_5to4_5 = 0;
  int numJetsRawPtGt30                = 0;
  int numJetsRawPtGt30AbsEtaLt2_5     = 0;
  int numJetsRawPtGt30AbsEta2_5to4_5  = 0;
  int numJetsCorrPtGt30               = 0;
  int numJetsCorrPtGt30AbsEtaLt2_5    = 0;
  int numJetsCorrPtGt30AbsEta2_5to4_5 = 0;
  for ( pat::JetCollection::const_iterator jet = jets->begin();
    jet != jets->end(); ++jet ) {

    reco::Candidate::LorentzVector rawJetP4 = jet->correctedP4("Uncorrected");
    double rawJetPt = rawJetP4.pt();
    double rawJetAbsEta = TMath::Abs(rawJetP4.eta());
    if ( rawJetAbsEta < 4.5 ) { // CV: do not consider any jet reconstructed outside eta range used in H -> tautau analysis
      histogramRawJetPt_->Fill(rawJetPt, muonRadCorrWeight*evtWeight);
      if      ( rawJetAbsEta < 2.5 ) histogramRawJetPtAbsEtaLt2_5_->Fill(rawJetPt, muonRadCorrWeight*evtWeight);
      else if ( rawJetAbsEta < 4.5 ) histogramRawJetPtAbsEta2_5to4_5_->Fill(rawJetPt, muonRadCorrWeight*evtWeight);
      if ( rawJetPt > 20. ) {
    histogramRawJetEtaPtGt20_->Fill(rawJetP4.eta(), muonRadCorrWeight*evtWeight);
    ++numJetsRawPtGt20;
    if      ( rawJetAbsEta < 2.5 ) ++numJetsRawPtGt20AbsEtaLt2_5;
    else if ( rawJetAbsEta < 4.5 ) ++numJetsRawPtGt20AbsEta2_5to4_5;
      }
      if ( rawJetPt > 30. ) {
    histogramRawJetEtaPtGt20_->Fill(rawJetP4.eta(), muonRadCorrWeight*evtWeight);
    ++numJetsRawPtGt30;
    if      ( rawJetAbsEta < 2.5 ) ++numJetsRawPtGt30AbsEtaLt2_5;
    else if ( rawJetAbsEta < 4.5 ) ++numJetsRawPtGt30AbsEta2_5to4_5;
      }
    }

    reco::Candidate::LorentzVector corrJetP4 = jet->p4();
    double corrJetPt = corrJetP4.pt();
    double corrJetAbsEta = TMath::Abs(corrJetP4.eta());
    if ( corrJetAbsEta < 4.5 ) { // CV: do not consider any jet reconstructed outside eta range used in H -> tautau analysis
      histogramCorrJetPt_->Fill(corrJetPt, muonRadCorrWeight*evtWeight);
      if      ( corrJetAbsEta < 2.5 ) histogramCorrJetPtAbsEtaLt2_5_->Fill(corrJetPt, muonRadCorrWeight*evtWeight);
      else if ( corrJetAbsEta < 4.5 ) histogramCorrJetPtAbsEta2_5to4_5_->Fill(corrJetPt, muonRadCorrWeight*evtWeight);
      if ( corrJetPt > 20. ) {
    histogramCorrJetEtaPtGt20_->Fill(corrJetP4.eta(), muonRadCorrWeight*evtWeight);
    ++numJetsCorrPtGt20;
    if      ( corrJetAbsEta < 2.5 ) ++numJetsCorrPtGt20AbsEtaLt2_5;
    else if ( corrJetAbsEta < 4.5 ) ++numJetsCorrPtGt20AbsEta2_5to4_5;
      }
      if ( corrJetPt > 30. ) {
    histogramCorrJetEtaPtGt20_->Fill(corrJetP4.eta(), muonRadCorrWeight*evtWeight);
    ++numJetsCorrPtGt30;
    if      ( corrJetAbsEta < 2.5 ) ++numJetsCorrPtGt30AbsEtaLt2_5;
    else if ( corrJetAbsEta < 4.5 ) ++numJetsCorrPtGt30AbsEta2_5to4_5;
      }
    }    
  }
  histogramNumJetsRawPtGt20_->Fill(numJetsRawPtGt20, muonRadCorrWeight*evtWeight);
  histogramNumJetsRawPtGt20AbsEtaLt2_5_->Fill(numJetsRawPtGt20AbsEtaLt2_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsRawPtGt20AbsEta2_5to4_5_->Fill(numJetsRawPtGt20AbsEta2_5to4_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt20_->Fill(numJetsCorrPtGt20, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt20AbsEtaLt2_5_->Fill(numJetsCorrPtGt20AbsEtaLt2_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt20AbsEta2_5to4_5_->Fill(numJetsCorrPtGt20AbsEta2_5to4_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsRawPtGt30_->Fill(numJetsRawPtGt30, muonRadCorrWeight*evtWeight);
  histogramNumJetsRawPtGt30AbsEtaLt2_5_->Fill(numJetsRawPtGt30AbsEtaLt2_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsRawPtGt30AbsEta2_5to4_5_->Fill(numJetsRawPtGt30AbsEta2_5to4_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt30_->Fill(numJetsCorrPtGt30, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt30AbsEtaLt2_5_->Fill(numJetsCorrPtGt30AbsEtaLt2_5, muonRadCorrWeight*evtWeight);
  histogramNumJetsCorrPtGt30AbsEta2_5to4_5_->Fill(numJetsCorrPtGt30AbsEta2_5to4_5, muonRadCorrWeight*evtWeight);

  edm::Handle<reco::MuonCollection> muons;
  evt.getByLabel(srcRecMuons_, muons);
  int numGlobalMuons     = 0;
  int numStandAloneMuons = 0;
  int numPFMuons         = 0;
  for ( reco::MuonCollection::const_iterator muon = muons->begin();
    muon != muons->end(); ++muon ) {
    if ( muon->isGlobalMuon()     ) ++numGlobalMuons;
    if ( muon->isStandAloneMuon() ) ++numStandAloneMuons;
    if ( muon->isPFMuon()         ) ++numPFMuons;
  }
  histogramNumGlobalMuons_->Fill(numGlobalMuons, muonRadCorrWeight*evtWeight);
  histogramNumStandAloneMuons_->Fill(numStandAloneMuons, muonRadCorrWeight*evtWeight);
  histogramNumPFMuons_->Fill(numPFMuons, muonRadCorrWeight*evtWeight);

  edm::Handle<reco::VertexCollection> theVertex;
  evt.getByLabel(srcTheRecVertex_, theVertex);
  if ( theVertex->size() >= 1 ) {
    const reco::Vertex::Point& theVertexPosition = theVertex->front().position();
    histogramTheRecVertexX_->Fill(theVertexPosition.x(), muonRadCorrWeight*evtWeight);
    histogramTheRecVertexY_->Fill(theVertexPosition.y(), muonRadCorrWeight*evtWeight);
    histogramTheRecVertexZ_->Fill(theVertexPosition.z(), muonRadCorrWeight*evtWeight);
  }
  edm::Handle<reco::VertexCollection> vertices;
  evt.getByLabel(srcRecVertices_, vertices);
  for ( reco::VertexCollection::const_iterator vertex = vertices->begin();
    vertex != vertices->end(); ++vertex ) {
    histogramRecVertexX_->Fill(vertex->position().x(), muonRadCorrWeight*evtWeight);
    histogramRecVertexY_->Fill(vertex->position().y(), muonRadCorrWeight*evtWeight);
    histogramRecVertexZ_->Fill(vertex->position().z(), muonRadCorrWeight*evtWeight);
  }
  histogramNumRecVertices_->Fill(vertices->size(), muonRadCorrWeight*evtWeight);
  edm::Handle<reco::VertexCollection> verticesWithBS;
  evt.getByLabel(srcRecVerticesWithBS_, verticesWithBS);
  for ( reco::VertexCollection::const_iterator vertex = verticesWithBS->begin();
    vertex != verticesWithBS->end(); ++vertex ) {
    histogramRecVertexWithBSx_->Fill(vertex->position().x(), muonRadCorrWeight*evtWeight);
    histogramRecVertexWithBSy_->Fill(vertex->position().y(), muonRadCorrWeight*evtWeight);
    histogramRecVertexWithBSz_->Fill(vertex->position().z(), muonRadCorrWeight*evtWeight);
  }
  histogramNumRecVerticesWithBS_->Fill(verticesWithBS->size(), muonRadCorrWeight*evtWeight);
  
  edm::Handle<reco::BeamSpot> beamSpot;
  evt.getByLabel(srcBeamSpot_, beamSpot);
  if ( beamSpot.isValid() ) { 
    histogramBeamSpotX_->Fill(beamSpot->position().x(), muonRadCorrWeight*evtWeight);
    histogramBeamSpotY_->Fill(beamSpot->position().y(), muonRadCorrWeight*evtWeight);
  }
      
  for ( CandidateView::const_iterator genDiTau = genDiTaus->begin();
    genDiTau != genDiTaus->end(); ++genDiTau ) {
    histogramGenDiTauPt_->Fill(genDiTau->pt(), muonRadCorrWeight*evtWeight);
    histogramGenDiTauEta_->Fill(genDiTau->eta(), muonRadCorrWeight*evtWeight);
    histogramGenDiTauPhi_->Fill(genDiTau->phi(), muonRadCorrWeight*evtWeight);
    histogramGenDiTauMass_->Fill(genDiTau->mass(), muonRadCorrWeight*evtWeight);
  }

  bool passesCutsBeforeRotation = false;
  if ( (replacedMuonPlus  && replacedMuonPlus->pt()  > replacedMuonPtThresholdHigh_ && replacedMuonMinus && replacedMuonMinus->pt() > replacedMuonPtThresholdLow_) ||
       (replacedMuonMinus && replacedMuonMinus->pt() > replacedMuonPtThresholdHigh_ && replacedMuonPlus  && replacedMuonPlus->pt()  > replacedMuonPtThresholdLow_) ) passesCutsBeforeRotation = true;
  bool passesCutsAfterRotation = false;
  for ( CandidateView::const_iterator genDiTau = genDiTaus->begin();
    genDiTau != genDiTaus->end(); ++genDiTau ) {
    const reco::CompositeCandidate* genDiTau_composite = dynamic_cast<const reco::CompositeCandidate*>(&(*genDiTau));
    if ( !(genDiTau_composite && genDiTau_composite->numberOfDaughters() == 2) ) continue;
    const reco::Candidate* genTau1 = genDiTau_composite->daughter(0);
    const reco::Candidate* genTau2 = genDiTau_composite->daughter(1);
    if ( !(genTau1 && genTau2) ) continue;
    if ( (genTau1->pt() > replacedMuonPtThresholdHigh_ && genTau2->pt() > replacedMuonPtThresholdLow_ ) ||
     (genTau1->pt() > replacedMuonPtThresholdLow_  && genTau2->pt() > replacedMuonPtThresholdHigh_) ) {
      passesCutsAfterRotation = true;
      break;
    }
    histogramRotationAngleMatrix_->Fill(passesCutsBeforeRotation, passesCutsAfterRotation, muonRadCorrWeight*evtWeight);   
  }
  if ( genTauPlus && genTauMinus ) {
    histogramGenDeltaPhiLeg1Leg2_->Fill(normalizedPhi(genTauPlus->phi() - genTauMinus->phi()), muonRadCorrWeight*evtWeight);
    reco::Particle::LorentzVector diTauP4_lab = genTauPlus->p4() + genTauMinus->p4();
    ROOT::Math::Boost boost_to_rf(diTauP4_lab.BoostToCM());
    reco::Particle::LorentzVector diTauP4_rf = boost_to_rf(diTauP4_lab);
    reco::Particle::LorentzVector tauPlusP4_rf = boost_to_rf(genTauPlus->p4());
    if ( (diTauP4_rf.P()*tauPlusP4_rf.P()) > 0. ) {
      double cosGjAngle = (diTauP4_rf.px()*tauPlusP4_rf.px() + diTauP4_rf.py()*tauPlusP4_rf.py() + diTauP4_rf.pz()*tauPlusP4_rf.pz())/(diTauP4_rf.P()*tauPlusP4_rf.P());
      double gjAngle = TMath::ACos(cosGjAngle);
      histogramGenDiTauDecayAngle_->Fill(gjAngle, muonRadCorrWeight*evtWeight);
    }
  }
  if ( replacedMuonPlus && genTauPlus && replacedMuonMinus && genTauMinus ) {
    histogramRotationLegPlusDeltaR_->Fill(deltaR(genTauPlus->p4(), replacedMuonPlus->p4()), muonRadCorrWeight*evtWeight);
    histogramRotationLegMinusDeltaR_->Fill(deltaR(genTauMinus->p4(), replacedMuonMinus->p4()), muonRadCorrWeight*evtWeight);
    
    reco::Particle::LorentzVector diTauP4_lab = genTauPlus->p4() + genTauMinus->p4();
    histogramPhiRotLegPlus_->Fill(SVfit_namespace::phiLabFromLabMomenta(replacedMuonPlus->p4(), diTauP4_lab), muonRadCorrWeight*evtWeight);
    histogramPhiRotLegMinus_->Fill(SVfit_namespace::phiLabFromLabMomenta(replacedMuonMinus->p4(), diTauP4_lab), muonRadCorrWeight*evtWeight);
  }
  
  fillVisPtEtaPhiMassDistributions(evt, srcGenLeg1_, srcGenLeg2_, histogramGenVisDiTauPt_, histogramGenVisDiTauEta_, histogramGenVisDiTauPhi_, histogramGenVisDiTauMass_, muonRadCorrWeight*evtWeight);
  fillVisPtEtaPhiMassDistributions(evt, srcRecLeg1_, srcRecLeg2_, histogramRecVisDiTauPt_, histogramRecVisDiTauEta_, histogramRecVisDiTauPhi_, histogramRecVisDiTauMass_, muonRadCorrWeight*evtWeight);

  typedef edm::View<reco::MET> METView;
  edm::Handle<METView> genCaloMETs;
  evt.getByLabel(srcGenCaloMEt_, genCaloMETs);
  const reco::Candidate::LorentzVector& genCaloMEtP4 = genCaloMETs->front().p4();
  histogramGenCaloMEt_->Fill(genCaloMEtP4.pt(), muonRadCorrWeight*evtWeight);
  edm::Handle<METView> genPFMETs;
  evt.getByLabel(srcGenPFMEt_, genPFMETs);
  const reco::Candidate::LorentzVector& genPFMEtP4 = genPFMETs->front().p4();
  histogramGenPFMEt_->Fill(genPFMEtP4.pt(), muonRadCorrWeight*evtWeight);

  fillX1andX2Distributions(evt, srcGenDiTaus_, srcGenLeg1_, srcGenLeg2_, genPFMEtP4,
               histogramGenTau1Pt_, histogramGenTau1Eta_, histogramGenTau1Phi_,
               histogramGenLeg1Pt_, histogramGenLeg1Eta_, histogramGenLeg1Phi_, histogramGenLeg1X_, 
               histogramGenLeg1XforGenLeg2X0_00to0_25_, histogramGenLeg1XforGenLeg2X0_25to0_50_, histogramGenLeg1XforGenLeg2X0_50to0_75_, histogramGenLeg1XforGenLeg2X0_75to1_00_,
               histogramGenLeg1Mt_, 
               histogramGenTau2Pt_, histogramGenTau2Eta_, histogramGenTau2Phi_,
               histogramGenLeg2Pt_, histogramGenLeg2Eta_, histogramGenLeg2Phi_, histogramGenLeg2X_, 
               histogramGenLeg2XforGenLeg1X0_00to0_25_, histogramGenLeg2XforGenLeg1X0_25to0_50_, histogramGenLeg2XforGenLeg1X0_50to0_75_, histogramGenLeg2XforGenLeg1X0_75to1_00_,
               histogramGenLeg2Mt_, 
               histogramGenVisDeltaPhiLeg1Leg2_, muonRadCorrWeight*evtWeight);

  edm::Handle<METView> recPFMETs;
  evt.getByLabel(srcRecPFMEt_, recPFMETs);
  const reco::Candidate::LorentzVector& recPFMEtP4 = recPFMETs->front().p4();
  
  fillX1andX2Distributions(evt, srcGenDiTaus_, srcRecLeg1_, srcRecLeg2_, recPFMEtP4,
               0, 0, 0, 
               0, 0, 0, histogramRecLeg1X_, 
               0, 0, 0, 0, 
               histogramRecLeg1PFMt_, 
               0, 0, 0, 
               0, 0, 0, histogramRecLeg2X_, 
               0, 0, 0, 0, histogramRecLeg2PFMt_, 
               histogramRecVisDeltaPhiLeg1Leg2_, muonRadCorrWeight*evtWeight);
  
  edm::Handle<reco::GenParticleCollection> genParticles;
  evt.getByLabel(srcGenParticles_, genParticles);     
  reco::Candidate::LorentzVector sumGenParticleP4;
  reco::Candidate::LorentzVector sumGenParticleP4_charged;
  for ( reco::GenParticleCollection::const_iterator genParticle = genParticles->begin();
    genParticle != genParticles->end(); ++genParticle ) {
    if ( genParticle->status() == 1 ) {
      int absPdgId = TMath::Abs(genParticle->pdgId());    
      if ( absPdgId == 12 || absPdgId == 14 || absPdgId == 16 ) continue;
      sumGenParticleP4 += genParticle->p4();
      if ( TMath::Abs(genParticle->charge()) > 0.5 ) sumGenParticleP4_charged += genParticle->p4();
    }
  }
  histogramSumGenParticlePt_->Fill(sumGenParticleP4.pt(), muonRadCorrWeight*evtWeight);
  
  edm::Handle<CaloTowerCollection> caloTowers;
  evt.getByLabel(srcCaloTowers_, caloTowers);
  reco::Candidate::LorentzVector sumCaloTowerP4_ecal;
  double sumEtCaloTowersECAL = 0.;
  reco::Candidate::LorentzVector sumCaloTowerP4_hcal;
  double sumEtCaloTowersHCAL = 0.;
  reco::Candidate::LorentzVector sumCaloTowerP4_hf;
  double sumEtCaloTowersHF   = 0.;
  reco::Candidate::LorentzVector sumCaloTowerP4_ho;
  double sumEtCaloTowersHO   = 0.;
  for ( CaloTowerCollection::const_iterator caloTower = caloTowers->begin();
    caloTower != caloTowers->end(); ++caloTower ) {
    if ( caloTower->energy() != 0. ) {
      double emFrac_ecal = caloTower->emEnergy()/caloTower->energy();
      double emFrac_hcal = (caloTower->energyInHB() + caloTower->energyInHE())/caloTower->energy();
      double emFrac_hf   = caloTower->energyInHF()/caloTower->energy();
      double emFrac_ho   = caloTower->energyInHO()/caloTower->energy();
      sumCaloTowerP4_ecal += (emFrac_ecal*caloTower->p4());
      sumEtCaloTowersECAL += (emFrac_ecal*caloTower->et());
      sumCaloTowerP4_hcal += (emFrac_hcal*caloTower->p4());
      sumEtCaloTowersHCAL += (emFrac_hcal*caloTower->et());
      sumCaloTowerP4_hf   += (emFrac_hf*caloTower->p4());
      sumEtCaloTowersHF   += (emFrac_hf*caloTower->et());
      sumCaloTowerP4_ho   += (emFrac_ho*caloTower->p4());      
      sumEtCaloTowersHO   += (emFrac_ho*caloTower->et());  
    }
  }
  histogramRecCaloMEtECAL_->Fill(sumCaloTowerP4_ecal.pt(), muonRadCorrWeight*evtWeight);
  histogramRecCaloSumEtECAL_->Fill(sumEtCaloTowersECAL, muonRadCorrWeight*evtWeight);
  histogramRecCaloMEtHCAL_->Fill(sumCaloTowerP4_hcal.pt(), muonRadCorrWeight*evtWeight);
  histogramRecCaloSumEtHCAL_->Fill(sumEtCaloTowersHCAL, muonRadCorrWeight*evtWeight);
  histogramRecCaloMEtHF_->Fill(sumCaloTowerP4_hf.pt(), muonRadCorrWeight*evtWeight);
  histogramRecCaloSumEtHF_->Fill(sumEtCaloTowersHF, muonRadCorrWeight*evtWeight);
  histogramRecCaloMEtHO_->Fill(sumCaloTowerP4_ho.pt(), muonRadCorrWeight*evtWeight);
  histogramRecCaloSumEtHO_->Fill(sumEtCaloTowersHO, muonRadCorrWeight*evtWeight);
  
  if ( srcMuonsBeforeRad_.label() != "" && srcMuonsAfterRad_.label() != "" ) {
    reco::Candidate::LorentzVector genMuonPlusP4_beforeRad;
    bool genMuonPlus_beforeRad_found = false;
    reco::Candidate::LorentzVector genMuonMinusP4_beforeRad;
    bool genMuonMinus_beforeRad_found = false;
    reco::Candidate::LorentzVector genMuonPlusP4_afterRad;
    bool genMuonPlus_afterRad_found = false;
    reco::Candidate::LorentzVector genMuonMinusP4_afterRad;
    bool genMuonMinus_afterRad_found = false;
    
    findMuons(evt, srcMuonsBeforeRad_, genMuonPlusP4_beforeRad, genMuonPlus_beforeRad_found, genMuonMinusP4_beforeRad, genMuonMinus_beforeRad_found);
    findMuons(evt, srcMuonsAfterRad_, genMuonPlusP4_afterRad, genMuonPlus_afterRad_found, genMuonMinusP4_afterRad, genMuonMinus_afterRad_found);
 
    bool genMuonPlus_found = (genMuonPlus_beforeRad_found && genMuonPlus_afterRad_found);
    bool genMuonMinus_found = (genMuonMinus_beforeRad_found && genMuonMinus_afterRad_found);

    if ( genTauPlus && genMuonPlus_found && genTauMinus && genMuonMinus_found ) {
      for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator muonRadCorrUncertaintyPlotEntry = muonRadCorrUncertaintyPlotEntries_beforeRad_.begin();
        muonRadCorrUncertaintyPlotEntry != muonRadCorrUncertaintyPlotEntries_beforeRad_.end(); ++muonRadCorrUncertaintyPlotEntry ) {
    (*muonRadCorrUncertaintyPlotEntry)->fillHistograms(numJetsCorrPtGt30, genMuonPlusP4_beforeRad, genMuonMinusP4_beforeRad, evtWeight, muonRadCorrWeight, muonRadCorrWeightUp, muonRadCorrWeightDown);
      }
      for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator muonRadCorrUncertaintyPlotEntry = muonRadCorrUncertaintyPlotEntries_afterRad_.begin();
        muonRadCorrUncertaintyPlotEntry != muonRadCorrUncertaintyPlotEntries_afterRad_.end(); ++muonRadCorrUncertaintyPlotEntry ) {
    (*muonRadCorrUncertaintyPlotEntry)->fillHistograms(numJetsCorrPtGt30, genMuonPlusP4_afterRad, genMuonMinusP4_afterRad, evtWeight, muonRadCorrWeight, muonRadCorrWeightUp, muonRadCorrWeightDown);
      }
      for ( std::vector<plotEntryTypeMuonRadCorrUncertainty*>::iterator muonRadCorrUncertaintyPlotEntry = muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_.begin();
        muonRadCorrUncertaintyPlotEntry != muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_.end(); ++muonRadCorrUncertaintyPlotEntry ) {
    (*muonRadCorrUncertaintyPlotEntry)->fillHistograms(numJetsCorrPtGt30, genTauPlus->p4(), genTauMinus->p4(), evtWeight, muonRadCorrWeight, muonRadCorrWeightUp, muonRadCorrWeightDown);
      }
    } else {
      if ( muonRadCorrUncertainty_numWarnings_ < muonRadCorrUncertainty_maxWarnings_ ) {
    edm::LogWarning ("<MCEmbeddingValidationAnalyzer::analyze>")
      << " (" << runLumiEventNumbers_to_string(evt) << ")" << std::endl
      << " Failed to match muons before and after radiation to embedded tau leptons !!" << std::endl;
    std::cout << "genTauPlus: ";
    if ( genTauPlus ) std::cout << "Pt = " << genTauPlus->pt() << ", eta = " << genTauPlus->eta() << ", phi = " << genTauPlus->phi() << std::endl;
    else std::cout << "NA" << std::endl;
    std::cout << "genMuonPlus (before Rad.): ";
    if ( genMuonPlus_found ) std::cout << "Pt = " << genMuonPlusP4_beforeRad.pt() << ", eta = " << genMuonPlusP4_beforeRad.eta() << ", phi = " << genMuonPlusP4_beforeRad.phi() << std::endl;
    else std::cout << "NA" << std::endl;
    std::cout << "genTauMinus: ";
    if ( genTauMinus ) std::cout << "Pt = " << genTauMinus->pt() << ", eta = " << genTauMinus->eta() << ", phi = " << genTauMinus->phi() << std::endl;
    else std::cout << "NA" << std::endl;
    std::cout << "genMuonMinus (before Rad.): ";
    if ( genMuonMinus_found ) std::cout << "Pt = " << genMuonMinusP4_beforeRad.pt() << ", eta = " << genMuonMinusP4_beforeRad.eta() << ", phi = " << genMuonMinusP4_beforeRad.phi() << std::endl;
    else std::cout << "NA" << std::endl;
    ++muonRadCorrUncertainty_numWarnings_;
      }
    }
  }

  edm::Handle<l1extra::L1EtMissParticleCollection> l1METs;
  evt.getByLabel(srcL1ETM_, l1METs);
  const reco::Candidate::LorentzVector& l1MEtP4 = l1METs->front().p4();
  //std::cout << "L1MEt: Pt = " << l1MEtP4.pt() << ", phi = " << l1MEtP4.phi() << " (Et = " << l1METs->front().etMiss() << ", Px = " << l1MEtP4.px() << ", Py = " << l1MEtP4.py() << ")" << std::endl;
  //if ( l1MEtP4.pt() > 75. ) std::cout << "--> CHECK !!" << std::endl;
  typedef edm::View<reco::MET> METView;
  edm::Handle<METView> recCaloMETs;
  evt.getByLabel(srcRecCaloMEt_, recCaloMETs);
  const reco::Candidate::LorentzVector& recCaloMEtP4 = recCaloMETs->front().p4();
  //std::cout << "recCaloMEt: Pt = " << recCaloMEtP4.pt() << ", phi = " << recCaloMEtP4.phi() << " (Px = " << recCaloMEtP4.px() << ", Py = " << recCaloMEtP4.py() << ")" << std::endl;
  for ( CandidateView::const_iterator genDiTau = genDiTaus->begin();
    genDiTau != genDiTaus->end(); ++genDiTau ) {
    const reco::CompositeCandidate* genDiTau_composite = dynamic_cast<const reco::CompositeCandidate*>(&(*genDiTau));
    if ( !(genDiTau_composite && genDiTau_composite->numberOfDaughters() == 2) ) continue;
    const reco::Candidate* genDaughter1 = genDiTau_composite->daughter(0);
    const reco::GenParticle* genTau1 = ( genDaughter1->hasMasterClone() ) ?
      dynamic_cast<const reco::GenParticle*>(&(*genDaughter1->masterClone())) : dynamic_cast<const reco::GenParticle*>(genDaughter1);
    const reco::Candidate* genDaughter2 = genDiTau_composite->daughter(1);
    const reco::GenParticle* genTau2 = ( genDaughter2->hasMasterClone() ) ?
      dynamic_cast<const reco::GenParticle*>(&(*genDaughter2->masterClone())) : dynamic_cast<const reco::GenParticle*>(genDaughter2);
    if ( !(genTau1 && genTau2) ) continue;
    std::string genTauDecayMode1 = getGenTauDecayMode(genTau1);
    std::string genTauDecayMode2 = getGenTauDecayMode(genTau2);
    std::string genTauDecayMode_ref;
    if      ( genTauDecayMode1 == "electron" || genTauDecayMode1 == "muon" ) genTauDecayMode_ref = genTauDecayMode2;
    else if ( genTauDecayMode2 == "electron" || genTauDecayMode2 == "muon" ) genTauDecayMode_ref = genTauDecayMode1;
    for ( std::vector<plotEntryTypeL1ETM*>::iterator l1ETMplotEntry = l1ETMplotEntries_.begin();
      l1ETMplotEntry != l1ETMplotEntries_.end(); ++l1ETMplotEntry ) {
      (*l1ETMplotEntry)->fillHistograms(genTauDecayMode_ref, l1MEtP4, genCaloMEtP4, recCaloMEtP4, genDiTau->p4(), muonRadCorrWeight*evtWeight);
    }
  }

  fillHistograms(electronDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(electronDistributionsExtra_, numJetsCorrPtGt30, evt, es, muonRadCorrWeight*evtWeight);
  fillHistograms(electronEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(gsfElectronEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(electronL1TriggerEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(muonDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(muonEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(muonL1TriggerEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(tauDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(tauDistributionsExtra_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(tauEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(l1ElectronDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(l1MuonDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(l1TauDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(l1CentralJetDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(l1ForwardJetDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(metDistributions_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);
  fillHistograms(metL1TriggerEfficiencies_, numJetsCorrPtGt30, evt, muonRadCorrWeight*evtWeight);

  // CV: Check for presence in the embedded event of high Pt tracks, charged PFCandidates and muons 
  //     reconstructed near the direction of the replaced muons
  //    (indicating that maybe not all muon signals were removed in the embedding)
  std::vector<const reco::Candidate*> replacedMuons;
  if ( replacedMuonPlus  ) replacedMuons.push_back(replacedMuonPlus);
  if ( replacedMuonMinus ) replacedMuons.push_back(replacedMuonMinus);
  for ( std::vector<const reco::Candidate*>::const_iterator replacedMuon = replacedMuons.begin();
    replacedMuon != replacedMuons.end(); ++replacedMuon ) {
    for ( reco::TrackCollection::const_iterator track = tracks->begin();
      track != tracks->end(); ++track ) {
      if ( track->pt() > 10. ) {
    double dR = deltaR(track->eta(), track->phi(), (*replacedMuon)->eta(), (*replacedMuon)->phi());
    if ( dR < 0.1 ) {
      edm::LogWarning("MCEmbeddingValidationAnalyzer") 
        << " (" << runLumiEventNumbers_to_string(evt) << ")" << std::endl
        << " Found track: Pt = " << track->pt() << ", eta = " << track->eta() << ", phi = " << track->phi() << ", charge = " << track->charge() << " in direction of"
        << " a replaced muon: Pt = " << (*replacedMuon)->pt()<< ", eta = " << (*replacedMuon)->eta() << ", phi = " << (*replacedMuon)->phi() << ", charge = " << (*replacedMuon)->charge()
        << " (dR = " << dR << "). This may point to a problem in removing all muon signals in the Embedding." << std::endl;
      histogramWarning_recTrackNearReplacedMuon_->Fill(track->charge()*(*replacedMuon)->charge(), muonRadCorrWeight*evtWeight);
    }
      }
    }
    for ( reco::PFCandidateCollection::const_iterator pfCandidate = pfCandidates->begin();
      pfCandidate != pfCandidates->end(); ++pfCandidate ) {
      if ( pfCandidate->pt() > 10. && TMath::Abs(pfCandidate->charge()) > 0.5 ) {
    double dR = deltaR(pfCandidate->eta(), pfCandidate->phi(), (*replacedMuon)->eta(), (*replacedMuon)->phi());
    if ( dR < 0.1 ) {
      edm::LogWarning("MCEmbeddingValidationAnalyzer") 
        << " (" << runLumiEventNumbers_to_string(evt) << ")" << std::endl
        << " Found charged PFCandidate: Pt = " << pfCandidate->pt() << ", eta = " << pfCandidate->eta() << ", phi = " << pfCandidate->phi() << ", charge = " << pfCandidate->charge() << " in direction of"
        << " a replaced muon: Pt = " << (*replacedMuon)->pt()<< ", eta = " << (*replacedMuon)->eta() << ", phi = " << (*replacedMuon)->phi() << ", charge = " << (*replacedMuon)->charge()
        << " (dR = " << dR << "). This may point to a problem in removing all muon signals in the Embedding." << std::endl;
      histogramWarning_recPFCandNearReplacedMuon_->Fill(pfCandidate->charge()*(*replacedMuon)->charge(), muonRadCorrWeight*evtWeight);
    }
      }
    }
    for ( reco::MuonCollection::const_iterator muon = muons->begin();
      muon != muons->end(); ++muon ) {
      if ( muon->pt() > 10. ) {
    double dR = deltaR(muon->eta(), muon->phi(), (*replacedMuon)->eta(), (*replacedMuon)->phi());
    if ( dR < 0.1 ) {
      edm::LogWarning("MCEmbeddingValidationAnalyzer") 
        << " (" << runLumiEventNumbers_to_string(evt) << ")" << std::endl
        << " Found track: Pt = " << muon->pt() << ", eta = " << muon->eta() << ", phi = " << muon->phi() << ", charge = " << muon->charge() << " in direction of"
        << " a replaced muon: Pt = " << (*replacedMuon)->pt()<< ", eta = " << (*replacedMuon)->eta() << ", phi = " << (*replacedMuon)->phi() << ", charge = " << (*replacedMuon)->charge()
        << " (dR = " << dR << "). This may point to a problem in removing all muon signals in the Embedding." << std::endl;
      histogramWarning_recMuonNearReplacedMuon_->Fill(muon->charge()*(*replacedMuon)->charge(), muonRadCorrWeight*evtWeight);
    }
      }
    }
  }
}

void MCEmbeddingValidationAnalyzer::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 375 of file MCEmbeddingValidationAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), MCEmbeddingValidationAnalyzer::plotEntryTypeEvtWeight::bookHistograms(), MCEmbeddingValidationAnalyzer::plotEntryTypeMuonRadCorrUncertainty::bookHistograms(), MCEmbeddingValidationAnalyzer::plotEntryTypeL1ETM::bookHistograms(), bookHistograms(), dqmDirectory_, electronDistributions_, electronDistributionsExtra_, electronEfficiencies_, electronL1TriggerEfficiencies_, evtWeightPlotEntries_, edm::hlt::Exception, JetMCTagUtils::genTauDecayMode(), gsfElectronEfficiencies_, histogramBeamSpotX_, histogramBeamSpotY_, histogramCorrJetEtaPtGt20_, histogramCorrJetEtaPtGt30_, histogramCorrJetPt_, histogramCorrJetPtAbsEta2_5to4_5_, histogramCorrJetPtAbsEtaLt2_5_, histogramEventCounter_, histogramGenCaloMEt_, histogramGenDeltaPhiLeg1Leg2_, histogramGenDiTauDecayAngle_, histogramGenDiTauEta_, histogramGenDiTauMass_, histogramGenDiTauPhi_, histogramGenDiTauPt_, histogramGenFilterEfficiency_, histogramGenLeg1Eta_, histogramGenLeg1Mt_, histogramGenLeg1Phi_, histogramGenLeg1Pt_, histogramGenLeg1X_, histogramGenLeg1XforGenLeg2X0_00to0_25_, histogramGenLeg1XforGenLeg2X0_25to0_50_, histogramGenLeg1XforGenLeg2X0_50to0_75_, histogramGenLeg1XforGenLeg2X0_75to1_00_, histogramGenLeg2Eta_, histogramGenLeg2Mt_, histogramGenLeg2Phi_, histogramGenLeg2Pt_, histogramGenLeg2X_, histogramGenLeg2XforGenLeg1X0_00to0_25_, histogramGenLeg2XforGenLeg1X0_25to0_50_, histogramGenLeg2XforGenLeg1X0_50to0_75_, histogramGenLeg2XforGenLeg1X0_75to1_00_, histogramGenPFMEt_, histogramGenTau1Eta_, histogramGenTau1Phi_, histogramGenTau1Pt_, histogramGenTau2Eta_, histogramGenTau2Phi_, histogramGenTau2Pt_, histogramGenVisDeltaPhiLeg1Leg2_, histogramGenVisDiTauEta_, histogramGenVisDiTauMass_, histogramGenVisDiTauPhi_, histogramGenVisDiTauPt_, histogramNumChargedPFCandsPtGt10_, histogramNumChargedPFCandsPtGt20_, histogramNumChargedPFCandsPtGt30_, histogramNumChargedPFCandsPtGt40_, histogramNumChargedPFCandsPtGt5_, histogramNumGlobalMuons_, histogramNumJetsCorrPtGt20_, histogramNumJetsCorrPtGt20AbsEta2_5to4_5_, histogramNumJetsCorrPtGt20AbsEtaLt2_5_, histogramNumJetsCorrPtGt30_, histogramNumJetsCorrPtGt30AbsEta2_5to4_5_, histogramNumJetsCorrPtGt30AbsEtaLt2_5_, histogramNumJetsRawPtGt20_, histogramNumJetsRawPtGt20AbsEta2_5to4_5_, histogramNumJetsRawPtGt20AbsEtaLt2_5_, histogramNumJetsRawPtGt30_, histogramNumJetsRawPtGt30AbsEta2_5to4_5_, histogramNumJetsRawPtGt30AbsEtaLt2_5_, histogramNumNeutralPFCandsPtGt10_, histogramNumNeutralPFCandsPtGt20_, histogramNumNeutralPFCandsPtGt30_, histogramNumNeutralPFCandsPtGt40_, histogramNumNeutralPFCandsPtGt5_, histogramNumPFMuons_, histogramNumRecVertices_, histogramNumRecVerticesWithBS_, histogramNumStandAloneMuons_, histogramNumTracksPtGt10_, histogramNumTracksPtGt20_, histogramNumTracksPtGt30_, histogramNumTracksPtGt40_, histogramNumTracksPtGt5_, histogramPhiRotLegMinus_, histogramPhiRotLegPlus_, histogramRawJetEtaPtGt20_, histogramRawJetEtaPtGt30_, histogramRawJetPt_, histogramRawJetPtAbsEta2_5to4_5_, histogramRawJetPtAbsEtaLt2_5_, histogramRecCaloMEtECAL_, histogramRecCaloMEtHCAL_, histogramRecCaloMEtHF_, histogramRecCaloMEtHO_, histogramRecCaloSumEtECAL_, histogramRecCaloSumEtHCAL_, histogramRecCaloSumEtHF_, histogramRecCaloSumEtHO_, histogramRecLeg1PFMt_, histogramRecLeg1X_, histogramRecLeg2PFMt_, histogramRecLeg2X_, histogramRecVertexWithBSx_, histogramRecVertexWithBSy_, histogramRecVertexWithBSz_, histogramRecVertexX_, histogramRecVertexY_, histogramRecVertexZ_, histogramRecVisDeltaPhiLeg1Leg2_, histogramRecVisDiTauEta_, histogramRecVisDiTauMass_, histogramRecVisDiTauPhi_, histogramRecVisDiTauPt_, histogramRotationAngleMatrix_, histogramRotationLegMinusDeltaR_, histogramRotationLegPlusDeltaR_, histogramSumGenParticlePt_, histogramSumGenParticlePt_charged_, histogramTheRecVertexX_, histogramTheRecVertexY_, histogramTheRecVertexZ_, histogramWarning_recMuonNearReplacedMuon_, histogramWarning_recPFCandNearReplacedMuon_, histogramWarning_recTrackNearReplacedMuon_, l1CentralJetDistributions_, l1ElectronDistributions_, l1ETMplotEntries_, l1ForwardJetDistributions_, l1MuonDistributions_, l1TauDistributions_, edm::InputTag::label(), metDistributions_, metL1TriggerEfficiencies_, muonDistributions_, muonEfficiencies_, muonL1TriggerEfficiencies_, muonRadCorrUncertainty_maxWarnings_, muonRadCorrUncertainty_numWarnings_, muonRadCorrUncertaintyPlotEntries_afterRad_, muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_, muonRadCorrUncertaintyPlotEntries_beforeRad_, Pi, DQMStore::setCurrentFolder(), srcMuonRadCorrWeight_, srcOtherWeights_, AlCaHLTBitMon_QueryRunRegistry::string, tauDistributions_, tauDistributionsExtra_, and tauEfficiencies_.

{
  if ( !edm::Service<DQMStore>().isAvailable() ) 
    throw cms::Exception("MuonRadiationAnalyzer::beginJob")
      << "Failed to access dqmStore !!\n";

  DQMStore& dqmStore = (*edm::Service<DQMStore>());
  dqmStore.setCurrentFolder(dqmDirectory_.data());

//--- book all histograms
  histogramEventCounter_                       = dqmStore.book1D("EventCounter",                       "EventCounter",                              1,     -0.5,          1.5);
 
  histogramGenFilterEfficiency_                = dqmStore.book1D("genFilterEfficiency",                "genFilterEfficiency",                     102,     -0.01,         1.01);

  histogramRotationAngleMatrix_                = dqmStore.book2D("rfRotationAngleMatrix",              "rfRotationAngleMatrix",                     2,     -0.5,          1.5, 2, -0.5, 1.5);
  histogramRotationLegPlusDeltaR_              = dqmStore.book1D("rfRotationLegPlusDeltaR",            "rfRotationLegPlusDeltaR",                 101,     -0.05,        10.05);
  histogramRotationLegMinusDeltaR_             = dqmStore.book1D("rfRotationLegMinusDeltaR",           "rfRotationLegMinusDeltaR",                101,     -0.05,        10.05);
  histogramPhiRotLegPlus_                      = dqmStore.book1D("rfPhiRotLegPlus",                    "rfPhiRotLegPlus",                          72, -TMath::Pi(), +TMath::Pi());
  histogramPhiRotLegMinus_                     = dqmStore.book1D("rfPhiRotLegMinus",                   "rfPhiRotLegMinus",                         72, -TMath::Pi(), +TMath::Pi());

  histogramNumTracksPtGt5_                     = dqmStore.book1D("numTracksPtGt5",                     "numTracksPtGt5",                           50,     -0.5,         49.5);
  histogramNumTracksPtGt10_                    = dqmStore.book1D("numTracksPtGt10",                    "numTracksPtGt10",                          50,     -0.5,         49.5);
  histogramNumTracksPtGt20_                    = dqmStore.book1D("numTracksPtGt20",                    "numTracksPtGt20",                          50,     -0.5,         49.5);
  histogramNumTracksPtGt30_                    = dqmStore.book1D("numTracksPtGt30",                    "numTracksPtGt30",                          50,     -0.5,         49.5);
  histogramNumTracksPtGt40_                    = dqmStore.book1D("numTracksPtGt40",                    "numTracksPtGt40",                          50,     -0.5,         49.5);
      
  histogramNumGlobalMuons_                     = dqmStore.book1D("numGlobalMuons",                     "numGlobalMuons",                           20,     -0.5,         19.5);
  histogramNumStandAloneMuons_                 = dqmStore.book1D("numStandAloneMuons",                 "numStandAloneMuons",                       20,     -0.5,         19.5);
  histogramNumPFMuons_                         = dqmStore.book1D("numPFMuons",                         "numPFMuons",                               20,     -0.5,         19.5);

  histogramNumChargedPFCandsPtGt5_             = dqmStore.book1D("numChargedPFCandsPtGt5",             "numChargedPFCandsPtGt5",                   50,     -0.5,         49.5);
  histogramNumChargedPFCandsPtGt10_            = dqmStore.book1D("numChargedPFCandsPtGt10",            "numChargedPFCandsPtGt10",                  50,     -0.5,         49.5);
  histogramNumChargedPFCandsPtGt20_            = dqmStore.book1D("numChargedPFCandsPtGt20",            "numChargedPFCandsPtGt20",                  50,     -0.5,         49.5);
  histogramNumChargedPFCandsPtGt30_            = dqmStore.book1D("numChargedPFCandsPtGt30",            "numChargedPFCandsPtGt30",                  50,     -0.5,         49.5);
  histogramNumChargedPFCandsPtGt40_            = dqmStore.book1D("numChargedPFCandsPtGt40",            "numChargedPFCandsPtGt40",                  50,     -0.5,         49.5);

  histogramNumNeutralPFCandsPtGt5_             = dqmStore.book1D("numNeutralPFCandsPtGt5",             "numNeutralPFCandsPtGt5",                   50,     -0.5,         49.5);
  histogramNumNeutralPFCandsPtGt10_            = dqmStore.book1D("numNeutralPFCandsPtGt10",            "numNeutralPFCandsPtGt10",                  50,     -0.5,         49.5);
  histogramNumNeutralPFCandsPtGt20_            = dqmStore.book1D("numNeutralPFCandsPtGt20",            "numNeutralPFCandsPtGt20",                  50,     -0.5,         49.5);
  histogramNumNeutralPFCandsPtGt30_            = dqmStore.book1D("numNeutralPFCandsPtGt30",            "numNeutralPFCandsPtGt30",                  50,     -0.5,         49.5);
  histogramNumNeutralPFCandsPtGt40_            = dqmStore.book1D("numNeutralPFCandsPtGt40",            "numNeutralPFCandsPtGt40",                  50,     -0.5,         49.5);
    
  histogramRawJetPt_                           = dqmStore.book1D("rawJetPt",                           "rawJetPt",                                250,      0.,         250.);
  histogramRawJetPtAbsEtaLt2_5_                = dqmStore.book1D("rawJetPtAbsEtaLt2_5",                "rawJetPtAbsEtaLt2_5",                     250,      0.,         250.);  
  histogramRawJetPtAbsEta2_5to4_5_             = dqmStore.book1D("rawJetPtAbsEta2_5to4_5",             "rawJetPtAbsEta2_5to4_5",                  250,      0.,         250.);
  histogramRawJetEtaPtGt20_                    = dqmStore.book1D("rawJetEtaPtGt20",                    "rawJetEtaPtGt20",                         198,     -9.9,         +9.9);
  histogramRawJetEtaPtGt30_                    = dqmStore.book1D("rawJetEtaPtGt30",                    "rawJetEtaPtGt30",                         198,     -9.9,         +9.9);
  histogramNumJetsRawPtGt20_                   = dqmStore.book1D("numJetsRawPtGt20",                   "numJetsRawPtGt20",                         50,     -0.5,         49.5);
  histogramNumJetsRawPtGt20AbsEtaLt2_5_        = dqmStore.book1D("numJetsRawPtGt20AbsEtaLt2_5",        "numJetsRawPtGt20AbsEtaLt2_5",              50,     -0.5,         49.5);
  histogramNumJetsRawPtGt20AbsEta2_5to4_5_     = dqmStore.book1D("numJetsRawPtGt20AbsEta2_5to4_5",     "numJetsRawPtGt20AbsEta2_5to4_5",           50,     -0.5,         49.5);
  histogramNumJetsRawPtGt30_                   = dqmStore.book1D("numJetsRawPtGt30",                   "numJetsRawPtGt30",                         50,     -0.5,         49.5);
  histogramNumJetsRawPtGt30AbsEtaLt2_5_        = dqmStore.book1D("numJetsRawPtGt30AbsEtaLt2_5",        "numJetsRawPtGt30AbsEtaLt2_5",              50,     -0.5,         49.5);
  histogramNumJetsRawPtGt30AbsEta2_5to4_5_     = dqmStore.book1D("numJetsRawPtGt30AbsEta2_5to4_5",     "numJetsRawPtGt30AbsEta2_5to4_5",           50,     -0.5,         49.5);
  histogramCorrJetPt_                          = dqmStore.book1D("corrJetPt",                          "corrJetPt",                               250,      0.,         250.);
  histogramCorrJetPtAbsEtaLt2_5_               = dqmStore.book1D("corrJetPtAbsEtaLt2_5",               "corrJetPtAbsEtaLt2_5",                    250,      0.,         250.);  
  histogramCorrJetPtAbsEta2_5to4_5_            = dqmStore.book1D("corrJetPtAbsEta2_5to4_5",            "corrJetPtAbsEta2_5to4_5",                 250,      0.,         250.);
  histogramCorrJetEtaPtGt20_                   = dqmStore.book1D("corrJetEtaPtGt20",                   "corrJetEtaPtGt20",                        198,     -9.9,         +9.9);
  histogramCorrJetEtaPtGt30_                   = dqmStore.book1D("corrJetEtaPtGt30",                   "corrJetEtaPtGt30",                        198,     -9.9,         +9.9);
  histogramNumJetsCorrPtGt20_                  = dqmStore.book1D("numJetsCorrPtGt20",                  "numJetsCorrPtGt20",                        20,     -0.5,         19.5);
  histogramNumJetsCorrPtGt20AbsEtaLt2_5_       = dqmStore.book1D("numJetsCorrPtGt20AbsEtaLt2_5",       "numJetsCorrPtGt20AbsEtaLt2_5",             20,     -0.5,         19.5);
  histogramNumJetsCorrPtGt20AbsEta2_5to4_5_    = dqmStore.book1D("numJetsCorrPtGt20AbsEta2_5to4_5",    "numJetsCorrPtGt20AbsEta2_5to4_5",          20,     -0.5,         19.5);  
  histogramNumJetsCorrPtGt30_                  = dqmStore.book1D("numJetsCorrPtGt30",                  "numJetsCorrPtGt30",                        20,     -0.5,         19.5);
  histogramNumJetsCorrPtGt30AbsEtaLt2_5_       = dqmStore.book1D("numJetsCorrPtGt30AbsEtaLt2_5",       "numJetsCorrPtGt30AbsEtaLt2_5",             20,     -0.5,         19.5);
  histogramNumJetsCorrPtGt30AbsEta2_5to4_5_    = dqmStore.book1D("numJetsCorrPtGt30AbsEta2_5to4_5",    "numJetsCorrPtGt30AbsEta2_5to4_5",          20,     -0.5,         19.5);
    
  histogramTheRecVertexX_                      = dqmStore.book1D("theRecVertexX",                      "theRecVertexX",                          2000,     -1.,          +1.);
  histogramTheRecVertexY_                      = dqmStore.book1D("theRecVertexY",                      "theRecVertexY",                          2000,     -1.,          +1.);
  histogramTheRecVertexZ_                      = dqmStore.book1D("theRecVertexZ",                      "theRecVertexZ",                           500,    -25.,         +25.);
  histogramRecVertexX_                         = dqmStore.book1D("recVertexX",                         "recVertexX",                             2000,     -1.,          +1.);
  histogramRecVertexY_                         = dqmStore.book1D("recVertexY",                         "recVertexY",                             2000,     -1.,          +1.);
  histogramRecVertexZ_                         = dqmStore.book1D("recVertexZ",                         "recVertexZ",                              500,    -25.,         +25.);
  histogramNumRecVertices_                     = dqmStore.book1D("numRecVertices",                     "numRecVertices",                           50,     -0.5,        +49.5);
  histogramRecVertexWithBSx_                   = dqmStore.book1D("recVertexWithBSx",                   "recVertexWithBSx",                       2000,     -1.,          +1.);
  histogramRecVertexWithBSy_                   = dqmStore.book1D("recVertexWithBSy",                   "recVertexWithBSy",                       2000,     -1.,          +1.);
  histogramRecVertexWithBSz_                   = dqmStore.book1D("recVertexWithBSz",                   "recVertexWithBSz",                        500,    -25.,         +25.);
  histogramNumRecVerticesWithBS_               = dqmStore.book1D("numRecVerticesWithBS",               "numRecVerticesWithBS",                     50,     -0.5,        +49.5);

  histogramBeamSpotX_                          = dqmStore.book1D("beamSpotX",                          "beamSpotX",                              2000,     -1.,          +1.);
  histogramBeamSpotY_                          = dqmStore.book1D("beamSpotY",                          "beamSpotY",                              2000,     -1.,          +1.);
  
  histogramGenDiTauPt_                         = dqmStore.book1D("genDiTauPt",                         "genDiTauPt",                              250,      0.,         250.);
  histogramGenDiTauEta_                        = dqmStore.book1D("genDiTauEta",                        "genDiTauEta",                             198,     -9.9,         +9.9);
  histogramGenDiTauPhi_                        = dqmStore.book1D("genDiTauPhi",                        "genDiTauPhi",                              72, -TMath::Pi(), +TMath::Pi());
  histogramGenDiTauMass_                       = dqmStore.book1D("genDiTauMass",                       "genDiTauMass",                            500,      0.,         500.);
  histogramGenDeltaPhiLeg1Leg2_                = dqmStore.book1D("genDeltaPhiLeg1Leg2",                "genDeltaPhiLeg1Leg2",                     180,      0.,      +TMath::Pi());
  histogramGenDiTauDecayAngle_                 = dqmStore.book1D("genDiTauDecayAngle",                 "genDiTauDecayAngle",                      180,      0.,      +TMath::Pi());

  histogramGenVisDiTauPt_                      = dqmStore.book1D("genVisDiTauPt",                      "genVisDiTauPt",                           250,      0.,         250.);
  histogramGenVisDiTauEta_                     = dqmStore.book1D("genVisDiTauEta",                     "genVisDiTauEta",                          198,     -9.9,         +9.9);
  histogramGenVisDiTauPhi_                     = dqmStore.book1D("genVisDiTauPhi",                     "genVisDiTauPhi",                           72, -TMath::Pi(), +TMath::Pi());
  histogramGenVisDiTauMass_                    = dqmStore.book1D("genVisDiTauMass",                    "genVisDiTauMass",                         500,      0.,         500.);
  histogramGenVisDeltaPhiLeg1Leg2_             = dqmStore.book1D("genVisDeltaPhiLeg1Leg2",             "genVisDeltaPhiLeg1Leg2",                  180,      0.,      +TMath::Pi());

  histogramRecVisDiTauPt_                      = dqmStore.book1D("recVisDiTauPt",                      "recVisDiTauPt",                           250,      0.,         250.);
  histogramRecVisDiTauEta_                     = dqmStore.book1D("recVisDiTauEta",                     "recVisDiTauEta",                          198,     -9.9,         +9.9);
  histogramRecVisDiTauPhi_                     = dqmStore.book1D("recVisDiTauPhi",                     "recVisDiTauPhi",                           72, -TMath::Pi(), +TMath::Pi());
  histogramRecVisDiTauMass_                    = dqmStore.book1D("recVisDiTauMass",                    "recVisDiTauMass",                         500,      0.,         500.);
  histogramRecVisDeltaPhiLeg1Leg2_             = dqmStore.book1D("recVisDeltaPhiLeg1Leg2",             "recVisDeltaPhiLeg1Leg2",                  180,      0.,      +TMath::Pi());

  histogramGenTau1Pt_                          = dqmStore.book1D("genTau1Pt",                          "genTau1Pt",                               250,      0.,         250.);
  histogramGenTau1Eta_                         = dqmStore.book1D("genTau1Eta",                         "genTau1Eta",                              198,     -9.9,         +9.9);
  histogramGenTau1Phi_                         = dqmStore.book1D("genTau1Phi",                         "genTau1Phi",                               72, -TMath::Pi(), +TMath::Pi());
  histogramGenLeg1Pt_                          = dqmStore.book1D("genLeg1Pt",                          "genLeg1Pt",                               250,      0.,         250.);
  histogramGenLeg1Eta_                         = dqmStore.book1D("genLeg1Eta",                         "genLeg1Eta",                              198,     -9.9,         +9.9);
  histogramGenLeg1Phi_                         = dqmStore.book1D("genLeg1Phi",                         "genLeg1Phi",                               72, -TMath::Pi(), +TMath::Pi());
  histogramGenLeg1X_                           = dqmStore.book1D("genLeg1X",                           "X_{1}^{gen}",                             102,     -0.01,         1.01);
  histogramGenLeg1XforGenLeg2X0_00to0_25_      = dqmStore.book1D("genLeg1XforGenLeg2X0_00to0_25",      "X_{1}^{gen} (0.00 < X_{2}^{gen} < 0.25)", 102,     -0.01,         1.01);
  histogramGenLeg1XforGenLeg2X0_25to0_50_      = dqmStore.book1D("genLeg1XforGenLeg2X0_25to0_50",      "X_{1}^{gen} (0.25 < X_{2}^{gen} < 0.50)", 102,     -0.01,         1.01);
  histogramGenLeg1XforGenLeg2X0_50to0_75_      = dqmStore.book1D("genLeg1XforGenLeg2X0_50to0_75",      "X_{1}^{gen} (0.50 < X_{2}^{gen} < 0.75)", 102,     -0.01,         1.01);
  histogramGenLeg1XforGenLeg2X0_75to1_00_      = dqmStore.book1D("genLeg1XforGenLeg2X0_75to1_00",      "X_{1}^{gen} (0.75 < X_{2}^{gen} < 1.00)", 102,     -0.01,         1.01);
  histogramGenLeg1Mt_                          = dqmStore.book1D("genLeg1Mt",                          "genLeg1Mt",                               250,      0.,         250.);
  histogramRecLeg1X_                           = dqmStore.book1D("recLeg1X",                           "recLeg1X",                                102,     -0.01,         1.01);
  histogramRecLeg1PFMt_                        = dqmStore.book1D("recLeg1PFMt",                        "recLeg1PFMt",                             250,      0.,         250.);
  histogramGenTau2Pt_                          = dqmStore.book1D("genTau2Pt",                          "genTau2Pt",                               250,      0.,         250.);
  histogramGenTau2Eta_                         = dqmStore.book1D("genTau2Eta",                         "genTau2Eta",                              198,     -9.9,         +9.9);
  histogramGenTau2Phi_                         = dqmStore.book1D("genTau2Phi",                         "genTau2Phi",                               72, -TMath::Pi(), +TMath::Pi());
  histogramGenLeg2Pt_                          = dqmStore.book1D("genLeg2Pt",                          "genLeg2Pt",                               250,      0.,         250.);
  histogramGenLeg2Eta_                         = dqmStore.book1D("genLeg2Eta",                         "genLeg2Eta",                              198,     -9.9,         +9.9);
  histogramGenLeg2Phi_                         = dqmStore.book1D("genLeg2Phi",                         "genLeg2Phi",                               72, -TMath::Pi(), +TMath::Pi());
  histogramGenLeg2X_                           = dqmStore.book1D("genLeg2X",                           "X_{2}^{gen}",                             102,     -0.01,         1.01);
  histogramGenLeg2XforGenLeg1X0_00to0_25_      = dqmStore.book1D("genLeg2XforGenLeg1X0_00to0_25",      "X_{2}^{gen} (0.00 < X_{1}^{gen} < 0.25)", 102,     -0.01,         1.01);
  histogramGenLeg2XforGenLeg1X0_25to0_50_      = dqmStore.book1D("genLeg2XforGenLeg1X0_25to0_50",      "X_{2}^{gen} (0.25 < X_{1}^{gen} < 0.50)", 102,     -0.01,         1.01);
  histogramGenLeg2XforGenLeg1X0_50to0_75_      = dqmStore.book1D("genLeg2XforGenLeg1X0_50to0_75",      "X_{2}^{gen} (0.50 < X_{1}^{gen} < 0.75)", 102,     -0.01,         1.01);
  histogramGenLeg2XforGenLeg1X0_75to1_00_      = dqmStore.book1D("genLeg2XforGenLeg1X0_75to1_00",      "X_{2}^{gen} (0.75 < X_{1}^{gen} < 1.00)", 102,     -0.01,         1.01);
  histogramGenLeg2Mt_                          = dqmStore.book1D("genLeg2Mt",                          "genLeg2Mt",                               250,      0.,         250.);
  histogramRecLeg2X_                           = dqmStore.book1D("recLeg2X",                           "recLeg2X",                                102,     -0.01,         1.01);
  histogramRecLeg2PFMt_                        = dqmStore.book1D("recLeg2PFMt",                        "recLeg2PFMt",                             250,      0.,         250.);

  histogramSumGenParticlePt_                   = dqmStore.book1D("sumGenParticlePt",                   "sumGenParticlePt",                        250,      0.,         250.);
  histogramSumGenParticlePt_charged_           = dqmStore.book1D("sumGenParticlePt_charged",           "sumGenParticlePt_charged",                250,      0.,         250.);
  histogramGenCaloMEt_                         = dqmStore.book1D("genCaloMEt",                         "genCaloMEt",                              250,      0.,         250.);
  histogramGenPFMEt_                           = dqmStore.book1D("genPFMEt",                           "genPFMEt",                                250,      0.,         250.);

  histogramRecCaloMEtECAL_                     = dqmStore.book1D("recCaloMEtECAL",                     "recCaloMEtECAL",                          250,      0.,         250.);
  histogramRecCaloSumEtECAL_                   = dqmStore.book1D("recCaloSumEtECAL",                   "recCaloSumEtECAL",                       2500,      0.,        2500.);
  histogramRecCaloMEtHCAL_                     = dqmStore.book1D("recCaloMEtHCAL",                     "recCaloMEtHCAL",                          250,      0.,         250.);
  histogramRecCaloSumEtHCAL_                   = dqmStore.book1D("recCaloSumEtHCAL",                   "recCaloSumEtHCAL",                       2500,      0.,        2500.);
  histogramRecCaloMEtHF_                       = dqmStore.book1D("recCaloMEtHF",                       "recCaloMEtHF",                            250,      0.,         250.);
  histogramRecCaloSumEtHF_                     = dqmStore.book1D("recCaloSumEtHF",                     "recCaloSumEtHF",                         2500,      0.,        2500.);
  histogramRecCaloMEtHO_                       = dqmStore.book1D("recCaloMEtHO",                       "recCaloMEtHO",                            250,      0.,         250.);  
  histogramRecCaloSumEtHO_                     = dqmStore.book1D("recCaloSumEtHO",                     "recCaloSumEtHO",                         2500,      0.,        2500.);

  // CV: Record presence in the embedded event of high Pt tracks, PFCandidates and muons 
  //     reconstructed near the direction of the replaced muons
  //    (indicating that maybe not all muon signals were removed in the embedding).
  //     Fill product of reconstructed track/PFCandidate/muon charge * charge of replaced muon into histogram.
  //     In case all matches happen just by chance, expect equal number of entries in positive and negative bins
  histogramWarning_recTrackNearReplacedMuon_   = dqmStore.book1D("Warning_recTrackNearReplacedMuon",   "Warning_recTrackNearReplacedMuon",          3,    -1.5,        +1.5);
  histogramWarning_recPFCandNearReplacedMuon_  = dqmStore.book1D("Warning_recPFCandNearReplacedMuon",  "Warning_recPFCandNearReplacedMuon",         3,    -1.5,        +1.5);
  histogramWarning_recMuonNearReplacedMuon_    = dqmStore.book1D("Warning_recMuonNearReplacedMuon",    "Warning_recMuonNearReplacedMuon",           3,    -1.5,        +1.5);

  for ( vInputTag::const_iterator srcWeight = srcOtherWeights_.begin();
    srcWeight != srcOtherWeights_.end(); ++srcWeight ) {
    plotEntryTypeEvtWeight* evtWeightPlotEntry = new plotEntryTypeEvtWeight(*srcWeight, dqmDirectory_);
    evtWeightPlotEntry->bookHistograms(dqmStore);
    evtWeightPlotEntries_.push_back(evtWeightPlotEntry);
  }
  if ( srcMuonRadCorrWeight_.label() != "" ) {
    plotEntryTypeEvtWeight* evtWeightPlotEntry = new plotEntryTypeEvtWeight(srcMuonRadCorrWeight_, dqmDirectory_);
    evtWeightPlotEntry->bookHistograms(dqmStore);
    evtWeightPlotEntries_.push_back(evtWeightPlotEntry);
  }

  typedef std::pair<int, int> pint;
  std::vector<pint> jetBins;
  jetBins.push_back(pint(-1, -1));
  jetBins.push_back(pint(0, 0));
  jetBins.push_back(pint(1, 1));
  jetBins.push_back(pint(2, 2));
  jetBins.push_back(pint(3, 1000));
  for ( std::vector<pint>::const_iterator jetBin = jetBins.begin();
    jetBin != jetBins.end(); ++jetBin ) {
    TString dqmDirectory_beforeRad = dqmDirectory_;
    if ( !dqmDirectory_beforeRad.EndsWith("/") ) dqmDirectory_beforeRad.Append("/");
    dqmDirectory_beforeRad.Append("beforeRad");
    plotEntryTypeMuonRadCorrUncertainty* muonRadCorrUncertaintyPlotEntry_beforeRad = new plotEntryTypeMuonRadCorrUncertainty(jetBin->first, jetBin->second, dqmDirectory_beforeRad.Data());
    muonRadCorrUncertaintyPlotEntry_beforeRad->bookHistograms(dqmStore);
    muonRadCorrUncertaintyPlotEntries_beforeRad_.push_back(muonRadCorrUncertaintyPlotEntry_beforeRad);
    TString dqmDirectory_afterRad = dqmDirectory_;
    if ( !dqmDirectory_afterRad.EndsWith("/") ) dqmDirectory_afterRad.Append("/");
    dqmDirectory_afterRad.Append("afterRad");
    plotEntryTypeMuonRadCorrUncertainty* muonRadCorrUncertaintyPlotEntry_afterRad = new plotEntryTypeMuonRadCorrUncertainty(jetBin->first, jetBin->second, dqmDirectory_afterRad.Data());
    muonRadCorrUncertaintyPlotEntry_afterRad->bookHistograms(dqmStore);
    muonRadCorrUncertaintyPlotEntries_afterRad_.push_back(muonRadCorrUncertaintyPlotEntry_afterRad);
    TString dqmDirectory_afterRadAndCorr = dqmDirectory_;
    if ( !dqmDirectory_afterRadAndCorr.EndsWith("/") ) dqmDirectory_afterRadAndCorr.Append("/");
    dqmDirectory_afterRadAndCorr.Append("afterRadAndCorr");
    plotEntryTypeMuonRadCorrUncertainty* muonRadCorrUncertaintyPlotEntry_afterRadAndCorr = new plotEntryTypeMuonRadCorrUncertainty(jetBin->first, jetBin->second, dqmDirectory_afterRadAndCorr.Data());
    muonRadCorrUncertaintyPlotEntry_afterRadAndCorr->bookHistograms(dqmStore);
    muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_.push_back(muonRadCorrUncertaintyPlotEntry_afterRadAndCorr);
  }
  muonRadCorrUncertainty_numWarnings_ = 0;
  muonRadCorrUncertainty_maxWarnings_ = 3;

  std::vector<std::string> genTauDecayModes;
  genTauDecayModes.push_back(std::string("")); // all tau decay modes
  genTauDecayModes.push_back(std::string("oneProng0Pi0"));
  genTauDecayModes.push_back(std::string("oneProng1Pi0"));
  genTauDecayModes.push_back(std::string("oneProng2Pi0"));
  genTauDecayModes.push_back(std::string("threeProng0Pi0"));
  genTauDecayModes.push_back(std::string("threeProng1Pi0"));
  for ( std::vector<std::string>::const_iterator genTauDecayMode = genTauDecayModes.begin();
    genTauDecayMode != genTauDecayModes.end(); ++genTauDecayMode ) {
    plotEntryTypeL1ETM* l1ETMplotEntry = new plotEntryTypeL1ETM(*genTauDecayMode, dqmDirectory_);
    l1ETMplotEntry->bookHistograms(dqmStore);
    l1ETMplotEntries_.push_back(l1ETMplotEntry);
  }

  bookHistograms(electronDistributions_, dqmStore);
  bookHistograms(electronDistributionsExtra_, dqmStore);
  bookHistograms(electronEfficiencies_, dqmStore);
  bookHistograms(gsfElectronEfficiencies_, dqmStore);
  bookHistograms(electronL1TriggerEfficiencies_, dqmStore);
  bookHistograms(muonDistributions_, dqmStore);
  bookHistograms(muonEfficiencies_, dqmStore);
  bookHistograms(muonL1TriggerEfficiencies_, dqmStore);
  bookHistograms(tauDistributions_, dqmStore);
  bookHistograms(tauDistributionsExtra_, dqmStore);
  bookHistograms(tauEfficiencies_, dqmStore);
  bookHistograms(l1ElectronDistributions_, dqmStore);
  bookHistograms(l1MuonDistributions_, dqmStore);
  bookHistograms(l1TauDistributions_, dqmStore);
  bookHistograms(l1CentralJetDistributions_, dqmStore);
  bookHistograms(l1ForwardJetDistributions_, dqmStore);
  bookHistograms(metDistributions_, dqmStore);
  bookHistograms(metL1TriggerEfficiencies_, dqmStore);
}

template<typename T >

void MCEmbeddingValidationAnalyzer::bookHistograms ( std::vector< T * >  collection, DQMStore &  dqmStore  )

inlineprivate

Definition at line 1493 of file MCEmbeddingValidationAnalyzer.h.

References dbtoconf::object.

Referenced by beginJob().

  {
    for ( typename std::vector<T*>::iterator object = collection.begin();
      object != collection.end(); ++object ) {
      (*object)->bookHistograms(dqmStore);
    }
  } 

template<typename T >

void MCEmbeddingValidationAnalyzer::cleanCollection ( std::vector< T * >  collection )

inlineprivate

Definition at line 1484 of file MCEmbeddingValidationAnalyzer.h.

References dbtoconf::object.

Referenced by ~MCEmbeddingValidationAnalyzer().

  {
    for ( typename std::vector<T*>::iterator object = collection.begin();
      object != collection.end(); ++object ) {
      delete (*object);
    }
  }

std::string MCEmbeddingValidationAnalyzer::dqmDirectory_full ( const std::string &  dqmSubDirectory )

inlineprivate

Definition at line 209 of file MCEmbeddingValidationAnalyzer.h.

References dqmDirectory_.

Referenced by setupElectronDistributionExtra(), setupL1ExtraObjectDistribution(), setupLeptonDistribution(), setupLeptonEfficiency(), setupLeptonL1TriggerEfficiency(), setupMEtDistribution(), setupMEtL1TriggerEfficiency(), and setupTauDistributionExtra().

  {
    TString dqmDirectory_full = dqmDirectory_.data();
    if ( !dqmDirectory_full.EndsWith("/") ) dqmDirectory_full.Append("/");
    dqmDirectory_full.Append(dqmSubDirectory);
    return dqmDirectory_full.Data();
  }

template<typename T >

void MCEmbeddingValidationAnalyzer::fillHistograms ( std::vector< T * >  collection, int  numJets, const edm::Event &  evt, double  evtWeight  )

inlineprivate

Definition at line 1502 of file MCEmbeddingValidationAnalyzer.h.

References dbtoconf::object.

Referenced by analyze().

  {
    for ( typename std::vector<T*>::iterator object = collection.begin();
      object != collection.end(); ++object ) {
      (*object)->fillHistograms(numJets, evt, evtWeight);
    }
  } 

template<typename T >

void MCEmbeddingValidationAnalyzer::fillHistograms ( std::vector< T * >  collection, int  numJets, const edm::Event &  evt, const edm::EventSetup &  es, double  evtWeight  )

inlineprivate

Definition at line 1510 of file MCEmbeddingValidationAnalyzer.h.

References dbtoconf::object.

  {
    for ( typename std::vector<T*>::iterator object = collection.begin();
      object != collection.end(); ++object ) {
      (*object)->fillHistograms(numJets, evt, es, evtWeight);
    }
  } 

void MCEmbeddingValidationAnalyzer::setupElectronDistributionExtra ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< electronDistributionExtra * > &  electronDistributionsExtra  )

private

Definition at line 233 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), ebRHToken_, eeRHToken_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), fastjetJetProducer_validation_cfg::srcGen, srcTheRecVertex_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgLeptonDistributions = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgLeptonDistribution = cfgLeptonDistributions.begin();
      cfgLeptonDistribution != cfgLeptonDistributions.end(); ++cfgLeptonDistribution ) {
      edm::InputTag srcGen = cfgLeptonDistribution->getParameter<edm::InputTag>("srcGen");
      std::string cutGen = cfgLeptonDistribution->exists("cutGen") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutGen") : "";
      edm::InputTag srcRec = cfgLeptonDistribution->getParameter<edm::InputTag>("srcRec");
      std::string cutRec = cfgLeptonDistribution->exists("cutRec") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutRec") : "";
      double dRmatch = cfgLeptonDistribution->exists("dRmatch") ? 
    cfgLeptonDistribution->getParameter<double>("dRmatch") : 0.3;
      std::string dqmDirectory = dqmDirectory_full(cfgLeptonDistribution->getParameter<std::string>("dqmDirectory"));
      electronDistributionExtra* electronDistribution = new electronDistributionExtra(minJets, maxJets, srcGen, cutGen, srcRec, cutRec, dRmatch, dqmDirectory, srcTheRecVertex_, ebRHToken_, eeRHToken_);
      electronDistributionsExtra.push_back(electronDistribution);
    }
  }
}

template<typename T >

void MCEmbeddingValidationAnalyzer::setupL1ExtraObjectDistribution ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< l1ExtraObjectDistributionT< T > * > &  l1ExtraObjectDistributions  )

private

Definition at line 323 of file MCEmbeddingValidationAnalyzer.cc.

References GOODCOLL_filter_cfg::cut, dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), alcazmumu_cfi::src, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgL1ExtraObjectDistributions = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgL1ExtraObjectDistribution = cfgL1ExtraObjectDistributions.begin();
      cfgL1ExtraObjectDistribution != cfgL1ExtraObjectDistributions.end(); ++cfgL1ExtraObjectDistribution ) {
      edm::InputTag src = cfgL1ExtraObjectDistribution->getParameter<edm::InputTag>("src");
      std::string cut = cfgL1ExtraObjectDistribution->exists("cut") ? 
    cfgL1ExtraObjectDistribution->getParameter<std::string>("cut") : "";
      std::string dqmDirectory = dqmDirectory_full(cfgL1ExtraObjectDistribution->getParameter<std::string>("dqmDirectory"));
      l1ExtraObjectDistributionT<T>* l1ExtraObjectDistribution = new l1ExtraObjectDistributionT<T>(minJets, maxJets, src, cut, dqmDirectory);
      l1ExtraObjectDistributions.push_back(l1ExtraObjectDistribution);
    }
  }
}

template<typename T >

void MCEmbeddingValidationAnalyzer::setupLeptonDistribution ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< leptonDistributionT< T > * > &  leptonDistributions  )

private

Definition at line 211 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), fastjetJetProducer_validation_cfg::srcGen, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgLeptonDistributions = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgLeptonDistribution = cfgLeptonDistributions.begin();
      cfgLeptonDistribution != cfgLeptonDistributions.end(); ++cfgLeptonDistribution ) {
      edm::InputTag srcGen = cfgLeptonDistribution->getParameter<edm::InputTag>("srcGen");
      std::string cutGen = cfgLeptonDistribution->exists("cutGen") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutGen") : "";
      edm::InputTag srcRec = cfgLeptonDistribution->getParameter<edm::InputTag>("srcRec");
      std::string cutRec = cfgLeptonDistribution->exists("cutRec") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutRec") : "";
      double dRmatch = cfgLeptonDistribution->exists("dRmatch") ? 
    cfgLeptonDistribution->getParameter<double>("dRmatch") : 0.3;
      std::string dqmDirectory = dqmDirectory_full(cfgLeptonDistribution->getParameter<std::string>("dqmDirectory"));
      leptonDistributionT<T>* leptonDistribution = new leptonDistributionT<T>(minJets, maxJets, srcGen, cutGen, srcRec, cutRec, dRmatch, dqmDirectory);
      leptonDistributions.push_back(leptonDistribution);
    }
  }
}

template<typename T >

void MCEmbeddingValidationAnalyzer::setupLeptonEfficiency ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< leptonEfficiencyT< T > * > &  leptonEfficiencies  )

private

Definition at line 278 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), fastjetJetProducer_validation_cfg::srcGen, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgLeptonEfficiencies = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgLeptonEfficiency = cfgLeptonEfficiencies.begin();
      cfgLeptonEfficiency != cfgLeptonEfficiencies.end(); ++cfgLeptonEfficiency ) {
      edm::InputTag srcGen = cfgLeptonEfficiency->getParameter<edm::InputTag>("srcGen");
      std::string cutGen = cfgLeptonEfficiency->exists("cutGen") ? 
    cfgLeptonEfficiency->getParameter<std::string>("cutGen") : "";
      edm::InputTag srcRec = cfgLeptonEfficiency->getParameter<edm::InputTag>("srcRec");
      std::string cutRec = cfgLeptonEfficiency->exists("cutRec") ? 
    cfgLeptonEfficiency->getParameter<std::string>("cutRec") : "";
      double dRmatch = cfgLeptonEfficiency->exists("dRmatch") ? 
    cfgLeptonEfficiency->getParameter<double>("dRmatch") : 0.3;
      std::string dqmDirectory = dqmDirectory_full(cfgLeptonEfficiency->getParameter<std::string>("dqmDirectory"));
      leptonEfficiencyT<T>* leptonEfficiency = new leptonEfficiencyT<T>(minJets, maxJets, srcGen, cutGen, srcRec, cutRec, dRmatch, dqmDirectory);
      leptonEfficiencies.push_back(leptonEfficiency);
    }
  }
}

template<typename T1 , typename T2 >

void MCEmbeddingValidationAnalyzer::setupLeptonL1TriggerEfficiency ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< leptonL1TriggerEfficiencyT1T2< T1, T2 > * > &  leptonL1TriggerEfficiencies  )

private

Definition at line 301 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgLeptonL1TriggerEfficiencies = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgLeptonL1TriggerEfficiency = cfgLeptonL1TriggerEfficiencies.begin();
      cfgLeptonL1TriggerEfficiency != cfgLeptonL1TriggerEfficiencies.end(); ++cfgLeptonL1TriggerEfficiency ) {
      edm::InputTag srcRef = cfgLeptonL1TriggerEfficiency->getParameter<edm::InputTag>("srcRef");
      std::string cutRef = cfgLeptonL1TriggerEfficiency->exists("cutRef") ? 
    cfgLeptonL1TriggerEfficiency->getParameter<std::string>("cutRef") : "";
      edm::InputTag srcL1 = cfgLeptonL1TriggerEfficiency->getParameter<edm::InputTag>("srcL1");
      std::string cutL1 = cfgLeptonL1TriggerEfficiency->getParameter<std::string>("cutL1");
      double dRmatch = cfgLeptonL1TriggerEfficiency->exists("dRmatch") ? 
    cfgLeptonL1TriggerEfficiency->getParameter<double>("dRmatch") : 0.3;
      std::string dqmDirectory = dqmDirectory_full(cfgLeptonL1TriggerEfficiency->getParameter<std::string>("dqmDirectory"));
      leptonL1TriggerEfficiencyT1T2<T1,T2>* leptonL1TriggerEfficiency = new leptonL1TriggerEfficiencyT1T2<T1,T2>(minJets, maxJets, srcRef, cutRef, srcL1, cutL1, dRmatch, dqmDirectory);
      leptonL1TriggerEfficiencies.push_back(leptonL1TriggerEfficiency);
    }
  }
}

void MCEmbeddingValidationAnalyzer::setupMEtDistribution ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< metDistributionType * > &  metDistributions  )

private

Definition at line 340 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), fastjetJetProducer_validation_cfg::srcGen, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgMEtDistributions = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgMEtDistribution = cfgMEtDistributions.begin();
      cfgMEtDistribution != cfgMEtDistributions.end(); ++cfgMEtDistribution ) {
      edm::InputTag srcGen = cfgMEtDistribution->getParameter<edm::InputTag>("srcGen");
      edm::InputTag srcRec = cfgMEtDistribution->getParameter<edm::InputTag>("srcRec");
      edm::InputTag srcGenZs = cfgMEtDistribution->getParameter<edm::InputTag>("srcGenZs");
      std::string dqmDirectory = dqmDirectory_full(cfgMEtDistribution->getParameter<std::string>("dqmDirectory"));
      metDistributionType* metDistribution = new metDistributionType(minJets, maxJets, srcGen, srcRec, srcGenZs, dqmDirectory);
      metDistributions.push_back(metDistribution);
    }
  }
}

void MCEmbeddingValidationAnalyzer::setupMEtL1TriggerEfficiency ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< metL1TriggerEfficiencyType * > &  metL1TriggerEfficiencies  )

private

Definition at line 357 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgMEtL1TriggerEfficiencies = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgMEtL1TriggerEfficiency = cfgMEtL1TriggerEfficiencies.begin();
      cfgMEtL1TriggerEfficiency != cfgMEtL1TriggerEfficiencies.end(); ++cfgMEtL1TriggerEfficiency ) {
      edm::InputTag srcRef = cfgMEtL1TriggerEfficiency->getParameter<edm::InputTag>("srcRef");
      edm::InputTag srcL1 = cfgMEtL1TriggerEfficiency->getParameter<edm::InputTag>("srcL1");
      double cutL1Et = cfgMEtL1TriggerEfficiency->getParameter<double>("cutL1Et");
      double cutL1Pt = cfgMEtL1TriggerEfficiency->getParameter<double>("cutL1Pt");
      std::string dqmDirectory = dqmDirectory_full(cfgMEtL1TriggerEfficiency->getParameter<std::string>("dqmDirectory"));
      metL1TriggerEfficiencyType* metL1TriggerEfficiency = new metL1TriggerEfficiencyType(minJets, maxJets, srcRef, srcL1, cutL1Et, cutL1Pt, dqmDirectory);
      metL1TriggerEfficiencies.push_back(metL1TriggerEfficiency);
    }
  }
}

void MCEmbeddingValidationAnalyzer::setupTauDistributionExtra ( int  minJets, int  maxJets, const edm::ParameterSet &  cfg, const std::string &  keyword, std::vector< tauDistributionExtra * > &  tauDistributionsExtra  )

private

Definition at line 255 of file MCEmbeddingValidationAnalyzer.cc.

References dqmDirectory_full(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), fastjetJetProducer_validation_cfg::srcGen, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MCEmbeddingValidationAnalyzer().

{
  if ( cfg.exists(keyword) ) {
    edm::VParameterSet cfgLeptonDistributions = cfg.getParameter<edm::VParameterSet>(keyword);
    for ( edm::VParameterSet::const_iterator cfgLeptonDistribution = cfgLeptonDistributions.begin();
      cfgLeptonDistribution != cfgLeptonDistributions.end(); ++cfgLeptonDistribution ) {
      edm::InputTag srcGen = cfgLeptonDistribution->getParameter<edm::InputTag>("srcGen");
      std::string cutGen = cfgLeptonDistribution->exists("cutGen") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutGen") : "";
      edm::InputTag srcRec = cfgLeptonDistribution->getParameter<edm::InputTag>("srcRec");
      std::string cutRec = cfgLeptonDistribution->exists("cutRec") ? 
    cfgLeptonDistribution->getParameter<std::string>("cutRec") : "";
      double dRmatch = cfgLeptonDistribution->exists("dRmatch") ? 
    cfgLeptonDistribution->getParameter<double>("dRmatch") : 0.3;
      std::string dqmDirectory = dqmDirectory_full(cfgLeptonDistribution->getParameter<std::string>("dqmDirectory"));
      tauDistributionExtra* tauDistribution = new tauDistributionExtra(minJets, maxJets, srcGen, cutGen, srcRec, cutRec, dRmatch, dqmDirectory);
      tauDistributionsExtra.push_back(tauDistribution);
    }
  }
}

Member Data Documentation

std::string MCEmbeddingValidationAnalyzer::dqmDirectory_

private

Definition at line 250 of file MCEmbeddingValidationAnalyzer.h.

Referenced by beginJob(), and dqmDirectory_full().

double MCEmbeddingValidationAnalyzer::dRminSeparation_

private

Definition at line 228 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<EcalRecHitCollection> MCEmbeddingValidationAnalyzer::ebRHToken_

private

Definition at line 248 of file MCEmbeddingValidationAnalyzer.h.

Referenced by MCEmbeddingValidationAnalyzer(), and setupElectronDistributionExtra().

edm::EDGetTokenT<EcalRecHitCollection> MCEmbeddingValidationAnalyzer::eeRHToken_

private

Definition at line 248 of file MCEmbeddingValidationAnalyzer.h.

Referenced by MCEmbeddingValidationAnalyzer(), and setupElectronDistributionExtra().

std::vector<leptonDistributionT<pat::Electron>*> MCEmbeddingValidationAnalyzer::electronDistributions_

private

Definition at line 1065 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<electronDistributionExtra*> MCEmbeddingValidationAnalyzer::electronDistributionsExtra_

private

Definition at line 1066 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonEfficiencyT<pat::Electron>*> MCEmbeddingValidationAnalyzer::electronEfficiencies_

private

Definition at line 1166 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonL1TriggerEfficiencyT1T2<pat::Electron, l1extra::L1EmParticle>*> MCEmbeddingValidationAnalyzer::electronL1TriggerEfficiencies_

private

Definition at line 1258 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<plotEntryTypeEvtWeight*> MCEmbeddingValidationAnalyzer::evtWeightPlotEntries_

private

Definition at line 423 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonEfficiencyT<reco::GsfElectron>*> MCEmbeddingValidationAnalyzer::gsfElectronEfficiencies_

private

Definition at line 1165 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramBeamSpotX_

private

Definition at line 321 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramBeamSpotY_

private

Definition at line 322 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramCorrJetEtaPtGt20_

private

Definition at line 300 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramCorrJetEtaPtGt30_

private

Definition at line 301 of file MCEmbeddingValidationAnalyzer.h.

Referenced by beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramCorrJetPt_

private

Definition at line 297 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramCorrJetPtAbsEta2_5to4_5_

private

Definition at line 299 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramCorrJetPtAbsEtaLt2_5_

private

Definition at line 298 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramEventCounter_

private

Definition at line 252 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenCaloMEt_

private

Definition at line 374 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDeltaPhiLeg1Leg2_

private

Definition at line 328 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDiTauDecayAngle_

private

Definition at line 329 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDiTauEta_

private

Definition at line 325 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDiTauMass_

private

Definition at line 327 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDiTauPhi_

private

Definition at line 326 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenDiTauPt_

private

Definition at line 324 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenFilterEfficiency_

private

Definition at line 254 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1Eta_

private

Definition at line 347 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1Mt_

private

Definition at line 354 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1Phi_

private

Definition at line 348 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1Pt_

private

Definition at line 346 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1X_

private

Definition at line 349 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1XforGenLeg2X0_00to0_25_

private

Definition at line 350 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1XforGenLeg2X0_25to0_50_

private

Definition at line 351 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1XforGenLeg2X0_50to0_75_

private

Definition at line 352 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg1XforGenLeg2X0_75to1_00_

private

Definition at line 353 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2Eta_

private

Definition at line 361 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2Mt_

private

Definition at line 368 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2Phi_

private

Definition at line 362 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2Pt_

private

Definition at line 360 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2X_

private

Definition at line 363 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2XforGenLeg1X0_00to0_25_

private

Definition at line 364 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2XforGenLeg1X0_25to0_50_

private

Definition at line 365 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2XforGenLeg1X0_50to0_75_

private

Definition at line 366 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenLeg2XforGenLeg1X0_75to1_00_

private

Definition at line 367 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenPFMEt_

private

Definition at line 375 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau1Eta_

private

Definition at line 344 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau1Phi_

private

Definition at line 345 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau1Pt_

private

Definition at line 343 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau2Eta_

private

Definition at line 358 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau2Phi_

private

Definition at line 359 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenTau2Pt_

private

Definition at line 357 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenVisDeltaPhiLeg1Leg2_

private

Definition at line 335 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenVisDiTauEta_

private

Definition at line 332 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenVisDiTauMass_

private

Definition at line 334 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenVisDiTauPhi_

private

Definition at line 333 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramGenVisDiTauPt_

private

Definition at line 331 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumChargedPFCandsPtGt10_

private

Definition at line 275 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumChargedPFCandsPtGt20_

private

Definition at line 276 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumChargedPFCandsPtGt30_

private

Definition at line 277 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumChargedPFCandsPtGt40_

private

Definition at line 278 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumChargedPFCandsPtGt5_

private

Definition at line 274 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumGlobalMuons_

private

Definition at line 270 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt20_

private

Definition at line 302 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt20AbsEta2_5to4_5_

private

Definition at line 304 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt20AbsEtaLt2_5_

private

Definition at line 303 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt30_

private

Definition at line 305 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt30AbsEta2_5to4_5_

private

Definition at line 307 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsCorrPtGt30AbsEtaLt2_5_

private

Definition at line 306 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt20_

private

Definition at line 291 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt20AbsEta2_5to4_5_

private

Definition at line 293 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt20AbsEtaLt2_5_

private

Definition at line 292 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt30_

private

Definition at line 294 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt30AbsEta2_5to4_5_

private

Definition at line 296 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumJetsRawPtGt30AbsEtaLt2_5_

private

Definition at line 295 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumNeutralPFCandsPtGt10_

private

Definition at line 281 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumNeutralPFCandsPtGt20_

private

Definition at line 282 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumNeutralPFCandsPtGt30_

private

Definition at line 283 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumNeutralPFCandsPtGt40_

private

Definition at line 284 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumNeutralPFCandsPtGt5_

private

Definition at line 280 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumPFMuons_

private

Definition at line 272 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumRecVertices_

private

Definition at line 315 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumRecVerticesWithBS_

private

Definition at line 319 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumStandAloneMuons_

private

Definition at line 271 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumTracksPtGt10_

private

Definition at line 265 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumTracksPtGt20_

private

Definition at line 266 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumTracksPtGt30_

private

Definition at line 267 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumTracksPtGt40_

private

Definition at line 268 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramNumTracksPtGt5_

private

Definition at line 264 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramPhiRotLegMinus_

private

Definition at line 262 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramPhiRotLegPlus_

private

Definition at line 261 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRawJetEtaPtGt20_

private

Definition at line 289 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRawJetEtaPtGt30_

private

Definition at line 290 of file MCEmbeddingValidationAnalyzer.h.

Referenced by beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRawJetPt_

private

Definition at line 286 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRawJetPtAbsEta2_5to4_5_

private

Definition at line 288 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRawJetPtAbsEtaLt2_5_

private

Definition at line 287 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloMEtECAL_

private

Definition at line 377 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloMEtHCAL_

private

Definition at line 379 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloMEtHF_

private

Definition at line 381 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloMEtHO_

private

Definition at line 383 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloSumEtECAL_

private

Definition at line 378 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloSumEtHCAL_

private

Definition at line 380 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloSumEtHF_

private

Definition at line 382 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecCaloSumEtHO_

private

Definition at line 384 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecLeg1PFMt_

private

Definition at line 356 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecLeg1X_

private

Definition at line 355 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecLeg2PFMt_

private

Definition at line 370 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecLeg2X_

private

Definition at line 369 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexWithBSx_

private

Definition at line 316 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexWithBSy_

private

Definition at line 317 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexWithBSz_

private

Definition at line 318 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexX_

private

Definition at line 312 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexY_

private

Definition at line 313 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVertexZ_

private

Definition at line 314 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVisDeltaPhiLeg1Leg2_

private

Definition at line 341 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVisDiTauEta_

private

Definition at line 338 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVisDiTauMass_

private

Definition at line 340 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVisDiTauPhi_

private

Definition at line 339 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRecVisDiTauPt_

private

Definition at line 337 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRotationAngleMatrix_

private

Definition at line 256 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRotationLegMinusDeltaR_

private

Definition at line 260 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramRotationLegPlusDeltaR_

private

Definition at line 259 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramSumGenParticlePt_

private

Definition at line 372 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramSumGenParticlePt_charged_

private

Definition at line 373 of file MCEmbeddingValidationAnalyzer.h.

Referenced by beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramTheRecVertexX_

private

Definition at line 309 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramTheRecVertexY_

private

Definition at line 310 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramTheRecVertexZ_

private

Definition at line 311 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramWarning_recMuonNearReplacedMuon_

private

Definition at line 388 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramWarning_recPFCandNearReplacedMuon_

private

Definition at line 387 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* MCEmbeddingValidationAnalyzer::histogramWarning_recTrackNearReplacedMuon_

private

Definition at line 386 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

std::vector<l1ExtraObjectDistributionT<l1extra::L1JetParticle>*> MCEmbeddingValidationAnalyzer::l1CentralJetDistributions_

private

Definition at line 1328 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<l1ExtraObjectDistributionT<l1extra::L1EmParticle>*> MCEmbeddingValidationAnalyzer::l1ElectronDistributions_

private

Definition at line 1325 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<plotEntryTypeL1ETM*> MCEmbeddingValidationAnalyzer::l1ETMplotEntries_

private

Definition at line 635 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCEmbeddingValidationAnalyzer().

std::vector<l1ExtraObjectDistributionT<l1extra::L1JetParticle>*> MCEmbeddingValidationAnalyzer::l1ForwardJetDistributions_

private

Definition at line 1329 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<l1ExtraObjectDistributionT<l1extra::L1MuonParticle>*> MCEmbeddingValidationAnalyzer::l1MuonDistributions_

private

Definition at line 1326 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<l1ExtraObjectDistributionT<l1extra::L1JetParticle>*> MCEmbeddingValidationAnalyzer::l1TauDistributions_

private

Definition at line 1327 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<metDistributionType*> MCEmbeddingValidationAnalyzer::metDistributions_

private

Definition at line 1412 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<metL1TriggerEfficiencyType*> MCEmbeddingValidationAnalyzer::metL1TriggerEfficiencies_

private

Definition at line 1481 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::string MCEmbeddingValidationAnalyzer::moduleLabel_

private

Definition at line 207 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

std::vector<leptonDistributionT<pat::Muon>*> MCEmbeddingValidationAnalyzer::muonDistributions_

private

Definition at line 1067 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonEfficiencyT<pat::Muon>*> MCEmbeddingValidationAnalyzer::muonEfficiencies_

private

Definition at line 1167 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonL1TriggerEfficiencyT1T2<pat::Muon, l1extra::L1MuonParticle>*> MCEmbeddingValidationAnalyzer::muonL1TriggerEfficiencies_

private

Definition at line 1259 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

int MCEmbeddingValidationAnalyzer::muonRadCorrUncertainty_maxWarnings_

private

Definition at line 531 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

int MCEmbeddingValidationAnalyzer::muonRadCorrUncertainty_numWarnings_

private

Definition at line 530 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

std::vector<plotEntryTypeMuonRadCorrUncertainty*> MCEmbeddingValidationAnalyzer::muonRadCorrUncertaintyPlotEntries_afterRad_

private

Definition at line 528 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCEmbeddingValidationAnalyzer().

std::vector<plotEntryTypeMuonRadCorrUncertainty*> MCEmbeddingValidationAnalyzer::muonRadCorrUncertaintyPlotEntries_afterRadAndCorr_

private

Definition at line 529 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCEmbeddingValidationAnalyzer().

std::vector<plotEntryTypeMuonRadCorrUncertainty*> MCEmbeddingValidationAnalyzer::muonRadCorrUncertaintyPlotEntries_beforeRad_

private

Definition at line 527 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCEmbeddingValidationAnalyzer().

double MCEmbeddingValidationAnalyzer::replacedMuonPtThresholdHigh_

private

Definition at line 257 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

double MCEmbeddingValidationAnalyzer::replacedMuonPtThresholdLow_

private

Definition at line 258 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcBeamSpot_

private

Definition at line 226 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcCaloTowers_

private

Definition at line 220 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenCaloMEt_

private

Definition at line 235 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenDiTaus_

private

Definition at line 227 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenFilterInfo_

private

Definition at line 246 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenLeg1_

private

Definition at line 229 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenLeg2_

private

Definition at line 231 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenParticles_

private

Definition at line 233 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcGenPFMEt_

private

Definition at line 236 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcL1ETM_

private

Definition at line 234 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcMuonRadCorrWeight_

private

Definition at line 241 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

edm::InputTag MCEmbeddingValidationAnalyzer::srcMuonRadCorrWeightDown_

private

Definition at line 243 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcMuonRadCorrWeightUp_

private

Definition at line 242 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcMuonsAfterRad_

private

Definition at line 240 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcMuonsBeforeRad_

private

Definition at line 239 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

vInputTag MCEmbeddingValidationAnalyzer::srcOtherWeights_

private

Definition at line 245 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and beginJob().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecCaloMEt_

private

Definition at line 237 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecJets_

private

Definition at line 222 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecLeg1_

private

Definition at line 230 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecLeg2_

private

Definition at line 232 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecMuons_

private

Definition at line 218 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecPFCandidates_

private

Definition at line 221 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecPFMEt_

private

Definition at line 238 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecTracks_

private

Definition at line 219 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecVertices_

private

Definition at line 224 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcRecVerticesWithBS_

private

Definition at line 225 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcReplacedMuons_

private

Definition at line 217 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze().

edm::InputTag MCEmbeddingValidationAnalyzer::srcTheRecVertex_

private

Definition at line 223 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and setupElectronDistributionExtra().

std::vector<leptonDistributionT<pat::Tau>*> MCEmbeddingValidationAnalyzer::tauDistributions_

private

Definition at line 1068 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<tauDistributionExtra*> MCEmbeddingValidationAnalyzer::tauDistributionsExtra_

private

Definition at line 1069 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

std::vector<leptonEfficiencyT<pat::Tau>*> MCEmbeddingValidationAnalyzer::tauEfficiencies_

private

Definition at line 1168 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), beginJob(), MCEmbeddingValidationAnalyzer(), and ~MCEmbeddingValidationAnalyzer().

int MCEmbeddingValidationAnalyzer::verbosity_

private

Definition at line 1518 of file MCEmbeddingValidationAnalyzer.h.

Referenced by analyze(), and MCEmbeddingValidationAnalyzer().