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.3 by yilmaz, Thu Jun 3 14:29:34 2010 UTC

# Line 1 | Line 1
1   import FWCore.ParameterSet.VarParsing as VarParsing
2  
3   ivars = VarParsing.VarParsing('standard')
4 + ivars.files = 'rfio:/castor/cern.ch/user/y/yilmaz/pat/CMSSW_3_6_0/UnquenchedDijet80to120_runs1to500.root'
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  
7 < #ivars.files = 'file:PyquenDijet80to120_d20100325_runs1to500.root'
8 < #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
7 > ivars.output = 'jets_pat_80to120.root'
8 > ivars.maxEvents = 3
9  
10   ivars.parseArguments()
11  
12   import FWCore.ParameterSet.Config as cms
13  
14 < process = cms.Process('MATCH')
14 > process = cms.Process('HIJETS')
15  
16   process.source = cms.Source("PoolSource",
17                              fileNames = cms.untracked.vstring(ivars.files)
18                              )
19  
20   process.maxEvents = cms.untracked.PSet(
21 <        input = cms.untracked.int32(ivars.maxEvents)
22 <        )
21 >            input = cms.untracked.int32(ivars.maxEvents)
22 >                    )
23  
24 + process.load('Configuration.StandardSequences.Geometry_cff')
25 + process.load('RecoJets.Configuration.GenJetParticles_cff')
26 + process.load('RecoHI.HiJetAlgos.HiGenJets_cff')
27 + process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
28   process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
29   from PhysicsTools.PatAlgos.tools.heavyIonTools import *
30   configureHeavyIons(process)
31  
32 < from PhysicsTools.PatAlgos.tools.jetTools import *
33 < switchJECSet( process, "Summer09_7TeV_ReReco332")
32 > #from PhysicsTools.PatAlgos.tools.jetTools import *
33 > #switchJECSet( process, "Summer09_7TeV_ReReco332")
34  
35 < process.ak5corr = process.patJetCorrFactors.clone()
36 < process.ak5corr.jetSource = cms.InputTag("akFastPu5CaloJets")
37 < process.ak5corr.corrLevels.L2Relative = "L2Relative_AK5Calo"
38 < process.ak5corr.corrLevels.L3Absolute = "L3Absolute_AK5Calo"
39 <
40 < process.ak7corr = process.patJetCorrFactors.clone()
41 < process.ak7corr.jetSource = cms.InputTag("akFastPu7CaloJets")
42 < process.ak7corr.corrLevels.L2Relative = cms.string('L2Relative_AK7Calo')
43 < process.ak7corr.corrLevels.L3Absolute = cms.string('L3Absolute_AK7Calo')
44 <
45 < process.akPu5corr = process.ak5corr.clone()
46 < process.akPu5corr.jetSource = cms.InputTag("akPu5CaloJets")
47 < process.akPu7corr = process.ak7corr.clone()
48 < process.akPu7corr.jetSource = cms.InputTag("akPu7CaloJets")
49 <
50 < process.kt4corr = process.patJetCorrFactors.clone()
51 < process.kt4corr.jetSource = cms.InputTag("ktFastPu4CaloJets")
52 < process.kt4corr.corrLevels.L2Relative = cms.string('L2Relative_KT4Calo')
53 < process.kt4corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT4Calo')
54 <
55 < process.kt6corr = process.patJetCorrFactors.clone()
56 < process.kt6corr.jetSource = cms.InputTag("ktFastPu6CaloJets")
57 < process.kt6corr.corrLevels.L2Relative = cms.string('L2Relative_KT6Calo')
58 < process.kt6corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT6Calo')
59 <
60 < process.ktPu4corr = process.patJetCorrFactors.clone()
61 < process.ktPu4corr.jetSource = cms.InputTag("ktPu4CaloJets")
62 < process.ktPu4corr.corrLevels.L2Relative = cms.string('L2Relative_KT4Calo')
63 < process.ktPu4corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT4Calo')
64 <
65 < process.ktPu6corr = process.patJetCorrFactors.clone()
66 < process.ktPu6corr.jetSource = cms.InputTag("ktPu6CaloJets")
67 < process.ktPu6corr.corrLevels.L2Relative = cms.string('L2Relative_KT6Calo')
68 < process.ktPu6corr.corrLevels.L3Absolute = cms.string('L3Absolute_KT6Calo')
69 <      
70 < process.icPu5corr = process.patJetCorrFactors.clone()
71 < process.icPu5corr.jetSource = cms.InputTag("iterativeConePu5CaloJets")
72 < process.icPu5corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo')
73 < process.icPu5corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo')
74 <
75 < process.icPu7corr = process.patJetCorrFactors.clone()
76 < process.icPu7corr.jetSource = cms.InputTag("iterativeConePu7CaloJets")
77 < process.icPu7corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo') # No IC7
78 < process.icPu7corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo') # No IC7
79 <
80 < process.ak5clean = process.heavyIonCleanedGenJets.clone()
81 < process.ak5clean.src = cms.untracked.string('ak5HiGenJets')
82 < process.ak7clean = process.heavyIonCleanedGenJets.clone()
83 < process.ak7clean.src = cms.untracked.string('ak7HiGenJets')
84 < process.kt4clean = process.heavyIonCleanedGenJets.clone()
85 < process.kt4clean.src = cms.untracked.string('kt4HiGenJets')
86 < process.kt6clean = process.heavyIonCleanedGenJets.clone()
87 < 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") )
35 > process.ak5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak5CaloJets"),
36 >                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
37 >                                                                        L3Absolute = cms.string("L3Absolute_AK5Calo")))
38 > process.ak5clean = process.heavyIonCleanedGenJets.clone(src = cms.InputTag('ak5HiGenJets'))
39 > process.ak5match = process.patJetGenJetMatch.clone(src = cms.InputTag("ak5CaloJets"),
40 >                                                   matched = cms.InputTag("ak5clean"))
41 > process.ak5parton = process.patJetPartonMatch.clone(src = cms.InputTag("ak5CaloJets"))
42 > process.ak5patJets = process.patJets.clone(jetSource  = cms.InputTag("ak5CaloJets"),
43 >                                           genJetMatch = cms.InputTag("ak5match"),
44 >                                           genPartonMatch= cms.InputTag("ak5parton"),
45 >                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr")))
46 > process.ak5patSequence = cms.Sequence(process.ak5CaloJets*process.ak5corr*process.ak5clean*process.ak5match*process.ak5parton*process.ak5patJets)
47 >
48 > process.ak7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak7CaloJets"),
49 >                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
50 >                                                                        L3Absolute = cms.string("L3Absolute_AK7Calo")))
51 > process.ak7clean = process.heavyIonCleanedGenJets.clone(src = cms.InputTag('ak7HiGenJets'))
52 > process.ak7match = process.patJetGenJetMatch.clone(src = cms.InputTag("ak7CaloJets"),
53 >                                                   matched = cms.InputTag("ak7clean"))
54 > process.ak7parton = process.patJetPartonMatch.clone(src = cms.InputTag("ak7CaloJets"))
55 > process.ak7patJets = process.patJets.clone(jetSource  = cms.InputTag("ak7CaloJets"),
56 >                                           genJetMatch = cms.InputTag("ak7match"),
57 >                                           genPartonMatch= cms.InputTag("ak7parton"),
58 >                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr")))
59 > process.ak7patSequence = cms.Sequence(process.ak7CaloJets*process.ak7corr*process.ak7clean*process.ak7match*process.ak7parton*process.ak7patJets)
60 >
61 > process.akPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu5CaloJets"),
62 >                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
63 >                                                                          L3Absolute = cms.string("L3Absolute_AK5Calo")))
64 > process.akPu5clean = process.heavyIonCleanedGenJets.clone(src = cms.InputTag('ak5HiGenJets'))
65 > process.akPu5match = process.patJetGenJetMatch.clone(src = cms.InputTag("akPu5CaloJets"),
66 >                                                     matched = cms.InputTag("akPu5clean"))
67 > process.akPu5parton = process.patJetPartonMatch.clone(src = cms.InputTag("akPu5CaloJets"))
68 > process.akPu5patJets = process.patJets.clone(jetSource  = cms.InputTag("akPu5CaloJets"),
69 >                                             genJetMatch = cms.InputTag("akPu5match"),
70 >                                             genPartonMatch= cms.InputTag("akPu5parton"),
71 >                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr")))
72 > process.akPu5patSequence = cms.Sequence(process.akPu5CaloJets*process.akPu5corr*process.akPu5clean*process.akPu5match*process.akPu5parton*process.akPu5patJets)
73 >
74 > process.akPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu7CaloJets"),
75 >                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
76 >                                                                          L3Absolute = cms.string("L3Absolute_AK7Calo")))
77 > process.akPu7clean = process.heavyIonCleanedGenJets.clone(src = cms.InputTag('ak7HiGenJets'))
78 > process.akPu7match = process.patJetGenJetMatch.clone(src = cms.InputTag("akPu7CaloJets"),
79 >                                                     matched = cms.InputTag("akPu7clean"))
80 > process.akPu7parton = process.patJetPartonMatch.clone(src = cms.InputTag("akPu7CaloJets"))
81 > process.akPu7patJets = process.patJets.clone(jetSource  = cms.InputTag("akPu7CaloJets"),
82 >                                             genJetMatch = cms.InputTag("akPu7match"),
83 >                                             genPartonMatch= cms.InputTag("akPu7parton"),
84 >                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr")))
85 > process.akPu7patSequence = cms.Sequence(process.akPu7CaloJets*process.akPu7corr*process.akPu7clean*process.akPu7match*process.akPu7parton*process.akPu7patJets)
86 >
87 > process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")
88  
89   process.output = cms.OutputModule("PoolOutputModule",
90 <                                  fileName = cms.untracked.string(ivars.output),
91 <                                  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 <    )
90 >                                  process.jetTrkSkimContent,
91 >                                  fileName = cms.untracked.string(ivars.output)
92                                    )
93  
94 < process.matches = cms.Path(
95 <    process.ak5clean +
96 <    process.ak7clean +
97 <    process.kt4clean +
98 <    process.kt6clean +
99 <    process.ic5clean +
100 <    process.ic7clean +
101 <
102 <    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 <    )
94 > process.runAllJets = cms.Path(process.hiGenParticlesForJets *
95 >                              process.hiRecoGenJets +
96 >                              process.iterativeConePu5CaloJets +
97 >                              process.makeHeavyIonJets +
98 >                              process.akPu5patSequence +
99 >                              process.akPu7patSequence +
100 >                              process.ak5patSequence +
101 >                              process.ak7patSequence
102 >                              )
103  
104   process.out_step = cms.EndPath(process.output)
105 +

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines