EgammaIsolationAlgos/test/egammaIsoTest_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("eleIso")
0004 process.load("Configuration.StandardSequences.GeometryDB_cff")
0005
0006 process.load("Configuration.EventContent.EventContent_cff")
0007 process.source = cms.Source("PoolSource",
0008     fileNames = cms.untracked.vstring(
0009         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/16EE7689-5EF4-DD11-97EF-001D09F2514F.root',
0010         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/66B68FBF-56F4-DD11-86DD-001617C3B710.root',
0011         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/8A6092F3-51F4-DD11-A213-001617DBD5AC.root',
0012         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/AE032B6E-54F4-DD11-ACA5-000423D9A212.root',
0013         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/C6DDF78F-55F4-DD11-8FF4-001617C3B66C.root',
0014         '/store/relval/CMSSW_3_1_0_pre1/RelValTTbar/GEN-SIM-RECO/IDEAL_30X_v1/0001/E83F680E-56F4-DD11-9976-001617C3B6CE.root'
0015
0016     )
0017 )
0018
0019 process.out = cms.OutputModule("PoolOutputModule",
0020     process.FEVTSIMEventContent,
0021     fileName = cms.untracked.string('file:eleIso.root')
0022 )
0023
0024 process.out.outputCommands.append('drop *_*_*_*')
0025 process.out.outputCommands.append('keep *_gsfElectrons_*_*')
0026 process.out.outputCommands.append('keep *_photons_*_*')
0027 process.out.outputCommands.append('keep *_*_*_eleIso')
0028
0029 process.maxEvents = cms.untracked.PSet(
0030     input = cms.untracked.int32(200)
0031 )
0032
0033 process.load("RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequence_cff")
0034 process.load("RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequencePAT_cff")
0035
0036 #Couple extra modules to test that the new abs(energy) was working:
0037 #Get IsoDeposits with no min energy cut, all RecHits stored
0038 process.eleIsoDepositEcalFromHitsNoCut = process.eleIsoDepositEcalFromHits.clone()
0039 process.eleIsoDepositEcalFromHitsNoCut.ExtractorPSet.energyMin = cms.double(0.0)
0040
0041 #Default Vetos on new IsoDeposit
0042 process.eleIsoFromDepsEcalFromHitsAbs = process.eleIsoFromDepsEcalFromHits.clone()
0043 process.eleIsoFromDepsEcalFromHitsAbs.deposits[0].src = 'eleIsoDepositEcalFromHitsNoCut'
0044
0045 #Vetos without the abs(energy) cut, using old single sided cut
0046 process.eleIsoFromDepsEcalFromHitsNoAbs = process.eleIsoFromDepsEcalFromHits.clone()
0047 process.eleIsoFromDepsEcalFromHitsNoAbs.deposits[0].src = 'eleIsoDepositEcalFromHitsNoCut'
0048 process.eleIsoFromDepsEcalFromHitsNoAbs.deposits[0].vetos = cms.vstring(
0049     'EcalBarrel:0.045',
0050     'EcalBarrel:RectangularEtaPhiVeto(-0.02,0.02,-0.5,0.5)',
0051     'EcalBarrel:ThresholdFromTransverse(0.08)',
0052     'EcalEndcaps:ThresholdFromTransverse(0.3)',
0053     'EcalEndcaps:0.070',
0054     'EcalEndcaps:RectangularEtaPhiVeto(-0.02,0.02,-0.5,0.5)'
0055 )
0056
0057 #Rediculous inner cone to make sure the cuts were working correctly:
0058 process.eleIsoFromDepsEcalFromHitsNoAbs2 = process.eleIsoFromDepsEcalFromHits.clone()
0059 process.eleIsoFromDepsEcalFromHitsNoAbs2.deposits[0].src = 'eleIsoDepositEcalFromHitsNoCut'
0060 process.eleIsoFromDepsEcalFromHitsNoAbs2.deposits[0].vetos = cms.vstring(
0061     'EcalBarrel:0.3',
0062     'EcalBarrel:RectangularEtaPhiVeto(-0.02,0.02,-0.5,0.5)',
0063     'EcalBarrel:ThresholdFromTransverse(0.08)',
0064     'EcalEndcaps:ThresholdFromTransverse(0.3)',
0065     'EcalEndcaps:0.3',
0066     'EcalEndcaps:RectangularEtaPhiVeto(-0.02,0.02,-0.5,0.5)'
0067 )
0068
0069
0070 process.p1 = cms.Path(
0071     process.egammaIsolationSequence +
0072     process.egammaIsolationSequencePAT +
0073     (
0074         process.eleIsoDepositEcalFromHitsNoCut *
0075         process.eleIsoFromDepsEcalFromHitsAbs *
0076         process.eleIsoFromDepsEcalFromHitsNoAbs *
0077         process.eleIsoFromDepsEcalFromHitsNoAbs2
0078     )
0079 )
0080
0081 process.outpath = cms.EndPath(process.out)