ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/makePatJets.py
(Generate patch)

Comparing UserCode/CmsHi/JetAnalysis/python/makePatJets.py (file contents):
Revision 1.1 by yilmaz, Thu Apr 22 16:41:09 2010 UTC vs.
Revision 1.2 by yilmaz, Tue May 4 17:35:19 2010 UTC

# Line 2 | Line 2 | import FWCore.ParameterSet.VarParsing as
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 +

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines