Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ 07b03c7

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 07b03c7 was 93b4a8c, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

removed leptons from isolation
  • Property mode set to 100644
File size: 25.2 KB
RevLine 
[5348ed9]1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Heather Gray (CERN)
5# Filip Moortgat (CERN)
6# Michele Selvaggi (CERN)
7#
8# Released on: Oct 02th, 2016
9#
10# Configuration: FCC-hh generic detector without dipole in the forward region, no pile-up
11#
12#######################################
13# Order of execution of various modules
14#######################################
15
16set ExecutionPath {
17 ParticlePropagator
18
19 ChargedHadronTrackingEfficiency
20 ElectronTrackingEfficiency
21 MuonTrackingEfficiency
22
23 ChargedHadronMomentumSmearing
24 ElectronMomentumSmearing
25 MuonMomentumSmearing
26
27 TrackMerger
28
29 ECal
30 HCal
31
32 Calorimeter
33 EFlowMerger
[93b4a8c]34 EFlowFilter
35
[5348ed9]36 PhotonEfficiency
37 PhotonIsolation
38
39 ElectronFilter
40 ElectronIsolation
41
[bdcf74c]42 ChargedHadronFilter
43
[5348ed9]44 MuonIsolation
45
46 NeutrinoFilter
47
48 MissingET
49 GenMissingET
50
51 GenJetFinder
52 FastJetFinder
53
54 JetEnergyScale
55
56 JetFlavorAssociation
57
58 BTagging
59 CTagging
60 TauTagging
61
62 ScalarHT
63
64 UniqueObjectFinder
65
66 TreeWriter
67}
68
69#################################
70# Propagate particles in cylinder
71#################################
72
73module ParticlePropagator ParticlePropagator {
74 set InputArray Delphes/stableParticles
75
76 set OutputArray stableParticles
77 set ChargedHadronOutputArray chargedHadrons
78 set ElectronOutputArray electrons
79 set MuonOutputArray muons
80
81 # radius of the magnetic field coverage, in m
82 set Radius 1.5
83 # half-length of the magnetic field coverage, in m
84 set HalfLength 5
85
86 # magnetic field
87 set Bz 4.0
88}
89
90####################################
91# Charged hadron tracking efficiency
92####################################
93
94module Efficiency ChargedHadronTrackingEfficiency {
95 set InputArray ParticlePropagator/chargedHadrons
96 set OutputArray chargedHadrons
97
98 # TBC (which eta_max ? which pT min?)
99
100 # tracking efficiency formula for charged hadrons
101
102 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
[82b7a80]103(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
104(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
105(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
106(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
[5348ed9]107(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
108(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
109(abs(eta) > 6.0) * (0.00)}
110
111}
112
113##############################
114# Electron tracking efficiency
115##############################
116
117module Efficiency ElectronTrackingEfficiency {
118 set InputArray ParticlePropagator/electrons
119 set OutputArray electrons
120
121# TBC (which eta_max ?)
122# putting same as charged hadrons for now...
123
124 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
[82b7a80]125 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
126 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
127 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
128 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
[5348ed9]129 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
130 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
131 (abs(eta) > 6.0) * (0.00)}
132
133}
134##########################
135# Muon tracking efficiency
136##########################
137
138module Efficiency MuonTrackingEfficiency {
139 set InputArray ParticlePropagator/muons
140 set OutputArray muons
141
142# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
143# what about high pT ?
144 # tracking efficiency formula for muons
145 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
146 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
147 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
148 (abs(eta) > 6.0) * (0.00)}
149
150}
151
152########################################
153# Momentum resolution for charged tracks
154########################################
155
156module MomentumSmearing ChargedHadronMomentumSmearing {
157 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
158 set OutputArray chargedHadrons
159
160 source momentumResolutionVsP.tcl
161}
162
163
164###################################
165# Momentum resolution for electrons
166###################################
167
168module MomentumSmearing ElectronMomentumSmearing {
169 set InputArray ElectronTrackingEfficiency/electrons
170 set OutputArray electrons
171
172 source momentumResolutionVsP.tcl
173}
174
175
176###############################
177# Momentum resolution for muons
178###############################
179
180module MomentumSmearing MuonMomentumSmearing {
181 set InputArray MuonTrackingEfficiency/muons
182 set OutputArray muons
183
184 # TBC for just putting tracker resolution/ need to add improvement at high pT
185
186 source muonMomentumResolutionVsP.tcl
187}
188
189##############
190# Track merger
191##############
192
193module Merger TrackMerger {
194# add InputArray InputArray
195 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
196 add InputArray ElectronMomentumSmearing/electrons
197 add InputArray MuonMomentumSmearing/muons
198 set OutputArray tracks
199}
200
201#############
202# ECAL
203#############
204
205# TBC : calos seems ok, check eta max value though.
206
207module SimpleCalorimeter ECal {
208 set ParticleInputArray ParticlePropagator/stableParticles
209 set TrackInputArray TrackMerger/tracks
210
211 set TowerOutputArray ecalTowers
212 set EFlowTrackOutputArray eflowTracks
213 set EFlowTowerOutputArray eflowPhotons
214
215 set IsEcal true
216
217 set EnergyMin 0.5
218 set EnergySignificanceMin 2.0
219
220 set SmearTowerCenter true
221
222 set pi [expr {acos(-1)}]
223
224 # lists of the edges of each tower in eta and phi
225 # each list starts with the lower edge of the first tower
226 # the list ends with the higher edged of the last tower
227
228 # 0.5 degree towers (5x5 mm^2)
229 set PhiBins {}
230 for {set i -256} {$i <= 256} {incr i} {
231 add PhiBins [expr {$i * $pi/256.0}]
232 }
233
234 # TBC
235 foreach eta {-2.4875 -2.475 -2.4625 -2.45 -2.4375 -2.425 -2.4125 -2.4 -2.3875 -2.375 -2.3625 -2.35 -2.3375 -2.325 -2.3125 -2.3 -2.2875 -2.275 -2.2625 -2.25 -2.2375 -2.225 -2.2125 -2.2 -2.1875 -2.175 -2.1625 -2.15 -2.1375 -2.125 -2.1125 -2.1 -2.0875 -2.075 -2.0625 -2.05 -2.0375 -2.025 -2.0125 -2 -1.9875 -1.975 -1.9625 -1.95 -1.9375 -1.925 -1.9125 -1.9 -1.8875 -1.875 -1.8625 -1.85 -1.8375 -1.825 -1.8125 -1.8 -1.7875 -1.775 -1.7625 -1.75 -1.7375 -1.725 -1.7125 -1.7 -1.6875 -1.675 -1.6625 -1.65 -1.6375 -1.625 -1.6125 -1.6 -1.5875 -1.575 -1.5625 -1.55 -1.5375 -1.525 -1.5125 -1.5 -1.4875 -1.475 -1.4625 -1.45 -1.4375 -1.425 -1.4125 -1.4 -1.3875 -1.375 -1.3625 -1.35 -1.3375 -1.325 -1.3125 -1.3 -1.2875 -1.275 -1.2625 -1.25 -1.2375 -1.225 -1.2125 -1.2 -1.1875 -1.175 -1.1625 -1.15 -1.1375 -1.125 -1.1125 -1.1 -1.0875 -1.075 -1.0625 -1.05 -1.0375 -1.025 -1.0125 -1.0 -0.9875 -0.975 -0.9625 -0.95 -0.9375 -0.925 -0.9125 -0.9 -0.8875 -0.875 -0.8625 -0.85 -0.8375 -0.825 -0.8125 -0.8 -0.7875 -0.775 -0.7625 -0.75 -0.7375 -0.725 -0.7125 -0.7 -0.6875 -0.675 -0.6625 -0.65 -0.6375 -0.625 -0.6125 -0.6 -0.5875 -0.575 -0.5625 -0.55 -0.5375 -0.525 -0.5125 -0.5 -0.4875 -0.475 -0.4625 -0.45 -0.4375 -0.425 -0.4125 -0.4 -0.3875 -0.375 -0.3625 -0.35 -0.3375 -0.325 -0.3125 -0.3 -0.2875 -0.275 -0.2625 -0.25 -0.2375 -0.225 -0.2125 -0.2 -0.1875 -0.175 -0.1625 -0.15 -0.1375 -0.125 -0.1125 -0.1 -0.0875 -0.075 -0.0625 -0.05 -0.0375 -0.025 -0.0125 0 0.0125 0.025 0.0375 0.05 0.0625 0.075 0.0875 0.1 0.1125 0.125 0.1375 0.15 0.1625 0.175 0.1875 0.2 0.2125 0.225 0.2375 0.25 0.2625 0.275 0.2875 0.3 0.3125 0.325 0.3375 0.35 0.3625 0.375 0.3875 0.4 0.4125 0.425 0.4375 0.45 0.4625 0.475 0.4875 0.5 0.5125 0.525 0.5375 0.55 0.5625 0.575 0.5875 0.6 0.6125 0.625 0.6375 0.65 0.6625 0.675 0.6875 0.7 0.7125 0.725 0.7375 0.75 0.7625 0.775 0.7875 0.8 0.8125 0.825 0.8375 0.85 0.8625 0.875 0.8875 0.9 0.9125 0.925 0.9375 0.95 0.9625 0.975 0.9875 1.0 1.0125 1.025 1.0375 1.05 1.0625 1.075 1.0875 1.1 1.1125 1.125 1.1375 1.15 1.1625 1.175 1.1875 1.2 1.2125 1.225 1.2375 1.25 1.2625 1.275 1.2875 1.3 1.3125 1.325 1.3375 1.35 1.3625 1.375 1.3875 1.4 1.4125 1.425 1.4375 1.45 1.4625 1.475 1.4875 1.5 1.5125 1.525 1.5375 1.55 1.5625 1.575 1.5875 1.6 1.6125 1.625 1.6375 1.65 1.6625 1.675 1.6875 1.7 1.7125 1.725 1.7375 1.75 1.7625 1.775 1.7875 1.8 1.8125 1.825 1.8375 1.85 1.8625 1.875 1.8875 1.9 1.9125 1.925 1.9375 1.95 1.9625 1.975 1.9875 2 2.0125 2.025 2.0375 2.05 2.0625 2.075 2.0875 2.1 2.1125 2.125 2.1375 2.15 2.1625 2.175 2.1875 2.2 2.2125 2.225 2.2375 2.25 2.2625 2.275 2.2875 2.3 2.3125 2.325 2.3375 2.35 2.3625 2.375 2.3875 2.4 2.4125 2.425 2.4375 2.45 2.4625 2.475 2.4875 2.5} {
236 add EtaPhiBins $eta $PhiBins
237 }
238 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
239 set PhiBins {}
240 for {set i -128} {$i <= 128} {incr i} {
241 add PhiBins [expr {$i * $pi/128.0}]
242 }
243 foreach eta {-3.975 -3.95 -3.925 -3.9 -3.875 -3.85 -3.825 -3.8 -3.775 -3.75 -3.725 -3.7 -3.675 -3.65 -3.625 -3.6 -3.575 -3.55 -3.525 -3.5 -3.475 -3.45 -3.425 -3.4 -3.375 -3.35 -3.325 -3.3 -3.275 -3.25 -3.225 -3.2 -3.175 -3.15 -3.125 -3.1 -3.075 -3.05 -3.025 -3.0 -2.975 -2.95 -2.925 -2.9 -2.875 -2.85 -2.825 -2.8 -2.775 -2.75 -2.725 -2.7 -2.675 -2.65 -2.625 -2.6 -2.575 -2.55 -2.525 -2.5 2.525 2.55 2.575 2.6 2.625 2.65 2.675 2.7 2.725 2.75 2.775 2.8 2.825 2.85 2.875 2.9 2.925 2.95 2.975 3.0 3.025 3.05 3.075 3.1 3.125 3.15 3.175 3.2 3.225 3.25 3.275 3.3 3.325 3.35 3.375 3.4 3.425 3.45 3.475 3.5 3.525 3.55 3.575 3.6 3.625 3.65 3.675 3.7 3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 3.925 3.95 3.975 4.0} {
244 add EtaPhiBins $eta $PhiBins
245 }
246
247 # 0.05 x 0.05 for eta between 4.0 and 6.0
248 set PhiBins {}
249 for {set i -64} {$i <= 64} {incr i} {
250 add PhiBins [expr {$i * $pi/64.0}]
251 }
252 foreach eta {-6.0 -5.95 -5.9 -5.85 -5.8 -5.75 -5.7 -5.65 -5.6 -5.55 -5.5 -5.45 -5.4 -5.35 -5.3 -5.25 -5.2 -5.15 -5.1 -5.05 -5.0 -4.95 -4.9 -4.85 -4.8 -4.75 -4.7 -4.65 -4.6 -4.55 -4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 -4 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 4.95 5.0 5.05 5.1 5.15 5.2 5.25 5.3 5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6.0} {
253 add EtaPhiBins $eta $PhiBins
254 }
255 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
256
257 add EnergyFraction {0} {0.0}
258 # energy fractions for e, gamma and pi0
259 add EnergyFraction {11} {1.0}
260 add EnergyFraction {22} {1.0}
261 add EnergyFraction {111} {1.0}
262 # energy fractions for muon, neutrinos and neutralinos
263 add EnergyFraction {12} {0.0}
264 add EnergyFraction {13} {0.0}
265 add EnergyFraction {14} {0.0}
266 add EnergyFraction {16} {0.0}
267 add EnergyFraction {1000022} {0.0}
268 add EnergyFraction {1000023} {0.0}
269 add EnergyFraction {1000025} {0.0}
270 add EnergyFraction {1000035} {0.0}
271 add EnergyFraction {1000045} {0.0}
272 # energy fractions for K0short and Lambda
273 # add EnergyFraction {310} {0.3}
274 # add EnergyFraction {3122} {0.3}
275
276 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
277 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
278 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
279
280
281}
282
283#############
284# HCAL
285#############
286
287module SimpleCalorimeter HCal {
288 set ParticleInputArray ParticlePropagator/stableParticles
289 set TrackInputArray ECal/eflowTracks
290
291 set TowerOutputArray hcalTowers
292 set EFlowTrackOutputArray eflowTracks
293 set EFlowTowerOutputArray eflowNeutralHadrons
294
295 set IsEcal false
296
297 set EnergyMin 1.0
298 set EnergySignificanceMin 2.0
299
300 set SmearTowerCenter true
301
302 set pi [expr {acos(-1)}]
303
304 # lists of the edges of each tower in eta and phi
305 # each list starts with the lower edge of the first tower
306 # the list ends with the higher edged of the last tower
307
308 # 6 degree towers
309 set PhiBins {}
310 for {set i -64} {$i <= 64} {incr i} {
311 add PhiBins [expr {$i * $pi/64.0}]
312 }
313
314 # TBC
315 foreach eta {-2.45 -2.4 -2.35 -2.3 -2.25 -2.2 -2.15 -2.1 -2.05 -2.0 -1.95 -1.9 -1.85 -1.8 -1.75 -1.7 -1.65 -1.6 -1.55 -1.5 -1.45 -1.4 -1.35 -1.3 -1.25 -1.2 -1.15 -1.1 -1.05 -1.0 -0.95 -0.9 -0.85 -0.8 -0.75 -0.7 -0.65 -0.6 -0.55 -0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5} {
316 add EtaPhiBins $eta $PhiBins
317 }
318
319 # 0.1 x 0.1 between eta 2.5 and 4.0
320 set PhiBins {}
321 for {set i -32} {$i <= 32} {incr i} {
322 add PhiBins [expr {$i * $pi/32.0}]
323 }
324 foreach eta {-3.9 -3.8 -3.7 -3.6 -3.5 -3.4 -3.3 -3.2 -3.1 -3.0 -2.9 -2.8 -2.7 -2.6 -2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0} {
325 add EtaPhiBins $eta $PhiBins
326 }
327
328 # 0.2 x 0.2 between eta 4.0 and 6.0
329 set PhiBins {}
330 for {set i -16} {$i <= 16} {incr i} {
331 add PhiBins [expr {$i * $pi/16.0}]
332 }
333 foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 -4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
334 add EtaPhiBins $eta $PhiBins
335 }
336
337 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
338 add EnergyFraction {0} {1.0}
339 # energy fractions for e, gamma and pi0
340 add EnergyFraction {11} {0.0}
341 add EnergyFraction {22} {0.0}
342 add EnergyFraction {111} {0.0}
343 # energy fractions for muon, neutrinos and neutralinos
344 add EnergyFraction {12} {0.0}
345 add EnergyFraction {13} {0.0}
346 add EnergyFraction {14} {0.0}
347 add EnergyFraction {16} {0.0}
348 add EnergyFraction {1000022} {0.0}
349 add EnergyFraction {1000023} {0.0}
350 add EnergyFraction {1000025} {0.0}
351 add EnergyFraction {1000035} {0.0}
352 add EnergyFraction {1000045} {0.0}
353 # energy fractions for K0short and Lambda
354 # add EnergyFraction {310} {0.7}
355 # add EnergyFraction {3122} {0.7}
356
357 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
358 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
359 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
360}
361
362#################
363# Electron filter
364#################
365
366module PdgCodeFilter ElectronFilter {
367 set InputArray HCal/eflowTracks
368 set OutputArray electrons
369 set Invert true
370 add PdgCode {11}
371 add PdgCode {-11}
372}
373
[bdcf74c]374
375######################
376# ChargedHadronFilter
377######################
378
379module PdgCodeFilter ChargedHadronFilter {
380 set InputArray HCal/eflowTracks
381 set OutputArray chargedHadrons
382
383 add PdgCode {11}
384 add PdgCode {-11}
385 add PdgCode {13}
386 add PdgCode {-13}
387}
388
389
[5348ed9]390###################################################
391# Tower Merger (in case not using e-flow algorithm)
392###################################################
393
394module Merger Calorimeter {
395# add InputArray InputArray
396 add InputArray ECal/ecalTowers
397 add InputArray HCal/hcalTowers
398 set OutputArray towers
399}
400
401####################
402# Energy flow merger
403####################
404
405module Merger EFlowMerger {
406# add InputArray InputArray
407 add InputArray HCal/eflowTracks
408 add InputArray ECal/eflowPhotons
409 add InputArray HCal/eflowNeutralHadrons
410 set OutputArray eflow
411}
412
[93b4a8c]413######################
414# EFlowFilter
415######################
416
417module PdgCodeFilter EFlowFilter {
418 set InputArray EFlowMerger/eflow
419 set OutputArray eflow
420
421 add PdgCode {11}
422 add PdgCode {-11}
423 add PdgCode {13}
424 add PdgCode {-13}
425}
426
[5348ed9]427
428###################
429# Missing ET merger
430###################
431
432module Merger MissingET {
433# add InputArray InputArray
434 add InputArray EFlowMerger/eflow
435 set MomentumOutputArray momentum
436}
437
438
439
440##################
441# Scalar HT merger
442##################
443
444module Merger ScalarHT {
445# add InputArray InputArray
446 add InputArray EFlowMerger/eflow
447 set EnergyOutputArray energy
448}
449
450#################
451# Neutrino Filter
452#################
453
454module PdgCodeFilter NeutrinoFilter {
455
456 set InputArray Delphes/stableParticles
457 set OutputArray filteredParticles
458
459 set PTMin 0.0
460
461 add PdgCode {12}
462 add PdgCode {14}
463 add PdgCode {16}
464 add PdgCode {-12}
465 add PdgCode {-14}
466 add PdgCode {-16}
467
468}
469
470
471#####################
472# MC truth jet finder
473#####################
474
475# TBC: is jet radius fine?
476
477module FastJetFinder GenJetFinder {
478# set InputArray NeutrinoFilter/filteredParticles
479 set InputArray Delphes/stableParticles
480
481 set OutputArray jets
482
483 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
484 set JetAlgorithm 6
485 set ParameterR 0.4
486
487 set JetPTMin 5.0
488}
489
490#########################
491# Gen Missing ET merger
492########################
493
494module Merger GenMissingET {
495
496# add InputArray InputArray
497 add InputArray NeutrinoFilter/filteredParticles
498 set MomentumOutputArray momentum
499}
500
501
502
503############
504# Jet finder
505############
506
507# TBC need to include jet substructure variables
508# TBC is jet radius fine?
509
510module FastJetFinder FastJetFinder {
511# set InputArray Calorimeter/towers
512 set InputArray EFlowMerger/eflow
513
514 set OutputArray jets
515
516 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
517 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
518
519 set JetAlgorithm 6
520 set ParameterR 0.4
521
522 set JetPTMin 30.0
523}
524
525##################
526# Jet Energy Scale
527##################
528
529module EnergyScale JetEnergyScale {
530 set InputArray FastJetFinder/jets
531 set OutputArray jets
532
533 # scale formula for jets
534 set ScaleFormula {1.00}
535}
536
537
538########################
539# Jet Flavor Association
540########################
541
542module JetFlavorAssociation JetFlavorAssociation {
543
544 set PartonInputArray Delphes/partons
545 set ParticleInputArray Delphes/allParticles
546 set ParticleLHEFInputArray Delphes/allParticlesLHEF
547 set JetInputArray JetEnergyScale/jets
548
549 set DeltaR 0.5
550 set PartonPTMin 5.0
551 set PartonEtaMax 6.0
552
553}
554
555###################
556# Photon efficiency
557###################
558
559module Efficiency PhotonEfficiency {
560 set InputArray ECal/eflowPhotons
561 set OutputArray photons
562
563 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
564
565 set EfficiencyFormula {
[2b1e147]566 (pt <= 1.0) * (0.00) + \
[2dfc52a]567 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
568 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
569 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
[2b1e147]570
[2dfc52a]571 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
572 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
573 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
[2b1e147]574
575 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
576 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
577 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
[5348ed9]578 (abs(eta) > 6.0) * (0.00)}
579
580}
581
582##################
583# Photon isolation
584##################
585
586# TBC: check values for iso cuts
587
588module Isolation PhotonIsolation {
589 set CandidateInputArray PhotonEfficiency/photons
[93b4a8c]590 set IsolationInputArray EFlowFilter/eflow
[5348ed9]591
592 set OutputArray photons
593
594 set DeltaRMax 0.3
595
596 set PTMin 0.5
597
598 set PTRatioMax 0.1
599}
600
601
602####################
603# Electron isolation
604####################
605
606# TBC: check values for iso cuts
607
608module Isolation ElectronIsolation {
609 set CandidateInputArray ElectronFilter/electrons
[93b4a8c]610 set IsolationInputArray EFlowFilter/eflow
[5348ed9]611
612 set OutputArray electrons
613
614 set DeltaRMax 0.3
615
616 set PTMin 0.5
617
618 set PTRatioMax 0.1
619}
620
621
622################
623# Muon isolation
624################
625
626# TBC: check values for iso cuts
627
628module Isolation MuonIsolation {
629 set CandidateInputArray MuonMomentumSmearing/muons
[93b4a8c]630 set IsolationInputArray EFlowFilter/eflow
[5348ed9]631
632 set OutputArray muons
633
634 set DeltaRMax 0.3
635
636 set PTMin 0.5
637
638 set PTRatioMax 0.2
639}
640
641
642###########
643# b-tagging
644###########
645
646module BTagging BTagging {
647 set JetInputArray JetEnergyScale/jets
648
649 set BitNumber 0
650
651 add EfficiencyFormula {0} {
652
[82b7a80]653 (pt <= 10.0) * (0.00) +
[44b481d]654 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
655 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
656 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
657 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
658 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
659 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
[5348ed9]660 (abs(eta) > 4.0) * (0.00)}
661
662 add EfficiencyFormula {4} {
663
[82b7a80]664 (pt <= 10.0) * (0.00) +
[44b481d]665 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
666 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
667 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
668 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
669 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
670 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
[5348ed9]671 (abs(eta) > 4.0) * (0.00)}
672
673 add EfficiencyFormula {5} {
674
[44b481d]675 (pt <= 10.0) * (0.00) +
676 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
677 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
678 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
679 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
680 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
681 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
682 (abs(eta) >= 4.0) * (0.00)}
[5348ed9]683
684}
685
686###########
687# c-tagging
688###########
689
690module BTagging CTagging {
691 set JetInputArray JetEnergyScale/jets
692
693 set BitNumber 1
694
695 add EfficiencyFormula {0} {
696
697 (pt <= 10.0) * (0.00) +
698 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
699 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
700
701 add EfficiencyFormula {4} {
702
703 (pt <= 10.0) * (0.00) +
704 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
705 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
706
707 add EfficiencyFormula {5} {
708
709 (pt <= 10.0) * (0.00) +
710 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
711 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
712
713}
714
715
716#############
717# tau-tagging
718#############
719
720
721module TauTagging TauTagging {
722 set ParticleInputArray Delphes/allParticles
723 set PartonInputArray Delphes/partons
724 set JetInputArray JetEnergyScale/jets
725
726 set DeltaR 0.5
727
728 set TauPTMin 1.0
729
730 set TauEtaMax 4.0
731
732 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
733 add EfficiencyFormula {0} {
734
[44b481d]735 (pt <= 10.0) * (0.00) +
736 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
737 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
738 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
739 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
740 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
741 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
[82b7a80]742 (abs(eta) > 4.0) * (0.00)}
[5348ed9]743
744 add EfficiencyFormula {11} {
745
[44b481d]746 (pt <= 10.0) * (0.00) +
747 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
748 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
749 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
750 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
751 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
752 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
[82b7a80]753 (abs(eta) > 4.0) * (0.00)}
[5348ed9]754
755 add EfficiencyFormula {15} {
756
[44b481d]757 (pt <= 10.0) * (0.00) +
758 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
759 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
760 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
761 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
762 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
763 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
764 (abs(eta) > 4.0) * (0.00)}
[5348ed9]765
766}
767
768#####################################################
769# Find uniquely identified photons/electrons/tau/jets
770#####################################################
771
772module UniqueObjectFinder UniqueObjectFinder {
773# earlier arrays take precedence over later ones
774# add InputArray InputArray OutputArray
775 add InputArray PhotonIsolation/photons photons
776 add InputArray ElectronIsolation/electrons electrons
777 add InputArray MuonIsolation/muons muons
778 add InputArray JetEnergyScale/jets jets
779
780}
781
782
783##################
784# ROOT tree writer
785##################
786
787module TreeWriter TreeWriter {
788# add Branch InputArray BranchName BranchClass
789 add Branch Delphes/allParticles Particle GenParticle
790
791 add Branch GenJetFinder/jets GenJet Jet
792 add Branch GenMissingET/momentum GenMissingET MissingET
793
794 add Branch TrackMerger/tracks Track Track
795 add Branch Calorimeter/towers Tower Tower
796
797 add Branch HCal/eflowTracks EFlowTrack Track
798 add Branch ECal/eflowPhotons EFlowPhoton Tower
799 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
800
801 add Branch UniqueObjectFinder/photons Photon Photon
802 add Branch UniqueObjectFinder/electrons Electron Electron
803 add Branch UniqueObjectFinder/muons Muon Muon
804 add Branch UniqueObjectFinder/jets Jet Jet
805
806 add Branch MissingET/momentum MissingET MissingET
807 add Branch ScalarHT/energy ScalarHT ScalarHT
808}
809
Note: See TracBrowser for help on using the repository browser.