-
Notifications
You must be signed in to change notification settings - Fork 0
/
ttjets.py
152 lines (131 loc) · 6.47 KB
/
ttjets.py
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAny_cfi")
process.load("Configuration.StandardSequences.Geometry_cff")
process.load("Geometry.CommonDetUnit.globalTrackingGeometry_cfi")
process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = 'IDEAL_V12::All'
# for match muon in a jet by using trackAssociator
#from TrackingTools.TrackAssociator.default_cfi import *
process.source = cms.Source("PoolSource",
debugFlag = cms.untracked.bool(False),
debugVebosity = cms.untracked.uint32(10),
skipEvents = cms.untracked.uint32(0),
fileNames = cms.untracked.vstring(
'file:/home/cms/sckao/Top/CMSSW_2_2_13/src/TopPhysics/TtAnalysis/test/TMass/patTest_SCJ_TC.root'
#'file:/home/cms/sckao/Top/CMSSW_2_2_13/src/TopPhysics/TtAnalysis/test/TMass/patTest_ICJ_TC.root'
#'file:/home/cms/sckao/Top/CMSSW_2_2_13/src/TopPhysics/TtAnalysis/test/TMass/patTest_oldMET_skim.root'
#'file:/home/cms/sckao/Data/FastPat2213/WJets_PAT_Test.root'
)
)
# replace the source files from a file list
#import TopPhysics.TtAnalysis.ttjetslist_skim2 as fileList
#process.source.fileNames = fileList.fileNames
process.maxEvents = cms.untracked.PSet(
# 100/pb => 9191
input = cms.untracked.int32(9191)
#input = cms.untracked.int32(24000)
)
process.MessageLogger = cms.Service("MessageLogger")
#from PhysicsTools.PatAlgos.tools.trigTools import switchOffTriggerMatchingOld
#switchOffTriggerMatchingOld( process )
process.ttAna = cms.EDAnalyzer("TtAnalysis",
# for matching Muon in a Jet
#TrackAssociatorParameterBlock,
#TrackAssociatorParameters,
debug = cms.untracked.bool(False),
btag = cms.untracked.bool(True),
bTagCut = cms.untracked.double(2),
#bTagAlgo = cms.untracked.string('softMuonBJetTags'),
bTagAlgo = cms.untracked.string("trackCountingHighEffBJetTags"),
needTree = cms.untracked.bool(False),
trigOn = cms.untracked.bool(False),
rootFileName = cms.untracked.string('pat2_SCJ_TC_btag_tt.root'),
genParticles = cms.InputTag("genParticles"),
#genJetSource = cms.InputTag("iterativeCone5GenJets"),
genJetSource = cms.InputTag("sisCone5GenJets"),
electronSource = cms.InputTag("cleanLayer1Electrons"),
photonSource = cms.InputTag("cleanLayer1Photons"),
jetSource = cms.InputTag("cleanLayer1Jets"),
jptSource = cms.InputTag("ZSPJetCorJetIcone5"),
tcMetSource = cms.InputTag("tcMet"),
metSource = cms.InputTag("layer1METsTC"),
genmetSource = cms.InputTag("genMet"),
muonSource = cms.InputTag("cleanLayer1Muons"),
caloSource = cms.InputTag("towerMaker"),
triggerSource = cms.InputTag("TriggerResults","","HLT"),
#triggerSource = cms.InputTag("TriggerResults","","PAT"),
#recoMuons = cms.untracked.string('paramMuons'),
recoMuons = cms.untracked.string('muons'),
recoAlgo = cms.untracked.string('zero'),
)
process.ttAna1 = cms.EDAnalyzer("TtAnalysis",
# for matching Muon in a Jet
#TrackAssociatorParameterBlock,
#TrackAssociatorParameters,
debug = cms.untracked.bool(False),
btag = cms.untracked.bool(False),
bTagCut = cms.untracked.double(2),
bTagAlgo = cms.untracked.string("trackCountingHighEffBJetTags"),
needTree = cms.untracked.bool(False),
trigOn = cms.untracked.bool(False),
rootFileName = cms.untracked.string('pat2_SCJ_TC_tt.root'),
genParticles = cms.InputTag("genParticles"),
genJetSource = cms.InputTag("sisCone5GenJets"),
electronSource = cms.InputTag("cleanLayer1Electrons"),
photonSource = cms.InputTag("cleanLayer1Photons"),
jetSource = cms.InputTag("cleanLayer1Jets"),
jptSource = cms.InputTag("ZSPJetCorJetIcone5"),
tcMetSource = cms.InputTag("tcMet"),
metSource = cms.InputTag("layer1METsTC"),
genmetSource = cms.InputTag("genMet"),
muonSource = cms.InputTag("cleanLayer1Muons"),
caloSource = cms.InputTag("towerMaker"),
triggerSource = cms.InputTag("TriggerResults","","HLT"),
recoMuons = cms.untracked.string('muons'),
recoAlgo = cms.untracked.string('zero'),
)
process.jetAna = cms.EDAnalyzer("JetAnalysis",
debug = cms.untracked.bool(False),
bTagCut = cms.untracked.double(2),
#bTagAlgo = cms.untracked.string('softMuonBJetTags'),
bTagAlgo = cms.untracked.string("trackCountingHighEffBJetTags"),
#bTagAlgo = cms.untracked.string('jetProbabilityBJetTags'),
rootFileName = cms.untracked.string('pat2_SCJ_TC_btag_jm.root'),
genParticles = cms.InputTag("genParticles"),
#genJetSource = cms.InputTag("iterativeCone5GenJets"),
genJetSource = cms.InputTag("sisCone5GenJets"),
electronSource = cms.InputTag("cleanLayer1Electrons"),
jetSource = cms.InputTag("cleanLayer1Jets"),
metSource = cms.InputTag("layer1METsTC"),
genmetSource = cms.InputTag("genMet"),
jptSource = cms.InputTag("ZSPJetCorJetIcone5"),
tcMetSource = cms.InputTag("tcMet"),
muonSource = cms.InputTag("cleanLayer1Muons"),
caloSource = cms.InputTag("towerMaker"),
#recoMuons = cms.untracked.string('paramMuons')
recoMuons = cms.untracked.string('muons')
)
process.muAna = cms.EDAnalyzer("MuonAnalysis",
debug = cms.untracked.bool(False),
rootFileName = cms.untracked.string('ttj_IsoMuAnalysis.root'),
genParticles = cms.InputTag("genParticles"),
genJetSource = cms.InputTag("iterativeCone5GenJets"),
electronSource = cms.InputTag("cleanLayer1Electrons"),
jetSource = cms.InputTag("cleanLayer1Jets"),
muonSource = cms.InputTag("cleanLayer1Muons"),
metSource = cms.InputTag("layer1METs"),
caloSource = cms.InputTag("towerMaker"),
recoMuons = cms.untracked.string('paramMuons')
)
#process.p = cms.Path( process.ttAna + process.jetAna + process.muAna )
#process.p = cms.Path( process.ttAna + process.ttAna1 + process.ttAna2 + process.ttAna3 )
process.p = cms.Path( process.ttAna + process.ttAna1 + process.jetAna )
#process.p = cms.Path( process.ttAna )
#process.ttAna.TrackAssociatorParameters.useEcal = False
#process.ttAna.TrackAssociatorParameters.useHcal = False
#process.ttAna.TrackAssociatorParameters.useCalo = True
#process.ttAna.TrackAssociatorParameters.useHO = False
#process.ttAna.TrackAssociatorParameters.useMuon = False