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:28] – taarre | exovvrun2:miniaod [2015/02/05 19:39] (current) – [Jet reclustering] taarre | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== miniAOD Workbook ====== | ====== miniAOD Workbook ====== | ||
- | First of all, read through the Twiki [[https:// | + | First of all, read through the Twiki [[https:// |
===== Jets ===== | ===== Jets ===== | ||
+ | |||
There are two jet collections stored. They are | There are two jet collections stored. They are | ||
- | * slimmedJets | + | * "slimmedJets" |
- | * slimmedJetsAK8 = AK8PFJetsCHS | + | * "slimmedJetsAK8" |
+ | |||
+ | For slimmedJets, | ||
+ | | ||
+ | The AK8 collection contains userfloats for the jets pruned, trimmed and filtered mass. Access via | ||
+ | | ||
+ | | ||
+ | | ||
+ | Flavour is stored for slimmedJets via | ||
+ | | ||
+ | not available for AK8. | ||
+ | |||
+ | ==== 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) |