2 |
|
|
3 |
|
ivars = VarParsing.VarParsing('standard') |
4 |
|
|
5 |
< |
#ivars.files = 'file:PyquenDijet80to120_d20100325_runs1to500.root' |
6 |
< |
#ivars.output = 'jets_pat_80to120.root' |
7 |
< |
#ivars.files = 'file:PyquenDijet120to170_d20100325_runs1to500.root' |
8 |
< |
#ivars.output = 'jets_pat_120to170.root' |
9 |
< |
|
10 |
< |
ivars.files = 'file:PyquenDijet50to80_d20100325_runs1to500.root' |
11 |
< |
ivars.output = 'jets_pat_50to80.root' |
12 |
< |
ivars.maxEvents = -1 |
5 |
> |
ivars.files = 'rfio:/castor/cern.ch/cms/store/relval/CMSSW_3_7_0_pre2/RelValPyquen_DiJet_pt80to120_2760GeV/GEN-SIM-RECO/MC_37Y_V1-v1/0018/1420F252-8453-DF11-9903-002618943882.root' |
6 |
> |
ivars.output = 'jets_pat_80to120.root' |
7 |
> |
ivars.maxEvents = 3 |
8 |
|
|
9 |
|
ivars.parseArguments() |
10 |
|
|
11 |
|
import FWCore.ParameterSet.Config as cms |
12 |
|
|
13 |
< |
process = cms.Process('MATCH') |
13 |
> |
process = cms.Process('HIJETS') |
14 |
|
|
15 |
|
process.source = cms.Source("PoolSource", |
16 |
|
fileNames = cms.untracked.vstring(ivars.files) |
17 |
|
) |
18 |
|
|
19 |
|
process.maxEvents = cms.untracked.PSet( |
20 |
< |
input = cms.untracked.int32(ivars.maxEvents) |
21 |
< |
) |
20 |
> |
input = cms.untracked.int32(ivars.maxEvents) |
21 |
> |
) |
22 |
|
|
23 |
+ |
process.load('Configuration.StandardSequences.Geometry_cff') |
24 |
+ |
process.load('RecoJets.Configuration.GenJetParticles_cff') |
25 |
+ |
process.load('RecoHI.HiJetAlgos.HiGenJets_cff') |
26 |
+ |
process.load('RecoHI.HiJetAlgos.HiRecoJets_cff') |
27 |
|
process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff') |
28 |
|
from PhysicsTools.PatAlgos.tools.heavyIonTools import * |
29 |
|
configureHeavyIons(process) |
30 |
|
|
31 |
< |
from PhysicsTools.PatAlgos.tools.jetTools import * |
32 |
< |
switchJECSet( process, "Summer09_7TeV_ReReco332") |
31 |
> |
#from PhysicsTools.PatAlgos.tools.jetTools import * |
32 |
> |
#switchJECSet( process, "Summer09_7TeV_ReReco332") |
33 |
|
|
34 |
< |
process.ak5corr = process.patJetCorrFactors.clone() |
35 |
< |
process.ak5corr.jetSource = cms.InputTag("akFastPu5CaloJets") |
36 |
< |
process.ak5corr.corrLevels.L2Relative = "L2Relative_AK5Calo" |
37 |
< |
process.ak5corr.corrLevels.L3Absolute = "L3Absolute_AK5Calo" |
38 |
< |
|
39 |
< |
process.ak7corr = process.patJetCorrFactors.clone() |
40 |
< |
process.ak7corr.jetSource = cms.InputTag("akFastPu7CaloJets") |
41 |
< |
process.ak7corr.corrLevels.L2Relative = cms.string('L2Relative_AK7Calo') |
42 |
< |
process.ak7corr.corrLevels.L3Absolute = cms.string('L3Absolute_AK7Calo') |
43 |
< |
|
44 |
< |
process.akPu5corr = process.ak5corr.clone() |
45 |
< |
process.akPu5corr.jetSource = cms.InputTag("akPu5CaloJets") |
46 |
< |
process.akPu7corr = process.ak7corr.clone() |
47 |
< |
process.akPu7corr.jetSource = cms.InputTag("akPu7CaloJets") |
48 |
< |
|
49 |
< |
process.kt4corr = process.patJetCorrFactors.clone() |
50 |
< |
process.kt4corr.jetSource = cms.InputTag("ktFastPu4CaloJets") |
51 |
< |
process.kt4corr.corrLevels.L2Relative = cms.string('L2Relative_KT4Calo') |
52 |
< |
process.kt4corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT4Calo') |
53 |
< |
|
54 |
< |
process.kt6corr = process.patJetCorrFactors.clone() |
55 |
< |
process.kt6corr.jetSource = cms.InputTag("ktFastPu6CaloJets") |
56 |
< |
process.kt6corr.corrLevels.L2Relative = cms.string('L2Relative_KT6Calo') |
57 |
< |
process.kt6corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT6Calo') |
58 |
< |
|
59 |
< |
process.ktPu4corr = process.patJetCorrFactors.clone() |
60 |
< |
process.ktPu4corr.jetSource = cms.InputTag("ktPu4CaloJets") |
61 |
< |
process.ktPu4corr.corrLevels.L2Relative = cms.string('L2Relative_KT4Calo') |
62 |
< |
process.ktPu4corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT4Calo') |
63 |
< |
|
64 |
< |
process.ktPu6corr = process.patJetCorrFactors.clone() |
65 |
< |
process.ktPu6corr.jetSource = cms.InputTag("ktPu6CaloJets") |
66 |
< |
process.ktPu6corr.corrLevels.L2Relative = cms.string('L2Relative_KT6Calo') |
67 |
< |
process.ktPu6corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT6Calo') |
68 |
< |
|
69 |
< |
process.icPu5corr = process.patJetCorrFactors.clone() |
70 |
< |
process.icPu5corr.jetSource = cms.InputTag("iterativeConePu5CaloJets") |
71 |
< |
process.icPu5corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo') |
72 |
< |
process.icPu5corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo') |
73 |
< |
|
74 |
< |
process.icPu7corr = process.patJetCorrFactors.clone() |
75 |
< |
process.icPu7corr.jetSource = cms.InputTag("iterativeConePu7CaloJets") |
76 |
< |
process.icPu7corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo') # No IC7 |
77 |
< |
process.icPu7corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo') # No IC7 |
78 |
< |
|
79 |
< |
process.ak5clean = process.heavyIonCleanedGenJets.clone() |
80 |
< |
process.ak5clean.src = cms.untracked.string('ak5HiGenJets') |
81 |
< |
process.ak7clean = process.heavyIonCleanedGenJets.clone() |
82 |
< |
process.ak7clean.src = cms.untracked.string('ak7HiGenJets') |
83 |
< |
process.kt4clean = process.heavyIonCleanedGenJets.clone() |
84 |
< |
process.kt4clean.src = cms.untracked.string('kt4HiGenJets') |
85 |
< |
process.kt6clean = process.heavyIonCleanedGenJets.clone() |
86 |
< |
process.kt6clean.src = cms.untracked.string('kt6HiGenJets') |
88 |
< |
process.ic5clean = process.heavyIonCleanedGenJets.clone() |
89 |
< |
process.ic5clean.src = cms.untracked.string('iterativeCone5HiGenJets') |
90 |
< |
process.ic7clean = process.heavyIonCleanedGenJets.clone() |
91 |
< |
process.ic7clean.src = cms.untracked.string('iterativeCone7HiGenJets') |
92 |
< |
|
93 |
< |
process.ak5match = process.patJetGenJetMatch.clone() |
94 |
< |
process.ak5match.src = cms.InputTag("akFastPu5CaloJets") |
95 |
< |
process.ak5match.matched = cms.InputTag("ak5clean") |
96 |
< |
process.ak7match = process.patJetGenJetMatch.clone() |
97 |
< |
process.ak7match.src = cms.InputTag("akFastPu7CaloJets") |
98 |
< |
process.ak7match.matched = cms.InputTag("ak7clean") |
99 |
< |
process.akPu5match = process.patJetGenJetMatch.clone() |
100 |
< |
process.akPu5match.src = cms.InputTag("akPu5CaloJets") |
101 |
< |
process.akPu5match.matched = cms.InputTag("ak5clean") |
102 |
< |
process.akPu7match = process.patJetGenJetMatch.clone() |
103 |
< |
process.akPu7match.src = cms.InputTag("akPu7CaloJets") |
104 |
< |
process.akPu7match.matched = cms.InputTag("ak7clean") |
105 |
< |
|
106 |
< |
process.kt4match = process.patJetGenJetMatch.clone() |
107 |
< |
process.kt4match.src = cms.InputTag("ktFastPu4CaloJets") |
108 |
< |
process.kt4match.matched = cms.InputTag("kt4clean") |
109 |
< |
process.kt6match = process.patJetGenJetMatch.clone() |
110 |
< |
process.kt6match.src = cms.InputTag("ktFastPu6CaloJets") |
111 |
< |
process.kt6match.matched = cms.InputTag("kt6clean") |
112 |
< |
process.ktPu4match = process.patJetGenJetMatch.clone() |
113 |
< |
process.ktPu4match.src = cms.InputTag("ktPu4CaloJets") |
114 |
< |
process.ktPu4match.matched = cms.InputTag("kt4clean") |
115 |
< |
process.ktPu6match = process.patJetGenJetMatch.clone() |
116 |
< |
process.ktPu6match.src = cms.InputTag("ktPu6CaloJets") |
117 |
< |
process.ktPu6match.matched = cms.InputTag("kt6clean") |
118 |
< |
|
119 |
< |
process.icPu5match = process.patJetGenJetMatch.clone() |
120 |
< |
process.icPu5match.src = cms.InputTag("iterativeConePu5CaloJets") |
121 |
< |
process.icPu5match.matched = cms.InputTag("ic5clean") |
122 |
< |
process.icPu7match = process.patJetGenJetMatch.clone() |
123 |
< |
process.icPu7match.src = cms.InputTag("iterativeConePu7CaloJets") |
124 |
< |
process.icPu7match.matched = cms.InputTag("ic7clean") |
125 |
< |
|
126 |
< |
process.patJets.addJetCorrFactors = True |
127 |
< |
process.patJets.addGenPartonMatch = False |
128 |
< |
process.patJets.addJetID = False |
129 |
< |
process.patJets.addGenJetMatch = True |
130 |
< |
process.patJets.embedGenJetMatch = True |
131 |
< |
|
132 |
< |
process.ak5patJets = process.patJets.clone() |
133 |
< |
process.ak5patJets.jetSource = cms.InputTag("akFastPu5CaloJets") |
134 |
< |
process.ak5patJets.genJetMatch = cms.InputTag("ak5match") |
135 |
< |
process.ak5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr") ) |
136 |
< |
|
137 |
< |
process.ak7patJets = process.patJets.clone() |
138 |
< |
process.ak7patJets.jetSource = cms.InputTag("akFastPu7CaloJets") |
139 |
< |
process.ak7patJets.genJetMatch = cms.InputTag("ak7match") |
140 |
< |
process.ak7patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr") ) |
141 |
< |
|
142 |
< |
process.akPu5patJets = process.patJets.clone() |
143 |
< |
process.akPu5patJets.jetSource = cms.InputTag("akPu5CaloJets") |
144 |
< |
process.akPu5patJets.genJetMatch = cms.InputTag("akPu5match") |
145 |
< |
process.akPu5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr") ) |
146 |
< |
|
147 |
< |
process.akPu7patJets = process.patJets.clone() |
148 |
< |
process.akPu7patJets.jetSource = cms.InputTag("akPu7CaloJets") |
149 |
< |
process.akPu7patJets.genJetMatch = cms.InputTag("akPu7match") |
150 |
< |
process.akPu7patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr") ) |
151 |
< |
|
152 |
< |
process.kt4patJets = process.patJets.clone() |
153 |
< |
process.kt4patJets.jetSource = cms.InputTag("ktFastPu4CaloJets") |
154 |
< |
process.kt4patJets.genJetMatch = cms.InputTag("kt4match") |
155 |
< |
process.kt4patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt4corr") ) |
156 |
< |
|
157 |
< |
process.kt6patJets = process.patJets.clone() |
158 |
< |
process.kt6patJets.jetSource = cms.InputTag("ktFastPu6CaloJets") |
159 |
< |
process.kt6patJets.genJetMatch = cms.InputTag("kt6match") |
160 |
< |
process.kt6patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt6corr") ) |
161 |
< |
|
162 |
< |
process.ktPu4patJets = process.patJets.clone() |
163 |
< |
process.ktPu4patJets.jetSource = cms.InputTag("ktPu4CaloJets") |
164 |
< |
process.ktPu4patJets.genJetMatch = cms.InputTag("ktPu4match") |
165 |
< |
process.ktPu4patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu4corr") ) |
166 |
< |
|
167 |
< |
process.ktPu6patJets = process.patJets.clone() |
168 |
< |
process.ktPu6patJets.jetSource = cms.InputTag("ktPu6CaloJets") |
169 |
< |
process.ktPu6patJets.genJetMatch = cms.InputTag("ktPu6match") |
170 |
< |
process.ktPu6patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu6corr") ) |
171 |
< |
|
172 |
< |
process.icPu5patJets = process.patJets.clone() |
173 |
< |
process.icPu5patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets") |
174 |
< |
process.icPu5patJets.genJetMatch = cms.InputTag("icPu5match") |
175 |
< |
process.icPu5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr") ) |
176 |
< |
|
177 |
< |
process.icPu7patJets = process.patJets.clone() |
178 |
< |
process.icPu7patJets.jetSource = cms.InputTag("iterativeConePu7CaloJets") |
179 |
< |
process.icPu7patJets.genJetMatch = cms.InputTag("icPu7match") |
180 |
< |
process.icPu7patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu7corr") ) |
34 |
> |
process.ak5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak5CaloJets"), |
35 |
> |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"), |
36 |
> |
L3Absolute = cms.string("L3Absolute_AK5Calo"))) |
37 |
> |
process.ak5clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('ak5HiGenJets')) |
38 |
> |
process.ak5match = process.patJetGenJetMatch.clone(src = cms.InputTag("ak5CaloJets"), |
39 |
> |
matched = cms.InputTag("ak5clean")) |
40 |
> |
process.ak5parton = process.patJetPartonMatch.clone(src = cms.InputTag("ak5CaloJets")) |
41 |
> |
process.ak5patJets = process.patJets.clone(jetSource = cms.InputTag("ak5CaloJets"), |
42 |
> |
genJetMatch = cms.InputTag("ak5match"), |
43 |
> |
genPartonMatch= cms.InputTag("ak5parton"), |
44 |
> |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr"))) |
45 |
> |
process.ak5patSequence = cms.Sequence(process.ak5CaloJets*process.ak5corr*process.ak5clean*process.ak5match*process.ak5parton*process.ak5patJets) |
46 |
> |
|
47 |
> |
process.ak7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak7CaloJets"), |
48 |
> |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"), |
49 |
> |
L3Absolute = cms.string("L3Absolute_AK7Calo"))) |
50 |
> |
process.ak7clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('ak7HiGenJets')) |
51 |
> |
process.ak7match = process.patJetGenJetMatch.clone(src = cms.InputTag("ak7CaloJets"), |
52 |
> |
matched = cms.InputTag("ak7clean")) |
53 |
> |
process.ak7parton = process.patJetPartonMatch.clone(src = cms.InputTag("ak7CaloJets")) |
54 |
> |
process.ak7patJets = process.patJets.clone(jetSource = cms.InputTag("ak7CaloJets"), |
55 |
> |
genJetMatch = cms.InputTag("ak7match"), |
56 |
> |
genPartonMatch= cms.InputTag("ak7parton"), |
57 |
> |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr"))) |
58 |
> |
process.ak7patSequence = cms.Sequence(process.ak7CaloJets*process.ak7corr*process.ak7clean*process.ak7match*process.ak7parton*process.ak7patJets) |
59 |
> |
|
60 |
> |
process.akPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu5CaloJets"), |
61 |
> |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"), |
62 |
> |
L3Absolute = cms.string("L3Absolute_AK5Calo"))) |
63 |
> |
process.akPu5clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('ak5HiGenJets')) |
64 |
> |
process.akPu5match = process.patJetGenJetMatch.clone(src = cms.InputTag("akPu5CaloJets"), |
65 |
> |
matched = cms.InputTag("akPu5clean")) |
66 |
> |
process.akPu5parton = process.patJetPartonMatch.clone(src = cms.InputTag("akPu5CaloJets")) |
67 |
> |
process.akPu5patJets = process.patJets.clone(jetSource = cms.InputTag("akPu5CaloJets"), |
68 |
> |
genJetMatch = cms.InputTag("akPu5match"), |
69 |
> |
genPartonMatch= cms.InputTag("akPu5parton"), |
70 |
> |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr"))) |
71 |
> |
process.akPu5patSequence = cms.Sequence(process.akPu5CaloJets*process.akPu5corr*process.akPu5clean*process.akPu5match*process.akPu5parton*process.akPu5patJets) |
72 |
> |
|
73 |
> |
process.akPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu7CaloJets"), |
74 |
> |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"), |
75 |
> |
L3Absolute = cms.string("L3Absolute_AK7Calo"))) |
76 |
> |
process.akPu7clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('ak7HiGenJets')) |
77 |
> |
process.akPu7match = process.patJetGenJetMatch.clone(src = cms.InputTag("akPu7CaloJets"), |
78 |
> |
matched = cms.InputTag("akPu7clean")) |
79 |
> |
process.akPu7parton = process.patJetPartonMatch.clone(src = cms.InputTag("akPu7CaloJets")) |
80 |
> |
process.akPu7patJets = process.patJets.clone(jetSource = cms.InputTag("akPu7CaloJets"), |
81 |
> |
genJetMatch = cms.InputTag("akPu7match"), |
82 |
> |
genPartonMatch= cms.InputTag("akPu7parton"), |
83 |
> |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr"))) |
84 |
> |
process.akPu7patSequence = cms.Sequence(process.akPu7CaloJets*process.akPu7corr*process.akPu7clean*process.akPu7match*process.akPu7parton*process.akPu7patJets) |
85 |
> |
|
86 |
> |
process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff") |
87 |
|
|
88 |
|
process.output = cms.OutputModule("PoolOutputModule", |
89 |
< |
fileName = cms.untracked.string(ivars.output), |
90 |
< |
outputCommands = cms.untracked.vstring( |
185 |
< |
'drop *', |
186 |
< |
'keep *_*FastPu*_*_JETS', |
187 |
< |
'drop recoCaloJets_*_*_*', |
188 |
< |
'keep recoGenJets_*_*_HISIGNAL', |
189 |
< |
'keep *_heavyIon_*_*', |
190 |
< |
'keep *_hiCentrality_*_*', |
191 |
< |
'keep *_towerMaker_*_JETS', |
192 |
< |
'keep patJets_*_*_*' |
193 |
< |
) |
89 |
> |
process.jetTrkSkimContent, |
90 |
> |
fileName = cms.untracked.string(ivars.output) |
91 |
|
) |
92 |
|
|
93 |
< |
process.matches = cms.Path( |
94 |
< |
process.ak5clean + |
95 |
< |
process.ak7clean + |
96 |
< |
process.kt4clean + |
97 |
< |
process.kt6clean + |
98 |
< |
process.ic5clean + |
99 |
< |
process.ic7clean + |
100 |
< |
|
204 |
< |
process.ak5match + |
205 |
< |
process.ak7match + |
206 |
< |
process.akPu5match + |
207 |
< |
process.akPu7match + |
208 |
< |
|
209 |
< |
process.kt4match + |
210 |
< |
process.kt6match + |
211 |
< |
process.ktPu4match + |
212 |
< |
process.ktPu6match + |
213 |
< |
|
214 |
< |
process.icPu5match + |
215 |
< |
process.icPu7match |
216 |
< |
) |
217 |
< |
|
218 |
< |
process.corrections = cms.Path( |
219 |
< |
process.ak5corr+ |
220 |
< |
process.ak7corr+ |
221 |
< |
process.akPu5corr+ |
222 |
< |
process.akPu7corr+ |
223 |
< |
|
224 |
< |
process.kt4corr+ |
225 |
< |
process.kt6corr+ |
226 |
< |
process.ktPu4corr+ |
227 |
< |
process.ktPu6corr+ |
228 |
< |
|
229 |
< |
process.icPu5corr+ |
230 |
< |
process.icPu7corr |
231 |
< |
) |
232 |
< |
|
233 |
< |
process.pats = cms.Path( |
234 |
< |
process.ak5patJets + |
235 |
< |
process.ak7patJets + |
236 |
< |
process.akPu5patJets + |
237 |
< |
process.akPu7patJets + |
238 |
< |
|
239 |
< |
process.kt4patJets + |
240 |
< |
process.kt6patJets + |
241 |
< |
process.ktPu4patJets + |
242 |
< |
process.ktPu6patJets + |
243 |
< |
|
244 |
< |
process.icPu5patJets + |
245 |
< |
process.icPu7patJets |
246 |
< |
) |
93 |
> |
process.runAllJets = cms.Path(process.hiGenParticlesForJets * |
94 |
> |
process.hiRecoGenJets + |
95 |
> |
process.makeHeavyIonJets + |
96 |
> |
process.akPu5patSequence + |
97 |
> |
process.akPu7patSequence + |
98 |
> |
process.ak5patSequence + |
99 |
> |
process.ak7patSequence |
100 |
> |
) |
101 |
|
|
102 |
|
process.out_step = cms.EndPath(process.output) |
103 |
+ |
|