EgammaCoreTools. More...

#include <ZToMuMuGammaAnalyzer.h>

Inheritance diagram for ZToMuMuGammaAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	ZToMuMuGammaAnalyzer (const edm::ParameterSet &)

virtual	~ZToMuMuGammaAnalyzer ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
bool	basicMuonSelection (const reco::Muon &m)

float	mumuGammaInvMass (const reco::Muon &mu1, const reco::Muon &mu2, const reco::PhotonRef &pho)

float	mumuInvMass (const reco::Muon &m1, const reco::Muon &m2)

bool	muonSelection (const reco::Muon &m, const reco::BeamSpot &bs)

bool	photonSelection (const reco::PhotonRef &p)

Private Attributes
edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > >	barrelRecHit_token_

edm::EDGetTokenT< reco::BeamSpot >	beamSpot_token_

std::stringstream	currentFolder_

int	eBin_

double	eMax_

double	eMin_

edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > >	endcapRecHit_token_

int	etaBin_

double	etaMax_

double	etaMin_

int	etBin_

double	etMax_

double	etMin_

float	farMuonEcalIso_

float	farMuonMinPt_

float	farMuonTrackIso_

std::string	fName_

MonitorElement *	h1_mumuGammaInvMass_ [3]

MonitorElement *	h1_mumuInvMass_ [3]
	photon histos More...

MonitorElement *	h2_e1x5VsEt_ [3]

MonitorElement *	h2_e1x5VsEta_ [3]

MonitorElement *	h2_e2x5VsEt_ [3]

MonitorElement *	h2_e2x5VsEta_ [3]

MonitorElement *	h2_ecalSumVsEt_ [3]

MonitorElement *	h2_ecalSumVsEta_ [3]

MonitorElement *	h2_hcalSumVsEt_ [3]

MonitorElement *	h2_hcalSumVsEta_ [3]

MonitorElement *	h2_nTrackIsolHollowVsEt_ [3]

MonitorElement *	h2_nTrackIsolHollowVsEta_ [3]

MonitorElement *	h2_nTrackIsolSolidVsEt_ [3]

MonitorElement *	h2_nTrackIsolSolidVsEta_ [3]

MonitorElement *	h2_r1x5VsEt_ [3]

MonitorElement *	h2_r1x5VsEta_ [3]

MonitorElement *	h2_r2x5VsEt_ [3]

MonitorElement *	h2_r2x5VsEta_ [3]

MonitorElement *	h2_r9VsEt_ [3]

MonitorElement *	h2_r9VsEta_ [3]

MonitorElement *	h2_sigmaIetaIetaVsEta_ [3]

MonitorElement *	h2_trackPtSumHollowVsEt_ [3]

MonitorElement *	h2_trackPtSumHollowVsEta_ [3]

MonitorElement *	h2_trackPtSumSolidVsEt_ [3]

MonitorElement *	h2_trackPtSumSolidVsEta_ [3]

MonitorElement *	h_chHadIso_ [3]

MonitorElement *	h_dRPhoPFcand_ChHad_Cleaned_ [3]

MonitorElement *	h_dRPhoPFcand_ChHad_unCleaned_ [3]

MonitorElement *	h_dRPhoPFcand_NeuHad_Cleaned_ [3]

MonitorElement *	h_dRPhoPFcand_NeuHad_unCleaned_ [3]

MonitorElement *	h_dRPhoPFcand_Pho_Cleaned_ [3]

MonitorElement *	h_dRPhoPFcand_Pho_unCleaned_ [3]

MonitorElement *	h_e1x5_ [3]

MonitorElement *	h_e2x5_ [3]

MonitorElement *	h_ecalSum_ [3]

MonitorElement *	h_etOutsideMustache_ [3]

MonitorElement *	h_h1OverE_ [3]

MonitorElement *	h_h2OverE_ [3]

MonitorElement *	h_hcalSum_ [3]

MonitorElement *	h_hOverE_ [3]

MonitorElement *	h_nCluOutsideMustache_ [3]

MonitorElement *	h_newhOverE_ [3]

MonitorElement *	h_nHadIso_ [3]

MonitorElement *	h_nPho_ [3]

MonitorElement *	h_nRecoVtx_

MonitorElement *	h_nTrackIsolHollow_ [3]

MonitorElement *	h_nTrackIsolSolid_ [3]

MonitorElement *	h_pfMva_ [3]

MonitorElement *	h_phoE_ [3]

MonitorElement *	h_phoEt_ [3]

MonitorElement *	h_phoEta_ [3]

MonitorElement *	h_phoIso_ [3]

MonitorElement *	h_phoPhi_ [3]

MonitorElement *	h_phoSigmaEoverE_ [3]

MonitorElement *	h_phoSigmaIetaIeta_ [3]

MonitorElement *	h_r1x5_ [3]

MonitorElement *	h_r2x5_ [3]

MonitorElement *	h_r9_ [3]

MonitorElement *	h_scEta_ [3]

MonitorElement *	h_scPhi_ [3]

MonitorElement *	h_SumPtOverPhoPt_ChHad_Cleaned_ [3]

MonitorElement *	h_SumPtOverPhoPt_ChHad_unCleaned_ [3]

MonitorElement *	h_SumPtOverPhoPt_NeuHad_Cleaned_ [3]

MonitorElement *	h_SumPtOverPhoPt_NeuHad_unCleaned_ [3]

MonitorElement *	h_SumPtOverPhoPt_Pho_Cleaned_ [3]

MonitorElement *	h_SumPtOverPhoPt_Pho_unCleaned_ [3]

MonitorElement *	h_trackPtSumHollow_ [3]

MonitorElement *	h_trackPtSumSolid_ [3]

int	hOverEBin_

double	hOverEMax_

double	hOverEMin_

bool	makeProfiles_

float	maxMumuGammaInvMass_

float	maxMumuInvMass_

float	minMumuGammaInvMass_

float	minMumuInvMass_

int	minPixStripHits_

edm::EDGetTokenT< std::vector< reco::Muon > >	muon_token_

int	muonMatches_

float	muonMaxChi2_

float	muonMaxDxy_

float	muonMinPt_

float	muonTightEta_

float	muonTrackIso_

float	nearMuonDr_

float	nearMuonHcalIso_

int	nEvt_

int	numberBin_

double	numberMax_

double	numberMin_

edm::EDGetTokenT< reco::VertexCollection >	offline_pvToken_

MonitorElement *	p_e1x5VsEt_ [3]

MonitorElement *	p_e1x5VsEta_ [3]

MonitorElement *	p_e2x5VsEt_ [3]

MonitorElement *	p_e2x5VsEta_ [3]

MonitorElement *	p_ecalSumVsEt_ [3]

MonitorElement *	p_ecalSumVsEta_ [3]

MonitorElement *	p_hcalSumVsEt_ [3]

MonitorElement *	p_hcalSumVsEta_ [3]

MonitorElement *	p_hOverEVsEt_ [3]

MonitorElement *	p_hOverEVsEta_ [3]

MonitorElement *	p_newhOverEVsEt_ [3]

MonitorElement *	p_newhOverEVsEta_ [3]

MonitorElement *	p_nTrackIsolHollowVsEt_ [3]

MonitorElement *	p_nTrackIsolHollowVsEta_ [3]

MonitorElement *	p_nTrackIsolSolidVsEt_ [3]

MonitorElement *	p_nTrackIsolSolidVsEta_ [3]

MonitorElement *	p_phoSigmaEoverEVsNVtx_ [3]

MonitorElement *	p_r1x5VsEt_ [3]

MonitorElement *	p_r1x5VsEta_ [3]

MonitorElement *	p_r2x5VsEt_ [3]

MonitorElement *	p_r2x5VsEta_ [3]

MonitorElement *	p_r9VsEt_ [3]

MonitorElement *	p_r9VsEta_ [3]

MonitorElement *	p_sigmaIetaIetaVsEta_ [3]

MonitorElement *	p_trackPtSumHollowVsEt_ [3]

MonitorElement *	p_trackPtSumHollowVsEta_ [3]

MonitorElement *	p_trackPtSumSolidVsEt_ [3]

MonitorElement *	p_trackPtSumSolidVsEta_ [3]

edm::EDGetTokenT< reco::PFCandidateCollection >	pfCandidates_

int	phiBin_

double	phiMax_

double	phiMin_

edm::EDGetTokenT< std::vector< reco::Photon > >	photon_token_

edm::EDGetTokenT< edm::ValueMap< bool > >	PhotonIDLoose_token_

edm::EDGetTokenT< edm::ValueMap< bool > >	PhotonIDTight_token_

edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > >	photonIsoValmap_token_

float	photonMaxEta_

float	photonMinEt_

float	photonTrackIso_

unsigned int	prescaleFactor_

int	r9Bin_

double	r9Max_

double	r9Min_

int	reducedEtaBin_

int	reducedEtBin_

int	reducedR9Bin_

int	reducedSumBin_

int	sigmaIetaBin_

double	sigmaIetaMax_

double	sigmaIetaMin_

int	sumBin_

double	sumMax_

double	sumMin_

edm::EDGetTokenT< trigger::TriggerEvent >	triggerEvent_token_

bool	use2DHistos_

int	validMuonHits_

int	validPixHits_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

EgammaCoreTools.

$Id: ZToMuMuGammaAnalyzer authors: Nancy Marinelli, U. of Notre Dame, US Nathan Kellams, U. of Notre Dame, US

$Id: ZToMuMuGammaAnalyzer authors: Nancy Marinelli, U. of Notre Dame, US Jamie Antonelli, U. of Notre Dame, US Nathan Kellams, U. of Notre Dame, US

Definition at line 100 of file ZToMuMuGammaAnalyzer.h.

Constructor & Destructor Documentation

ZToMuMuGammaAnalyzer::ZToMuMuGammaAnalyzer ( const edm::ParameterSet & pset )

explicit

Definition at line 18 of file ZToMuMuGammaAnalyzer.cc.

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

 {
   fName_                 = pset.getParameter<std::string>("analyzerName");
   prescaleFactor_        = pset.getUntrackedParameter<int>("prescaleFactor",1);
   use2DHistos_           = pset.getParameter<bool>("use2DHistos");
   makeProfiles_          = pset.getParameter<bool>("makeProfiles");
   
   triggerEvent_token_    = consumes<trigger::TriggerEvent>(pset.getParameter<edm::InputTag>("triggerEvent"));
   offline_pvToken_       = consumes<reco::VertexCollection>(pset.getUntrackedParameter<edm::InputTag>("offlinePV", edm::InputTag("offlinePrimaryVertices")));
   photon_token_          = consumes<vector<reco::Photon> >(pset.getParameter<edm::InputTag>("phoProducer"));
   muon_token_            = consumes<vector<reco::Muon> >(pset.getParameter<edm::InputTag>("muonProducer"));
   pfCandidates_          = consumes<reco::PFCandidateCollection>(pset.getParameter<edm::InputTag>("pfCandidates"));
   photonIsoValmap_token_ = consumes<edm::ValueMap<std::vector<reco::PFCandidateRef> > >(pset.getParameter<edm::InputTag>("particleBasedIso"));
   barrelRecHit_token_    = consumes<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > >(pset.getParameter<edm::InputTag>("barrelRecHitProducer"));
   endcapRecHit_token_    = consumes<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > >(pset.getParameter<edm::InputTag>("endcapRecHitProducer"));
   beamSpot_token_        = consumes<reco::BeamSpot>(pset.getParameter<edm::InputTag>("beamSpot"));
   
   nEvt_=0;
   
   // Muon selection
   muonMinPt_             = pset.getParameter<double>("muonMinPt");
   minPixStripHits_       = pset.getParameter<int>("minPixStripHits");
   muonMaxChi2_           = pset.getParameter<double>("muonMaxChi2");
   muonMaxDxy_            = pset.getParameter<double>("muonMaxDxy");
   muonMatches_           = pset.getParameter<int>("muonMatches");
   validPixHits_          = pset.getParameter<int>("validPixHits");
   validMuonHits_         = pset.getParameter<int>("validMuonHits");
   muonTrackIso_          = pset.getParameter<double>("muonTrackIso");
   muonTightEta_          = pset.getParameter<double>("muonTightEta");
   // Dimuon selection
   minMumuInvMass_        = pset.getParameter<double>("minMumuInvMass");
   maxMumuInvMass_        = pset.getParameter<double>("maxMumuInvMass");
   // Photon selection
   photonMinEt_           = pset.getParameter<double>("photonMinEt");
   photonMaxEta_          = pset.getParameter<double>("photonMaxEta");
   photonTrackIso_        = pset.getParameter<double>("photonTrackIso");
   // mumuGamma selection
   nearMuonDr_            = pset.getParameter<double>("nearMuonDr");
   nearMuonHcalIso_       = pset.getParameter<double>("nearMuonHcalIso");
   farMuonEcalIso_        = pset.getParameter<double>("farMuonEcalIso");
   farMuonTrackIso_       = pset.getParameter<double>("farMuonTrackIso");
   farMuonMinPt_          = pset.getParameter<double>("farMuonMinPt");
   minMumuGammaInvMass_   = pset.getParameter<double>("minMumuGammaInvMass");
   maxMumuGammaInvMass_   = pset.getParameter<double>("maxMumuGammaInvMass");
 
   // Histogram parameters 
   eMin_                  = pset.getParameter<double>("eMin");
   eMax_                  = pset.getParameter<double>("eMax");
   eBin_                  = pset.getParameter<int>("eBin");
 
   etMin_                 = pset.getParameter<double>("etMin");
   etMax_                 = pset.getParameter<double>("etMax");
   etBin_                 = pset.getParameter<int>("etBin");
 
   sumMin_                = pset.getParameter<double>("sumMin");
   sumMax_                = pset.getParameter<double>("sumMax");
   sumBin_                = pset.getParameter<int>("sumBin");
 
   etaMin_                = pset.getParameter<double>("etaMin");
   etaMax_                = pset.getParameter<double>("etaMax");
   etaBin_                = pset.getParameter<int>("etaBin");
 
   phiMin_                = pset.getParameter<double>("phiMin");
   phiMax_                = pset.getParameter<double>("phiMax");
   phiBin_                = pset.getParameter<int>("phiBin");
 
   r9Min_                 = pset.getParameter<double>("r9Min");
   r9Max_                 = pset.getParameter<double>("r9Max");
   r9Bin_                 = pset.getParameter<int>("r9Bin");
 
   hOverEMin_             = pset.getParameter<double>("hOverEMin");
   hOverEMax_             = pset.getParameter<double>("hOverEMax");
   hOverEBin_             = pset.getParameter<int>("hOverEBin");
 
   numberMin_             = pset.getParameter<double>("numberMin");
   numberMax_             = pset.getParameter<double>("numberMax");
   numberBin_             = pset.getParameter<int>("numberBin");
 
   sigmaIetaMin_          = pset.getParameter<double>("sigmaIetaMin");
   sigmaIetaMax_          = pset.getParameter<double>("sigmaIetaMax");
   sigmaIetaBin_          = pset.getParameter<int>("sigmaIetaBin");
 
   reducedEtBin_          = etBin_/4;
   reducedEtaBin_         = etaBin_/4;
   reducedSumBin_         = sumBin_/4;
   reducedR9Bin_          = r9Bin_/4;
 }

ZToMuMuGammaAnalyzer::~ZToMuMuGammaAnalyzer ( )

virtual

Definition at line 106 of file ZToMuMuGammaAnalyzer.cc.

107 {

108 }

Member Function Documentation

void ZToMuMuGammaAnalyzer::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	esup
	)

overridevirtual

uncleaned

Definition at line 421 of file ZToMuMuGammaAnalyzer.cc.

References begin, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::PFCandidate::e, reco::MuonIsolation::emEt, edm::ValueMap< T >::end(), HcalObjRepresent::Fill(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), reco::PFCandidate::gamma, edm::Event::getByToken(), reco::PFCandidate::h, reco::PFCandidate::h0, reco::MuonIsolation::hadEt, mps_fire::i, edm::EventBase::id(), reco::Muon::isolationR03(), edm::HandleBase::isValid(), edm::InputTag::label(), diffTwoXMLs::label, HiRecoMuon_cff::muonCollection, softMuonTagInfos_cfi::muonSelection, correctedPhotons_cfi::photonCollection, edm::Handle< T >::product(), reco::LeafCandidate::pt(), trigger::TriggerEvent::sizeFilters(), createPayload::skip, and reco::MuonIsolation::sumPt.

 {
   using namespace edm;
 
   if (nEvt_% prescaleFactor_ ) return;
   nEvt_++;
   LogInfo("ZToMuMuGammaAnalyzer") << "ZToMuMuGammaAnalyzer Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
 
   // Get the trigger results
   bool validTriggerEvent=true;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   trigger::TriggerEvent triggerEvent;
   e.getByToken(triggerEvent_token_,triggerEventHandle);
   if(!triggerEventHandle.isValid()) {
     edm::LogInfo("PhotonAnalyzer") << "Error! Can't get the product: triggerEvent_token_" << endl;
     validTriggerEvent=false;
   }
   if(validTriggerEvent) triggerEvent = *(triggerEventHandle.product());
 
   // Get the reconstructed photons
   //  bool validPhotons=true;
   Handle<reco::PhotonCollection> photonHandle;
   reco::PhotonCollection photonCollection;
   e.getByToken(photon_token_ , photonHandle);
   if ( !photonHandle.isValid()) {
     edm::LogInfo("ZToMuMuGammaAnalyzer") << "Error! Can't get the product: photon_token_" << endl;
     //validPhotons=false;
   }
   //  if(validPhotons) photonCollection = *(photonHandle.product());
 
   // Get the  PF refined cluster  collection
   Handle<reco::PFCandidateCollection> pfCandidateHandle;
   e.getByToken(pfCandidates_,pfCandidateHandle);
   if (!pfCandidateHandle.isValid()) {
     edm::LogError("PhotonValidator") << "Error! Can't get the product pfCandidates "<< std::endl ;
   }
   
   edm::Handle<edm::ValueMap<std::vector<reco::PFCandidateRef> > > phoToParticleBasedIsoMapHandle;
   edm::ValueMap<std::vector<reco::PFCandidateRef> > phoToParticleBasedIsoMap;
   if ( fName_ == "zmumugammaGedValidation") {
    e.getByToken(photonIsoValmap_token_,phoToParticleBasedIsoMapHandle);
     //   e.getByLabel("particleBasedIsolation",valueMapPhoPFCandIso_,phoToParticleBasedIsoMapHandle);
     if ( ! phoToParticleBasedIsoMapHandle.isValid()) {
       edm::LogInfo("PhotonValidator") << "Error! Can't get the product: valueMap photons to particle based iso " << std::endl;
       
     }
     phoToParticleBasedIsoMap = *(phoToParticleBasedIsoMapHandle.product());
   }
 
   // Get the reconstructed muons
   bool validMuons=true;
   Handle<reco::MuonCollection> muonHandle;
   reco::MuonCollection muonCollection;
   e.getByToken(muon_token_, muonHandle);
   if ( !muonHandle.isValid()) {
     edm::LogInfo("ZToMuMuGammaAnalyzer") << "Error! Can't get the product: muon_token_" << endl;
     validMuons=false;
   }
   if(validMuons) muonCollection = *(muonHandle.product());
 
   // Get the beam spot
   edm::Handle<reco::BeamSpot> bsHandle;
   e.getByToken(beamSpot_token_, bsHandle);
   if (!bsHandle.isValid()) {
       edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "<< "\n";
       return;
   }
   const reco::BeamSpot &thebs = *bsHandle.product();
 
   //Prepare list of photon-related HLT filter names
   vector<int> Keys;
   for(uint filterIndex=0;filterIndex<triggerEvent.sizeFilters();++filterIndex){  //loop over all trigger filters in event (i.e. filters passed)
     string label = triggerEvent.filterTag(filterIndex).label();
     if(label.find( "Photon" ) != string::npos ) {  //get photon-related filters
       for(uint filterKeyIndex=0;filterKeyIndex<triggerEvent.filterKeys(filterIndex).size();++filterKeyIndex){  //loop over keys to objects passing this filter
         Keys.push_back(triggerEvent.filterKeys(filterIndex)[filterKeyIndex]);  //add keys to a vector for later reference
       }
     }
   }
   
   // sort Keys vector in ascending order
   // and erases duplicate entries from the vector
   sort(Keys.begin(),Keys.end());
   for ( uint i=0; i<Keys.size(); ) {
     if (i!=(Keys.size()-1)) {
       if (Keys[i]==Keys[i+1]) {
         Keys.erase(Keys.begin()+i+1);
       } else {
         ++i;
       }
     } else {
       ++i;
     }
   }
 
   edm::Handle<reco::VertexCollection> vtxH;
   e.getByToken(offline_pvToken_, vtxH);
   h_nRecoVtx_ ->Fill (float(vtxH->size()));
 
   //photon counters
   int nPho = 0;
   int nPhoBarrel = 0;
   int nPhoEndcap = 0;
 
   if ( muonCollection.size() < 2 ) return;
 
   for( reco::MuonCollection::const_iterator  iMu = muonCollection.begin(); iMu != muonCollection.end(); iMu++) {
     if ( !basicMuonSelection (*iMu) ) continue;
  
     for( reco::MuonCollection::const_iterator  iMu2 = iMu+1; iMu2 != muonCollection.end(); iMu2++) {
       if ( !basicMuonSelection (*iMu2) ) continue;
       if ( iMu->charge()*iMu2->charge() > 0) continue;
 
       if ( !muonSelection(*iMu,thebs) && !muonSelection(*iMu2,thebs) ) continue;
     
       float mumuMass = mumuInvMass(*iMu,*iMu2) ;
       if ( mumuMass <  minMumuInvMass_  ||  mumuMass >  maxMumuInvMass_ ) continue;
 
       h1_mumuInvMass_[0] -> Fill (mumuMass);      
 
       if (   photonHandle->size() < 1 ) continue;
 
       reco::Muon nearMuon;
       reco::Muon farMuon;
       for(unsigned int iPho=0; iPho < photonHandle->size(); iPho++) {
     reco::PhotonRef aPho(reco::PhotonRef(photonHandle, iPho));
     //
         double dr1 = deltaR((*iMu).eta(), aPho->eta(), (*iMu).phi(), aPho->phi());
         double dr2 = deltaR((*iMu2).eta(), aPho->eta(), (*iMu2).phi(), aPho->phi());
     double drNear = dr1;
         if (dr1 < dr2) {
       nearMuon =*iMu ; farMuon  = *iMu2; drNear = dr1;
     } else {
       nearMuon = *iMu2; farMuon  = *iMu; drNear = dr2;
     }
     //        
     if ( nearMuon.isolationR03().hadEt > nearMuonHcalIso_ )  continue;
         if ( farMuon.isolationR03().sumPt > farMuonTrackIso_ )  continue;
         if ( farMuon.isolationR03().emEt  > farMuonEcalIso_ )  continue;
         if ( farMuon.pt() < farMuonMinPt_ )       continue;
         if ( drNear > nearMuonDr_)                continue;
     //
         if ( !photonSelection (aPho) ) continue;        
         float mumuGammaMass = mumuGammaInvMass(*iMu,*iMu2,aPho) ;
         if ( mumuGammaMass < minMumuGammaInvMass_ || mumuGammaMass > maxMumuGammaInvMass_ ) continue;
     //
         //counter: number of photons
     int iDet=0;
         if ( aPho->isEB() || aPho->isEE() ) {
       nPho++;
     }
     if ( aPho->isEB()  ) {
       iDet=1; 
       nPhoBarrel++;
     }
         if ( aPho->isEE() ) {
       iDet=2;       
       nPhoEndcap++;           
     }
         
         //PHOTON RELATED HISTOGRAMS
         h1_mumuGammaInvMass_[0] ->Fill (mumuGammaMass);
         h1_mumuGammaInvMass_[iDet] ->Fill (mumuGammaMass);
     //ENERGY        
         h_phoE_[0]  ->Fill (aPho->energy());
     h_phoSigmaEoverE_[0] ->Fill( aPho->getCorrectedEnergyError(aPho->getCandidateP4type())/aPho->energy() ); 
         h_phoEt_[0] ->Fill (aPho->et());
         h_phoE_[iDet]  ->Fill (aPho->energy());
     h_phoSigmaEoverE_[iDet] ->Fill( aPho->getCorrectedEnergyError(aPho->getCandidateP4type())/aPho->energy() ); 
     p_phoSigmaEoverEVsNVtx_[iDet] ->Fill( float(vtxH->size()), aPho->getCorrectedEnergyError(aPho->getCandidateP4type())/aPho->energy() ); 
         h_phoEt_[iDet] ->Fill (aPho->et());
         //GEOMETRICAL
         h_phoEta_[0] ->Fill (aPho->eta());
         h_phoPhi_[0] ->Fill (aPho->phi());
         h_scEta_[0]  ->Fill (aPho->superCluster()->eta());
         h_scPhi_[0]  ->Fill (aPho->superCluster()->phi());
         h_phoEta_[iDet] ->Fill (aPho->eta());
         h_phoPhi_[iDet] ->Fill (aPho->phi());
         h_scEta_[iDet]  ->Fill (aPho->superCluster()->eta());
         h_scPhi_[iDet]  ->Fill (aPho->superCluster()->phi());
         //SHOWER SHAPE
         h_r9_[0]     ->Fill (aPho->r9());
     h_e1x5_[0]->Fill(aPho->e1x5());
     h_e2x5_[0]->Fill(aPho->e2x5());
     h_r1x5_[0]->Fill(aPho->r1x5());
     h_r2x5_[0]->Fill(aPho->r2x5());
         h_phoSigmaIetaIeta_[0]    ->Fill(aPho->sigmaIetaIeta());
     //
         h_r9_[iDet]     ->Fill (aPho->r9());    
     h_e1x5_[iDet] ->Fill(aPho->e1x5()); 
     h_e2x5_[iDet] ->Fill(aPho->e2x5());   
     h_r1x5_[iDet] ->Fill( aPho->r1x5());     
     h_r2x5_[iDet] ->Fill(aPho->r2x5());   
         h_phoSigmaIetaIeta_[iDet]    ->Fill(aPho->sigmaIetaIeta());
         //TRACK ISOLATION
     h_nTrackIsolSolid_[0]      ->Fill(aPho->nTrkSolidConeDR04());
         h_nTrackIsolHollow_[0]      ->Fill(aPho->nTrkHollowConeDR04());    
     h_trackPtSumSolid_[0]       ->Fill(aPho->trkSumPtSolidConeDR04());        
     h_trackPtSumHollow_[0]      ->Fill(aPho->trkSumPtSolidConeDR04());
     h_nTrackIsolSolid_[iDet]      ->Fill(aPho->nTrkSolidConeDR04());
         h_nTrackIsolHollow_[iDet]      ->Fill(aPho->nTrkHollowConeDR04());    
     h_trackPtSumSolid_[iDet]       ->Fill(aPho->trkSumPtSolidConeDR04());        
     h_trackPtSumHollow_[iDet]      ->Fill(aPho->trkSumPtSolidConeDR04());
         //CALORIMETER ISOLATION
     h_ecalSum_[0]      ->Fill(aPho->ecalRecHitSumEtConeDR04());
     h_hcalSum_[0]      ->Fill(aPho->hcalTowerSumEtConeDR04());
     h_hOverE_[0]       ->Fill(aPho->hadTowOverEm());
         h_h1OverE_[0]      ->Fill(aPho->hadTowDepth1OverEm());
         h_h2OverE_[0]      ->Fill(aPho->hadTowDepth2OverEm());
     h_newhOverE_[0]->Fill( aPho->hadTowOverEm());
     h_ecalSum_[iDet]      ->Fill(aPho->ecalRecHitSumEtConeDR04());
     h_hcalSum_[iDet]      ->Fill(aPho->hcalTowerSumEtConeDR04());
     h_hOverE_[iDet]       ->Fill(aPho->hadTowOverEm());
         h_h1OverE_[iDet]      ->Fill(aPho->hadTowDepth1OverEm());
         h_h2OverE_[iDet]      ->Fill(aPho->hadTowDepth2OverEm());
     h_newhOverE_[iDet]->Fill( aPho->hadTowOverEm());
     // Isolation from particle flow
     h_chHadIso_[0]-> Fill (aPho->chargedHadronIso());
     h_nHadIso_[0]-> Fill (aPho->neutralHadronIso());
     h_phoIso_[0]-> Fill (aPho->photonIso());
     h_nCluOutsideMustache_[0]->Fill(float(aPho->nClusterOutsideMustache()));
     h_etOutsideMustache_[0]->Fill(aPho->etOutsideMustache());
     h_pfMva_[0]->Fill(aPho->pfMVA());
     h_chHadIso_[iDet]-> Fill (aPho->chargedHadronIso());
     h_nHadIso_[iDet]-> Fill (aPho->neutralHadronIso());
     h_phoIso_[iDet]-> Fill (aPho->photonIso());
     h_nCluOutsideMustache_[iDet]->Fill(float(aPho->nClusterOutsideMustache()));
     h_etOutsideMustache_[iDet]->Fill(aPho->etOutsideMustache());
     h_pfMva_[iDet]->Fill(aPho->pfMVA());
 
     if ( fName_ == "zmumugammaGedValidation") {
       
       float SumPtIsoValCh = 0.; 
       float SumPtIsoValNh = 0.;
       float SumPtIsoValPh = 0.;
       
       float SumPtIsoValCleanCh = 0.;    
       float SumPtIsoValCleanNh = 0.;
       float SumPtIsoValCleanPh = 0.;
       
       for(unsigned int lCand=0; lCand < pfCandidateHandle->size(); lCand++) {
         reco::PFCandidateRef pfCandRef(reco::PFCandidateRef(pfCandidateHandle,lCand));
         float dR= deltaR(aPho->eta(),  aPho->phi(),pfCandRef->eta(),  pfCandRef->phi()); 
         if ( dR<0.4) {
           reco::PFCandidate::ParticleType type = pfCandRef->particleId();
           if ( type == reco::PFCandidate::e ) continue; 
           if ( type == reco::PFCandidate::gamma && pfCandRef->mva_nothing_gamma() > 0.) continue;
           
           if( type == reco::PFCandidate::h ) {
         SumPtIsoValCh += pfCandRef->pt();
         h_dRPhoPFcand_ChHad_unCleaned_[0]->Fill(dR);
         h_dRPhoPFcand_ChHad_unCleaned_[iDet]->Fill(dR);
           }
           if( type == reco::PFCandidate::h0 ) {
         SumPtIsoValNh += pfCandRef->pt();
         h_dRPhoPFcand_NeuHad_unCleaned_[0]->Fill(dR);
         h_dRPhoPFcand_NeuHad_unCleaned_[iDet]->Fill(dR);
           }
           if( type == reco::PFCandidate::gamma ) {
         SumPtIsoValPh += pfCandRef->pt();
         h_dRPhoPFcand_Pho_unCleaned_[0]->Fill(dR);
         h_dRPhoPFcand_Pho_unCleaned_[iDet]->Fill(dR);
           }
           bool skip=false;
           for( std::vector<reco::PFCandidateRef>::const_iterator i = phoToParticleBasedIsoMap[aPho].begin(); i != phoToParticleBasedIsoMap[aPho].end(); ++i ) {
         //        std::cout << " PhotonValidator PfCand pt " << pfCandRef->pt() << " id " <<pfCandRef->particleId() <<  " and in the map " << (*i)->pt() << " type " << (*i)->particleId() << std::endl;
         if ( (*i) == pfCandRef ) {
           skip=true;
         }
           } // loop over the PFCandidates flagged as overlapping with the photon
           
           if ( skip ) continue;
           if( type == reco::PFCandidate::h ) {
         SumPtIsoValCleanCh += pfCandRef->pt();
         h_dRPhoPFcand_ChHad_Cleaned_[0]->Fill(dR);
         h_dRPhoPFcand_ChHad_Cleaned_[iDet]->Fill(dR);
           }
           if( type == reco::PFCandidate::h0 ) {
         SumPtIsoValCleanNh += pfCandRef->pt();
         h_dRPhoPFcand_NeuHad_Cleaned_[0]->Fill(dR);
         h_dRPhoPFcand_NeuHad_Cleaned_[iDet]->Fill(dR);
           }
           if( type == reco::PFCandidate::gamma ) {
         SumPtIsoValCleanPh += pfCandRef->pt();
         h_dRPhoPFcand_Pho_Cleaned_[0]->Fill(dR);
         h_dRPhoPFcand_Pho_Cleaned_[iDet]->Fill(dR);
           }
         }  // dr=0.4          
       }  // loop over all PF Candidates
       
       h_SumPtOverPhoPt_ChHad_Cleaned_[0]->Fill(SumPtIsoValCleanCh/aPho->pt());
       h_SumPtOverPhoPt_NeuHad_Cleaned_[0]->Fill(SumPtIsoValCleanNh/aPho->pt());
       h_SumPtOverPhoPt_Pho_Cleaned_[0]->Fill(SumPtIsoValCleanPh/aPho->pt());
       h_SumPtOverPhoPt_ChHad_unCleaned_[0]->Fill(SumPtIsoValCh/aPho->pt());
       h_SumPtOverPhoPt_NeuHad_unCleaned_[0]->Fill(SumPtIsoValNh/aPho->pt());
       h_SumPtOverPhoPt_Pho_unCleaned_[0]->Fill(SumPtIsoValPh/aPho->pt());
       //
       h_SumPtOverPhoPt_ChHad_Cleaned_[iDet]->Fill(SumPtIsoValCleanCh/aPho->pt());
       h_SumPtOverPhoPt_NeuHad_Cleaned_[iDet]->Fill(SumPtIsoValCleanNh/aPho->pt());
       h_SumPtOverPhoPt_Pho_Cleaned_[iDet]->Fill(SumPtIsoValCleanPh/aPho->pt());
       h_SumPtOverPhoPt_ChHad_unCleaned_[iDet]->Fill(SumPtIsoValCh/aPho->pt());
       h_SumPtOverPhoPt_NeuHad_unCleaned_[iDet]->Fill(SumPtIsoValNh/aPho->pt());
       h_SumPtOverPhoPt_Pho_unCleaned_[iDet]->Fill(SumPtIsoValPh/aPho->pt());
     } // only for zmumugammaGedValidation
         
     if ( makeProfiles_ ) {
       p_r9VsEt_[iDet] ->Fill (aPho->et(),aPho->r9());
           p_r9VsEta_[0]->Fill (aPho->eta(),aPho->r9());
       p_e1x5VsEt_[iDet] ->Fill(aPho->et(), aPho->e1x5()); 
       p_e1x5VsEta_[0]->Fill(aPho->eta(),aPho->e1x5());
       p_e2x5VsEt_[iDet] ->Fill(aPho->et(), aPho->e2x5());     
       p_e2x5VsEta_[0]->Fill(aPho->eta(),aPho->e2x5());
       p_r1x5VsEt_[iDet] ->Fill(aPho->et(), aPho->r1x5());    
       p_r1x5VsEta_[0]->Fill(aPho->eta(),aPho->r1x5());
       p_r2x5VsEt_[iDet] ->Fill(aPho->et(), aPho->r2x5());     
       p_r2x5VsEta_[0]->Fill(aPho->eta(),aPho->r2x5());
       //
       p_sigmaIetaIetaVsEta_[0]  ->Fill(aPho->eta(),aPho->sigmaIetaIeta());
           p_nTrackIsolSolidVsEt_[iDet]  ->Fill(aPho->et(), aPho->nTrkSolidConeDR04());
           p_nTrackIsolSolidVsEta_[0] ->Fill(aPho->eta(),aPho->nTrkSolidConeDR04());
           p_nTrackIsolHollowVsEt_[iDet]  ->Fill(aPho->et(), aPho->nTrkHollowConeDR04());
           p_nTrackIsolHollowVsEta_[0] ->Fill(aPho->eta(),aPho->nTrkHollowConeDR04());
           p_trackPtSumSolidVsEt_[iDet]   ->Fill(aPho->et(), aPho->trkSumPtSolidConeDR04());
           p_trackPtSumSolidVsEta_[0]  ->Fill(aPho->eta(),aPho->trkSumPtSolidConeDR04());
           p_trackPtSumHollowVsEt_[iDet]  ->Fill(aPho->et(), aPho->trkSumPtHollowConeDR04());
           p_trackPtSumHollowVsEta_[0] ->Fill(aPho->eta(),aPho->trkSumPtHollowConeDR04());
       //
           p_ecalSumVsEt_[iDet]  ->Fill(aPho->et(), aPho->ecalRecHitSumEtConeDR04());
       p_ecalSumVsEta_[0] ->Fill(aPho->eta(),aPho->ecalRecHitSumEtConeDR04());
       p_hcalSumVsEt_[iDet]  ->Fill(aPho->et(), aPho->hcalTowerSumEtConeDR04());
           p_hcalSumVsEta_[0] ->Fill(aPho->eta(),aPho->hcalTowerSumEtConeDR04());
           p_hOverEVsEt_[iDet]   ->Fill(aPho->et(), aPho->hadTowOverEm());    
           p_hOverEVsEta_[0]  ->Fill(aPho->eta(),aPho->hadTowOverEm());
       
     }
 
         if(use2DHistos_){
           //SHOWER SHAPE
           h2_r9VsEt_[iDet] ->Fill (aPho->et(),aPho->r9());
           h2_r9VsEta_[0]->Fill (aPho->eta(),aPho->r9());
       h2_e1x5VsEt_[iDet] ->Fill(aPho->et(), aPho->e1x5()); 
       h2_e1x5VsEta_[0]->Fill(aPho->eta(),aPho->e1x5());
       h2_e2x5VsEta_[0]->Fill(aPho->eta(),aPho->e2x5());
           h2_e2x5VsEt_[iDet] ->Fill(aPho->et(), aPho->e2x5());
           h2_r1x5VsEta_[0]->Fill(aPho->eta(),aPho->r1x5());
           h2_r1x5VsEt_[iDet] ->Fill(aPho->et(), aPho->r1x5());
       h2_r2x5VsEt_[iDet] ->Fill(aPho->et(), aPho->r2x5()); 
       h2_r2x5VsEta_[0]->Fill(aPho->eta(),aPho->r2x5());
       h2_sigmaIetaIetaVsEta_[0]  ->Fill(aPho->eta(),aPho->sigmaIetaIeta());
           //TRACK ISOLATION
           h2_nTrackIsolSolidVsEt_[0]  ->Fill(aPho->et(), aPho->nTrkSolidConeDR04());
           h2_nTrackIsolSolidVsEta_[0] ->Fill(aPho->eta(),aPho->nTrkSolidConeDR04());
           h2_nTrackIsolHollowVsEt_[0]  ->Fill(aPho->et(), aPho->nTrkHollowConeDR04());
           h2_nTrackIsolHollowVsEta_[0] ->Fill(aPho->eta(),aPho->nTrkHollowConeDR04());
           h2_trackPtSumSolidVsEt_[0]   ->Fill(aPho->et(), aPho->trkSumPtSolidConeDR04());
           h2_trackPtSumSolidVsEta_[0]  ->Fill(aPho->eta(),aPho->trkSumPtSolidConeDR04());
           h2_trackPtSumHollowVsEt_[0]  ->Fill(aPho->et(), aPho->trkSumPtHollowConeDR04());
           h2_trackPtSumHollowVsEta_[0] ->Fill(aPho->eta(),aPho->trkSumPtHollowConeDR04());
           //CALORIMETER ISOLATION
           h2_ecalSumVsEt_[iDet]  ->Fill(aPho->et(), aPho->ecalRecHitSumEtConeDR04());
           h2_ecalSumVsEta_[0] ->Fill(aPho->eta(),aPho->ecalRecHitSumEtConeDR04());
       h2_hcalSumVsEt_[iDet]  ->Fill(aPho->et(), aPho->hcalTowerSumEtConeDR04());
           h2_hcalSumVsEta_[0] ->Fill(aPho->eta(),aPho->hcalTowerSumEtConeDR04());
         }
       } //end photon loop
 
       h_nPho_[0] ->Fill (float(nPho));
       h_nPho_[1] ->Fill (float(nPhoBarrel));
       h_nPho_[2] ->Fill (float(nPhoEndcap));
     } //end inner muon loop
   } //end outer muon loop
 }//End of Analyze method

bool ZToMuMuGammaAnalyzer::basicMuonSelection ( const reco::Muon & m )

private

Definition at line 799 of file ZToMuMuGammaAnalyzer.cc.

References reco::LeafCandidate::eta(), reco::Muon::globalTrack(), reco::Muon::innerTrack(), reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), reco::LeafCandidate::pt(), and mps_fire::result.

                                                                    {
   bool result=true;
   if (!mu.innerTrack().isNonnull())    result=false;
   if (!mu.globalTrack().isNonnull())   result=false;
   if ( !mu.isGlobalMuon() )            result=false; 
   if ( mu.pt() < muonMinPt_ )          result=false;
   if ( fabs(mu.eta())>2.4 )            result=false;
 
   int pixHits=0;
   int tkHits=0;
   if ( mu.innerTrack().isNonnull() ) {
     pixHits=mu.innerTrack()->hitPattern().numberOfValidPixelHits();
     tkHits=mu.innerTrack()->hitPattern().numberOfValidStripHits();
   }
 
   if ( pixHits+tkHits < minPixStripHits_ ) result=false;
   
   return result;  
 }

void ZToMuMuGammaAnalyzer::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 110 of file ZToMuMuGammaAnalyzer.cc.

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

 {
   iBooker.setCurrentFolder("Egamma/"+fName_+"/ZToMuMuGamma");
   
   h1_mumuInvMass_[0]      = iBooker.book1D("mumuInvMass","Two muon invariant mass: M (GeV)",etBin_,etMin_,etMax_);
   h1_mumuGammaInvMass_[0] = iBooker.book1D("mumuGammaInvMass","Two-muon plus gamma invariant mass: M (GeV)",etBin_,etMin_,etMax_);
   h1_mumuGammaInvMass_[1] = iBooker.book1D("mumuGammaInvMassBarrel","Two-muon plus gamma invariant mass: M (GeV)",etBin_,etMin_,etMax_);
   h1_mumuGammaInvMass_[2] = iBooker.book1D("mumuGammaInvMassEndcap","Two-muon plus gamma invariant mass: M (GeV)",etBin_,etMin_,etMax_);
 
   h_nRecoVtx_ =  iBooker.book1D("nOfflineVtx","# of Offline Vertices",80, -0.5, 79.5);  
   
   //ENERGY
   h_phoE_[0]  = iBooker.book1D("phoE","Energy;E (GeV)",eBin_,eMin_,eMax_);
   h_phoSigmaEoverE_[0]  = iBooker.book1D("phoSigmaEoverE","All Ecal: #sigma_{E}/E;#sigma_{E}/E",eBin_,eMin_,eMax_);
   h_phoEt_[0] = iBooker.book1D("phoEt","E_{T};E_{T} (GeV)", etBin_,etMin_,etMax_);
 
   //NUMBER OF PHOTONS
   h_nPho_[0] = iBooker.book1D("nPho", "Number of Photons per Event;# #gamma", numberBin_,numberMin_,numberMax_);
 
   //GEOMETRICAL
   h_phoEta_[0] = iBooker.book1D("phoEta", "#eta;#eta",etaBin_,etaMin_,etaMax_);
   h_phoPhi_[0] = iBooker.book1D("phoPhi", "#phi;#phi",phiBin_,phiMin_,phiMax_);
 
   h_scEta_[0]  = iBooker.book1D("scEta", "SuperCluster #eta;#eta",etaBin_,etaMin_,etaMax_);
   h_scPhi_[0]  = iBooker.book1D("scPhi", "SuperCluster #phi;#phi",phiBin_,phiMin_,phiMax_);
 
   //SHOWER SHAPE
   h_r9_[0]      = iBooker.book1D("r9","R9;R9",r9Bin_,r9Min_, r9Max_);
   h_e1x5_[0]  = iBooker.book1D("e1x5","E1x5;E1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_e2x5_[0]  = iBooker.book1D("e2x5","E2x5;E2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r1x5_[0]  = iBooker.book1D("r1x5","r1x5;r1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r2x5_[0]  = iBooker.book1D("r2x5","r2x5;r2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_phoSigmaIetaIeta_[0]   = iBooker.book1D("phoSigmaIetaIeta","#sigma_{i#etai#eta};#sigma_{i#etai#eta}",sigmaIetaBin_,sigmaIetaMin_,sigmaIetaMax_);
   //TRACK ISOLATION
   h_nTrackIsolSolid_[0]       = iBooker.book1D("nIsoTracksSolid","Number Of Tracks in the Solid Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_nTrackIsolHollow_[0]      = iBooker.book1D("nIsoTracksHollow","Number Of Tracks in the Hollow Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_trackPtSumSolid_[0]       = iBooker.book1D("isoPtSumSolid","Track P_{T} Sum in the Solid Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   h_trackPtSumHollow_[0]      = iBooker.book1D("isoPtSumHollow","Track P_{T} Sum in the Hollow Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   //CALORIMETER ISOLATION VARIABLES
   h_ecalSum_[0]      = iBooker.book1D("ecalSum","Ecal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   h_hcalSum_[0]      = iBooker.book1D("hcalSum","Hcal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   h_hOverE_[0]       = iBooker.book1D("hOverE","H/E;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h1OverE_[0]      = iBooker.book1D("h1OverE","H/E for Depth 1;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h2OverE_[0]      = iBooker.book1D("h2OverE","H/E for Depth 2;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   string histname = "newhOverE";
   h_newhOverE_[0] = iBooker.book1D(histname+"All",   "new H/E: All Ecal",100,0., 0.1) ;
   //  Information from Particle Flow 
   histname = "chargedHadIso";
   h_chHadIso_[0]=  iBooker.book1D(histname+"All",   "PF chargedHadIso:  All Ecal",etBin_,etMin_,20.);
   histname = "neutralHadIso";
   h_nHadIso_[0]=  iBooker.book1D(histname+"All",   "PF neutralHadIso:  All Ecal",etBin_,etMin_,20.);
   histname = "photonIso";
   h_phoIso_[0]=  iBooker.book1D(histname+"All",   "PF photonIso:  All Ecal",etBin_,etMin_,20.);
   histname = "nCluOutMustache";
   h_nCluOutsideMustache_[0]= iBooker.book1D(histname+"All",   "PF number of clusters outside Mustache:  All Ecal",50,0.,50.);
   histname = "etOutMustache";
   h_etOutsideMustache_[0]= iBooker.book1D(histname+"All",   "PF et outside Mustache:  All Ecal",etBin_,etMin_,20.);
   histname = "pfMVA";
   h_pfMva_[0]= iBooker.book1D(histname+"All",   "PF MVA output:  All Ecal",50,-1.,2.);
   histname = "SumPtOverPhoPt_ChHad_Cleaned";
   h_SumPtOverPhoPt_ChHad_Cleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Charged Hadrons:  All Ecal",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_NeuHad_Cleaned";
   h_SumPtOverPhoPt_NeuHad_Cleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Neutral Hadrons:  All Ecal",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_Pho_Cleaned";
   h_SumPtOverPhoPt_Pho_Cleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Photons Hadrons:  All Ecal",etBin_,etMin_,2.);
   histname = "dRPhoPFcand_ChHad_Cleaned";
   h_dRPhoPFcand_ChHad_Cleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Charged Hadrons : All Ecal",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_NeuHad_Cleaned";
   h_dRPhoPFcand_NeuHad_Cleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Neutral Hadrons : All Ecal",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_Pho_Cleaned";
   h_dRPhoPFcand_Pho_Cleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Photons : All Ecal",etBin_,etMin_,0.7);
   //
   histname = "SumPtOverPhoPt_ChHad_unCleaned";
   h_SumPtOverPhoPt_ChHad_unCleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Charged Hadrons :  All Ecal",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_NeuHad_unCleaned";
   h_SumPtOverPhoPt_NeuHad_unCleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Neutral Hadrons :  All Ecal",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_Pho_unCleaned";
   h_SumPtOverPhoPt_Pho_unCleaned_[0]=  iBooker.book1D(histname+"All",   "Pf Cand Sum Pt Over photon pt Photons:  All Ecal",etBin_,etMin_,2.);
   histname = "dRPhoPFcand_ChHad_unCleaned";
   h_dRPhoPFcand_ChHad_unCleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Charged Hadrons :  All Ecal",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_NeuHad_unCleaned";
   h_dRPhoPFcand_NeuHad_unCleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Neutral Hadrons :  All Ecal",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_Pho_unCleaned";
   h_dRPhoPFcand_Pho_unCleaned_[0]=  iBooker.book1D(histname+"All",   "dR(pho,cand) Photons:  All Ecal",etBin_,etMin_,0.7);
 
   // NUMBER OF PHOTONS
   h_nPho_[1]  = iBooker.book1D("nPhoBarrel","Number of Photons per Event;# #gamma", numberBin_,numberMin_,numberMax_);
   h_nPho_[2]  = iBooker.book1D("nPhoEndcap","Number of Photons per Event;# #gamma", numberBin_,numberMin_,numberMax_);
   //EB ENERGY
   h_phoE_[1]  = iBooker.book1D("phoEBarrel","Energy for Barrel;E (GeV)",eBin_,eMin_,eMax_);
   h_phoSigmaEoverE_[1]  = iBooker.book1D("phoSigmaEoverEBarrel","Barrel: #sigma_E/E;#sigma_{E}/E",eBin_,eMin_,eMax_);
   h_phoEt_[1] = iBooker.book1D("phoEtBarrel","E_{T};E_{T} (GeV)", etBin_,etMin_,etMax_);
   //EE ENERGY
   h_phoEt_[2] = iBooker.book1D("phoEtEndcap","E_{T};E_{T} (GeV)", etBin_,etMin_,etMax_);
   h_phoE_[2]  = iBooker.book1D("phoEEndcap","Energy for Endcap;E (GeV)",eBin_,eMin_,eMax_);
   h_phoSigmaEoverE_[2]  = iBooker.book1D("phoSigmaEoverEEndcap","Endcap: #sigma_{E}/E;#sigma_{E}/E",eBin_,eMin_,eMax_);
   //EB GEOMETRICAL
   h_phoEta_[1] = iBooker.book1D("phoEtaBarrel","#eta;#eta",etaBin_,etaMin_,etaMax_);
   h_phoPhi_[1] = iBooker.book1D("phoPhiBarrel","#phi;#phi",phiBin_,phiMin_,phiMax_);
   h_scEta_[1]  = iBooker.book1D("scEtaBarrel","SuperCluster #eta;#eta",etaBin_,etaMin_,etaMax_);
   h_scPhi_[1]  = iBooker.book1D("scPhiBarrel","SuperCluster #phi;#phi",phiBin_,phiMin_,phiMax_);
   //EE GEOMETRICAL
   h_phoEta_[2] = iBooker.book1D("phoEtaEndcap","#eta;#eta",etaBin_,etaMin_,etaMax_);
   h_phoPhi_[2] = iBooker.book1D("phoPhiEndcap","#phi;#phi",phiBin_,phiMin_,phiMax_);
   h_scEta_[2]  = iBooker.book1D("scEtaEndcap","SuperCluster #eta;#eta",etaBin_,etaMin_,etaMax_);
   h_scPhi_[2]  = iBooker.book1D("scPhiEndcap","SuperCluster #phi;#phi",phiBin_,phiMin_,phiMax_);
   //SHOWER SHAPES
   h_r9_[1]      = iBooker.book1D("r9Barrel","R9;R9",r9Bin_,r9Min_, r9Max_);
   h_r9_[2]      = iBooker.book1D("r9Endcap","R9;R9",r9Bin_,r9Min_, r9Max_);
   h_e1x5_[1]  = iBooker.book1D("e1x5Barrel","E1x5;E1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_e1x5_[2]  = iBooker.book1D("e1x5Endcap","E1x5;E1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_e2x5_[1]  = iBooker.book1D("e2x5Barrel","E2x5;E2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_e2x5_[2]  = iBooker.book1D("e2x5Endcap","E2x5;E2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r1x5_[1]  = iBooker.book1D("r1x5Barrel","r1x5;r1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r1x5_[2]  = iBooker.book1D("r1x5Endcap","r1x5;r1X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r2x5_[1]  = iBooker.book1D("r2x5Barrel","r2x5;r2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_r2x5_[2]  = iBooker.book1D("r2x5Endcap","r2x5;r2X5 (GeV)",reducedEtBin_,etMin_,etMax_);
   h_phoSigmaIetaIeta_[1]   = iBooker.book1D("phoSigmaIetaIetaBarrel","#sigma_{i#etai#eta};#sigma_{i#etai#eta}",sigmaIetaBin_,sigmaIetaMin_,sigmaIetaMax_);
   h_phoSigmaIetaIeta_[2]   = iBooker.book1D("phoSigmaIetaIetaEndcap","#sigma_{i#etai#eta};#sigma_{i#etai#eta}",sigmaIetaBin_,sigmaIetaMin_,sigmaIetaMax_);
   // TRACK ISOLATION
   h_nTrackIsolSolid_[1]       = iBooker.book1D("nIsoTracksSolidBarrel","Number Of Tracks in the Solid Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_nTrackIsolSolid_[2]       = iBooker.book1D("nIsoTracksSolidEndcap","Number Of Tracks in the Solid Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_nTrackIsolHollow_[1]      = iBooker.book1D("nIsoTracksHollowBarrel","Number Of Tracks in the Hollow Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_nTrackIsolHollow_[2]      = iBooker.book1D("nIsoTracksHollowEndcap","Number Of Tracks in the Hollow Iso Cone;# tracks",numberBin_,numberMin_,numberMax_);
   h_trackPtSumSolid_[1]       = iBooker.book1D("isoPtSumSolidBarrel","Track P_{T} Sum in the Solid Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   h_trackPtSumSolid_[2]       = iBooker.book1D("isoPtSumSolidEndcap","Track P_{T} Sum in the Solid Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   h_trackPtSumHollow_[1]      = iBooker.book1D("isoPtSumHollowBarrel","Track P_{T} Sum in the Hollow Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   h_trackPtSumHollow_[2]      = iBooker.book1D("isoPtSumHollowEndcap","Track P_{T} Sum in the Hollow Iso Cone;P_{T} (GeV)",sumBin_,sumMin_,sumMax_);
   // CALORIMETER ISOLATION VARIABLES
   h_ecalSum_[1]      = iBooker.book1D("ecalSumBarrel","Ecal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   h_ecalSum_[2]      = iBooker.book1D("ecalSumEndcap","Ecal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   h_hcalSum_[1]      = iBooker.book1D("hcalSumBarrel","Hcal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   h_hcalSum_[2]      = iBooker.book1D("hcalSumEndcap","Hcal Sum in the Iso Cone;E (GeV)",sumBin_,sumMin_,sumMax_);
   //H/E
   // EB
   h_hOverE_[1]       = iBooker.book1D("hOverEBarrel","H/E;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h1OverE_[1]      = iBooker.book1D("h1OverEBarrel","H/E for Depth 1;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h2OverE_[1]      = iBooker.book1D("h2OverEBarrel","H/E for Depth 2;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   histname = "newhOverE";
   h_newhOverE_[1] = iBooker.book1D(histname+"Barrel",   "new H/E: Barrel",100,0., 0.1) ;
   //EE 
   h_hOverE_[2]       = iBooker.book1D("hOverEEndcap","H/E;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h1OverE_[2]      = iBooker.book1D("h1OverEEndcap","H/E for Depth 1;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   h_h2OverE_[2]      = iBooker.book1D("h2OverEEndcap","H/E for Depth 2;H/E",hOverEBin_,hOverEMin_,hOverEMax_);
   histname = "newhOverE";
   h_newhOverE_[2] = iBooker.book1D(histname+"Endcap",   "new H/E: Endcap",100,0., 0.1) ;
   // Information from Particle Flow
   histname = "chargedHadIso";
   h_chHadIso_[1]=  iBooker.book1D(histname+"Barrel",   "PF chargedHadIso:  Barrel",etBin_,etMin_,20.);
   h_chHadIso_[2]=  iBooker.book1D(histname+"Endcap",   "PF chargedHadIso:  Endcap",etBin_,etMin_,20.);
   histname = "neutralHadIso";
   h_nHadIso_[1]=  iBooker.book1D(histname+"Barrel",   "PF neutralHadIso:  Barrel",etBin_,etMin_,20.);
   h_nHadIso_[2]=  iBooker.book1D(histname+"Endcap",   "PF neutralHadIso:  Endcap",etBin_,etMin_,20.);
   histname = "photonIso";
   h_phoIso_[1]=  iBooker.book1D(histname+"Barrel",   "PF photonIso:  Barrel",etBin_,etMin_,20.);
   h_phoIso_[2]=  iBooker.book1D(histname+"Endcap",   "PF photonIso:  Endcap",etBin_,etMin_,20.);
   histname = "nCluOutMustache";
   h_nCluOutsideMustache_[1]= iBooker.book1D(histname+"Barrel",   "PF number of clusters outside Mustache:  Barrel",50,0.,50.);
   h_nCluOutsideMustache_[2]= iBooker.book1D(histname+"Endcap",   "PF number of clusters outside Mustache:  Endcap",50,0.,50.);
   histname = "etOutMustache";
   h_etOutsideMustache_[1]= iBooker.book1D(histname+"Barrel",   "PF et outside Mustache:  Barrel",etBin_,etMin_,20.);
   h_etOutsideMustache_[2]= iBooker.book1D(histname+"Endcap",   "PF et outside Mustache:  Endcap",etBin_,etMin_,20.);
   histname = "pfMVA";
   h_pfMva_[1]= iBooker.book1D(histname+"Barrel",   "PF MVA output:  Barrel",50,-1.,2.);
   h_pfMva_[2]= iBooker.book1D(histname+"Endcap",   "PF MVA output:  Endcap",50,-1,2.);
   histname = "SumPtOverPhoPt_ChHad_Cleaned";
   h_SumPtOverPhoPt_ChHad_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Charged Hadrons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_ChHad_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Charged Hadrons:  Endcap",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_NeuHad_Cleaned";
   h_SumPtOverPhoPt_NeuHad_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Neutral Hadrons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_NeuHad_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Neutral Hadrons:  Endcap",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_Pho_Cleaned";
   h_SumPtOverPhoPt_Pho_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Photons Hadrons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_Pho_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Photons Hadrons:  Endcap",etBin_,etMin_,2.);
   histname = "dRPhoPFcand_ChHad_Cleaned";
   h_dRPhoPFcand_ChHad_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Charged Hadrons :  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_ChHad_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Charged Hadrons :  Endcap",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_NeuHad_Cleaned";
   h_dRPhoPFcand_NeuHad_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Neutral Hadrons :  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_NeuHad_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Neutral Hadrons :  Endcap",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_Pho_Cleaned";
   h_dRPhoPFcand_Pho_Cleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Photons :  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_Pho_Cleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Photons :  Endcap",etBin_,etMin_,0.7);
   //
   histname = "SumPtOverPhoPt_ChHad_unCleaned";
   h_SumPtOverPhoPt_ChHad_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Charged Hadrons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_ChHad_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Charged Hadrons:  Endcap",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_NeuHad_unCleaned";
   h_SumPtOverPhoPt_NeuHad_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Neutral Hadrons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_NeuHad_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Neutral Hadrons:  Endcap",etBin_,etMin_,2.);
   histname = "SumPtOverPhoPt_Pho_unCleaned";
   h_SumPtOverPhoPt_Pho_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","PF Cand Sum Pt Over photon pt Photons:  Barrel",etBin_,etMin_,2.);
   h_SumPtOverPhoPt_Pho_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","PF Cand Sum Pt Over photon pt Photons:  Endcap",etBin_,etMin_,2.);
   histname = "dRPhoPFcand_ChHad_unCleaned";
   h_dRPhoPFcand_ChHad_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Charged Hadrons :  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_ChHad_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Charged Hadrons :  Endcap",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_NeuHad_unCleaned";
   h_dRPhoPFcand_NeuHad_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Neutral Hadrons :  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_NeuHad_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Neutral Hadrons :  Endcap",etBin_,etMin_,0.7);
   histname = "dRPhoPFcand_Pho_unCleaned";
   h_dRPhoPFcand_Pho_unCleaned_[1]=  iBooker.book1D(histname+"Barrel","dR(pho,cand) Photons:  Barrel",etBin_,etMin_,0.7);
   h_dRPhoPFcand_Pho_unCleaned_[2]=  iBooker.book1D(histname+"Endcap","dR(pho,cand) Photons:  Endcap",etBin_,etMin_,0.7);
 
   if(makeProfiles_){
     //     p_r9VsEt_[0]  = iBooker.bookProfile("r9VsEt","Avg R9 vs E_{T};E_{T} (GeV);R9",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r9VsEt_[1]  = iBooker.bookProfile("r9VsEtBarrel","Avg R9 vs E_{T};E_{T} (GeV);R9",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r9VsEt_[2]  = iBooker.bookProfile("r9VsEtEndcap","Avg R9 vs E_{T};E_{T} (GeV);R9",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r9VsEta_[0] = iBooker.bookProfile("r9VsEta","Avg R9 vs #eta;#eta;R9",etaBin_,etaMin_,etaMax_,r9Bin_,r9Min_,r9Max_);
     //
     p_sigmaIetaIetaVsEta_[0] = iBooker.bookProfile("sigmaIetaIetaVsEta","Avg #sigma_{i#etai#eta} vs #eta;#eta;#sigma_{i#etai#eta}",etaBin_,etaMin_,etaMax_,sigmaIetaBin_,sigmaIetaMin_,sigmaIetaMax_);
     p_e1x5VsEt_[1]  = iBooker.bookProfile("e1x5VsEtBarrel","Avg E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",etBin_,etMin_,etMax_,etBin_,etMin_,etMax_);
     p_e1x5VsEt_[2]  = iBooker.bookProfile("e1x5VsEtEndcap","Avg E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",etBin_,etMin_,etMax_,etBin_,etMin_,etMax_);
     p_e1x5VsEta_[0] = iBooker.bookProfile("e1x5VsEta","Avg E1x5 vs #eta;#eta;E1X5 (GeV)",etaBin_,etaMin_,etaMax_,etBin_,etMin_,etMax_);
     p_e2x5VsEt_[1]  = iBooker.bookProfile("e2x5VsEtBarrel","Avg E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",etBin_,etMin_,etMax_,etBin_,etMin_,etMax_);
     p_e2x5VsEt_[2]  = iBooker.bookProfile("e2x5VsEtEndcap","Avg E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",etBin_,etMin_,etMax_,etBin_,etMin_,etMax_);
     p_e2x5VsEta_[0] = iBooker.bookProfile("e2x5VsEta","Avg E2x5 vs #eta;#eta;E2X5 (GeV)",etaBin_,etaMin_,etaMax_,etBin_,etMin_,etMax_);
     p_r1x5VsEt_[1]  = iBooker.bookProfile("r1x5VsEtBarrel","Avg R1x5 vs E_{T};E_{T} (GeV);R1X5",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r1x5VsEt_[2]  = iBooker.bookProfile("r1x5VsEtEndcap","Avg R1x5 vs E_{T};E_{T} (GeV);R1X5",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r1x5VsEta_[0] = iBooker.bookProfile("r1x5VsEta","Avg R1x5 vs #eta;#eta;R1X5",etaBin_,etaMin_,etaMax_,r9Bin_,r9Min_,r9Max_);
     p_r2x5VsEt_[1]  = iBooker.bookProfile("r2x5VsEtBarrel","Avg R2x5 vs E_{T};E_{T} (GeV);R2X5",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r2x5VsEt_[2]  = iBooker.bookProfile("r2x5VsEtEndcap","Avg R2x5 vs E_{T};E_{T} (GeV);R2X5",etBin_,etMin_,etMax_,r9Bin_,r9Min_,r9Max_);
     p_r2x5VsEta_[0] = iBooker.bookProfile("r2x5VsEta","Avg R2x5 vs #eta;#eta;R2X5",etaBin_,etaMin_,etaMax_,r9Bin_,r9Min_,r9Max_);
     p_nTrackIsolSolidVsEt_[1]   = iBooker.bookProfile("nIsoTracksSolidVsEtBarrel","Avg Number Of Tracks in the Solid Iso Cone vs E_{T};E_{T};# tracks",etBin_,etMin_,etMax_,numberBin_,numberMin_,numberMax_);
     p_nTrackIsolSolidVsEt_[2]   = iBooker.bookProfile("nIsoTracksSolidVsEtEndcap","Avg Number Of Tracks in the Solid Iso Cone vs E_{T};E_{T};# tracks",etBin_,etMin_,etMax_,numberBin_,numberMin_,numberMax_);
     p_nTrackIsolSolidVsEta_[0]  = iBooker.bookProfile("nIsoTracksSolidVsEta","Avg Number Of Tracks in the Solid Iso Cone vs #eta;#eta;# tracks",etaBin_,etaMin_, etaMax_,numberBin_,numberMin_,numberMax_);
     p_nTrackIsolHollowVsEt_[1]  = iBooker.bookProfile("nIsoTracksHollowVsEtBarrel","Avg Number Of Tracks in the Hollow Iso Cone vs E_{T};E_{T};# tracks",etBin_,etMin_,etMax_,numberBin_,numberMin_,numberMax_);
     p_nTrackIsolHollowVsEt_[2]  = iBooker.bookProfile("nIsoTracksHollowVsEtEndcap","Avg Number Of Tracks in the Hollow Iso Cone vs E_{T};E_{T};# tracks",etBin_,etMin_,etMax_,numberBin_,numberMin_,numberMax_);
     p_nTrackIsolHollowVsEta_[0] = iBooker.bookProfile("nIsoTracksHollowVsEta","Avg Number Of Tracks in the Hollow Iso Cone vs #eta;#eta;# tracks",etaBin_,etaMin_, etaMax_,numberBin_,numberMin_,numberMax_);
     p_trackPtSumSolidVsEt_[1]   = iBooker.bookProfile("isoPtSumSolidVsEtBarrel","Avg Track P_{T} Sum in the Solid Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",etBin_,etMin_,etMax_,sumBin_,sumMin_,sumMax_);
     p_trackPtSumSolidVsEt_[2]   = iBooker.bookProfile("isoPtSumSolidVsEtEndcap","Avg Track P_{T} Sum in the Solid Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",etBin_,etMin_,etMax_,sumBin_,sumMin_,sumMax_);
     p_trackPtSumSolidVsEta_[0]  = iBooker.bookProfile("isoPtSumSolidVsEta","Avg Track P_{T} Sum in the Solid Iso Cone vs #eta;#eta;P_{T} (GeV)",etaBin_,etaMin_, etaMax_,sumBin_,sumMin_,sumMax_);
     p_trackPtSumHollowVsEt_[1]  = iBooker.bookProfile("isoPtSumHollowVsEtBarrel","Avg Track P_{T} Sum in the Hollow Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",etBin_,etMin_,etMax_,sumBin_,sumMin_,sumMax_);
     p_trackPtSumHollowVsEt_[2]  = iBooker.bookProfile("isoPtSumHollowVsEtEndcap","Avg Track P_{T} Sum in the Hollow Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",etBin_,etMin_,etMax_,sumBin_,sumMin_,sumMax_);
     p_trackPtSumHollowVsEta_[0] = iBooker.bookProfile("isoPtSumHollowVsEta","Avg Track P_{T} Sum in the Hollow Iso Cone vs #eta;#eta;P_{T} (GeV)",etaBin_,etaMin_, etaMax_,sumBin_,sumMin_,sumMax_);
     p_ecalSumVsEt_[1]  = iBooker.bookProfile("ecalSumVsEtBarrel","Avg Ecal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",etBin_,etMin_, etMax_,sumBin_,sumMin_,sumMax_);
     p_ecalSumVsEt_[2]  = iBooker.bookProfile("ecalSumVsEtEndcap","Avg Ecal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",etBin_,etMin_, etMax_,sumBin_,sumMin_,sumMax_);
     p_ecalSumVsEta_[0] = iBooker.bookProfile("ecalSumVsEta","Avg Ecal Sum in the Iso Cone vs #eta;#eta;E (GeV)",etaBin_,etaMin_, etaMax_,sumBin_,sumMin_,sumMax_);
     p_hcalSumVsEt_[1]  = iBooker.bookProfile("hcalSumVsEtBarrel","Avg Hcal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",etBin_,etMin_, etMax_,sumBin_,sumMin_,sumMax_);
     p_hcalSumVsEt_[2]  = iBooker.bookProfile("hcalSumVsEtEndcap","Avg Hcal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",etBin_,etMin_, etMax_,sumBin_,sumMin_,sumMax_);
     p_hcalSumVsEta_[0] = iBooker.bookProfile("hcalSumVsEta","Avg Hcal Sum in the Iso Cone vs #eta;#eta;E (GeV)",etaBin_,etaMin_, etaMax_,sumBin_,sumMin_,sumMax_);
     p_hOverEVsEt_[1]   = iBooker.bookProfile("p_hOverEVsEtBarrel","Avg H/E vs Et;E_{T} (GeV);H/E",etBin_,etMin_,etMax_,hOverEBin_,hOverEMin_,hOverEMax_);
     p_hOverEVsEt_[2]   = iBooker.bookProfile("p_hOverEVsEtEndcap","Avg H/E vs Et;E_{T} (GeV);H/E",etBin_,etMin_,etMax_,hOverEBin_,hOverEMin_,hOverEMax_);
     p_hOverEVsEta_[0]  = iBooker.bookProfile("p_hOverEVsEta","Avg H/E vs #eta;#eta;H/E",etaBin_,etaMin_,etaMax_,hOverEBin_,hOverEMin_,hOverEMax_);
 
     // sigmaE/E
     histname = "sigmaEoverEVsNVtx";
     p_phoSigmaEoverEVsNVtx_[1] = iBooker.bookProfile(histname+"Barrel","Photons #sigma_{E}/E vs N_{vtx}: Barrel; N_{vtx}; #sigma_{E}/E ",80, -0.5, 79.5, 100,0., 0.08, "");
     p_phoSigmaEoverEVsNVtx_[2] = iBooker.bookProfile(histname+"Endcap","Photons #sigma_{E}/E vs N_{vtx}: Endcap;  N_{vtx}; #sigma_{E}/E",80, -0.5, 79.5, 100,0., 0.08, "");
   }
 
   if(use2DHistos_){
     //SHOWER SHAPE
     //r9    
     h2_r9VsEt_[0]  = iBooker.book2D("r9VsEt2D","R9 vs E_{T};E_{T} (GeV);R9",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r9VsEt_[1]  = iBooker.book2D("r9VsEt2DBarrel","R9 vs E_{T};E_{T} (GeV);R9",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r9VsEt_[2]  = iBooker.book2D("r9VsEt2DEndcap","R9 vs E_{T};E_{T} (GeV);R9",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r9VsEta_[0] = iBooker.book2D("r9VsEta2D","R9 vs #eta;#eta;R9",reducedEtaBin_,etaMin_,etaMax_,reducedR9Bin_,r9Min_,r9Max_);
     //sigmaIetaIeta
     h2_sigmaIetaIetaVsEta_[0] = iBooker.book2D("sigmaIetaIetaVsEta2D","#sigma_{i#etai#eta} vs #eta;#eta;#sigma_{i#etai#eta}",reducedEtaBin_,etaMin_,etaMax_,sigmaIetaBin_,sigmaIetaMin_,sigmaIetaMax_);
     //e1x5
     h2_e1x5VsEt_[0]  = iBooker.book2D("e1x5VsEt2D","E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e1x5VsEt_[1]  = iBooker.book2D("e1x5VsEt2DBarrel","E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e1x5VsEt_[2]  = iBooker.book2D("e1x5VsEt2DEndcap","E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e1x5VsEta_[0] = iBooker.book2D("e1x5VsEta2D","E1x5 vs #eta;#eta;E1X5 (GeV)",reducedEtaBin_,etaMin_,etaMax_,reducedEtBin_,etMin_,etMax_);
     //e2x5
     h2_e2x5VsEt_[0]  = iBooker.book2D("e2x5VsEt2D","E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e2x5VsEt_[1]  = iBooker.book2D("e2x5VsEt2DBarrel","E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e2x5VsEt_[2]  = iBooker.book2D("e2x5VsEt2DEndcap","E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",reducedEtBin_,etMin_,etMax_,reducedEtBin_,etMin_,etMax_);
     h2_e2x5VsEta_[0] = iBooker.book2D("e2x5VsEta2D","E2x5 vs #eta;#eta;E2X5 (GeV)",reducedEtaBin_,etaMin_,etaMax_,reducedEtBin_,etMin_,etMax_);
     //r1x5
     h2_r1x5VsEt_[0]  = iBooker.book2D("r1x5VsEt2D","R1x5 vs E_{T};E_{T} (GeV);R1X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r1x5VsEt_[1]  = iBooker.book2D("r1x5VsEt2DBarrel","R1x5 vs E_{T};E_{T} (GeV);R1X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r1x5VsEt_[2]  = iBooker.book2D("r1x5VsEt2DEndcap","R1x5 vs E_{T};E_{T} (GeV);R1X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r1x5VsEta_[0] = iBooker.book2D("r1x5VsEta2D","R1x5 vs #eta;#eta;R1X5",reducedEtaBin_,etaMin_,etaMax_,reducedR9Bin_,r9Min_,r9Max_);
     //r2x5
     h2_r2x5VsEt_[0]  = iBooker.book2D("r2x5VsEt2D","R2x5 vs E_{T};E_{T} (GeV);R2X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r2x5VsEt_[1]  = iBooker.book2D("r2x5VsEt2DBarrel","R2x5 vs E_{T};E_{T} (GeV);R2X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r2x5VsEt_[2]  = iBooker.book2D("r2x5VsEt2DEndcap","R2x5 vs E_{T};E_{T} (GeV);R2X5",reducedEtBin_,etMin_,etMax_,reducedR9Bin_,r9Min_,r9Max_);
     h2_r2x5VsEta_[0] = iBooker.book2D("r2x5VsEta2D","R2x5 vs #eta;#eta;R2X5",reducedEtaBin_,etaMin_,etaMax_,reducedR9Bin_,r9Min_,r9Max_);
      //TRACK ISOLATION
     //nTrackIsolSolid
     h2_nTrackIsolSolidVsEt_[0]   = iBooker.book2D("nIsoTracksSolidVsEt2D","Number Of Tracks in the Solid Iso Cone vs E_{T};E_{T};# tracks",reducedEtBin_,etMin_, etMax_,numberBin_,numberMin_,numberMax_);
     h2_nTrackIsolSolidVsEta_[0]  = iBooker.book2D("nIsoTracksSolidVsEta2D","Number Of Tracks in the Solid Iso Cone vs #eta;#eta;# tracks",reducedEtaBin_,etaMin_, etaMax_,numberBin_,numberMin_,numberMax_);
     //nTrackIsolHollow
     h2_nTrackIsolHollowVsEt_[0]  = iBooker.book2D("nIsoTracksHollowVsEt2D","Number Of Tracks in the Hollow Iso Cone vs E_{T};E_{T};# tracks",reducedEtBin_,etMin_, etMax_,numberBin_,numberMin_,numberMax_);
     h2_nTrackIsolHollowVsEta_[0] = iBooker.book2D("nIsoTracksHollowVsEta2D","Number Of Tracks in the Hollow Iso Cone vs #eta;#eta;# tracks",reducedEtaBin_,etaMin_, etaMax_,numberBin_,numberMin_,numberMax_);
     //trackPtSumSolid
     h2_trackPtSumSolidVsEt_[0]   = iBooker.book2D("isoPtSumSolidVsEt2D","Track P_{T} Sum in the Solid Iso Cone;E_{T} (GeV);P_{T} (GeV)",reducedEtBin_,etMin_, etMax_,reducedSumBin_,sumMin_,sumMax_);
     h2_trackPtSumSolidVsEta_[0]  = iBooker.book2D("isoPtSumSolidVsEta2D","Track P_{T} Sum in the Solid Iso Cone;#eta;P_{T} (GeV)",reducedEtaBin_,etaMin_, etaMax_,reducedSumBin_,sumMin_,sumMax_);
     //trackPtSumHollow
     h2_trackPtSumHollowVsEt_[0]  = iBooker.book2D("isoPtSumHollowVsEt2D","Track P_{T} Sum in the Hollow Iso Cone;E_{T} (GeV);P_{T} (GeV)",reducedEtBin_,etMin_, etMax_,reducedSumBin_,sumMin_,sumMax_);
     h2_trackPtSumHollowVsEta_[0] = iBooker.book2D("isoPtSumHollowVsEta2D","Track P_{T} Sum in the Hollow Iso Cone;#eta;P_{T} (GeV)",reducedEtaBin_,etaMin_, etaMax_,reducedSumBin_,sumMin_,sumMax_);
     //CALORIMETER ISOLATION VARIABLES
     //ecal sum
     h2_ecalSumVsEt_[0]  = iBooker.book2D("ecalSumVsEt2D","Ecal Sum in the Iso Cone;E_{T} (GeV);E (GeV)",reducedEtBin_,etMin_, etMax_,reducedSumBin_,sumMin_,sumMax_);
     h2_ecalSumVsEta_[0] = iBooker.book2D("ecalSumVsEta2D","Ecal Sum in the Iso Cone;#eta;E (GeV)",reducedEtaBin_,etaMin_, etaMax_,reducedSumBin_,sumMin_,sumMax_);
     //hcal sum
     h2_hcalSumVsEt_[0]  = iBooker.book2D("hcalSumVsEt2D","Hcal Sum in the Iso Cone;E_{T} (GeV);E (GeV)",reducedEtBin_,etMin_, etMax_,reducedSumBin_,sumMin_,sumMax_);
     h2_hcalSumVsEta_[0] = iBooker.book2D("hcalSumVsEta2D","Hcal Sum in the Iso Cone;#eta;E (GeV)",reducedEtaBin_,etaMin_, etaMax_,reducedSumBin_,sumMin_,sumMax_);
   }
 }

float ZToMuMuGammaAnalyzer::mumuGammaInvMass	(	const reco::Muon &	mu1,
		const reco::Muon &	mu2,
		const reco::PhotonRef &	pho
	)

private

Definition at line 869 of file ZToMuMuGammaAnalyzer.cc.

References reco::LeafCandidate::p4(), and mathSSE::sqrt().

                                                                                                                  {
   math::XYZTLorentzVector p12 = mu1.p4()+mu2.p4()+pho->p4() ;
   float Mass2 = p12.Dot(p12) ;
   float invMass = sqrt(Mass2) ;
   return invMass ;
 }

float ZToMuMuGammaAnalyzer::mumuInvMass	(	const reco::Muon &	m1,
		const reco::Muon &	m2
	)

private

Definition at line 862 of file ZToMuMuGammaAnalyzer.cc.

References reco::LeafCandidate::p4(), and mathSSE::sqrt().

                                                                                  {
   math::XYZTLorentzVector p12 = mu1.p4()+mu2.p4() ;
   float mumuMass2 = p12.Dot(p12) ;
   float invMass = sqrt(mumuMass2) ;
   return invMass ;
 }

bool ZToMuMuGammaAnalyzer::muonSelection	(	const reco::Muon &	m,
		const reco::BeamSpot &	bs
	)

private

Definition at line 819 of file ZToMuMuGammaAnalyzer.cc.

References reco::LeafCandidate::eta(), reco::Muon::globalTrack(), reco::Muon::isolationR03(), reco::Muon::isTrackerMuon(), reco::Muon::numberOfMatches(), mps_fire::result, reco::MuonIsolation::sumPt, and reco::Muon::track().

                                                                                              {
   bool result=true;
   if ( mu.globalTrack()->normalizedChi2() > muonMaxChi2_ )          result=false;
   if ( fabs( mu.globalTrack()->dxy(beamSpot)) > muonMaxDxy_ )       result=false;
   if ( mu.numberOfMatches() < muonMatches_ )                                   result=false;
 
   if ( mu.track()-> hitPattern().numberOfValidPixelHits() <  validPixHits_ )     result=false;
   if ( mu.globalTrack()->hitPattern().numberOfValidMuonHits() < validMuonHits_ ) result=false;
   if ( !mu.isTrackerMuon() )                                        result=false;
   // track isolation 
   if ( mu.isolationR03().sumPt > muonTrackIso_ )                                result=false;
   if ( fabs(mu.eta())>  muonTightEta_ )                                         result=false;
  
   return result;  
 }

bool ZToMuMuGammaAnalyzer::photonSelection ( const reco::PhotonRef & p )

private

remove after moriond if (EtCorrEcalIso>4.0) result=false;

remove after moriond if (EtCorrEcalIso>50.0) result=false;

Definition at line 835 of file ZToMuMuGammaAnalyzer.cc.

References funct::false, and mps_fire::result.

                                                                       {
   bool result=true;
   if ( pho->pt() < photonMinEt_ )          result=false;
   if ( fabs(pho->eta())> photonMaxEta_ )   result=false;
   if ( pho->isEBEEGap() )       result=false;
 
   double EtCorrHcalIso = pho->hcalTowerSumEtConeDR03() - 0.005*pho->pt();
   double EtCorrTrkIso  = pho->trkSumPtHollowConeDR03() - 0.002*pho->pt();
 
   if (pho->r9() <=0.9) {
     if (pho->isEB() && (pho->hadTowOverEm()>0.075 || pho->sigmaIetaIeta() > 0.014)) result=false;
     if (pho->isEE() && (pho->hadTowOverEm()>0.075 || pho->sigmaIetaIeta() > 0.034)) result=false;
     if (EtCorrHcalIso>4.0) result=false;
     if (EtCorrTrkIso>4.0) result=false ;
     if ( pho->chargedHadronIso()  > 4 )  result=false;
   } else {
     if (pho->isEB() && (pho->hadTowOverEm()>0.082 || pho->sigmaIetaIeta() > 0.014)) result=false;
     if (pho->isEE() && (pho->hadTowOverEm()>0.075 || pho->sigmaIetaIeta() > 0.034)) result=false;
     if (EtCorrHcalIso>50.0) result=false;
     if (EtCorrTrkIso>50.0) result=false;
     if ( pho->chargedHadronIso()  > 4 )  result=false;
   }
   return result;  
 }

Member Data Documentation

edm::EDGetTokenT<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > ZToMuMuGammaAnalyzer::barrelRecHit_token_

private

Definition at line 113 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<reco::BeamSpot> ZToMuMuGammaAnalyzer::beamSpot_token_

private

Definition at line 116 of file ZToMuMuGammaAnalyzer.h.

std::stringstream ZToMuMuGammaAnalyzer::currentFolder_

private

Definition at line 125 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::eBin_

private

Definition at line 158 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::eMax_

private

Definition at line 157 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::eMin_

private

Definition at line 156 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > ZToMuMuGammaAnalyzer::endcapRecHit_token_

private

Definition at line 114 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::etaBin_

private

Definition at line 170 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::etaMax_

private

Definition at line 169 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::etaMin_

private

Definition at line 168 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::etBin_

private

Definition at line 162 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::etMax_

private

Definition at line 161 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::etMin_

private

Definition at line 160 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::farMuonEcalIso_

private

Definition at line 149 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::farMuonMinPt_

private

Definition at line 151 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::farMuonTrackIso_

private

Definition at line 150 of file ZToMuMuGammaAnalyzer.h.

std::string ZToMuMuGammaAnalyzer::fName_

private

Definition at line 121 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h1_mumuGammaInvMass_[3]

private

Definition at line 206 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h1_mumuInvMass_[3]

private

photon histos

Definition at line 205 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_e1x5VsEt_[3]

private

Definition at line 229 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_e1x5VsEta_[3]

private

Definition at line 227 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_e2x5VsEt_[3]

private

Definition at line 235 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_e2x5VsEta_[3]

private

Definition at line 233 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_ecalSumVsEt_[3]

private

Definition at line 279 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_ecalSumVsEta_[3]

private

Definition at line 281 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_hcalSumVsEt_[3]

private

Definition at line 285 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_hcalSumVsEta_[3]

private

Definition at line 287 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_nTrackIsolHollowVsEt_[3]

private

Definition at line 261 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_nTrackIsolHollowVsEta_[3]

private

Definition at line 263 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_nTrackIsolSolidVsEt_[3]

private

Definition at line 255 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_nTrackIsolSolidVsEta_[3]

private

Definition at line 257 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r1x5VsEt_[3]

private

Definition at line 241 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r1x5VsEta_[3]

private

Definition at line 239 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r2x5VsEt_[3]

private

Definition at line 247 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r2x5VsEta_[3]

private

Definition at line 245 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r9VsEt_[3]

private

Definition at line 221 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_r9VsEta_[3]

private

Definition at line 223 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_sigmaIetaIetaVsEta_[3]

private

Definition at line 251 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_trackPtSumHollowVsEt_[3]

private

Definition at line 273 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_trackPtSumHollowVsEta_[3]

private

Definition at line 275 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_trackPtSumSolidVsEt_[3]

private

Definition at line 267 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h2_trackPtSumSolidVsEta_[3]

private

Definition at line 269 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_chHadIso_[3]

private

Definition at line 301 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_ChHad_Cleaned_[3]

private

Definition at line 309 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_ChHad_unCleaned_[3]

private

Definition at line 312 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_NeuHad_Cleaned_[3]

private

Definition at line 310 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_NeuHad_unCleaned_[3]

private

Definition at line 313 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_Pho_Cleaned_[3]

private

Definition at line 311 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_dRPhoPFcand_Pho_unCleaned_[3]

private

Definition at line 314 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_e1x5_[3]

private

Definition at line 226 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_e2x5_[3]

private

Definition at line 232 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_ecalSum_[3]

private

Definition at line 278 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_etOutsideMustache_[3]

private

Definition at line 306 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_h1OverE_[3]

private

Definition at line 293 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_h2OverE_[3]

private

Definition at line 294 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_hcalSum_[3]

private

Definition at line 284 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_hOverE_[3]

private

Definition at line 290 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nCluOutsideMustache_[3]

private

Definition at line 305 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_newhOverE_[3]

private

Definition at line 296 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nHadIso_[3]

private

Definition at line 302 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nPho_[3]

private

Definition at line 213 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nRecoVtx_

private

Definition at line 203 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nTrackIsolHollow_[3]

private

Definition at line 260 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_nTrackIsolSolid_[3]

private

Definition at line 254 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_pfMva_[3]

private

Definition at line 307 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoE_[3]

private

Definition at line 208 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoEt_[3]

private

Definition at line 211 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoEta_[3]

private

Definition at line 215 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoIso_[3]

private

Definition at line 303 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoPhi_[3]

private

Definition at line 216 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoSigmaEoverE_[3]

private

Definition at line 209 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_phoSigmaIetaIeta_[3]

private

Definition at line 250 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_r1x5_[3]

private

Definition at line 238 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_r2x5_[3]

private

Definition at line 244 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_r9_[3]

private

Definition at line 220 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_scEta_[3]

private

Definition at line 217 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_scPhi_[3]

private

Definition at line 218 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_ChHad_Cleaned_[3]

private

Definition at line 315 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_ChHad_unCleaned_[3]

private

Definition at line 318 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_NeuHad_Cleaned_[3]

private

Definition at line 316 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_NeuHad_unCleaned_[3]

private

Definition at line 319 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_Pho_Cleaned_[3]

private

Definition at line 317 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_SumPtOverPhoPt_Pho_unCleaned_[3]

private

Definition at line 320 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_trackPtSumHollow_[3]

private

Definition at line 272 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::h_trackPtSumSolid_[3]

private

Definition at line 266 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::hOverEBin_

private

Definition at line 182 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::hOverEMax_

private

Definition at line 181 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::hOverEMin_

private

Definition at line 180 of file ZToMuMuGammaAnalyzer.h.

bool ZToMuMuGammaAnalyzer::makeProfiles_

private

Definition at line 123 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::maxMumuGammaInvMass_

private

Definition at line 153 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::maxMumuInvMass_

private

Definition at line 140 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::minMumuGammaInvMass_

private

Definition at line 152 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::minMumuInvMass_

private

Definition at line 139 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::minPixStripHits_

private

Definition at line 130 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<std::vector<reco::Muon> > ZToMuMuGammaAnalyzer::muon_token_

private

Definition at line 110 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::muonMatches_

private

Definition at line 133 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::muonMaxChi2_

private

Definition at line 131 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::muonMaxDxy_

private

Definition at line 132 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::muonMinPt_

private

Definition at line 129 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::muonTightEta_

private

Definition at line 137 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::muonTrackIso_

private

Definition at line 136 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::nearMuonDr_

private

Definition at line 147 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::nearMuonHcalIso_

private

Definition at line 148 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::nEvt_

private

Definition at line 126 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::numberBin_

private

Definition at line 186 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::numberMax_

private

Definition at line 185 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::numberMin_

private

Definition at line 184 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<reco::VertexCollection> ZToMuMuGammaAnalyzer::offline_pvToken_

private

Definition at line 118 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_e1x5VsEt_[3]

private

Definition at line 230 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_e1x5VsEta_[3]

private

Definition at line 228 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_e2x5VsEt_[3]

private

Definition at line 236 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_e2x5VsEta_[3]

private

Definition at line 234 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_ecalSumVsEt_[3]

private

Definition at line 280 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_ecalSumVsEta_[3]

private

Definition at line 282 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_hcalSumVsEt_[3]

private

Definition at line 286 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_hcalSumVsEta_[3]

private

Definition at line 288 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_hOverEVsEt_[3]

private

Definition at line 291 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_hOverEVsEta_[3]

private

Definition at line 292 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_newhOverEVsEt_[3]

private

Definition at line 298 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_newhOverEVsEta_[3]

private

Definition at line 297 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_nTrackIsolHollowVsEt_[3]

private

Definition at line 262 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_nTrackIsolHollowVsEta_[3]

private

Definition at line 264 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_nTrackIsolSolidVsEt_[3]

private

Definition at line 256 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_nTrackIsolSolidVsEta_[3]

private

Definition at line 258 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_phoSigmaEoverEVsNVtx_[3]

private

Definition at line 210 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r1x5VsEt_[3]

private

Definition at line 242 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r1x5VsEta_[3]

private

Definition at line 240 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r2x5VsEt_[3]

private

Definition at line 248 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r2x5VsEta_[3]

private

Definition at line 246 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r9VsEt_[3]

private

Definition at line 222 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_r9VsEta_[3]

private

Definition at line 224 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_sigmaIetaIetaVsEta_[3]

private

Definition at line 252 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_trackPtSumHollowVsEt_[3]

private

Definition at line 274 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_trackPtSumHollowVsEta_[3]

private

Definition at line 276 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_trackPtSumSolidVsEt_[3]

private

Definition at line 268 of file ZToMuMuGammaAnalyzer.h.

MonitorElement* ZToMuMuGammaAnalyzer::p_trackPtSumSolidVsEta_[3]

private

Definition at line 270 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<reco::PFCandidateCollection> ZToMuMuGammaAnalyzer::pfCandidates_

private

Definition at line 117 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::phiBin_

private

Definition at line 174 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::phiMax_

private

Definition at line 173 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::phiMin_

private

Definition at line 172 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<std::vector<reco::Photon> > ZToMuMuGammaAnalyzer::photon_token_

private

Definition at line 109 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<edm::ValueMap<bool> > ZToMuMuGammaAnalyzer::PhotonIDLoose_token_

private

Definition at line 111 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<edm::ValueMap<bool> > ZToMuMuGammaAnalyzer::PhotonIDTight_token_

private

Definition at line 112 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<edm::ValueMap<std::vector<reco::PFCandidateRef> > > ZToMuMuGammaAnalyzer::photonIsoValmap_token_

private

Definition at line 119 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::photonMaxEta_

private

Definition at line 143 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::photonMinEt_

private

Definition at line 142 of file ZToMuMuGammaAnalyzer.h.

float ZToMuMuGammaAnalyzer::photonTrackIso_

private

Definition at line 144 of file ZToMuMuGammaAnalyzer.h.

unsigned int ZToMuMuGammaAnalyzer::prescaleFactor_

private

Definition at line 124 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::r9Bin_

private

Definition at line 178 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::r9Max_

private

Definition at line 177 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::r9Min_

private

Definition at line 176 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::reducedEtaBin_

private

Definition at line 193 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::reducedEtBin_

private

Definition at line 192 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::reducedR9Bin_

private

Definition at line 195 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::reducedSumBin_

private

Definition at line 194 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::sigmaIetaBin_

private

Definition at line 190 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::sigmaIetaMax_

private

Definition at line 189 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::sigmaIetaMin_

private

Definition at line 188 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::sumBin_

private

Definition at line 166 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::sumMax_

private

Definition at line 165 of file ZToMuMuGammaAnalyzer.h.

double ZToMuMuGammaAnalyzer::sumMin_

private

Definition at line 164 of file ZToMuMuGammaAnalyzer.h.

edm::EDGetTokenT<trigger::TriggerEvent> ZToMuMuGammaAnalyzer::triggerEvent_token_

private

Definition at line 115 of file ZToMuMuGammaAnalyzer.h.

bool ZToMuMuGammaAnalyzer::use2DHistos_

private

Definition at line 122 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::validMuonHits_

private

Definition at line 135 of file ZToMuMuGammaAnalyzer.h.

int ZToMuMuGammaAnalyzer::validPixHits_

private

Definition at line 134 of file ZToMuMuGammaAnalyzer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation