JetAnalysis/python/makePatJetsData.py

import FWCore.ParameterSet.VarParsing as VarParsing

ivars = VarParsing.VarParsing('standard')
ivars.register('initialEvent',mult=ivars.multiplicity.singleton,info="for testing")
#ivars.files = 'dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/nart/MinimumBiasHI/HighPtMinBiasAnaSkimHardEnriched_01/6edc9684fb662a4038079f7b1d3fbd82/minbiashardenriched_2_1_wUe.root'
#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'

ivars.files = "dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/nart/MinimumBiasHI/HighPtHardEnrichedAnaSkimNZS_03/b651c142d334b6448e9d198ad759ddac/allenrichedNZS_4_1_OCO.root"

ivars.output = 'jets_pat_80to120.root'
ivars.maxEvents = -1
ivars.initialEvent = 1

ivars.parseArguments()

import FWCore.ParameterSet.Config as cms

process = cms.Process('HIJETS')

process.source = cms.Source("PoolSource",
                            fileNames = cms.untracked.vstring(ivars.files)
                            )

process.maxEvents = cms.untracked.PSet(
            input = cms.untracked.int32(ivars.maxEvents)
                    )

process.load('Configuration.StandardSequences.Geometry_cff')
process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
from PhysicsTools.PatAlgos.tools.heavyIonTools import *
configureHeavyIons(process)

#from PhysicsTools.PatAlgos.tools.jetTools import *
#switchJECSet( process, "Summer09_7TeV_ReReco332")

process.patJets.addGenPartonMatch   = False
process.patJets.addGenJetMatch      = False

process.icPu5CaloJets = process.iterativeConePu5CaloJets.clone()
process.icPu5corr = process.patJetCorrFactors.clone()
process.icPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu5CaloJets"),
                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
                                                                        L3Absolute = cms.string("L3Absolute_IC5Calo")))
process.icPu5patJets = process.patJets.clone(jetSource  = cms.InputTag("icPu5CaloJets"),
                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))

process.icPu5patSequence = cms.Sequence(process.icPu5CaloJets*process.icPu5corr*process.icPu5patJets)

process.icPu7CaloJets = process.iterativeConePu7CaloJets.clone()
process.icPu7corr = process.patJetCorrFactors.clone()
process.icPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu7CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_IC5Calo")))
process.icPu7patJets = process.patJets.clone(jetSource  = cms.InputTag("icPu7CaloJets"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu7corr")))

process.icPu7patSequence = cms.Sequence(process.icPu7CaloJets*process.icPu7corr*process.icPu7patJets)

process.ak5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak5CaloJets"),
                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
                                                                        L3Absolute = cms.string("L3Absolute_AK5Calo")))
process.ak5patJets = process.patJets.clone(jetSource  = cms.InputTag("ak5CaloJets"),
                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr")))
process.ak5patSequence = cms.Sequence(process.ak5CaloJets*process.ak5corr*process.ak5patJets)


process.ak7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak7CaloJets"),
                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
                                                                        L3Absolute = cms.string("L3Absolute_AK7Calo")))
process.ak7patJets = process.patJets.clone(jetSource  = cms.InputTag("ak7CaloJets"),
                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr")))
process.ak7patSequence = cms.Sequence(process.ak7CaloJets*process.ak7corr*process.ak7patJets)

process.akPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu5CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_AK5Calo")))
process.akPu5patJets = process.patJets.clone(jetSource  = cms.InputTag("akPu5CaloJets"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr")))
process.akPu5patSequence = cms.Sequence(process.akPu5CaloJets*process.akPu5corr*process.akPu5patJets)


process.akPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu7CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_AK7Calo")))
process.akPu7patJets = process.patJets.clone(jetSource  = cms.InputTag("akPu7CaloJets"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr")))
process.akPu7patSequence = cms.Sequence(process.akPu7CaloJets*process.akPu7corr*process.akPu7patJets)

process.kt4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt4CaloJets"),
                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
                                                                        L3Absolute = cms.string("L3Absolute_KT4Calo")))
process.kt4patJets = process.patJets.clone(jetSource  = cms.InputTag("kt4CaloJets"),
                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt4corr")))
process.kt4patSequence = cms.Sequence(process.kt4CaloJets*process.kt4corr*process.kt4patJets)

process.kt6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt6CaloJets"),
                                                  corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
                                                                        L3Absolute = cms.string("L3Absolute_KT6Calo")))
process.kt6patJets = process.patJets.clone(jetSource  = cms.InputTag("kt6CaloJets"),
                                           jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt6corr")))

process.kt6patSequence = cms.Sequence(process.kt6CaloJets*process.kt6corr*process.kt6patJets)

process.ktPu4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu4CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_KT4Calo")))
process.ktPu4patJets = process.patJets.clone(jetSource  = cms.InputTag("ktPu4CaloJets"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu4corr")))
process.ktPu4patSequence = cms.Sequence(process.ktPu4CaloJets*process.ktPu4corr*process.ktPu4patJets)

process.ktPu6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu6CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_KT6Calo")))
process.ktPu6patJets = process.patJets.clone(jetSource  = cms.InputTag("ktPu6CaloJets"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu6corr")))

process.ktPu6patSequence = cms.Sequence(process.ktPu6CaloJets*process.ktPu6corr*process.ktPu6patJets)

process.load("HeavyIonsAnalysis.Configuration.analysisProducers_cff")
process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")

process.output = cms.OutputModule("PoolOutputModule",
                                  process.jetTrkSkimContent,
                                  fileName = cms.untracked.string(ivars.output)
                                  )

process.output.outputCommands.extend(["drop *_towerMaker_*_*"])

process.runAllJets = cms.Path(
    process.allTracks +
    process.icPu5patSequence +
    process.icPu7patSequence +
#    process.ktPu4patSequence +
#    process.ktPu6patSequence +
    process.kt4patSequence +
    process.kt6patSequence +
#    process.akPu5patSequence +
#    process.akPu7patSequence +
    process.ak5patSequence +
    process.ak7patSequence 
    )

process.out_step = cms.EndPath(process.output)

Revision:	1.4
Committed:	Tue Jul 13 08:07:55 2010 UTC (14 years, 9 months ago) by yilmaz
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	HiForest_V02_85, HiForest_V02_84, HiForest_V02_83, HiForest_V02_82, HiForest_V02_81, HiForest_V02_80, HiForest_V02_79, HiForest_V02_78, HiForest_V02_77, HiForest_V02_76, HiForest_V02_75, HiForest_V02_74, HiForest_V02_73, HiForest_V02_72, HiForest_V02_71, HiForest_V02_70, HiForest_V02_69, HiForest_V02_68, HiForest_V02_67, HiForest_V02_66, HiForest_V02_65, HiForest_V02_64, HiForest_V02_63, HiForest_V02_62, HiForest_V02_61, HiForest_V02_60, HiForest_V02_59, HiForest_V02_58, HiForest_V02_57, HiForest_V02_56, HiForest_V02_55, HiForest_V02_54, HiForest_V02_53, HiForest_V02_52, HiForest_V02_51, HiForest_V02_50, HiForest_V02_49, HiForest_V02_48, HiForest_V02_47, HiForest_V02_46, HiForest_V02_45, HiForest_V02_44, HiForest_V02_43, HiForest_V02_42, HiForest_V02_41, HiForest_V02_40, HiForest_V02_39, HiForest_V02_38, HiForest_V02_37, HiForest_V02_36, HiForest_V02_35, HiForest_V02_34, HiForest_V02_33, HiForest_V02_32, HiForest_V02_31, HiForest_V02_30, HiForest_V02_27, HiForest_V02_26, QM_2012, HiForest_V02_25, HiForest_V02_24, HiForest_V02_23, HiForest_V02_22, HiForest_V02_21, HiForest_V02_20, HiForest_V02_19, HiForest_V02_18, HiForest_V02_17, HiForest_V02_16, HiForest_V02_15, HiForest_V02_14, HiForest_V02_13, HiForest_V02_12, HiForest_V02_11, HiForest_V02_10, HiForest_V02_09, HiForest_V02_08, HiForest_V02_07, HiForest_V02_06, HiForest_V02_05, HiForest_V02_04, HiForest_V02_03, HiForest_V02_02, HiForest_V02_01, HiForest_V02_00, hi44X_02, hi441_1, hi441_0, hi413_11, hi413_10, hi413_09, hi413_08, hi413_07, hi413_06, hi413_05, hi413_04, hi413_03, hi413_02, hi413_01, hi39X_01, tag_d20110915, cmssw39x_base, cmssw39X_base, HEAD
Branch point for:	branch_44x, cmssw39x_branch
Changes since 1.3:	+16 -5 lines
Error occurred while calculating annotation data.
Log Message:	remove time-taking algos
#	Content
1	import FWCore.ParameterSet.VarParsing as VarParsing
2
3	ivars = VarParsing.VarParsing('standard')
4	ivars.register('initialEvent',mult=ivars.multiplicity.singleton,info="for testing")
5	#ivars.files = 'dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/nart/MinimumBiasHI/HighPtMinBiasAnaSkimHardEnriched_01/6edc9684fb662a4038079f7b1d3fbd82/minbiashardenriched_2_1_wUe.root'
6	#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'
7
8	ivars.files = "dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/nart/MinimumBiasHI/HighPtHardEnrichedAnaSkimNZS_03/b651c142d334b6448e9d198ad759ddac/allenrichedNZS_4_1_OCO.root"
9
10	ivars.output = 'jets_pat_80to120.root'
11	ivars.maxEvents = -1
12	ivars.initialEvent = 1
13
14	ivars.parseArguments()
15
16	import FWCore.ParameterSet.Config as cms
17
18	process = cms.Process('HIJETS')
19
20	process.source = cms.Source("PoolSource",
21	fileNames = cms.untracked.vstring(ivars.files)
22	)
23
24	process.maxEvents = cms.untracked.PSet(
25	input = cms.untracked.int32(ivars.maxEvents)
26	)
27
28	process.load('Configuration.StandardSequences.Geometry_cff')
29	process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
30	process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
31	from PhysicsTools.PatAlgos.tools.heavyIonTools import *
32	configureHeavyIons(process)
33
34	#from PhysicsTools.PatAlgos.tools.jetTools import *
35	#switchJECSet( process, "Summer09_7TeV_ReReco332")
36
37	process.patJets.addGenPartonMatch = False
38	process.patJets.addGenJetMatch = False
39
40	process.icPu5CaloJets = process.iterativeConePu5CaloJets.clone()
41	process.icPu5corr = process.patJetCorrFactors.clone()
42	process.icPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu5CaloJets"),
43	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
44	L3Absolute = cms.string("L3Absolute_IC5Calo")))
45	process.icPu5patJets = process.patJets.clone(jetSource = cms.InputTag("icPu5CaloJets"),
46	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))
47
48	process.icPu5patSequence = cms.Sequence(process.icPu5CaloJetsprocess.icPu5corrprocess.icPu5patJets)
49
50	process.icPu7CaloJets = process.iterativeConePu7CaloJets.clone()
51	process.icPu7corr = process.patJetCorrFactors.clone()
52	process.icPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu7CaloJets"),
53	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
54	L3Absolute = cms.string("L3Absolute_IC5Calo")))
55	process.icPu7patJets = process.patJets.clone(jetSource = cms.InputTag("icPu7CaloJets"),
56	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu7corr")))
57
58	process.icPu7patSequence = cms.Sequence(process.icPu7CaloJetsprocess.icPu7corrprocess.icPu7patJets)
59
60	process.ak5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak5CaloJets"),
61	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
62	L3Absolute = cms.string("L3Absolute_AK5Calo")))
63	process.ak5patJets = process.patJets.clone(jetSource = cms.InputTag("ak5CaloJets"),
64	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr")))
65	process.ak5patSequence = cms.Sequence(process.ak5CaloJetsprocess.ak5corrprocess.ak5patJets)
66
67
68	process.ak7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak7CaloJets"),
69	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
70	L3Absolute = cms.string("L3Absolute_AK7Calo")))
71	process.ak7patJets = process.patJets.clone(jetSource = cms.InputTag("ak7CaloJets"),
72	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr")))
73	process.ak7patSequence = cms.Sequence(process.ak7CaloJetsprocess.ak7corrprocess.ak7patJets)
74
75	process.akPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu5CaloJets"),
76	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
77	L3Absolute = cms.string("L3Absolute_AK5Calo")))
78	process.akPu5patJets = process.patJets.clone(jetSource = cms.InputTag("akPu5CaloJets"),
79	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr")))
80	process.akPu5patSequence = cms.Sequence(process.akPu5CaloJetsprocess.akPu5corrprocess.akPu5patJets)
81
82
83	process.akPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu7CaloJets"),
84	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
85	L3Absolute = cms.string("L3Absolute_AK7Calo")))
86	process.akPu7patJets = process.patJets.clone(jetSource = cms.InputTag("akPu7CaloJets"),
87	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr")))
88	process.akPu7patSequence = cms.Sequence(process.akPu7CaloJetsprocess.akPu7corrprocess.akPu7patJets)
89
90	process.kt4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt4CaloJets"),
91	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
92	L3Absolute = cms.string("L3Absolute_KT4Calo")))
93	process.kt4patJets = process.patJets.clone(jetSource = cms.InputTag("kt4CaloJets"),
94	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt4corr")))
95	process.kt4patSequence = cms.Sequence(process.kt4CaloJetsprocess.kt4corrprocess.kt4patJets)
96
97	process.kt6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt6CaloJets"),
98	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
99	L3Absolute = cms.string("L3Absolute_KT6Calo")))
100	process.kt6patJets = process.patJets.clone(jetSource = cms.InputTag("kt6CaloJets"),
101	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt6corr")))
102
103	process.kt6patSequence = cms.Sequence(process.kt6CaloJetsprocess.kt6corrprocess.kt6patJets)
104
105	process.ktPu4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu4CaloJets"),
106	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
107	L3Absolute = cms.string("L3Absolute_KT4Calo")))
108	process.ktPu4patJets = process.patJets.clone(jetSource = cms.InputTag("ktPu4CaloJets"),
109	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu4corr")))
110	process.ktPu4patSequence = cms.Sequence(process.ktPu4CaloJetsprocess.ktPu4corrprocess.ktPu4patJets)
111
112	process.ktPu6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu6CaloJets"),
113	corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
114	L3Absolute = cms.string("L3Absolute_KT6Calo")))
115	process.ktPu6patJets = process.patJets.clone(jetSource = cms.InputTag("ktPu6CaloJets"),
116	jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu6corr")))
117
118	process.ktPu6patSequence = cms.Sequence(process.ktPu6CaloJetsprocess.ktPu6corrprocess.ktPu6patJets)
119
120	process.load("HeavyIonsAnalysis.Configuration.analysisProducers_cff")
121	process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")
122
123	process.output = cms.OutputModule("PoolOutputModule",
124	process.jetTrkSkimContent,
125	fileName = cms.untracked.string(ivars.output)
126	)
127
128	process.output.outputCommands.extend(["drop _towerMaker__*"])
129
130	process.runAllJets = cms.Path(
131	process.allTracks +
132	process.icPu5patSequence +
133	process.icPu7patSequence +
134	# process.ktPu4patSequence +
135	# process.ktPu6patSequence +
136	process.kt4patSequence +
137	process.kt6patSequence +
138	# process.akPu5patSequence +
139	# process.akPu7patSequence +
140	process.ak5patSequence +
141	process.ak7patSequence
142	)
143
144	process.out_step = cms.EndPath(process.output)
145