RecoJets/python/JetValidationHeavyIons_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 ### genjet cleaning for improved matching in HI environment
0004
0005 from RecoHI.HiJetAlgos.HiGenCleaner_cff import *
0006
0007 iterativeCone5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets'))
0008 #iterativeCone7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone7HiGenJets'))
0009 ak2HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = 'ak2HiGenJets')
0010 ak3HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = 'ak3HiGenJets')
0011 ak4HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = 'ak4HiGenJets')
0012 ak5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = 'ak5HiGenJets')
0013
0014 ak7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = 'ak7HiGenJets')
0015
0016 ### jet analyzer for several radii
0017 ### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU
0018
0019
0020 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0021 JetAnalyzerICPU5Calo = DQMEDAnalyzer('JetTester_HeavyIons',
0022                                       JetType = cms.untracked.string('calo'),
0023                                       UEAlgo = cms.untracked.string('Pu'),
0024                                       OutputFile = cms.untracked.string(''),
0025                                       src = cms.InputTag("iterativeConePu5CaloJets"),
0026                                       srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"),
0027                                       PFcands = cms.InputTag("particleFlowTmp"),
0028                                       Background = cms.InputTag("voronoiBackgroundCalo"),
0029                                       #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0030                                       centralitycollection = cms.InputTag("hiCentrality"),
0031                                       centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0032                                       JetCorrections = cms.string(""),
0033                                       recoJetPtThreshold = cms.double(10),
0034                                       genEnergyFractionThreshold = cms.double(0.05),
0035                                       genPtThreshold = cms.double(1.0),
0036                                       matchGenPtThreshold                 = cms.double(20.0),
0037                                       RThreshold = cms.double(0.3),
0038                                       reverseEnergyFractionThreshold = cms.double(0.5)
0039 )
0040 '''
0041 JetAnalyzerICPU7Calo = DQMEDAnalyzer('JetTester_HeavyIons',
0042                                       JetType = cms.untracked.string('calo'),
0043                                       UEAlgo = cms.untracked.string('Pu'),
0044                                       OutputFile = cms.untracked.string(''),
0045                                       src = cms.InputTag("iterativeConePu7CaloJets"),
0046                                       srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"),
0047                                       #srcRho = cms.InputTag("iterativeConePu7CaloJets","rho"),
0048                                       JetCorrections = cms.string(""),
0049                                       recoJetPtThreshold = cms.double(10),
0050                                       genEnergyFractionThreshold = cms.double(0.05),
0051                                       genPtThreshold = cms.double(1.0),
0052                                       matchGenPtThreshold                 = cms.double(20.0),
0053                                       RThreshold = cms.double(0.3),
0054                                       reverseEnergyFractionThreshold = cms.double(0.5)
0055 )
0056 '''
0057 JetAnalyzerAkPU3Calo = DQMEDAnalyzer('JetTester_HeavyIons',
0058                                       JetType = cms.untracked.string('calo'),
0059                                       UEAlgo = cms.untracked.string('Pu'),
0060                                       OutputFile = cms.untracked.string(''),
0061                                       src = cms.InputTag("akPu3CaloJets"),
0062                                       srcGen = cms.InputTag("ak3HiCleanedGenJets"),
0063                                       PFcands = cms.InputTag("particleFlowTmp"),
0064                                       Background = cms.InputTag("voronoiBackgroundCalo"),
0065                                       #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0066                                       centralitycollection = cms.InputTag("hiCentrality"),
0067                                       centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0068                                       JetCorrections = cms.string(""),
0069                                       recoJetPtThreshold = cms.double(10),
0070                                       genEnergyFractionThreshold = cms.double(0.05),
0071                                       genPtThreshold = cms.double(1.0),
0072                                       matchGenPtThreshold                 = cms.double(20.0),
0073                                       RThreshold = cms.double(0.3),
0074                                       reverseEnergyFractionThreshold = cms.double(0.5)
0075 )
0076
0077 JetAnalyzerAkPU4Calo = DQMEDAnalyzer('JetTester_HeavyIons',
0078                                       JetType = cms.untracked.string('calo'),
0079                                       UEAlgo = cms.untracked.string('Pu'),
0080                                       OutputFile = cms.untracked.string(''),
0081                                       src = cms.InputTag("akPu4CaloJets"),
0082                                       srcGen = cms.InputTag("ak4HiCleanedGenJets"),
0083                                       PFcands = cms.InputTag("particleFlowTmp"),
0084                                       Background = cms.InputTag("voronoiBackgroundCalo"),
0085                                       #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0086                                       centralitycollection = cms.InputTag("hiCentrality"),
0087                                       centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0088                                       JetCorrections = cms.string(""),
0089                                       recoJetPtThreshold = cms.double(10),
0090                                       genEnergyFractionThreshold = cms.double(0.05),
0091                                       genPtThreshold = cms.double(1.0),
0092                                       matchGenPtThreshold                 = cms.double(20.0),
0093                                       RThreshold = cms.double(0.3),
0094                                       reverseEnergyFractionThreshold = cms.double(0.5)
0095 )
0096
0097 JetAnalyzerAkPU5Calo = DQMEDAnalyzer('JetTester_HeavyIons',
0098                                       JetType = cms.untracked.string('calo'),
0099                                       UEAlgo = cms.untracked.string('Pu'),
0100                                       OutputFile = cms.untracked.string(''),
0101                                       src = cms.InputTag("akPu5CaloJets"),
0102                                       srcGen = cms.InputTag("ak5HiCleanedGenJets"),
0103                                       PFcands = cms.InputTag("particleFlowTmp"),
0104                                       Background = cms.InputTag("voronoiBackgroundCalo"),
0105                                       #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0106                                       centralitycollection = cms.InputTag("hiCentrality"),
0107                                       centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0108                                       JetCorrections = cms.string(""),
0109                                       recoJetPtThreshold = cms.double(10),
0110                                       genEnergyFractionThreshold = cms.double(0.05),
0111                                       genPtThreshold = cms.double(1.0),
0112                                       matchGenPtThreshold                 = cms.double(20.0),
0113                                       RThreshold = cms.double(0.3),
0114                                       reverseEnergyFractionThreshold = cms.double(0.5)
0115 )
0116
0117 JetAnalyzerAkPU3PF = DQMEDAnalyzer('JetTester_HeavyIons',
0118                                     JetType = cms.untracked.string('pf'),
0119                                     UEAlgo = cms.untracked.string('Pu'),
0120                                     OutputFile = cms.untracked.string(''),
0121                                     src = cms.InputTag("akPu3PFJets"),
0122                                     srcGen = cms.InputTag("ak3HiCleanedGenJets"),
0123                                     PFcands = cms.InputTag("particleFlowTmp"),
0124                                     Background = cms.InputTag("voronoiBackgroundPF"),
0125                                     #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0126                                     centralitycollection = cms.InputTag("hiCentrality"),
0127                                     centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0128                                     JetCorrections = cms.string(""),
0129                                     recoJetPtThreshold = cms.double(10),
0130                                     genEnergyFractionThreshold = cms.double(0.05),
0131                                     genPtThreshold = cms.double(1.0),
0132                                     matchGenPtThreshold                 = cms.double(20.0),
0133                                     RThreshold = cms.double(0.3),
0134                                     reverseEnergyFractionThreshold = cms.double(0.5)
0135 )
0136
0137 JetAnalyzerAkPU4PF = DQMEDAnalyzer('JetTester_HeavyIons',
0138                                     JetType = cms.untracked.string('pf'),
0139                                     UEAlgo = cms.untracked.string('Pu'),
0140                                     OutputFile = cms.untracked.string(''),
0141                                     src = cms.InputTag("akPu4PFJets"),
0142                                     srcGen = cms.InputTag("ak4HiCleanedGenJets"),
0143                                     PFcands = cms.InputTag("particleFlowTmp"),
0144                                     Background = cms.InputTag("voronoiBackgroundPF"),
0145                                     #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0146                                     centralitycollection = cms.InputTag("hiCentrality"),
0147                                     centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0148                                     JetCorrections = cms.string(""),
0149                                     recoJetPtThreshold = cms.double(10),
0150                                     genEnergyFractionThreshold = cms.double(0.05),
0151                                     genPtThreshold = cms.double(1.0),
0152                                     matchGenPtThreshold                 = cms.double(20.0),
0153                                     RThreshold = cms.double(0.3),
0154                                     reverseEnergyFractionThreshold = cms.double(0.5)
0155 )
0156
0157 JetAnalyzerAkPU5PF = DQMEDAnalyzer('JetTester_HeavyIons',
0158                                     JetType = cms.untracked.string('pf'),
0159                                     UEAlgo = cms.untracked.string('Pu'),
0160                                     OutputFile = cms.untracked.string(''),
0161                                     src = cms.InputTag("akPu5PFJets"),
0162                                     srcGen = cms.InputTag("ak5HiCleanedGenJets"),
0163                                     PFcands = cms.InputTag("particleFlowTmp"),
0164                                     Background = cms.InputTag("voronoiBackgroundPF"),
0165                                     #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0166                                     centralitycollection = cms.InputTag("hiCentrality"),
0167                                     centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0168                                     JetCorrections = cms.string(""),
0169                                     recoJetPtThreshold = cms.double(10),
0170                                     genEnergyFractionThreshold = cms.double(0.05),
0171                                     genPtThreshold = cms.double(1.0),
0172                                     matchGenPtThreshold                 = cms.double(20.0),
0173                                     RThreshold = cms.double(0.3),
0174                                     reverseEnergyFractionThreshold = cms.double(0.5)
0175 )
0176
0177
0178 hiJetValidation = cms.Sequence(
0179     ak3HiCleanedGenJets
0180     * ak4HiCleanedGenJets
0181     * ak5HiCleanedGenJets
0182     * JetAnalyzerAkPU3Calo
0183     * JetAnalyzerAkPU4Calo
0184     * JetAnalyzerAkPU5Calo
0185
0186     * JetAnalyzerAkPU3PF
0187     * JetAnalyzerAkPU4PF
0188     * JetAnalyzerAkPU5PF
0189
0190 )