ZToMuMuGammaAnalyzer Class Reference

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
	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 &	itemsToGetFromEvent () 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(), HLT_FULL_cff::InputTag, 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.

{
}

Member Function Documentation

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

overridevirtual

uncleaned

Implements edm::stream::EDAnalyzerBase.

Definition at line 421 of file ZToMuMuGammaAnalyzer.cc.

References begin, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::PFCandidate::e, reco::MuonIsolation::emEt, HcalObjRepresent::Fill(), reco::PFCandidate::gamma, edm::Event::getByToken(), reco::PFCandidate::h, reco::PFCandidate::h0, reco::MuonIsolation::hadEt, i, edm::EventBase::id(), reco::Muon::isolationR03(), diffTwoXMLs::label, HLT_FULL_cff::muonSelection, reco::LeafCandidate::pt(), 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 edm::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.