| 1 |
frankma |
1.1 |
import FWCore.ParameterSet.Config as cms
|
| 2 |
|
|
|
| 3 |
|
|
process = cms.Process('JetAna')
|
| 4 |
|
|
|
| 5 |
|
|
process.options = cms.untracked.PSet(
|
| 6 |
|
|
wantSummary = cms.untracked.bool(True)
|
| 7 |
|
|
)
|
| 8 |
|
|
|
| 9 |
|
|
# Input source
|
| 10 |
|
|
process.source = cms.Source("PoolSource",
|
| 11 |
|
|
fileNames = cms.untracked.vstring(
|
| 12 |
|
|
'/store/relval/CMSSW_3_9_9/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-RECO/START39_V7HI-v1/0000/40B4255D-5A3D-E011-9D80-001BFCDBD100.root'
|
| 13 |
|
|
))
|
| 14 |
|
|
|
| 15 |
|
|
process.maxEvents = cms.untracked.PSet(
|
| 16 |
|
|
input = cms.untracked.int32(10))
|
| 17 |
|
|
|
| 18 |
|
|
#load some general stuff
|
| 19 |
|
|
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
|
| 20 |
|
|
process.load('Configuration.StandardSequences.Services_cff')
|
| 21 |
|
|
process.load('Configuration.StandardSequences.GeometryExtended_cff')
|
| 22 |
|
|
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
|
| 23 |
|
|
process.load('FWCore.MessageService.MessageLogger_cfi')
|
| 24 |
|
|
#process.MessageLogger.cerr.FwkReport.reportEvery = 100
|
| 25 |
|
|
process.GlobalTag.globaltag = 'GR_R_39X_V6B::All'
|
| 26 |
|
|
|
| 27 |
|
|
# load centrality
|
| 28 |
|
|
from CmsHi.Analysis2010.CommonFunctions_cff import *
|
| 29 |
|
|
overrideCentrality(process)
|
| 30 |
|
|
process.HeavyIonGlobalParameters = cms.PSet(
|
| 31 |
|
|
centralityVariable = cms.string("HFhits"),
|
| 32 |
|
|
nonDefaultGlauberModel = cms.string(""),
|
| 33 |
|
|
centralitySrc = cms.InputTag("hiCentrality")
|
| 34 |
|
|
)
|
| 35 |
|
|
|
| 36 |
|
|
# tree output
|
| 37 |
|
|
process.TFileService = cms.Service("TFileService",
|
| 38 |
|
|
fileName=cms.string("JetAnaTrees.root"))
|
| 39 |
|
|
|
| 40 |
|
|
|
| 41 |
|
|
# Define Jet Algo parameters
|
| 42 |
|
|
process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
|
| 43 |
|
|
#process.load('RecoHI.HiJetAlgos.HiRecoPFJets_cff')
|
| 44 |
|
|
|
| 45 |
|
|
process.iterativeConePu5CaloJets.doPVCorrection = cms.bool(True)
|
| 46 |
|
|
process.iterativeConePu5CaloJets.srcPVs = 'hiSelectedVertex'
|
| 47 |
|
|
process.iterativeConePu5CaloJets.puPtMin = cms.double(10.0)
|
| 48 |
|
|
|
| 49 |
|
|
|
| 50 |
|
|
# Define Jet energy correction
|
| 51 |
|
|
process.jec = cms.ESSource("PoolDBESSource",
|
| 52 |
|
|
DBParameters = cms.PSet(messageLevel = cms.untracked.int32(0)),
|
| 53 |
|
|
timetype = cms.string('runnumber'),
|
| 54 |
|
|
toGet = cms.VPSet(
|
| 55 |
|
|
cms.PSet(record = cms.string("JetCorrectionsRecord"),
|
| 56 |
|
|
tag = cms.string("JetCorrectorParametersCollection_HI_hiGoodTightTracks_D6T_399_IC5Calo"),
|
| 57 |
|
|
label = cms.untracked.string("IC5Calo")),
|
| 58 |
|
|
cms.PSet(record = cms.string("JetCorrectionsRecord"),
|
| 59 |
|
|
tag = cms.string("JetCorrectorParametersCollection_HI_PFTowers_hiGoodTightTracks_D6T_399_AK3PF"),
|
| 60 |
|
|
label = cms.untracked.string("AK3PF")),
|
| 61 |
|
|
cms.PSet(record = cms.string("JetCorrectionsRecord"),
|
| 62 |
|
|
tag = cms.string("JetCorrectorParametersCollection_HI_PFTowers_hiGoodTightTracks_D6T_399_AK4PF"),
|
| 63 |
|
|
label = cms.untracked.string("AK4PF")),
|
| 64 |
|
|
cms.PSet(record = cms.string("JetCorrectionsRecord"),
|
| 65 |
|
|
tag = cms.string("JetCorrectorParametersCollection_HI_PFTowers_hiGoodTightTracks_D6T_399_AK5PF"),
|
| 66 |
|
|
label = cms.untracked.string("AK5PF")
|
| 67 |
|
|
),
|
| 68 |
|
|
),
|
| 69 |
|
|
connect = cms.string("sqlite_file:JEC_HI_PFTowers_hiGoodTightTracks_D6T_399.db"),
|
| 70 |
|
|
)
|
| 71 |
|
|
process.es_prefer_jec = cms.ESPrefer('PoolDBESSource','jec')
|
| 72 |
|
|
|
| 73 |
|
|
# Define Pat
|
| 74 |
|
|
process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
|
| 75 |
|
|
process.patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets")
|
| 76 |
|
|
process.patJets.addBTagInfo = False
|
| 77 |
|
|
process.patJets.addTagInfos = False
|
| 78 |
|
|
process.patJets.addDiscriminators = False
|
| 79 |
|
|
process.patJets.addAssociatedTracks = False
|
| 80 |
|
|
process.patJets.addJetCharge = False
|
| 81 |
|
|
process.patJets.addJetID = False
|
| 82 |
|
|
process.patJets.getJetMCFlavour = False
|
| 83 |
|
|
process.patJets.addGenPartonMatch = False
|
| 84 |
|
|
process.patJets.addGenJetMatch = False
|
| 85 |
|
|
process.patJets.embedGenJetMatch = False
|
| 86 |
|
|
process.patJets.embedGenPartonMatch = False
|
| 87 |
|
|
|
| 88 |
|
|
process.icPu5corr = process.patJetCorrFactors.clone(
|
| 89 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
| 90 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
| 91 |
|
|
payload = cms.string('IC5Calo')
|
| 92 |
|
|
)
|
| 93 |
|
|
process.icPu5patJets = process.patJets.clone(
|
| 94 |
|
|
jetSource = cms.InputTag("iterativeConePu5CaloJets"),
|
| 95 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))
|
| 96 |
|
|
|
| 97 |
|
|
process.icPu5patSequence = cms.Sequence(process.iterativeConePu5CaloJets*process.icPu5corr*process.icPu5patJets)
|
| 98 |
|
|
|
| 99 |
|
|
|
| 100 |
|
|
# Define the analyzer modules
|
| 101 |
|
|
process.load("HLTrigger.HLTanalyzers.HI_HLTAnalyser_cff")
|
| 102 |
|
|
process.hltanalysis.RunParameters.Debug = False
|
| 103 |
|
|
process.hltanalysis.RunParameters.UseTFileService = True
|
| 104 |
|
|
process.hltanalysis.RunParameters.Monte = (False)
|
| 105 |
|
|
process.hltanalysis.RunParameters.DoMC = (False)
|
| 106 |
|
|
process.hltanalysis.RunParameters.DoJets = True
|
| 107 |
|
|
process.hltanalysis.RunParameters.DoPhotons = True
|
| 108 |
|
|
process.hltanalysis.RunParameters.DoSuperClusters = True
|
| 109 |
|
|
process.hltanalysis.recjets = "iterativeConePu5CaloJets"
|
| 110 |
|
|
process.hltanalysis.BarrelPhoton = "correctedIslandBarrelSuperClusters"
|
| 111 |
|
|
process.hltanalysis.EndcapPhoton = "correctedIslandEndcapSuperClusters"
|
| 112 |
|
|
process.hltanalysis.l1GtReadoutRecord = cms.InputTag("gtDigis")
|
| 113 |
|
|
|
| 114 |
|
|
process.load("MNguyen.InclusiveJetAnalyzer.inclusiveJetAnalyzer_cff")
|
| 115 |
|
|
process.inclusiveJetAnalyzer.isMC = cms.untracked.bool(False)
|
| 116 |
|
|
process.icPu5JetAnalyzer = process.inclusiveJetAnalyzer.clone()
|
| 117 |
|
|
|
| 118 |
|
|
# =============== Final Paths =====================
|
| 119 |
|
|
process.reco_extra = cms.Path( process.icPu5patSequence )
|
| 120 |
|
|
process.ana_step = cms.Path( process.centralityBin * process.hltanalysis * process.icPu5JetAnalyzer)
|
| 121 |
|
|
process.schedule = cms.Schedule(process.reco_extra,process.ana_step)
|
