| 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 |
|
|
|
| 20 |
|
|
# full reco
|
| 21 |
|
|
icPu5corr = patJetCorrFactors.clone(
|
| 22 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
| 23 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
| 24 |
|
|
payload = cms.string('IC5Calo')
|
| 25 |
|
|
)
|
| 26 |
|
|
icPu5patJets = patJets.clone(
|
| 27 |
|
|
jetSource = cms.InputTag("iterativeConePu5CaloJets"),
|
| 28 |
|
|
genJetMatch = cms.InputTag("icPu5match"),
|
| 29 |
|
|
genPartonMatch = cms.InputTag("icPu5parton"),
|
| 30 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr"))
|
| 31 |
|
|
)
|
| 32 |
|
|
|
| 33 |
|
|
akPu5PFcorr = icPu5corr.clone(
|
| 34 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
| 35 |
|
|
payload = cms.string('AK5PF')
|
| 36 |
|
|
)
|
| 37 |
|
|
akPu5PFpatJets = patJets.clone(
|
| 38 |
|
|
jetSource = cms.InputTag("akPu5PFJets"),
|
| 39 |
|
|
genJetMatch = cms.InputTag("akPu5PFmatch"),
|
| 40 |
|
|
genPartonMatch = cms.InputTag("akPu5PFparton"),
|
| 41 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr"))
|
| 42 |
|
|
)
|
| 43 |
|
|
|
| 44 |
|
|
akPu3PFcorr = icPu5corr.clone(
|
| 45 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
| 46 |
|
|
payload = cms.string('AK5PF')
|
| 47 |
|
|
)
|
| 48 |
|
|
akPu3PFpatJets = patJets.clone(
|
| 49 |
|
|
jetSource = cms.InputTag("akPu3PFJets"),
|
| 50 |
|
|
genJetMatch = cms.InputTag("akPu3PFmatch"),
|
| 51 |
|
|
genPartonMatch = cms.InputTag("akPu3PFparton"),
|
| 52 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr"))
|
| 53 |
|
|
)
|
| 54 |
|
|
|
| 55 |
|
|
# === data sequences ===
|
| 56 |
|
|
# Note still need to use enableData function in cfg to remove mc dep of patjet
|
| 57 |
|
|
icPu5patSequence_data = cms.Sequence( icPu5corr * icPu5patJets )
|
| 58 |
|
|
akPu5PFpatSequence_data = cms.Sequence(akPu5PFcorr * akPu5PFpatJets )
|
| 59 |
|
|
akPu3PFpatSequence_data = cms.Sequence( akPu3PFcorr * akPu3PFpatJets )
|
| 60 |
|
|
|
| 61 |
|
|
# mc matching
|
| 62 |
|
|
patJetPartonMatch.matched = cms.InputTag("hiPartons")
|
| 63 |
|
|
|
| 64 |
|
|
icPu5clean = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets') )
|
| 65 |
|
|
icPu5match = patJetGenJetMatch.clone(
|
| 66 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
| 67 |
|
|
matched = cms.InputTag("icPu5clean")
|
| 68 |
|
|
)
|
| 69 |
|
|
icPu5parton = patJetPartonMatch.clone(
|
| 70 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets")
|
| 71 |
|
|
)
|
| 72 |
|
|
|
| 73 |
|
|
|
| 74 |
|
|
akPu5PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets') )
|
| 75 |
|
|
akPu5PFmatch = patJetGenJetMatch.clone(
|
| 76 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
| 77 |
|
|
matched = cms.InputTag("akPu5PFclean")
|
| 78 |
|
|
)
|
| 79 |
|
|
akPu5PFparton = patJetPartonMatch.clone(
|
| 80 |
|
|
src = cms.InputTag("akPu5PFJets")
|
| 81 |
|
|
)
|
| 82 |
|
|
|
| 83 |
|
|
akPu3PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets') )
|
| 84 |
|
|
akPu3PFmatch = patJetGenJetMatch.clone(
|
| 85 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
| 86 |
|
|
matched = cms.InputTag("akPu3PFclean")
|
| 87 |
|
|
)
|
| 88 |
|
|
akPu3PFparton = patJetPartonMatch.clone(
|
| 89 |
|
|
src = cms.InputTag("akPu3PFJets")
|
| 90 |
|
|
)
|
| 91 |
|
|
|
| 92 |
|
|
# === mc sequences ===
|
| 93 |
|
|
icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5patJets)
|
| 94 |
|
|
akPu5PFpatSequence = cms.Sequence(akPu5PFcorr * akPu5PFclean * akPu5PFmatch * akPu5PFparton * akPu5PFpatJets)
|
| 95 |
|
|
akPu3PFpatSequence = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFpatJets)
|
| 96 |
frankma |
1.1.2.1 |
|
| 97 |
|
|
# Pat Photon
|
| 98 |
|
|
from PhysicsTools.PatAlgos.producersHeavyIons.heavyIonPhotons_cff import *
|
| 99 |
|
|
from PhysicsTools.PatAlgos.producersLayer1.photonProducer_cff import *
|
| 100 |
|
|
from PhysicsTools.PatAlgos.selectionLayer1.photonSelector_cfi import selectedPatPhotons
|
| 101 |
|
|
from RecoHI.HiEgammaAlgos.HiEgammaIsolation_cff import * # (must be decleared after PAT sequence Yen-Jie)
|
| 102 |
|
|
|
| 103 |
|
|
# from PhysicsTools/PatAlgos/python/tools/heavyIonTools.py
|
| 104 |
|
|
photonMatch.matched = cms.InputTag("hiGenParticles")
|
| 105 |
|
|
patPhotons.addPhotonID = cms.bool(True)
|
| 106 |
|
|
patPhotons.addGenMatch = cms.bool(True)
|
| 107 |
|
|
patPhotons.embedGenMatch = cms.bool(True)
|
| 108 |
|
|
patPhotons.userData.userFloats.src = cms.VInputTag(
|
| 109 |
|
|
cms.InputTag( "isoCC1"),cms.InputTag( "isoCC2"),cms.InputTag( "isoCC3"),cms.InputTag( "isoCC4"),cms.InputTag("isoCC5"),
|
| 110 |
|
|
cms.InputTag( "isoCR1"),cms.InputTag( "isoCR2"),cms.InputTag( "isoCR3"),cms.InputTag( "isoCR4"),cms.InputTag("isoCR5"),
|
| 111 |
|
|
cms.InputTag( "isoT11"),cms.InputTag( "isoT12"),cms.InputTag( "isoT13"),cms.InputTag( "isoT14"),
|
| 112 |
|
|
cms.InputTag( "isoT21"),cms.InputTag( "isoT22"),cms.InputTag( "isoT23"),cms.InputTag( "isoT24"),
|
| 113 |
|
|
cms.InputTag( "isoT31"),cms.InputTag( "isoT32"),cms.InputTag( "isoT33"),cms.InputTag( "isoT34"),
|
| 114 |
|
|
cms.InputTag( "isoT41"),cms.InputTag( "isoT42"),cms.InputTag( "isoT43"),cms.InputTag( "isoT44"),
|
| 115 |
|
|
cms.InputTag("isoDR11"),cms.InputTag("isoDR12"),cms.InputTag("isoDR13"),cms.InputTag("isoDR14"),
|
| 116 |
|
|
cms.InputTag("isoDR21"),cms.InputTag("isoDR22"),cms.InputTag("isoDR23"),cms.InputTag("isoDR24"),
|
| 117 |
|
|
cms.InputTag("isoDR31"),cms.InputTag("isoDR32"),cms.InputTag("isoDR33"),cms.InputTag("isoDR34"),
|
| 118 |
|
|
cms.InputTag("isoDR41"),cms.InputTag("isoDR42"),cms.InputTag("isoDR43"),cms.InputTag("isoDR44")
|
| 119 |
|
|
)
|
| 120 |
|
|
del patPhotons.photonIDSources
|
| 121 |
|
|
|
| 122 |
|
|
# pat photon selection
|
| 123 |
|
|
selectedPatPhotons.cut = cms.string('pt > 0. & abs(eta) < 12.')
|
| 124 |
|
|
|
| 125 |
|
|
makeHeavyIonPhotons = cms.Sequence(
|
| 126 |
|
|
# reco pre-production
|
| 127 |
|
|
hiEgammaIsolationSequence *
|
| 128 |
|
|
patPhotonIsolation *
|
| 129 |
|
|
# pat and HI specifics
|
| 130 |
|
|
photonMatch *
|
| 131 |
|
|
# object production
|
| 132 |
|
|
patPhotons *
|
| 133 |
|
|
selectedPatPhotons
|
| 134 |
|
|
) |
