test/utilities/skimPhotonPhoton.py

0001 import FWCore.ParameterSet.Config as cms
0002 process = cms.Process("SKIM")
0003
0004 #------------------------------------------
0005 # Load standard sequences.
0006 #------------------------------------------
0007 process.load('FWCore.MessageService.MessageLogger_cfi')
0008
0009
0010 #  Photons!!! ################
0011 process.goodPhotons = cms.EDFilter("PhotonSelector",
0012                                    src = cms.InputTag("photons"),
0013                                    cut = cms.string("hadronicOverEm<0.15"
0014                                                     " && (superCluster.rawEnergy*sin(superCluster.position.theta)>15.0)")
0015                                    )
0016
0017
0018 process.ZToEE = cms.EDProducer("CandViewShallowCloneCombiner",
0019     cut = cms.string('40 < mass < 7000'),
0020     decay = cms.string('goodPhotons goodPhotons'),
0021     checkCharge = cms.bool(False),
0022
0023 )
0024
0025 process.ZToEEFilter = cms.EDFilter("CandViewCountFilter",
0026                                     src = cms.InputTag("ZToEE"),
0027                                     minNumber = cms.uint32(1)
0028                                     )
0029
0030 process.ZeeSequence = cms.Sequence(process.goodPhotons *
0031                                     process.ZToEE)
0032
0033
0034
0035
0036 process.mypath = cms.Path(
0037     process.ZeeSequence*
0038     process.ZToEEFilter)
0039
0040
0041
0042
0043 #############  output module if just want to skim by HLT path ##############
0044 process.out = cms.OutputModule("PoolOutputModule",
0045     outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0046     SelectEvents = cms.untracked.PSet(
0047        SelectEvents = cms.vstring('mypath')
0048        ),
0049     fileName = cms.untracked.string('mySkim.root')
0050 )
0051 process.p = cms.EndPath(process.out)
0052
0053
0054 #############   Format MessageLogger #################
0055 process.MessageLogger.cerr.FwkReport.reportEvery = 100
0056
0057
0058
0059 #############   Set the number of events #############
0060 process.maxEvents = cms.untracked.PSet(
0061     input = cms.untracked.int32(-1)
0062 )
0063 #############   Define the source file ###############
0064 process.source = cms.Source("PoolSource",
0065     fileNames = cms.untracked.vstring(
0066         '/store/data/Run2010A/EG/RECO/v4/000/139/195/FA470877-BE85-DF11-840C-003048F117B4.root',
0067         '/store/data/Run2010A/EG/RECO/v4/000/139/195/ACECA275-D885-DF11-931C-003048F11DE2.root',
0068         '/store/data/Run2010A/EG/RECO/v4/000/139/195/A6DB3770-D885-DF11-9A0B-003048F024DE.root',
0069
0070      )
0071 )
0072 process.source.inputCommands = cms.untracked.vstring("keep *","drop *_MEtoEDMConverter_*_*")
0073
0074
0075