|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("PFISO")
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
)
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring('/store/relval/CMSSW_7_1_0_pre1/RelValProdTTbar/GEN-SIM-RECO/START70_V5-v1/00000/14842A6B-2086-E311-B5CB-02163E00E8DA.root'))
# Tae Jeong, could you please remove these from the PF2PAT sequence?
# they are not used, and are creating problems in such kinds of manipulations:
# process.pfElectronIsolationSequence.remove( process.pfElectronIsoDepositsSequence )
# process.pfElectronIsolationSequence.remove( process.pfElectronIsolationFromDepositsSequence )
# process.load("CommonTools.ParticleFlow.PFBRECO_cff")
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso, setupPFMuonIso, setupPFPhotonIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.muIsoSequence = setupPFMuonIso(process, 'muons')
process.phoIsoSequence = setupPFPhotonIso(process, 'photons')
process.TFileService = cms.Service("TFileService", fileName = cms.string("histo.root") )
process.pfIsoReader = cms.EDAnalyzer("PFIsoReaderDemo",
Electrons = cms.InputTag('gsfElectrons'),
Photons = cms.InputTag('photons'),
PFCandidateMap = cms.InputTag('particleFlow:electrons'),
PrintElectrons = cms.bool(True),
PrintPhotons = cms.bool(True),
IsoDepElectron = cms.VInputTag(cms.InputTag('elPFIsoDepositChargedPFIso'),
cms.InputTag('elPFIsoDepositGammaPFIso'),
cms.InputTag('elPFIsoDepositNeutralPFIso')),
IsoValElectronPF = cms.VInputTag(cms.InputTag('elPFIsoValueCharged03PFIdPFIso'),
cms.InputTag('elPFIsoValueGamma03PFIdPFIso'),
cms.InputTag('elPFIsoValueNeutral03PFIdPFIso')),
IsoDepPhoton = cms.VInputTag(cms.InputTag('phPFIsoDepositChargedPFIso'),
cms.InputTag('phPFIsoDepositGammaPFIso'),
cms.InputTag('phPFIsoDepositNeutralPFIso')),
IsoValPhoton = cms.VInputTag(cms.InputTag('phPFIsoValueCharged03PFIdPFIso'),
cms.InputTag('phPFIsoValueGamma03PFIdPFIso'),
cms.InputTag('phPFIsoValueNeutral03PFIdPFIso'))
)
process.p = cms.Path(
# process.pfNoPileUpSequence +
process.pfParticleSelectionSequence +
process.eleIsoSequence +
process.muIsoSequence+
process.phoIsoSequence+
process.pfIsoReader
)
# output ------------------------------------------------------------
process.out = cms.OutputModule("PoolOutputModule",
outputCommands = cms.untracked.vstring('keep *'),
fileName = cms.untracked.string('pfIsolation.root')
)
process.outpath = cms.EndPath(
process.out
)
# other stuff
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 10
|