Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ c6667c0

ImprovedOutputFile Timing dual_readout llp
Last change on this file since c6667c0 was 44b481d, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

fixed b tag and tau tag formulae
  • Property mode set to 100644
File size: 24.9 KB
Line 
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
34
35 PhotonEfficiency
36 PhotonIsolation
37
38 ElectronFilter
39 ElectronIsolation
40
41 ChargedHadronFilter
42
43 MuonIsolation
44
45 NeutrinoFilter
46
47 MissingET
48 GenMissingET
49
50 GenJetFinder
51 FastJetFinder
52
53 JetEnergyScale
54
55 JetFlavorAssociation
56
57 BTagging
58 CTagging
59 TauTagging
60
61 ScalarHT
62
63 UniqueObjectFinder
64
65 TreeWriter
66}
67
68#################################
69# Propagate particles in cylinder
70#################################
71
72module ParticlePropagator ParticlePropagator {
73 set InputArray Delphes/stableParticles
74
75 set OutputArray stableParticles
76 set ChargedHadronOutputArray chargedHadrons
77 set ElectronOutputArray electrons
78 set MuonOutputArray muons
79
80 # radius of the magnetic field coverage, in m
81 set Radius 1.5
82 # half-length of the magnetic field coverage, in m
83 set HalfLength 5
84
85 # magnetic field
86 set Bz 4.0
87}
88
89####################################
90# Charged hadron tracking efficiency
91####################################
92
93module Efficiency ChargedHadronTrackingEfficiency {
94 set InputArray ParticlePropagator/chargedHadrons
95 set OutputArray chargedHadrons
96
97 # TBC (which eta_max ? which pT min?)
98
99 # tracking efficiency formula for charged hadrons
100
101 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
102(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
103(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
104(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
105(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
106(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
107(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
108(abs(eta) > 6.0) * (0.00)}
109
110}
111
112##############################
113# Electron tracking efficiency
114##############################
115
116module Efficiency ElectronTrackingEfficiency {
117 set InputArray ParticlePropagator/electrons
118 set OutputArray electrons
119
120# TBC (which eta_max ?)
121# putting same as charged hadrons for now...
122
123 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
124 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
125 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
126 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
127 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
128 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
129 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
130 (abs(eta) > 6.0) * (0.00)}
131
132}
133##########################
134# Muon tracking efficiency
135##########################
136
137module Efficiency MuonTrackingEfficiency {
138 set InputArray ParticlePropagator/muons
139 set OutputArray muons
140
141# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
142# what about high pT ?
143 # tracking efficiency formula for muons
144 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
145 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
146 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
147 (abs(eta) > 6.0) * (0.00)}
148
149}
150
151########################################
152# Momentum resolution for charged tracks
153########################################
154
155module MomentumSmearing ChargedHadronMomentumSmearing {
156 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
157 set OutputArray chargedHadrons
158
159 source momentumResolutionVsP.tcl
160}
161
162
163###################################
164# Momentum resolution for electrons
165###################################
166
167module MomentumSmearing ElectronMomentumSmearing {
168 set InputArray ElectronTrackingEfficiency/electrons
169 set OutputArray electrons
170
171 source momentumResolutionVsP.tcl
172}
173
174
175###############################
176# Momentum resolution for muons
177###############################
178
179module MomentumSmearing MuonMomentumSmearing {
180 set InputArray MuonTrackingEfficiency/muons
181 set OutputArray muons
182
183 # TBC for just putting tracker resolution/ need to add improvement at high pT
184
185 source muonMomentumResolutionVsP.tcl
186}
187
188##############
189# Track merger
190##############
191
192module Merger TrackMerger {
193# add InputArray InputArray
194 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
195 add InputArray ElectronMomentumSmearing/electrons
196 add InputArray MuonMomentumSmearing/muons
197 set OutputArray tracks
198}
199
200#############
201# ECAL
202#############
203
204# TBC : calos seems ok, check eta max value though.
205
206module SimpleCalorimeter ECal {
207 set ParticleInputArray ParticlePropagator/stableParticles
208 set TrackInputArray TrackMerger/tracks
209
210 set TowerOutputArray ecalTowers
211 set EFlowTrackOutputArray eflowTracks
212 set EFlowTowerOutputArray eflowPhotons
213
214 set IsEcal true
215
216 set EnergyMin 0.5
217 set EnergySignificanceMin 2.0
218
219 set SmearTowerCenter true
220
221 set pi [expr {acos(-1)}]
222
223 # lists of the edges of each tower in eta and phi
224 # each list starts with the lower edge of the first tower
225 # the list ends with the higher edged of the last tower
226
227 # 0.5 degree towers (5x5 mm^2)
228 set PhiBins {}
229 for {set i -256} {$i <= 256} {incr i} {
230 add PhiBins [expr {$i * $pi/256.0}]
231 }
232
233 # TBC
234 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} {
235 add EtaPhiBins $eta $PhiBins
236 }
237 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
238 set PhiBins {}
239 for {set i -128} {$i <= 128} {incr i} {
240 add PhiBins [expr {$i * $pi/128.0}]
241 }
242 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} {
243 add EtaPhiBins $eta $PhiBins
244 }
245
246 # 0.05 x 0.05 for eta between 4.0 and 6.0
247 set PhiBins {}
248 for {set i -64} {$i <= 64} {incr i} {
249 add PhiBins [expr {$i * $pi/64.0}]
250 }
251 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} {
252 add EtaPhiBins $eta $PhiBins
253 }
254 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
255
256 add EnergyFraction {0} {0.0}
257 # energy fractions for e, gamma and pi0
258 add EnergyFraction {11} {1.0}
259 add EnergyFraction {22} {1.0}
260 add EnergyFraction {111} {1.0}
261 # energy fractions for muon, neutrinos and neutralinos
262 add EnergyFraction {12} {0.0}
263 add EnergyFraction {13} {0.0}
264 add EnergyFraction {14} {0.0}
265 add EnergyFraction {16} {0.0}
266 add EnergyFraction {1000022} {0.0}
267 add EnergyFraction {1000023} {0.0}
268 add EnergyFraction {1000025} {0.0}
269 add EnergyFraction {1000035} {0.0}
270 add EnergyFraction {1000045} {0.0}
271 # energy fractions for K0short and Lambda
272 # add EnergyFraction {310} {0.3}
273 # add EnergyFraction {3122} {0.3}
274
275 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
276 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
277 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
278
279
280}
281
282#############
283# HCAL
284#############
285
286module SimpleCalorimeter HCal {
287 set ParticleInputArray ParticlePropagator/stableParticles
288 set TrackInputArray ECal/eflowTracks
289
290 set TowerOutputArray hcalTowers
291 set EFlowTrackOutputArray eflowTracks
292 set EFlowTowerOutputArray eflowNeutralHadrons
293
294 set IsEcal false
295
296 set EnergyMin 1.0
297 set EnergySignificanceMin 2.0
298
299 set SmearTowerCenter true
300
301 set pi [expr {acos(-1)}]
302
303 # lists of the edges of each tower in eta and phi
304 # each list starts with the lower edge of the first tower
305 # the list ends with the higher edged of the last tower
306
307 # 6 degree towers
308 set PhiBins {}
309 for {set i -64} {$i <= 64} {incr i} {
310 add PhiBins [expr {$i * $pi/64.0}]
311 }
312
313 # TBC
314 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} {
315 add EtaPhiBins $eta $PhiBins
316 }
317
318 # 0.1 x 0.1 between eta 2.5 and 4.0
319 set PhiBins {}
320 for {set i -32} {$i <= 32} {incr i} {
321 add PhiBins [expr {$i * $pi/32.0}]
322 }
323 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} {
324 add EtaPhiBins $eta $PhiBins
325 }
326
327 # 0.2 x 0.2 between eta 4.0 and 6.0
328 set PhiBins {}
329 for {set i -16} {$i <= 16} {incr i} {
330 add PhiBins [expr {$i * $pi/16.0}]
331 }
332 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} {
333 add EtaPhiBins $eta $PhiBins
334 }
335
336 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
337 add EnergyFraction {0} {1.0}
338 # energy fractions for e, gamma and pi0
339 add EnergyFraction {11} {0.0}
340 add EnergyFraction {22} {0.0}
341 add EnergyFraction {111} {0.0}
342 # energy fractions for muon, neutrinos and neutralinos
343 add EnergyFraction {12} {0.0}
344 add EnergyFraction {13} {0.0}
345 add EnergyFraction {14} {0.0}
346 add EnergyFraction {16} {0.0}
347 add EnergyFraction {1000022} {0.0}
348 add EnergyFraction {1000023} {0.0}
349 add EnergyFraction {1000025} {0.0}
350 add EnergyFraction {1000035} {0.0}
351 add EnergyFraction {1000045} {0.0}
352 # energy fractions for K0short and Lambda
353 # add EnergyFraction {310} {0.7}
354 # add EnergyFraction {3122} {0.7}
355
356 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
357 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
358 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
359}
360
361#################
362# Electron filter
363#################
364
365module PdgCodeFilter ElectronFilter {
366 set InputArray HCal/eflowTracks
367 set OutputArray electrons
368 set Invert true
369 add PdgCode {11}
370 add PdgCode {-11}
371}
372
373
374######################
375# ChargedHadronFilter
376######################
377
378module PdgCodeFilter ChargedHadronFilter {
379 set InputArray HCal/eflowTracks
380 set OutputArray chargedHadrons
381
382 add PdgCode {11}
383 add PdgCode {-11}
384 add PdgCode {13}
385 add PdgCode {-13}
386}
387
388
389###################################################
390# Tower Merger (in case not using e-flow algorithm)
391###################################################
392
393module Merger Calorimeter {
394# add InputArray InputArray
395 add InputArray ECal/ecalTowers
396 add InputArray HCal/hcalTowers
397 set OutputArray towers
398}
399
400####################
401# Energy flow merger
402####################
403
404module Merger EFlowMerger {
405# add InputArray InputArray
406 add InputArray HCal/eflowTracks
407 add InputArray ECal/eflowPhotons
408 add InputArray HCal/eflowNeutralHadrons
409 set OutputArray eflow
410}
411
412
413###################
414# Missing ET merger
415###################
416
417module Merger MissingET {
418# add InputArray InputArray
419 add InputArray EFlowMerger/eflow
420 set MomentumOutputArray momentum
421}
422
423
424
425##################
426# Scalar HT merger
427##################
428
429module Merger ScalarHT {
430# add InputArray InputArray
431 add InputArray EFlowMerger/eflow
432 set EnergyOutputArray energy
433}
434
435#################
436# Neutrino Filter
437#################
438
439module PdgCodeFilter NeutrinoFilter {
440
441 set InputArray Delphes/stableParticles
442 set OutputArray filteredParticles
443
444 set PTMin 0.0
445
446 add PdgCode {12}
447 add PdgCode {14}
448 add PdgCode {16}
449 add PdgCode {-12}
450 add PdgCode {-14}
451 add PdgCode {-16}
452
453}
454
455
456#####################
457# MC truth jet finder
458#####################
459
460# TBC: is jet radius fine?
461
462module FastJetFinder GenJetFinder {
463# set InputArray NeutrinoFilter/filteredParticles
464 set InputArray Delphes/stableParticles
465
466 set OutputArray jets
467
468 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
469 set JetAlgorithm 6
470 set ParameterR 0.4
471
472 set JetPTMin 5.0
473}
474
475#########################
476# Gen Missing ET merger
477########################
478
479module Merger GenMissingET {
480
481# add InputArray InputArray
482 add InputArray NeutrinoFilter/filteredParticles
483 set MomentumOutputArray momentum
484}
485
486
487
488############
489# Jet finder
490############
491
492# TBC need to include jet substructure variables
493# TBC is jet radius fine?
494
495module FastJetFinder FastJetFinder {
496# set InputArray Calorimeter/towers
497 set InputArray EFlowMerger/eflow
498
499 set OutputArray jets
500
501 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
502 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
503
504 set JetAlgorithm 6
505 set ParameterR 0.4
506
507 set JetPTMin 30.0
508}
509
510##################
511# Jet Energy Scale
512##################
513
514module EnergyScale JetEnergyScale {
515 set InputArray FastJetFinder/jets
516 set OutputArray jets
517
518 # scale formula for jets
519 set ScaleFormula {1.00}
520}
521
522
523########################
524# Jet Flavor Association
525########################
526
527module JetFlavorAssociation JetFlavorAssociation {
528
529 set PartonInputArray Delphes/partons
530 set ParticleInputArray Delphes/allParticles
531 set ParticleLHEFInputArray Delphes/allParticlesLHEF
532 set JetInputArray JetEnergyScale/jets
533
534 set DeltaR 0.5
535 set PartonPTMin 5.0
536 set PartonEtaMax 6.0
537
538}
539
540###################
541# Photon efficiency
542###################
543
544module Efficiency PhotonEfficiency {
545 set InputArray ECal/eflowPhotons
546 set OutputArray photons
547
548 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
549
550 set EfficiencyFormula {
551 (pt <= 1.0) * (0.00) + \
552 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
553 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
554 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
555
556 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
557 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
558 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
559
560 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
561 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
562 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
563 (abs(eta) > 6.0) * (0.00)}
564
565}
566
567##################
568# Photon isolation
569##################
570
571# TBC: check values for iso cuts
572
573module Isolation PhotonIsolation {
574 set CandidateInputArray PhotonEfficiency/photons
575 set IsolationInputArray EFlowMerger/eflow
576
577 set OutputArray photons
578
579 set DeltaRMax 0.3
580
581 set PTMin 0.5
582
583 set PTRatioMax 0.1
584}
585
586
587####################
588# Electron isolation
589####################
590
591# TBC: check values for iso cuts
592
593module Isolation ElectronIsolation {
594 set CandidateInputArray ElectronFilter/electrons
595 set IsolationInputArray EFlowMerger/eflow
596
597 set OutputArray electrons
598
599 set DeltaRMax 0.3
600
601 set PTMin 0.5
602
603 set PTRatioMax 0.1
604}
605
606
607################
608# Muon isolation
609################
610
611# TBC: check values for iso cuts
612
613module Isolation MuonIsolation {
614 set CandidateInputArray MuonMomentumSmearing/muons
615 set IsolationInputArray EFlowMerger/eflow
616
617 set OutputArray muons
618
619 set DeltaRMax 0.3
620
621 set PTMin 0.5
622
623 set PTRatioMax 0.2
624}
625
626
627###########
628# b-tagging
629###########
630
631module BTagging BTagging {
632 set JetInputArray JetEnergyScale/jets
633
634 set BitNumber 0
635
636 add EfficiencyFormula {0} {
637
638 (pt <= 10.0) * (0.00) +
639 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
640 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
641 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
642 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
643 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
644 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
645 (abs(eta) > 4.0) * (0.00)}
646
647 add EfficiencyFormula {4} {
648
649 (pt <= 10.0) * (0.00) +
650 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
651 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
652 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
653 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
654 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
655 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
656 (abs(eta) > 4.0) * (0.00)}
657
658 add EfficiencyFormula {5} {
659
660 (pt <= 10.0) * (0.00) +
661 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
662 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
663 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
664 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
665 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
666 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
667 (abs(eta) >= 4.0) * (0.00)}
668
669}
670
671###########
672# c-tagging
673###########
674
675module BTagging CTagging {
676 set JetInputArray JetEnergyScale/jets
677
678 set BitNumber 1
679
680 add EfficiencyFormula {0} {
681
682 (pt <= 10.0) * (0.00) +
683 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
684 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
685
686 add EfficiencyFormula {4} {
687
688 (pt <= 10.0) * (0.00) +
689 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
690 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
691
692 add EfficiencyFormula {5} {
693
694 (pt <= 10.0) * (0.00) +
695 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
696 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
697
698}
699
700
701#############
702# tau-tagging
703#############
704
705
706module TauTagging TauTagging {
707 set ParticleInputArray Delphes/allParticles
708 set PartonInputArray Delphes/partons
709 set JetInputArray JetEnergyScale/jets
710
711 set DeltaR 0.5
712
713 set TauPTMin 1.0
714
715 set TauEtaMax 4.0
716
717 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
718 add EfficiencyFormula {0} {
719
720 (pt <= 10.0) * (0.00) +
721 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
722 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
723 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
724 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
725 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
726 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
727 (abs(eta) > 4.0) * (0.00)}
728
729 add EfficiencyFormula {11} {
730
731 (pt <= 10.0) * (0.00) +
732 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
733 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
734 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
735 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
736 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
737 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
738 (abs(eta) > 4.0) * (0.00)}
739
740 add EfficiencyFormula {15} {
741
742 (pt <= 10.0) * (0.00) +
743 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
744 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
745 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
746 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
747 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
748 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
749 (abs(eta) > 4.0) * (0.00)}
750
751}
752
753#####################################################
754# Find uniquely identified photons/electrons/tau/jets
755#####################################################
756
757module UniqueObjectFinder UniqueObjectFinder {
758# earlier arrays take precedence over later ones
759# add InputArray InputArray OutputArray
760 add InputArray PhotonIsolation/photons photons
761 add InputArray ElectronIsolation/electrons electrons
762 add InputArray MuonIsolation/muons muons
763 add InputArray JetEnergyScale/jets jets
764
765}
766
767
768##################
769# ROOT tree writer
770##################
771
772module TreeWriter TreeWriter {
773# add Branch InputArray BranchName BranchClass
774 add Branch Delphes/allParticles Particle GenParticle
775
776 add Branch GenJetFinder/jets GenJet Jet
777 add Branch GenMissingET/momentum GenMissingET MissingET
778
779 add Branch TrackMerger/tracks Track Track
780 add Branch Calorimeter/towers Tower Tower
781
782 add Branch HCal/eflowTracks EFlowTrack Track
783 add Branch ECal/eflowPhotons EFlowPhoton Tower
784 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
785
786 add Branch UniqueObjectFinder/photons Photon Photon
787 add Branch UniqueObjectFinder/electrons Electron Electron
788 add Branch UniqueObjectFinder/muons Muon Muon
789 add Branch UniqueObjectFinder/jets Jet Jet
790
791 add Branch MissingET/momentum MissingET MissingET
792 add Branch ScalarHT/energy ScalarHT ScalarHT
793}
794
Note: See TracBrowser for help on using the repository browser.