| 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 = "file:/net/hisrv0001/home/nart/scratch/july/CMSSW_3_7_0_patch4/src/CmsHi/JetAnalysis/test/MinBias0711_runs11to20.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 |
skipEvents = cms.untracked.uint32(ivars.initialEvent)
|
| 23 |
)
|
| 24 |
|
| 25 |
process.maxEvents = cms.untracked.PSet(
|
| 26 |
input = cms.untracked.int32(ivars.maxEvents)
|
| 27 |
)
|
| 28 |
|
| 29 |
process.load('Configuration.StandardSequences.Geometry_cff')
|
| 30 |
process.load('RecoJets.Configuration.GenJetParticles_cff')
|
| 31 |
process.load('RecoHI.HiJetAlgos.HiGenJets_cff')
|
| 32 |
process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
|
| 33 |
process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
|
| 34 |
from PhysicsTools.PatAlgos.tools.heavyIonTools import *
|
| 35 |
configureHeavyIons(process)
|
| 36 |
|
| 37 |
#from PhysicsTools.PatAlgos.tools.jetTools import *
|
| 38 |
#switchJECSet( process, "Summer09_7TeV_ReReco332")
|
| 39 |
|
| 40 |
process.patJets.addGenPartonMatch = False
|
| 41 |
process.patJets.addGenJetMatch = False
|
| 42 |
|
| 43 |
process.icPu5CaloJets = process.iterativeConePu5CaloJets.clone()
|
| 44 |
process.icPu5corr = process.patJetCorrFactors.clone()
|
| 45 |
process.icPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu5CaloJets"),
|
| 46 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
|
| 47 |
L3Absolute = cms.string("L3Absolute_IC5Calo")))
|
| 48 |
process.icPu5patJets = process.patJets.clone(jetSource = cms.InputTag("icPu5CaloJets"),
|
| 49 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))
|
| 50 |
|
| 51 |
process.icPu5patSequence = cms.Sequence(process.icPu5CaloJets*process.icPu5corr*process.icPu5patJets)
|
| 52 |
|
| 53 |
process.ak5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak5CaloJets"),
|
| 54 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
|
| 55 |
L3Absolute = cms.string("L3Absolute_AK5Calo")))
|
| 56 |
process.ak5patJets = process.patJets.clone(jetSource = cms.InputTag("ak5CaloJets"),
|
| 57 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr")))
|
| 58 |
process.ak5patSequence = cms.Sequence(process.ak5CaloJets*process.ak5corr*process.ak5patJets)
|
| 59 |
|
| 60 |
|
| 61 |
process.ak7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ak7CaloJets"),
|
| 62 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
|
| 63 |
L3Absolute = cms.string("L3Absolute_AK7Calo")))
|
| 64 |
process.ak7patJets = process.patJets.clone(jetSource = cms.InputTag("ak7CaloJets"),
|
| 65 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak7corr")))
|
| 66 |
process.ak7patSequence = cms.Sequence(process.ak7CaloJets*process.ak7corr*process.ak7patJets)
|
| 67 |
|
| 68 |
process.akPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu5CaloJets"),
|
| 69 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK5Calo"),
|
| 70 |
L3Absolute = cms.string("L3Absolute_AK5Calo")))
|
| 71 |
process.akPu5patJets = process.patJets.clone(jetSource = cms.InputTag("akPu5CaloJets"),
|
| 72 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr")))
|
| 73 |
process.akPu5patSequence = cms.Sequence(process.akPu5CaloJets*process.akPu5corr*process.akPu5patJets)
|
| 74 |
|
| 75 |
|
| 76 |
process.akPu7corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("akPu7CaloJets"),
|
| 77 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_AK7Calo"),
|
| 78 |
L3Absolute = cms.string("L3Absolute_AK7Calo")))
|
| 79 |
process.akPu7patJets = process.patJets.clone(jetSource = cms.InputTag("akPu7CaloJets"),
|
| 80 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu7corr")))
|
| 81 |
process.akPu7patSequence = cms.Sequence(process.akPu7CaloJets*process.akPu7corr*process.akPu7patJets)
|
| 82 |
|
| 83 |
process.kt4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt4CaloJets"),
|
| 84 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
|
| 85 |
L3Absolute = cms.string("L3Absolute_KT4Calo")))
|
| 86 |
process.kt4patJets = process.patJets.clone(jetSource = cms.InputTag("kt4CaloJets"),
|
| 87 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt4corr")))
|
| 88 |
process.kt4patSequence = cms.Sequence(process.kt4CaloJets*process.kt4corr*process.kt4patJets)
|
| 89 |
|
| 90 |
process.kt6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("kt6CaloJets"),
|
| 91 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
|
| 92 |
L3Absolute = cms.string("L3Absolute_KT6Calo")))
|
| 93 |
process.kt6patJets = process.patJets.clone(jetSource = cms.InputTag("kt6CaloJets"),
|
| 94 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("kt6corr")))
|
| 95 |
|
| 96 |
process.kt6patSequence = cms.Sequence(process.kt6CaloJets*process.kt6corr*process.kt6patJets)
|
| 97 |
|
| 98 |
process.ktPu4corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu4CaloJets"),
|
| 99 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT4Calo"),
|
| 100 |
L3Absolute = cms.string("L3Absolute_KT4Calo")))
|
| 101 |
process.ktPu4patJets = process.patJets.clone(jetSource = cms.InputTag("ktPu4CaloJets"),
|
| 102 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu4corr")))
|
| 103 |
process.ktPu4patSequence = cms.Sequence(process.ktPu4CaloJets*process.ktPu4corr*process.ktPu4patJets)
|
| 104 |
|
| 105 |
process.ktPu6corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("ktPu6CaloJets"),
|
| 106 |
corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_KT6Calo"),
|
| 107 |
L3Absolute = cms.string("L3Absolute_KT6Calo")))
|
| 108 |
process.ktPu6patJets = process.patJets.clone(jetSource = cms.InputTag("ktPu6CaloJets"),
|
| 109 |
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ktPu6corr")))
|
| 110 |
|
| 111 |
process.ktPu6patSequence = cms.Sequence(process.ktPu6CaloJets*process.ktPu6corr*process.ktPu6patJets)
|
| 112 |
|
| 113 |
process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")
|
| 114 |
|
| 115 |
process.output = cms.OutputModule("PoolOutputModule",
|
| 116 |
process.jetTrkSkimContent,
|
| 117 |
fileName = cms.untracked.string(ivars.output)
|
| 118 |
)
|
| 119 |
|
| 120 |
process.runAllJets = cms.Path(#process.hiGenParticlesForJets *
|
| 121 |
#process.hiRecoGenJets +
|
| 122 |
process.iterativeConePu5CaloJets +
|
| 123 |
# process.makeHeavyIonJets +
|
| 124 |
process.ktPu4patSequence +
|
| 125 |
process.ktPu6patSequence +
|
| 126 |
process.kt4patSequence +
|
| 127 |
process.kt6patSequence +
|
| 128 |
|
| 129 |
process.akPu5patSequence +
|
| 130 |
process.akPu7patSequence +
|
| 131 |
process.ak5patSequence +
|
| 132 |
process.ak7patSequence
|
| 133 |
|
| 134 |
)
|
| 135 |
|
| 136 |
process.out_step = cms.EndPath(process.output)
|
| 137 |
|
