ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/makePatJets.py
Revision: 1.3
Committed: Thu Jun 3 14:29:34 2010 UTC (14 years, 11 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.2: +7 -5 lines
Log Message: 
updated cfg

File Contents

# Content
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.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('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 )
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")
34
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 process.jetTrkSkimContent,
91 fileName = cms.untracked.string(ivars.output)
92 )
93
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