Skims/python/IsoPhotonEBSkim_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002 from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
0003 IsoPhotonEBHLTFilter = hltHighLevel.clone(
0004      throw = False,
0005      HLTPaths = ["HLT_Photon110EB_TightID_TightIso_v*"]
0006 )
0007
0008 # run on MIONAOD
0009 RUN_ON_MINIAOD = False
0010
0011
0012 # cuts
0013 PHOTON_CUT=("pt > 110 && abs(eta)<1.4442")
0014
0015 # single lepton selectors
0016 if RUN_ON_MINIAOD:
0017     goodPhotons = cms.EDFilter("PATElectronRefSelector",
0018                                     src = cms.InputTag("slimmedPhotons"),
0019                                     cut = cms.string(PHOTON_CUT)
0020                                     )
0021 else:
0022     goodPhotons = cms.EDFilter("PhotonRefSelector",
0023                                     src = cms.InputTag("gedPhotons"),
0024                                     cut = cms.string(PHOTON_CUT)
0025                                     )
0026
0027 photonIDWP = cms.PSet( #first for barrel, second for endcap.
0028     full5x5_sigmaIEtaIEtaCut       = cms.vdouble(0.011 ,-1. )  , # full5x5_sigmaIEtaIEtaCut
0029     hOverECut                      = cms.vdouble(0.1  ,-1. )  , # hOverECut
0030     relCombIsolationWithEACut      = cms.vdouble(0.05  ,-1. )   # relCombIsolationWithEALowPtCut
0031 )
0032
0033
0034 identifiedPhotons = cms.EDFilter("IsoPhotonEBSelectorAndSkim",
0035                                    src    = cms.InputTag("goodPhotons"),
0036                                    phID = photonIDWP,
0037                                    absEtaMin=cms.vdouble( 0.0000, 1.0000, 1.4790, 2.0000, 2.2000, 2.3000, 2.4000),
0038                                    absEtaMax=cms.vdouble( 1.0000,  1.4790, 2.0000,  2.2000, 2.3000, 2.4000, 5.0000),
0039                                    effectiveAreaValues=cms.vdouble( 0.1703, 0.1715, 0.1213, 0.1230, 0.1635, 0.1937, 0.2393),
0040                                    rho = cms.InputTag("fixedGridRhoFastjetCentralCalo")
0041                          )
0042
0043 identifiedPhotonsCountFilter = cms.EDFilter("CandViewCountFilter",
0044                                     src = cms.InputTag("identifiedPhotons"),
0045                                     minNumber = cms.uint32(1)
0046                                     )
0047
0048
0049
0050 #sequences
0051 isoPhotonEBSequence = cms.Sequence(IsoPhotonEBHLTFilter*goodPhotons*identifiedPhotons*identifiedPhotonsCountFilter )