CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
JetValidationHeavyIons_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 ### genjet cleaning for improved matching in HI environment
4 
6 
7 iterativeCone5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets'))
8 #iterativeCone7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone7HiGenJets'))
9 ak2HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak2HiGenJets'))
10 ak3HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets'))
11 ak4HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak4HiGenJets'))
12 ak5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets'))
13 
14 ak7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak7HiGenJets'))
15 
16 ### jet analyzer for several radii
17 ### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU
18 
19 
20 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
21 JetAnalyzerICPU5Calo = DQMEDAnalyzer('JetTester_HeavyIons',
22  JetType = cms.untracked.string('calo'),
23  UEAlgo = cms.untracked.string('Pu'),
24  OutputFile = cms.untracked.string(''),
25  src = cms.InputTag("iterativeConePu5CaloJets"),
26  srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"),
27  PFcands = cms.InputTag("particleFlowTmp"),
28  Background = cms.InputTag("voronoiBackgroundCalo"),
29  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
30  centralitycollection = cms.InputTag("hiCentrality"),
31  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
32  JetCorrections = cms.string(""),
33  recoJetPtThreshold = cms.double(10),
34  genEnergyFractionThreshold = cms.double(0.05),
35  genPtThreshold = cms.double(1.0),
36  matchGenPtThreshold = cms.double(20.0),
37  RThreshold = cms.double(0.3),
38  reverseEnergyFractionThreshold = cms.double(0.5)
39 )
40 '''
41 JetAnalyzerICPU7Calo = DQMEDAnalyzer('JetTester_HeavyIons',
42  JetType = cms.untracked.string('calo'),
43  UEAlgo = cms.untracked.string('Pu'),
44  OutputFile = cms.untracked.string(''),
45  src = cms.InputTag("iterativeConePu7CaloJets"),
46  srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"),
47  #srcRho = cms.InputTag("iterativeConePu7CaloJets","rho"),
48  JetCorrections = cms.string(""),
49  recoJetPtThreshold = cms.double(10),
50  genEnergyFractionThreshold = cms.double(0.05),
51  genPtThreshold = cms.double(1.0),
52  matchGenPtThreshold = cms.double(20.0),
53  RThreshold = cms.double(0.3),
54  reverseEnergyFractionThreshold = cms.double(0.5)
55 )
56 '''
57 JetAnalyzerAkPU3Calo = DQMEDAnalyzer('JetTester_HeavyIons',
58  JetType = cms.untracked.string('calo'),
59  UEAlgo = cms.untracked.string('Pu'),
60  OutputFile = cms.untracked.string(''),
61  src = cms.InputTag("akPu3CaloJets"),
62  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
63  PFcands = cms.InputTag("particleFlowTmp"),
64  Background = cms.InputTag("voronoiBackgroundCalo"),
65  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
66  centralitycollection = cms.InputTag("hiCentrality"),
67  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
68  JetCorrections = cms.string(""),
69  recoJetPtThreshold = cms.double(10),
70  genEnergyFractionThreshold = cms.double(0.05),
71  genPtThreshold = cms.double(1.0),
72  matchGenPtThreshold = cms.double(20.0),
73  RThreshold = cms.double(0.3),
74  reverseEnergyFractionThreshold = cms.double(0.5)
75 )
76 
77 JetAnalyzerAkPU4Calo = DQMEDAnalyzer('JetTester_HeavyIons',
78  JetType = cms.untracked.string('calo'),
79  UEAlgo = cms.untracked.string('Pu'),
80  OutputFile = cms.untracked.string(''),
81  src = cms.InputTag("akPu4CaloJets"),
82  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
83  PFcands = cms.InputTag("particleFlowTmp"),
84  Background = cms.InputTag("voronoiBackgroundCalo"),
85  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
86  centralitycollection = cms.InputTag("hiCentrality"),
87  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
88  JetCorrections = cms.string(""),
89  recoJetPtThreshold = cms.double(10),
90  genEnergyFractionThreshold = cms.double(0.05),
91  genPtThreshold = cms.double(1.0),
92  matchGenPtThreshold = cms.double(20.0),
93  RThreshold = cms.double(0.3),
94  reverseEnergyFractionThreshold = cms.double(0.5)
95 )
96 
97 JetAnalyzerAkPU5Calo = DQMEDAnalyzer('JetTester_HeavyIons',
98  JetType = cms.untracked.string('calo'),
99  UEAlgo = cms.untracked.string('Pu'),
100  OutputFile = cms.untracked.string(''),
101  src = cms.InputTag("akPu5CaloJets"),
102  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
103  PFcands = cms.InputTag("particleFlowTmp"),
104  Background = cms.InputTag("voronoiBackgroundCalo"),
105  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
106  centralitycollection = cms.InputTag("hiCentrality"),
107  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
108  JetCorrections = cms.string(""),
109  recoJetPtThreshold = cms.double(10),
110  genEnergyFractionThreshold = cms.double(0.05),
111  genPtThreshold = cms.double(1.0),
112  matchGenPtThreshold = cms.double(20.0),
113  RThreshold = cms.double(0.3),
114  reverseEnergyFractionThreshold = cms.double(0.5)
115 )
116 
117 JetAnalyzerAkPU3PF = DQMEDAnalyzer('JetTester_HeavyIons',
118  JetType = cms.untracked.string('pf'),
119  UEAlgo = cms.untracked.string('Pu'),
120  OutputFile = cms.untracked.string(''),
121  src = cms.InputTag("akPu3PFJets"),
122  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
123  PFcands = cms.InputTag("particleFlowTmp"),
124  Background = cms.InputTag("voronoiBackgroundPF"),
125  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
126  centralitycollection = cms.InputTag("hiCentrality"),
127  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
128  JetCorrections = cms.string(""),
129  recoJetPtThreshold = cms.double(10),
130  genEnergyFractionThreshold = cms.double(0.05),
131  genPtThreshold = cms.double(1.0),
132  matchGenPtThreshold = cms.double(20.0),
133  RThreshold = cms.double(0.3),
134  reverseEnergyFractionThreshold = cms.double(0.5)
135 )
136 
137 JetAnalyzerAkPU4PF = DQMEDAnalyzer('JetTester_HeavyIons',
138  JetType = cms.untracked.string('pf'),
139  UEAlgo = cms.untracked.string('Pu'),
140  OutputFile = cms.untracked.string(''),
141  src = cms.InputTag("akPu4PFJets"),
142  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
143  PFcands = cms.InputTag("particleFlowTmp"),
144  Background = cms.InputTag("voronoiBackgroundPF"),
145  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
146  centralitycollection = cms.InputTag("hiCentrality"),
147  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
148  JetCorrections = cms.string(""),
149  recoJetPtThreshold = cms.double(10),
150  genEnergyFractionThreshold = cms.double(0.05),
151  genPtThreshold = cms.double(1.0),
152  matchGenPtThreshold = cms.double(20.0),
153  RThreshold = cms.double(0.3),
154  reverseEnergyFractionThreshold = cms.double(0.5)
155 )
156 
157 JetAnalyzerAkPU5PF = DQMEDAnalyzer('JetTester_HeavyIons',
158  JetType = cms.untracked.string('pf'),
159  UEAlgo = cms.untracked.string('Pu'),
160  OutputFile = cms.untracked.string(''),
161  src = cms.InputTag("akPu5PFJets"),
162  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
163  PFcands = cms.InputTag("particleFlowTmp"),
164  Background = cms.InputTag("voronoiBackgroundPF"),
165  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
166  centralitycollection = cms.InputTag("hiCentrality"),
167  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
168  JetCorrections = cms.string(""),
169  recoJetPtThreshold = cms.double(10),
170  genEnergyFractionThreshold = cms.double(0.05),
171  genPtThreshold = cms.double(1.0),
172  matchGenPtThreshold = cms.double(20.0),
173  RThreshold = cms.double(0.3),
174  reverseEnergyFractionThreshold = cms.double(0.5)
175 )
176 
177 
178 hiJetValidation = cms.Sequence(
179  ak3HiCleanedGenJets
180  * ak4HiCleanedGenJets
181  * ak5HiCleanedGenJets
182  * JetAnalyzerAkPU3Calo
183  * JetAnalyzerAkPU4Calo
184  * JetAnalyzerAkPU5Calo
185 
186  * JetAnalyzerAkPU3PF
187  * JetAnalyzerAkPU4PF
188  * JetAnalyzerAkPU5PF
189 
190 )
genjet cleaning for improved matching in HI environment