Skip to content

Commit

Permalink
Merge pull request cms-sw#103 from dgulhan/forest_JEC_pA_80X_8160GeV_v0
Browse files Browse the repository at this point in the history 
Forest jec p a 80 x 8160 ge v v0
  • Loading branch information
@kurtejung
kurtejung committed Dec 13, 2016
2 parents af73978 + 01647f2 commit c9bd685
Show file tree
Hide file tree
Showing 34 changed files with 3,556 additions and 47 deletions.
37 changes: 13 additions & 24 deletions HeavyIonsAnalysis/Configuration/python/CommonFunctions_cff.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -755,46 +755,35 @@ def overrideJEC_pPb8TeV(process):
process.GlobalTag.toGet.extend([
## no PU
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4Calo"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AK4Calo")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4PF"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AK4PF")
),
##adding PU = nonPU TEMPORARILY
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4Calo"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4Calo_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AK4Calo_offline")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4PF"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4PF_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AK4PF_offline")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4Calo"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AKPu4Calo")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4PF"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4PF_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AKPu4PF")
),
label = cms.untracked.string("AK3PF_offline")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4Calo"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4Calo_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AKPu4Calo_offline")
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_Spring16_25nsV6_MC_AK4PF"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4PF_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AKPu4PF_offline")
)
),
cms.PSet(record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_pA_80X_8160GeV_v0_AK4PF_offline"),
connect = cms.string("frontier://FrontierProd/CMS_CONDITIONS"),
label = cms.untracked.string("AKPu3PF_offline")
)
])
return process

Expand Down
251 changes: 251 additions & 0 deletions HeavyIonsAnalysis/JetAnalysis/python/jets/ak3CaloJetSequence_pPb_data_cff.py
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,251 @@


import FWCore.ParameterSet.Config as cms
from HeavyIonsAnalysis.JetAnalysis.patHeavyIonSequences_cff import patJetGenJetMatch, patJetPartonMatch, patJetCorrFactors, patJets
from HeavyIonsAnalysis.JetAnalysis.inclusiveJetAnalyzer_cff import *
from HeavyIonsAnalysis.JetAnalysis.bTaggers_cff import *
from RecoJets.JetProducers.JetIDParams_cfi import *
from RecoJets.JetProducers.nJettinessAdder_cfi import Njettiness

ak3Calomatch = patJetGenJetMatch.clone(
src = cms.InputTag("ak3CaloJets"),
matched = cms.InputTag("ak3GenJets"),
resolveByMatchQuality = cms.bool(False),
maxDeltaR = 0.3
)

ak3CalomatchGroomed = patJetGenJetMatch.clone(
src = cms.InputTag("ak3GenJets"),
matched = cms.InputTag("ak3GenJets"),
resolveByMatchQuality = cms.bool(False),
maxDeltaR = 0.3
)

ak3Caloparton = patJetPartonMatch.clone(src = cms.InputTag("ak3CaloJets")
)

ak3Calocorr = patJetCorrFactors.clone(
useNPV = cms.bool(False),
useRho = cms.bool(False),
# primaryVertices = cms.InputTag("hiSelectedVertex"),
levels = cms.vstring('L2Relative','L3Absolute'),
src = cms.InputTag("ak3CaloJets"),
payload = "AK3Calo_offline"
)

ak3CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak3CaloJets'))

#ak3Caloclean = heavyIonCleanedGenJets.clone(src = cms.InputTag('ak3GenJets'))

ak3CalobTagger = bTaggers("ak3Calo",0.3)

#create objects locally since they dont load properly otherwise
#ak3Calomatch = ak3CalobTagger.match
ak3Caloparton = patJetPartonMatch.clone(src = cms.InputTag("ak3CaloJets"), matched = cms.InputTag("genParticles"))
ak3CaloPatJetFlavourAssociationLegacy = ak3CalobTagger.PatJetFlavourAssociationLegacy
ak3CaloPatJetPartons = ak3CalobTagger.PatJetPartons
ak3CaloJetTracksAssociatorAtVertex = ak3CalobTagger.JetTracksAssociatorAtVertex
ak3CaloJetTracksAssociatorAtVertex.tracks = cms.InputTag("highPurityTracks")
ak3CaloSimpleSecondaryVertexHighEffBJetTags = ak3CalobTagger.SimpleSecondaryVertexHighEffBJetTags
ak3CaloSimpleSecondaryVertexHighPurBJetTags = ak3CalobTagger.SimpleSecondaryVertexHighPurBJetTags
ak3CaloCombinedSecondaryVertexBJetTags = ak3CalobTagger.CombinedSecondaryVertexBJetTags
ak3CaloCombinedSecondaryVertexV2BJetTags = ak3CalobTagger.CombinedSecondaryVertexV2BJetTags
ak3CaloJetBProbabilityBJetTags = ak3CalobTagger.JetBProbabilityBJetTags
ak3CaloSoftPFMuonByPtBJetTags = ak3CalobTagger.SoftPFMuonByPtBJetTags
ak3CaloSoftPFMuonByIP3dBJetTags = ak3CalobTagger.SoftPFMuonByIP3dBJetTags
ak3CaloTrackCountingHighEffBJetTags = ak3CalobTagger.TrackCountingHighEffBJetTags
ak3CaloTrackCountingHighPurBJetTags = ak3CalobTagger.TrackCountingHighPurBJetTags
ak3CaloPatJetPartonAssociationLegacy = ak3CalobTagger.PatJetPartonAssociationLegacy

ak3CaloImpactParameterTagInfos = ak3CalobTagger.ImpactParameterTagInfos
ak3CaloImpactParameterTagInfos.primaryVertex = cms.InputTag("offlinePrimaryVertices")
ak3CaloJetProbabilityBJetTags = ak3CalobTagger.JetProbabilityBJetTags

ak3CaloSecondaryVertexTagInfos = ak3CalobTagger.SecondaryVertexTagInfos
ak3CaloSimpleSecondaryVertexHighEffBJetTags = ak3CalobTagger.SimpleSecondaryVertexHighEffBJetTags
ak3CaloSimpleSecondaryVertexHighPurBJetTags = ak3CalobTagger.SimpleSecondaryVertexHighPurBJetTags
ak3CaloCombinedSecondaryVertexBJetTags = ak3CalobTagger.CombinedSecondaryVertexBJetTags
ak3CaloCombinedSecondaryVertexV2BJetTags = ak3CalobTagger.CombinedSecondaryVertexV2BJetTags

ak3CaloSecondaryVertexNegativeTagInfos = ak3CalobTagger.SecondaryVertexNegativeTagInfos
ak3CaloNegativeSimpleSecondaryVertexHighEffBJetTags = ak3CalobTagger.NegativeSimpleSecondaryVertexHighEffBJetTags
ak3CaloNegativeSimpleSecondaryVertexHighPurBJetTags = ak3CalobTagger.NegativeSimpleSecondaryVertexHighPurBJetTags
ak3CaloNegativeCombinedSecondaryVertexBJetTags = ak3CalobTagger.NegativeCombinedSecondaryVertexBJetTags
ak3CaloPositiveCombinedSecondaryVertexBJetTags = ak3CalobTagger.PositiveCombinedSecondaryVertexBJetTags
ak3CaloNegativeCombinedSecondaryVertexV2BJetTags = ak3CalobTagger.NegativeCombinedSecondaryVertexV2BJetTags
ak3CaloPositiveCombinedSecondaryVertexV2BJetTags = ak3CalobTagger.PositiveCombinedSecondaryVertexV2BJetTags

ak3CaloSoftPFMuonsTagInfos = ak3CalobTagger.SoftPFMuonsTagInfos
ak3CaloSoftPFMuonsTagInfos.primaryVertex = cms.InputTag("offlinePrimaryVertices")
ak3CaloSoftPFMuonBJetTags = ak3CalobTagger.SoftPFMuonBJetTags
ak3CaloSoftPFMuonByIP3dBJetTags = ak3CalobTagger.SoftPFMuonByIP3dBJetTags
ak3CaloSoftPFMuonByPtBJetTags = ak3CalobTagger.SoftPFMuonByPtBJetTags
ak3CaloNegativeSoftPFMuonByPtBJetTags = ak3CalobTagger.NegativeSoftPFMuonByPtBJetTags
ak3CaloPositiveSoftPFMuonByPtBJetTags = ak3CalobTagger.PositiveSoftPFMuonByPtBJetTags
ak3CaloPatJetFlavourIdLegacy = cms.Sequence(ak3CaloPatJetPartonAssociationLegacy*ak3CaloPatJetFlavourAssociationLegacy)
#Not working with our PU sub, but keep it here for reference
#ak3CaloPatJetFlavourAssociation = ak3CalobTagger.PatJetFlavourAssociation
#ak3CaloPatJetFlavourId = cms.Sequence(ak3CaloPatJetPartons*ak3CaloPatJetFlavourAssociation)

ak3CaloJetBtaggingIP = cms.Sequence(ak3CaloImpactParameterTagInfos *
(ak3CaloTrackCountingHighEffBJetTags +
ak3CaloTrackCountingHighPurBJetTags +
ak3CaloJetProbabilityBJetTags +
ak3CaloJetBProbabilityBJetTags
)
)

ak3CaloJetBtaggingSV = cms.Sequence(ak3CaloImpactParameterTagInfos
*
ak3CaloSecondaryVertexTagInfos
* (ak3CaloSimpleSecondaryVertexHighEffBJetTags+
ak3CaloSimpleSecondaryVertexHighPurBJetTags+
ak3CaloCombinedSecondaryVertexBJetTags+
ak3CaloCombinedSecondaryVertexV2BJetTags
)
)

ak3CaloJetBtaggingNegSV = cms.Sequence(ak3CaloImpactParameterTagInfos
*
ak3CaloSecondaryVertexNegativeTagInfos
* (ak3CaloNegativeSimpleSecondaryVertexHighEffBJetTags+
ak3CaloNegativeSimpleSecondaryVertexHighPurBJetTags+
ak3CaloNegativeCombinedSecondaryVertexBJetTags+
ak3CaloPositiveCombinedSecondaryVertexBJetTags+
ak3CaloNegativeCombinedSecondaryVertexV2BJetTags+
ak3CaloPositiveCombinedSecondaryVertexV2BJetTags
)
)

ak3CaloJetBtaggingMu = cms.Sequence(ak3CaloSoftPFMuonsTagInfos * (ak3CaloSoftPFMuonBJetTags
+
ak3CaloSoftPFMuonByIP3dBJetTags
+
ak3CaloSoftPFMuonByPtBJetTags
+
ak3CaloNegativeSoftPFMuonByPtBJetTags
+
ak3CaloPositiveSoftPFMuonByPtBJetTags
)
)

ak3CaloJetBtagging = cms.Sequence(ak3CaloJetBtaggingIP
*ak3CaloJetBtaggingSV
*ak3CaloJetBtaggingNegSV
# *ak3CaloJetBtaggingMu
)

ak3CalopatJetsWithBtagging = patJets.clone(jetSource = cms.InputTag("ak3CaloJets"),
genJetMatch = cms.InputTag("ak3Calomatch"),
genPartonMatch = cms.InputTag("ak3Caloparton"),
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3Calocorr")),
JetPartonMapSource = cms.InputTag("ak3CaloPatJetFlavourAssociationLegacy"),
JetFlavourInfoSource = cms.InputTag("ak3CaloPatJetFlavourAssociation"),
trackAssociationSource = cms.InputTag("ak3CaloJetTracksAssociatorAtVertex"),
useLegacyJetMCFlavour = True,
discriminatorSources = cms.VInputTag(cms.InputTag("ak3CaloSimpleSecondaryVertexHighEffBJetTags"),
cms.InputTag("ak3CaloSimpleSecondaryVertexHighPurBJetTags"),
cms.InputTag("ak3CaloCombinedSecondaryVertexBJetTags"),
cms.InputTag("ak3CaloCombinedSecondaryVertexV2BJetTags"),
cms.InputTag("ak3CaloJetBProbabilityBJetTags"),
cms.InputTag("ak3CaloJetProbabilityBJetTags"),
#cms.InputTag("ak3CaloSoftPFMuonByPtBJetTags"),
#cms.InputTag("ak3CaloSoftPFMuonByIP3dBJetTags"),
cms.InputTag("ak3CaloTrackCountingHighEffBJetTags"),
cms.InputTag("ak3CaloTrackCountingHighPurBJetTags"),
),
jetIDMap = cms.InputTag("ak3CaloJetID"),
addBTagInfo = True,
addTagInfos = True,
addDiscriminators = True,
addAssociatedTracks = True,
addJetCharge = False,
addJetID = False,
getJetMCFlavour = False,
addGenPartonMatch = False,
addGenJetMatch = False,
embedGenJetMatch = False,
embedGenPartonMatch = False,
# embedCaloTowers = False,
# embedPFCandidates = True
)

ak3CaloNjettiness = Njettiness.clone(
src = cms.InputTag("ak3CaloJets"),
R0 = cms.double( 0.3)
)
ak3CalopatJetsWithBtagging.userData.userFloats.src += ['ak3CaloNjettiness:tau1','ak3CaloNjettiness:tau2','ak3CaloNjettiness:tau3']

ak3CaloJetAnalyzer = inclusiveJetAnalyzer.clone(jetTag = cms.InputTag("ak3CalopatJetsWithBtagging"),
genjetTag = 'ak3GenJets',
rParam = 0.3,
matchJets = cms.untracked.bool(False),
matchTag = 'patJetsWithBtagging',
pfCandidateLabel = cms.untracked.InputTag('particleFlow'),
trackTag = cms.InputTag("generalTracks"),
fillGenJets = False,
isMC = False,
doSubEvent = False,
useHepMC = cms.untracked.bool(False),
genParticles = cms.untracked.InputTag("genParticles"),
eventInfoTag = cms.InputTag("generator"),
doLifeTimeTagging = cms.untracked.bool(True),
doLifeTimeTaggingExtras = cms.untracked.bool(False),
bTagJetName = cms.untracked.string("ak3Calo"),
jetName = cms.untracked.string("ak3Calo"),
genPtMin = cms.untracked.double(5),
hltTrgResults = cms.untracked.string('TriggerResults::'+'HISIGNAL'),
doTower = cms.untracked.bool(True),
doSubJets = cms.untracked.bool(False),
doGenSubJets = cms.untracked.bool(False),
subjetGenTag = cms.untracked.InputTag("ak3GenJets"),
doGenTaus = False
)

ak3CaloJetSequence_mc = cms.Sequence(
#ak3Caloclean
#*
ak3Calomatch
#*
#ak3CalomatchGroomed
*
ak3Caloparton
*
ak3Calocorr
*
#ak3CaloJetID
#*
ak3CaloPatJetFlavourIdLegacy
#*
#ak3CaloPatJetFlavourId # Use legacy algo till PU implemented
*
ak3CaloJetTracksAssociatorAtVertex
*
ak3CaloJetBtagging
*
ak3CaloNjettiness #No constituents for calo jets in pp. Must be removed for pp calo jets but I'm not sure how to do this transparently (Marta)
*
ak3CalopatJetsWithBtagging
*
ak3CaloJetAnalyzer
)

ak3CaloJetSequence_data = cms.Sequence(ak3Calocorr
*
#ak3CaloJetID
#*
ak3CaloJetTracksAssociatorAtVertex
*
ak3CaloJetBtagging
*
ak3CaloNjettiness
*
ak3CalopatJetsWithBtagging
*
ak3CaloJetAnalyzer
)

ak3CaloJetSequence_jec = cms.Sequence(ak3CaloJetSequence_mc)
ak3CaloJetSequence_mb = cms.Sequence(ak3CaloJetSequence_mc)

ak3CaloJetSequence = cms.Sequence(ak3CaloJetSequence_data)
Loading

0 comments on commit c9bd685

Please sign in to comment.