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
Log Message:	remove time-taking algos
#	User	Rev	Content
1	yilmaz	1.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	yilmaz	1.4	ivars.files = "dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/nart/MinimumBiasHI/HighPtHardEnrichedAnaSkimNZS_03/b651c142d334b6448e9d198ad759ddac/allenrichedNZS_4_1_OCO.root"
9	yilmaz	1.1
10			ivars.output = 'jets_pat_80to120.root'
11	yilmaz	1.3	ivars.maxEvents = -1
12	yilmaz	1.1	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	yilmaz	1.3	fileNames = cms.untracked.vstring(ivars.files)
22	yilmaz	1.1	)
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	yilmaz	1.4	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	yilmaz	1.1	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	yilmaz	1.2	process.load("HeavyIonsAnalysis.Configuration.analysisProducers_cff")
121	yilmaz	1.1	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	yilmaz	1.2	process.output.outputCommands.extend(["drop _towerMaker__*"])
129
130			process.runAllJets = cms.Path(
131			process.allTracks +
132			process.icPu5patSequence +
133	yilmaz	1.4	process.icPu7patSequence +
134			# process.ktPu4patSequence +
135			# process.ktPu6patSequence +
136	yilmaz	1.1	process.kt4patSequence +
137			process.kt6patSequence +
138	yilmaz	1.4	# process.akPu5patSequence +
139			# process.akPu7patSequence +
140	yilmaz	1.1	process.ak5patSequence +
141			process.ak7patSequence
142	yilmaz	1.2	)
143	yilmaz	1.1
144			process.out_step = cms.EndPath(process.output)
145