PhotonIdentification/test/testPhotonIDProducer_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PhotonIDProc")
0004
0005 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0006 process.load("RecoEgamma.PhotonIdentification.photonId_cff")
0007 process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
0008 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
0009 #process.load("MagneticField.Engine.volumeBasedMagneticField_cfi")
0010 process.load("Configuration.StandardSequences.MagneticField_cff")
0011 process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
0012 process.load("RecoEgamma.EgammaPhotonProducers.photonSequence_cff")
0013
0014 process.source = cms.Source("PoolSource",
0015     fileNames = cms.untracked.vstring('/store/relval/2008/6/22/RelVal-RelValSingleGammaPt35-1213986417-IDEAL_V2-2nd/0004/443BCAED-CB40-DD11-AB37-000423D6B48C.root')
0016     #fileNames = cms.untracked.vstring('dcap://cmsdcap.hep.wisc.edu:22125/pnfs/hep.wisc.edu/store/user/mbanderson/PhotonJet20-200/PhotonJet20-200-0000.root')
0017 )
0018
0019 process.maxEvents = cms.untracked.PSet(
0020     input = cms.untracked.int32(-1)
0021 )
0022
0023 process.Out = cms.OutputModule("PoolOutputModule",
0024     outputCommands = cms.untracked.vstring('drop *',
0025         'keep edmHepMCProduct_*_*_*',
0026         'keep recoBasicClusters_*_*_*',
0027         'keep recoSuperClusters_*_*_*',
0028         'keep *_PhotonIDProd_*_*',
0029         'keep *_PhotonIDProd_*_*',
0030         'keep recoPhotons_*_*_*'),
0031     fileName = cms.untracked.string('Photest.root')
0032 )
0033
0034 process.photonIDAna = cms.EDAnalyzer("PhotonIDSimpleAnalyzer",
0035     outputFile  = cms.string('PhoIDHists.root'),
0036
0037     # Variables that must be passed before a photon candidate (a SuperCluster)
0038     #  gets placed into the histograms.  Basic, simple cuts.
0039     # Minimum Et
0040     minPhotonEt     = cms.double(10.0),
0041     # Minimum and max abs(eta)
0042     minPhotonAbsEta = cms.double(0.0),
0043     maxPhotonAbsEta = cms.double(3.0),
0044     # Minimum R9 = E(3x3) / E(SuperCluster)
0045     minPhotonR9     = cms.double(0.3),
0046     # Maximum HCAL / ECAL Energy
0047     maxPhotonHoverE = cms.double(0.2),
0048
0049     # Optionally produce a TTree of photons (set to False or True).
0050     # This slows down the analyzer, and if running
0051     # over 100,000+ events, this can create a large ROOT file
0052     createPhotonTTree  = cms.bool(True)
0053 )
0054
0055 #process.p = cms.Path(process.photonSequence*process.photonIDSequence*process.photonIDAna)
0056 process.p = cms.Path(process.photonIDSequence*process.photonIDAna)
0057 process.e = cms.EndPath(process.Out)
0058
0059 process.PoolSource.fileNames = [
0060     '/store/relval/CMSSW_3_2_5/RelValZEE/GEN-SIM-RECO/MC_31X_V5-v1/0011/84E3AE4D-738E-DE11-A5E4-003048D374F2.root',
0061 '/store/relval/CMSSW_3_2_5/RelValZEE/GEN-SIM-RECO/MC_31X_V5-v1/0011/7886FE89-738E-DE11-B893-001D09F241F0.root',
0062 '/store/relval/CMSSW_3_2_5/RelValZEE/GEN-SIM-RECO/MC_31X_V5-v1/0011/60D71215-828E-DE11-8BFE-000423D98804.root',
0063 '/store/relval/CMSSW_3_2_5/RelValZEE/GEN-SIM-RECO/MC_31X_V5-v1/0011/4AA439F6-728E-DE11-AF2D-000423D98EA8.root'
0064 ]