| 1 |
frankma |
1.1 |
import FWCore.ParameterSet.Config as cms
|
| 2 |
|
|
|
| 3 |
|
|
# Pat Jets
|
| 4 |
|
|
from PhysicsTools.PatAlgos.patHeavyIonSequences_cff import *
|
| 5 |
|
|
patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets")
|
| 6 |
|
|
patJets.addBTagInfo = False
|
| 7 |
|
|
patJets.addTagInfos = False
|
| 8 |
|
|
patJets.addDiscriminators = False
|
| 9 |
|
|
patJets.addAssociatedTracks = False
|
| 10 |
|
|
patJets.addJetCharge = False
|
| 11 |
|
|
patJets.addJetID = False
|
| 12 |
|
|
patJets.getJetMCFlavour = False
|
| 13 |
|
|
patJets.addGenPartonMatch = True
|
| 14 |
|
|
patJets.addGenJetMatch = True
|
| 15 |
|
|
patJets.embedGenJetMatch = True
|
| 16 |
|
|
patJets.embedGenPartonMatch = True
|
| 17 |
|
|
patJets.embedCaloTowers = False
|
| 18 |
|
|
|
| 19 |
yilmaz |
1.5 |
patJetCorrFactors.useNPV = False
|
| 20 |
frankma |
1.1 |
# full reco
|
| 21 |
|
|
icPu5corr = patJetCorrFactors.clone(
|
| 22 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
| 23 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
| 24 |
yilmaz |
1.7 |
payload = cms.string('IC5Calo_2760GeV')
|
| 25 |
frankma |
1.1 |
)
|
| 26 |
yilmaz |
1.6 |
|
| 27 |
frankma |
1.1 |
icPu5patJets = patJets.clone(
|
| 28 |
|
|
jetSource = cms.InputTag("iterativeConePu5CaloJets"),
|
| 29 |
|
|
genJetMatch = cms.InputTag("icPu5match"),
|
| 30 |
|
|
genPartonMatch = cms.InputTag("icPu5parton"),
|
| 31 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr"))
|
| 32 |
|
|
)
|
| 33 |
|
|
|
| 34 |
|
|
akPu5PFcorr = icPu5corr.clone(
|
| 35 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
| 36 |
yilmaz |
1.7 |
payload = cms.string('AK5PFTowers_hiGoodTightTracks')
|
| 37 |
frankma |
1.1 |
)
|
| 38 |
|
|
akPu5PFpatJets = patJets.clone(
|
| 39 |
|
|
jetSource = cms.InputTag("akPu5PFJets"),
|
| 40 |
|
|
genJetMatch = cms.InputTag("akPu5PFmatch"),
|
| 41 |
|
|
genPartonMatch = cms.InputTag("akPu5PFparton"),
|
| 42 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr"))
|
| 43 |
|
|
)
|
| 44 |
|
|
|
| 45 |
|
|
akPu3PFcorr = icPu5corr.clone(
|
| 46 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
| 47 |
yilmaz |
1.7 |
payload = cms.string('AK3PFTowers_hiGoodTightTracks')
|
| 48 |
frankma |
1.1 |
)
|
| 49 |
yilmaz |
1.6 |
|
| 50 |
frankma |
1.1 |
akPu3PFpatJets = patJets.clone(
|
| 51 |
|
|
jetSource = cms.InputTag("akPu3PFJets"),
|
| 52 |
|
|
genJetMatch = cms.InputTag("akPu3PFmatch"),
|
| 53 |
|
|
genPartonMatch = cms.InputTag("akPu3PFparton"),
|
| 54 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr"))
|
| 55 |
|
|
)
|
| 56 |
|
|
|
| 57 |
yilmaz |
1.6 |
akPu5corr = icPu5corr.clone(
|
| 58 |
|
|
src = cms.InputTag("akPu5CaloJets"),
|
| 59 |
|
|
payload = cms.string('AK5Calo')
|
| 60 |
|
|
)
|
| 61 |
|
|
|
| 62 |
|
|
akPu5patJets = patJets.clone(
|
| 63 |
|
|
jetSource = cms.InputTag("akPu5CaloJets"),
|
| 64 |
|
|
genJetMatch = cms.InputTag("akPu5match"),
|
| 65 |
|
|
genPartonMatch = cms.InputTag("akPu5parton"),
|
| 66 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr"))
|
| 67 |
|
|
)
|
| 68 |
|
|
|
| 69 |
|
|
akPu3corr = icPu5corr.clone(
|
| 70 |
|
|
src = cms.InputTag("akPu3CaloJets"),
|
| 71 |
yilmaz |
1.7 |
payload = cms.string('AK5Calo')
|
| 72 |
yilmaz |
1.6 |
)
|
| 73 |
|
|
|
| 74 |
|
|
akPu3patJets = patJets.clone(
|
| 75 |
|
|
jetSource = cms.InputTag("akPu3CaloJets"),
|
| 76 |
|
|
genJetMatch = cms.InputTag("akPu3match"),
|
| 77 |
|
|
genPartonMatch = cms.InputTag("akPu3parton"),
|
| 78 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3corr"))
|
| 79 |
|
|
)
|
| 80 |
|
|
|
| 81 |
|
|
|
| 82 |
frankma |
1.1 |
# === data sequences ===
|
| 83 |
|
|
# Note still need to use enableData function in cfg to remove mc dep of patjet
|
| 84 |
|
|
icPu5patSequence_data = cms.Sequence( icPu5corr * icPu5patJets )
|
| 85 |
|
|
akPu5PFpatSequence_data = cms.Sequence(akPu5PFcorr * akPu5PFpatJets )
|
| 86 |
|
|
akPu3PFpatSequence_data = cms.Sequence( akPu3PFcorr * akPu3PFpatJets )
|
| 87 |
yilmaz |
1.6 |
akPu5patSequence_data = cms.Sequence( akPu5corr * akPu5patJets )
|
| 88 |
|
|
akPu3patSequence_data = cms.Sequence( akPu3corr * akPu3patJets )
|
| 89 |
frankma |
1.1 |
|
| 90 |
|
|
# mc matching
|
| 91 |
|
|
patJetPartonMatch.matched = cms.InputTag("hiPartons")
|
| 92 |
|
|
|
| 93 |
|
|
icPu5clean = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets') )
|
| 94 |
|
|
icPu5match = patJetGenJetMatch.clone(
|
| 95 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
| 96 |
|
|
matched = cms.InputTag("icPu5clean")
|
| 97 |
|
|
)
|
| 98 |
|
|
icPu5parton = patJetPartonMatch.clone(
|
| 99 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets")
|
| 100 |
|
|
)
|
| 101 |
|
|
|
| 102 |
yilmaz |
1.3 |
|
| 103 |
frankma |
1.1 |
akPu5PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets') )
|
| 104 |
|
|
akPu5PFmatch = patJetGenJetMatch.clone(
|
| 105 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
| 106 |
|
|
matched = cms.InputTag("akPu5PFclean")
|
| 107 |
|
|
)
|
| 108 |
|
|
akPu5PFparton = patJetPartonMatch.clone(
|
| 109 |
|
|
src = cms.InputTag("akPu5PFJets")
|
| 110 |
|
|
)
|
| 111 |
|
|
|
| 112 |
|
|
akPu3PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets') )
|
| 113 |
|
|
akPu3PFmatch = patJetGenJetMatch.clone(
|
| 114 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
| 115 |
|
|
matched = cms.InputTag("akPu3PFclean")
|
| 116 |
|
|
)
|
| 117 |
|
|
akPu3PFparton = patJetPartonMatch.clone(
|
| 118 |
|
|
src = cms.InputTag("akPu3PFJets")
|
| 119 |
|
|
)
|
| 120 |
|
|
|
| 121 |
yilmaz |
1.6 |
akPu5match = patJetGenJetMatch.clone(
|
| 122 |
|
|
src = cms.InputTag("akPu5CaloJets"),
|
| 123 |
|
|
matched = cms.InputTag("akPu5PFclean")
|
| 124 |
|
|
)
|
| 125 |
|
|
akPu5parton = patJetPartonMatch.clone(
|
| 126 |
|
|
src = cms.InputTag("akPu5CaloJets")
|
| 127 |
|
|
)
|
| 128 |
|
|
|
| 129 |
|
|
akPu3match = patJetGenJetMatch.clone(
|
| 130 |
|
|
src = cms.InputTag("akPu3CaloJets"),
|
| 131 |
|
|
matched = cms.InputTag("akPu3PFclean")
|
| 132 |
|
|
)
|
| 133 |
|
|
akPu3parton = patJetPartonMatch.clone(
|
| 134 |
|
|
src = cms.InputTag("akPu3CaloJets")
|
| 135 |
|
|
)
|
| 136 |
|
|
|
| 137 |
frankma |
1.1 |
# === mc sequences ===
|
| 138 |
|
|
icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5patJets)
|
| 139 |
|
|
akPu5PFpatSequence = cms.Sequence(akPu5PFcorr * akPu5PFclean * akPu5PFmatch * akPu5PFparton * akPu5PFpatJets)
|
| 140 |
|
|
akPu3PFpatSequence = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFpatJets)
|
| 141 |
yilmaz |
1.6 |
akPu5patSequence = cms.Sequence(akPu5corr * akPu5PFclean * akPu5match * akPu5parton * akPu5patJets)
|
| 142 |
|
|
akPu3patSequence = cms.Sequence(akPu3corr * akPu3PFclean * akPu3match * akPu3parton * akPu3patJets)
|
