Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add PUPPI JEC configurations and fix MiniAOD to use PUPPI JEC #13175

Merged
merged 1 commit into from Feb 4, 2016
Merged

Add PUPPI JEC configurations and fix MiniAOD to use PUPPI JEC #13175

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
50 changes: 50 additions & 0 deletions JetMETCorrections/Configuration/python/JetCorrectors_cff.py
View file
Open in desktop
Expand Up @@ -16,6 +16,7 @@
ak4PFL1OffsetCorrector = ak4CaloL1OffsetCorrector.clone(algorithm = 'AK4PF')
ak4PFCHSL1OffsetCorrector = ak4CaloL1OffsetCorrector.clone(algorithm = 'AK4PFchs')
ak4JPTL1OffsetCorrector = ak4CaloL1OffsetCorrector.clone(algorithm = 'AK4JPT')
ak4PFPuppiL1OffsetCorrector = ak4CaloL1OffsetCorrector.clone(algorithm = 'AK4PFPuppi')

# L1 (JPT Offset) CORRECTOR
ak4L1JPTOffsetCorrector = cms.EDProducer(
Expand Down Expand Up @@ -48,6 +49,12 @@
srcRho = cms.InputTag( 'fixedGridRhoFastjetAll' )
)
ak4JPTL1FastjetCorrector = ak4CaloL1FastjetCorrector.clone()
ak4PFPuppiL1FastjetCorrector = cms.EDProducer(
'L1FastjetCorrectorProducer',
level = cms.string('L1FastJet'),
algorithm = cms.string('AK4PFPuppi'),
srcRho = cms.InputTag( 'fixedGridRhoFastjetAll' )
)

# L2 (relative eta-conformity) CORRECTORS
ak4CaloL2RelativeCorrector = cms.EDProducer(
Expand All @@ -59,6 +66,7 @@
ak4PFCHSL2RelativeCorrector = ak4CaloL2RelativeCorrector.clone( algorithm = 'AK4PFchs' )
ak4JPTL2RelativeCorrector = ak4CaloL2RelativeCorrector.clone( algorithm = 'AK4JPT' )
ak4TrackL2RelativeCorrector = ak4CaloL2RelativeCorrector.clone( algorithm = 'AK4TRK' )
ak4PFPuppiL2RelativeCorrector = ak4CaloL2RelativeCorrector.clone( algorithm = 'AK4PFPuppi' )

# L3 (absolute) CORRECTORS
ak4CaloL3AbsoluteCorrector = cms.EDProducer(
Expand All @@ -70,6 +78,7 @@
ak4PFCHSL3AbsoluteCorrector = ak4CaloL3AbsoluteCorrector.clone( algorithm = 'AK4PFchs' )
ak4JPTL3AbsoluteCorrector = ak4CaloL3AbsoluteCorrector.clone( algorithm = 'AK4JPT' )
ak4TrackL3AbsoluteCorrector = ak4CaloL3AbsoluteCorrector.clone( algorithm = 'AK4TRK' )
ak4PFPuppiL3AbsoluteCorrector = ak4CaloL3AbsoluteCorrector.clone( algorithm = 'AK4PFPuppi' )

# Residual CORRECTORS
ak4CaloResidualCorrector = cms.EDProducer(
Expand All @@ -80,6 +89,7 @@
ak4PFResidualCorrector = ak4CaloResidualCorrector.clone( algorithm = 'AK4PF' )
ak4PFCHSResidualCorrector = ak4CaloResidualCorrector.clone( algorithm = 'AK4PFchs' )
ak4JPTResidualCorrector = ak4CaloResidualCorrector.clone( algorithm = 'AK4JPT' )
ak4PFPuppiResidualCorrector = ak4CaloResidualCorrector.clone( algorithm = 'AK4PFPuppi' )

# L6 (semileptonically decaying b-jet) Correction Services
ak4CaloL6SLBCorrector = cms.EDProducer(
Expand Down Expand Up @@ -142,6 +152,13 @@
ak4TrackL2L3CorrectorChain = cms.Sequence(
ak4TrackL2RelativeCorrector * ak4TrackL3AbsoluteCorrector * ak4TrackL2L3Corrector
)
ak4PFPuppiL2L3Corrector = cms.EDProducer(
'ChainedJetCorrectorProducer',
correctors = cms.VInputTag('ak4PFPuppiL2RelativeCorrector','ak4PFPuppiL3AbsoluteCorrector')
)
ak4PFPuppiL2L3CorrectorChain = cms.Sequence(
ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiL2L3Corrector
)

# L2L3Residual CORRECTORS
ak4CaloL2L3ResidualCorrector = cms.EDProducer(
Expand Down Expand Up @@ -174,6 +191,13 @@
ak4JPTL2L3ResidualCorrectorChain = cms.Sequence(
ak4L1JPTOffsetCorrectorChain * ak4JPTL2RelativeCorrector * ak4JPTL3AbsoluteCorrector * ak4JPTResidualCorrector * ak4JPTL2L3ResidualCorrector
)
ak4PFPuppiL2L3ResidualCorrector = cms.EDProducer(
'ChainedJetCorrectorProducer',
correctors = cms.VInputTag('ak4PFPuppiL2RelativeCorrector','ak4PFPuppiL3AbsoluteCorrector','ak4PFPuppiResidualCorrector')
)
ak4PFPuppiL2L3ResidualCorrectorChain = cms.Sequence(
ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiResidualCorrector * ak4PFPuppiL2L3ResidualCorrector
)

# L1L2L3 CORRECTORS
ak4CaloL1L2L3Corrector = cms.EDProducer(
Expand Down Expand Up @@ -206,6 +230,13 @@
ak4JPTL1L2L3CorrectorChain = cms.Sequence(
ak4L1JPTOffsetCorrectorChain * ak4JPTL2RelativeCorrector * ak4JPTL3AbsoluteCorrector * ak4JPTL1L2L3Corrector
)
ak4PFPuppiL1L2L3Corrector = cms.EDProducer(
'ChainedJetCorrectorProducer',
correctors = cms.VInputTag('ak4PFPuppiL1OffsetCorrector','ak4PFPuppiL2RelativeCorrector','ak4PFPuppiL3AbsoluteCorrector')
)
ak4PFPuppiL1L2L3CorrectorChain = cms.Sequence(
ak4PFPuppiL1OffsetCorrector * ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiL1L2L3Corrector
)

# L1L2L3Residual CORRECTORS
ak4CaloL1L2L3ResidualCorrector = cms.EDProducer(
Expand Down Expand Up @@ -238,6 +269,13 @@
ak4JPTL1L2L3ResidualCorrectorChain = cms.Sequence(
ak4L1JPTOffsetCorrectorChain * ak4JPTL2RelativeCorrector * ak4JPTL3AbsoluteCorrector * ak4JPTResidualCorrector * ak4JPTL1L2L3ResidualCorrector
)
ak4PFPuppiL1L2L3ResidualCorrector = cms.EDProducer(
'ChainedJetCorrectorProducer',
correctors = cms.VInputTag('ak4PFPuppiL1OffsetCorrector','ak4PFPuppiL2RelativeCorrector','ak4PFPuppiL3AbsoluteCorrector','ak4PFPuppiResidualCorrector')
)
ak4PFPuppiL1L2L3ResidualCorrectorChain = cms.Sequence(
ak4PFPuppiL1OffsetCorrector * ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiResidualCorrector * ak4PFPuppiL1L2L3ResidualCorrector
)

# L1L2L3 CORRECTORS WITH FASTJET
ak4CaloL1FastL2L3Corrector = ak4CaloL2L3Corrector.clone()
Expand All @@ -264,6 +302,11 @@
ak4JPTL1FastL2L3CorrectorChain = cms.Sequence(
ak4JPTL1FastjetCorrector * ak4JPTL2RelativeCorrector * ak4JPTL3AbsoluteCorrector * ak4JPTL1FastL2L3Corrector
)
ak4PFPuppiL1FastL2L3Corrector = ak4PFPuppiL2L3Corrector.clone()
ak4PFPuppiL1FastL2L3Corrector.correctors.insert(0,'ak4PFPuppiL1FastjetCorrector')
ak4PFPuppiL1FastL2L3CorrectorChain = cms.Sequence(
ak4PFPuppiL1FastjetCorrector * ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiL1FastL2L3Corrector
)

# L1L2L3Residual CORRECTORS WITH FASTJET
ak4CaloL1FastL2L3ResidualCorrector = cms.EDProducer(
Expand Down Expand Up @@ -296,6 +339,13 @@
ak4JPTL1FastL2L3ResidualCorrectorChain = cms.Sequence(
ak4JPTL1FastjetCorrector * ak4JPTL2RelativeCorrector * ak4JPTL3AbsoluteCorrector * ak4JPTResidualCorrector * ak4JPTL1FastL2L3ResidualCorrector
)
ak4PFPuppiL1FastL2L3ResidualCorrector = cms.EDProducer(
'ChainedJetCorrectorProducer',
correctors = cms.VInputTag('ak4PFPuppiL1FastjetCorrector','ak4PFPuppiL2RelativeCorrector','ak4PFPuppiL3AbsoluteCorrector','ak4PFPuppiResidualCorrector')
)
ak4PFPuppiL1FastL2L3ResidualCorrectorChain = cms.Sequence(
ak4PFPuppiL1FastjetCorrector * ak4PFPuppiL2RelativeCorrector * ak4PFPuppiL3AbsoluteCorrector * ak4PFPuppiResidualCorrector * ak4PFPuppiL1FastL2L3ResidualCorrector
)

# L2L3L6 CORRECTORS
ak4CaloL2L3L6Corrector = ak4CaloL2L3Corrector.clone()
Expand Down
4 changes: 2 additions & 2 deletions PhysicsTools/PatAlgos/python/slimming/miniAOD_tools.py
View file
Open in desktop
Expand Up @@ -237,7 +237,7 @@ def miniAOD_customizeCommon(process):
)

addJetCollection(process, postfix = "", labelName = 'Puppi', jetSource = cms.InputTag('ak4PFJetsPuppi'),
jetCorrections = ('AK4PFchs', ['L2Relative', 'L3Absolute'], ''),
jetCorrections = ('AK4PFPuppi', ['L2Relative', 'L3Absolute'], ''),
algo= 'AK', rParam = 0.4, btagDiscriminators = map(lambda x: x.value() ,process.patJets.discriminatorSources)
)

Expand All @@ -264,7 +264,7 @@ def miniAOD_customizeCommon(process):
# type1 correction, from puppi jets
process.corrPfMetType1Puppi = process.corrPfMetType1.clone(
src = 'ak4PFJetsPuppi',
jetCorrLabel = 'ak4PFCHSL2L3Corrector',
jetCorrLabel = 'ak4PFPuppiL2L3Corrector',
)
del process.corrPfMetType1Puppi.offsetCorrLabel # no L1 for PUPPI jets
process.pfMetT1Puppi = process.pfMetT1.clone(
Expand Down