This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exovvrun2:miniaod [2015/02/05 19:35] – taarre | exovvrun2:miniaod [2015/02/05 19:39] (current) – [Jet reclustering] taarre | ||
---|---|---|---|
Line 16: | Line 16: | ||
Flavour is stored for slimmedJets via | Flavour is stored for slimmedJets via | ||
| | ||
- | , not for AK8. | + | not available |
- | ==== Jet reslustering | + | ==== Jet reclustering |
- | | + | In config file |
+ | |||
+ | < | ||
+ | ####### Redo Jet clustering sequence ########## | ||
+ | |||
+ | from RecoJets.Configuration.RecoPFJets_cff import ak8PFJetsCHS, | ||
+ | |||
+ | process.chs = cms.EDFilter(" | ||
+ | src = cms.InputTag(' | ||
+ | cut = cms.string(' | ||
+ | ) | ||
+ | |||
+ | process.ak8PFJetsCHS = ak8PFJetsCHS.clone( src = ' | ||
+ | process.ak8PFJetsCHSPruned = ak8PFJetsCHSPruned.clone( src = ' | ||
+ | process.ak8PFJetsCHSSoftDrop = ak8PFJetsCHSSoftDrop.clone( src = ' | ||
+ | # process.ak8PFJetsCSTrimmed = ak8PFJetsCSTrimmed.clone( src = ' | ||
+ | # process.ak8PFJetsCSFiltered = ak8PFJetsCSFiltered.clone( src = ' | ||
+ | process.ak8PFJetsCHSPrunedLinks = ak8PFJetsCHSPrunedLinks.clone() | ||
+ | process.ak8PFJetsCHSSoftDropLinks = ak8PFJetsCHSSoftDropLinks.clone() | ||
+ | # process.ak8PFJetsCHSFilteredLinks = ak8PFJetsCHSFilteredLinks.clone() | ||
+ | # process.ak8PFJetsCHSTrimmedLinks = ak8PFJetsCHSTrimmedLinks.clone() | ||
+ | process.NjettinessAK8 = cms.EDProducer(" | ||
+ | src = cms.InputTag(" | ||
+ | Njets = cms.vuint32(1, | ||
+ | # variables for measure definition : | ||
+ | | ||
+ | beta = cms.double(1.0), | ||
+ | R0 = cms.double( 0.8 ), # CMS default is jet cone size | ||
+ | | ||
+ | # variables for axes definition : | ||
+ | | ||
+ | nPass = cms.int32(-999), | ||
+ | | ||
+ | ) | ||
+ | |||
+ | |||
+ | process.substructureSequence = cms.Sequence(process.chs + | ||
+ | process.ak8PFJetsCHS + | ||
+ | process.ak8PFJetsCHSPruned + | ||
+ | process.ak8PFJetsCHSSoftDrop + | ||
+ | # process.ak8PFJetsCSTrimmed + | ||
+ | # process.ak8PFJetsCSFiltered + | ||
+ | process.NjettinessAK8 + | ||
+ | process.ak8PFJetsCHSPrunedLinks + | ||
+ | process.ak8PFJetsCHSSoftDropLinks) # + | ||
+ | # | ||
+ | # | ||
+ | |||
+ | ####### Redo pat jets sequence ########## | ||
+ | |||
+ | from ExoDiBosonResonances.EDBRJets.redoPatJets_cff import patJetCorrFactorsAK8, | ||
+ | |||
+ | # Redo ak8PFJetsCHS | ||
+ | process.patJetCorrFactorsAK8 = patJetCorrFactorsAK8.clone( src = ' | ||
+ | process.patJetsAK8 = patJetsAK8.clone( jetSource = ' | ||
+ | process.patJetsAK8.jetCorrFactorsSource = cms.VInputTag( cms.InputTag(" | ||
+ | process.selectedPatJetsAK8 = selectedPatJetsAK8.clone( cut = cms.string(' | ||
+ | |||
+ | process.redoPatJets = cms.Sequence( process.patJetCorrFactorsAK8 + process.patJetsAK8 + process.selectedPatJetsAK8 ) | ||
+ | |||
+ | # Redo ak8PFJetsCHSPruned | ||
+ | process.patPrunedJetCorrFactorsAK8 = patJetCorrFactorsAK8.clone( src = ' | ||
+ | process.patPrunedJetsAK8 = patJetsAK8.clone( jetSource = ' | ||
+ | process.patPrunedJetsAK8.jetCorrFactorsSource = cms.VInputTag( cms.InputTag(" | ||
+ | process.patPrunedJetsAK8.userData.userFloats =cms.PSet(src = cms.VInputTag("" | ||
+ | process.selectedPrunedPatJetsAK8 = selectedPatJetsAK8.clone(cut = 'pt > 20', src = " | ||
+ | |||
+ | process.redoPrunedPatJets = cms.Sequence( process.patPrunedJetCorrFactorsAK8 + process.patPrunedJetsAK8 + process.selectedPrunedPatJetsAK8 ) | ||
+ | |||
+ | # Redo ak8PFJetsCHSSoftDrop | ||
+ | process.patSoftDropJetCorrFactorsAK8 = patJetCorrFactorsAK8.clone( src = ' | ||
+ | process.patSoftDropJetsAK8 = patJetsAK8.clone( jetSource = ' | ||
+ | process.patSoftDropJetsAK8.jetCorrFactorsSource = cms.VInputTag( cms.InputTag(" | ||
+ | process.patSoftDropJetsAK8.userData.userFloats =cms.PSet(src = cms.VInputTag("" | ||
+ | process.selectedSoftDropPatJetsAK8 = selectedPatJetsAK8.clone(cut = 'pt > 20', src = " | ||
+ | |||
+ | process.redoSoftDropPatJets = cms.Sequence( process.patSoftDropJetCorrFactorsAK8 + process.patSoftDropJetsAK8 + process.selectedSoftDropPatJetsAK8 ) | ||
+ | </ | ||
+ | |||
+ | Remember to add to your path: | ||
+ | process.p = cms.Path(process.substructureSequence*process.redoPatJets*process.redoPrunedPatJets*process.redoSoftDropPatJets*process.ntuplizer) |