Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage3.tcl@ 4692fd9

Last change on this file since 4692fd9 was 0e6b14e, checked in by Ulrike Schnoor <ulrike.schnoor@…>, 5 years ago

add flavor tagging (b, tau) to the smeared jets in the CLICdet card
  • Property mode set to 100644
File size: 40.8 KB
Line 
1#######################################
2# CLICdet model
3# based on arXiv:1812.07337 and
4# CLICdp-Note-2017-001
5#
6# Ulrike Schnoor ulrike.schnoor@cern.ch
7#
8# For the high energy stage of
9# CLIC: 3 TeV
10#
11# Jet finding with Valencia algorithm:
12# use exclusive clustering with njets
13# according to final state
14#
15# c-tagging capabilities of CLICdet are
16# not yet implemented here. Please
17# contact us if you want to use it.
18#######################################
19
20#######################################
21# Order of execution of various modules
22#######################################
23
24set ExecutionPath {
25 ParticlePropagator
26
27 ChargedHadronTrackingEfficiency
28 ElectronTrackingEfficiency
29 MuonTrackingEfficiency
30
31 ChargedHadronMomentumSmearing
32 ElectronMomentumSmearing
33 MuonMomentumSmearing
34
35 TrackMerger
36
37 ECal
38 HCal
39
40 Calorimeter
41 EFlowMerger
42
43
44 PhotonEfficiency
45 PhotonIsolation
46
47 ElectronFilter
48 ElectronEfficiency
49 ElectronIsolation
50
51 ChargedHadronFilter
52
53 MuonEfficiency
54 MuonIsolation
55
56 EFlowFilter
57
58 NeutrinoFilter
59 GenJetFinder
60 FastJetFinderKt
61 FastJetFinderVLC_R05_N2
62 FastJetFinderVLC_R05_N3
63 FastJetFinderVLC_R05_N4
64 FastJetFinderVLC_R05_N5
65 FastJetFinderVLC_R05_N6
66 FastJetFinderVLC_R07_N2
67 FastJetFinderVLC_R07_N3
68 FastJetFinderVLC_R07_N4
69 FastJetFinderVLC_R07_N5
70 FastJetFinderVLC_R07_N6
71 FastJetFinderVLC_R10_N2
72 FastJetFinderVLC_R10_N3
73 FastJetFinderVLC_R10_N4
74 FastJetFinderVLC_R10_N5
75 FastJetFinderVLC_R10_N6
76 FastJetFinderVLC_R12_N2
77 FastJetFinderVLC_R12_N3
78 FastJetFinderVLC_R12_N4
79 FastJetFinderVLC_R12_N5
80 FastJetFinderVLC_R12_N6
81 FastJetFinderVLC_R15_N2
82 FastJetFinderVLC_R15_N3
83 FastJetFinderVLC_R15_N4
84 FastJetFinderVLC_R15_N5
85 FastJetFinderVLC_R15_N6
86
87 FastJetFinderVLC_R05_inclusive
88 FastJetFinderVLC_R07_inclusive
89 FastJetFinderVLC_R10_inclusive
90 FastJetFinderVLC_R12_inclusive
91 FastJetFinderVLC_R15_inclusive
92
93
94 MissingET
95 GenMissingET
96
97
98 JetMomentumSmearing_VLCR05N2
99 JetMomentumSmearing_VLCR05N3
100 JetMomentumSmearing_VLCR05N4
101 JetMomentumSmearing_VLCR05N5
102 JetMomentumSmearing_VLCR05N6
103 JetMomentumSmearing_VLCR05_inclusive
104 JetMomentumSmearing_VLCR07N2
105 JetMomentumSmearing_VLCR07N3
106 JetMomentumSmearing_VLCR07N4
107 JetMomentumSmearing_VLCR07N5
108 JetMomentumSmearing_VLCR07N6
109 JetMomentumSmearing_VLCR07_inclusive
110 JetMomentumSmearing_VLCR10N2
111 JetMomentumSmearing_VLCR10N3
112 JetMomentumSmearing_VLCR10N4
113 JetMomentumSmearing_VLCR10N5
114 JetMomentumSmearing_VLCR10N6
115 JetMomentumSmearing_VLCR10_inclusive
116 JetMomentumSmearing_VLCR12N2
117 JetMomentumSmearing_VLCR12N3
118 JetMomentumSmearing_VLCR12N4
119 JetMomentumSmearing_VLCR12N5
120 JetMomentumSmearing_VLCR12N6
121 JetMomentumSmearing_VLCR12_inclusive
122 JetMomentumSmearing_VLCR15N2
123 JetMomentumSmearing_VLCR15N3
124 JetMomentumSmearing_VLCR15N4
125 JetMomentumSmearing_VLCR15N5
126 JetMomentumSmearing_VLCR15N6
127 JetMomentumSmearing_VLCR15_inclusive
128
129
130
131
132 JetFlavorAssociation_R05N2
133 JetFlavorAssociation_R05N3
134 JetFlavorAssociation_R05N4
135 JetFlavorAssociation_R05N5
136 JetFlavorAssociation_R05N6
137
138 JetFlavorAssociation_R07N2
139 JetFlavorAssociation_R07N3
140 JetFlavorAssociation_R07N4
141 JetFlavorAssociation_R07N5
142 JetFlavorAssociation_R07N6
143
144 JetFlavorAssociation_R10N2
145 JetFlavorAssociation_R10N3
146 JetFlavorAssociation_R10N4
147 JetFlavorAssociation_R10N5
148 JetFlavorAssociation_R10N6
149
150 JetFlavorAssociation_R12N2
151 JetFlavorAssociation_R12N3
152 JetFlavorAssociation_R12N4
153 JetFlavorAssociation_R12N5
154 JetFlavorAssociation_R12N6
155
156 JetFlavorAssociation_R15N2
157 JetFlavorAssociation_R15N3
158 JetFlavorAssociation_R15N4
159 JetFlavorAssociation_R15N5
160 JetFlavorAssociation_R15N6
161
162 JetFlavorAssociation_R05_inclusive
163 JetFlavorAssociation_R07_inclusive
164 JetFlavorAssociation_R10_inclusive
165 JetFlavorAssociation_R12_inclusive
166 JetFlavorAssociation_R15_inclusive
167
168 JetFlavorAssociation_JER_R05N2
169 JetFlavorAssociation_JER_R05N3
170 JetFlavorAssociation_JER_R05N4
171 JetFlavorAssociation_JER_R05N5
172 JetFlavorAssociation_JER_R05N6
173
174 JetFlavorAssociation_JER_R07N2
175 JetFlavorAssociation_JER_R07N3
176 JetFlavorAssociation_JER_R07N4
177 JetFlavorAssociation_JER_R07N5
178 JetFlavorAssociation_JER_R07N6
179
180 JetFlavorAssociation_JER_R10N2
181 JetFlavorAssociation_JER_R10N3
182 JetFlavorAssociation_JER_R10N4
183 JetFlavorAssociation_JER_R10N5
184 JetFlavorAssociation_JER_R10N6
185
186 JetFlavorAssociation_JER_R12N2
187 JetFlavorAssociation_JER_R12N3
188 JetFlavorAssociation_JER_R12N4
189 JetFlavorAssociation_JER_R12N5
190 JetFlavorAssociation_JER_R12N6
191
192 JetFlavorAssociation_JER_R15N2
193 JetFlavorAssociation_JER_R15N3
194 JetFlavorAssociation_JER_R15N4
195 JetFlavorAssociation_JER_R15N5
196 JetFlavorAssociation_JER_R15N6
197
198 JetFlavorAssociation_JER_R05_inclusive
199 JetFlavorAssociation_JER_R07_inclusive
200 JetFlavorAssociation_JER_R10_inclusive
201 JetFlavorAssociation_JER_R12_inclusive
202 JetFlavorAssociation_JER_R15_inclusive
203
204 BTaggingWP50_R05N2
205 BTaggingWP70_R05N2
206 BTaggingWP90_R05N2
207 BTaggingWP50_R05N3
208 BTaggingWP70_R05N3
209 BTaggingWP90_R05N3
210 BTaggingWP50_R05N4
211 BTaggingWP70_R05N4
212 BTaggingWP90_R05N4
213 BTaggingWP50_R05N5
214 BTaggingWP70_R05N5
215 BTaggingWP90_R05N5
216 BTaggingWP50_R05N6
217 BTaggingWP70_R05N6
218 BTaggingWP90_R05N6
219 BTaggingWP50_R07N2
220 BTaggingWP70_R07N2
221 BTaggingWP90_R07N2
222 BTaggingWP50_R07N3
223 BTaggingWP70_R07N3
224 BTaggingWP90_R07N3
225 BTaggingWP50_R07N4
226 BTaggingWP70_R07N4
227 BTaggingWP90_R07N4
228 BTaggingWP50_R07N5
229 BTaggingWP70_R07N5
230 BTaggingWP90_R07N5
231 BTaggingWP50_R07N6
232 BTaggingWP70_R07N6
233 BTaggingWP90_R07N6
234 BTaggingWP50_R10N2
235 BTaggingWP70_R10N2
236 BTaggingWP90_R10N2
237 BTaggingWP50_R10N3
238 BTaggingWP70_R10N3
239 BTaggingWP90_R10N3
240 BTaggingWP50_R10N4
241 BTaggingWP70_R10N4
242 BTaggingWP90_R10N4
243 BTaggingWP50_R10N5
244 BTaggingWP70_R10N5
245 BTaggingWP90_R10N5
246 BTaggingWP50_R10N6
247 BTaggingWP70_R10N6
248 BTaggingWP90_R10N6
249 BTaggingWP50_R12N2
250 BTaggingWP70_R12N2
251 BTaggingWP90_R12N2
252 BTaggingWP50_R12N3
253 BTaggingWP70_R12N3
254 BTaggingWP90_R12N3
255 BTaggingWP50_R12N4
256 BTaggingWP70_R12N4
257 BTaggingWP90_R12N4
258 BTaggingWP50_R12N5
259 BTaggingWP70_R12N5
260 BTaggingWP90_R12N5
261 BTaggingWP50_R12N6
262 BTaggingWP70_R12N6
263 BTaggingWP90_R12N6
264 BTaggingWP50_R15N2
265 BTaggingWP70_R15N2
266 BTaggingWP90_R15N2
267 BTaggingWP50_R15N3
268 BTaggingWP70_R15N3
269 BTaggingWP90_R15N3
270 BTaggingWP50_R15N4
271 BTaggingWP70_R15N4
272 BTaggingWP90_R15N4
273 BTaggingWP50_R15N5
274 BTaggingWP70_R15N5
275 BTaggingWP90_R15N5
276 BTaggingWP50_R15N6
277 BTaggingWP70_R15N6
278 BTaggingWP90_R15N6
279 BTaggingWP50_R05_inclusive
280 BTaggingWP70_R05_inclusive
281 BTaggingWP90_R05_inclusive
282 BTaggingWP50_R07_inclusive
283 BTaggingWP70_R07_inclusive
284 BTaggingWP90_R07_inclusive
285 BTaggingWP50_R10_inclusive
286 BTaggingWP70_R10_inclusive
287 BTaggingWP90_R10_inclusive
288 BTaggingWP50_R12_inclusive
289 BTaggingWP70_R12_inclusive
290 BTaggingWP90_R12_inclusive
291 BTaggingWP50_R15_inclusive
292 BTaggingWP70_R15_inclusive
293 BTaggingWP90_R15_inclusive
294
295 BTagging_JER_WP50_R05N2
296 BTagging_JER_WP70_R05N2
297 BTagging_JER_WP90_R05N2
298 BTagging_JER_WP50_R05N3
299 BTagging_JER_WP70_R05N3
300 BTagging_JER_WP90_R05N3
301 BTagging_JER_WP50_R05N4
302 BTagging_JER_WP70_R05N4
303 BTagging_JER_WP90_R05N4
304 BTagging_JER_WP50_R05N5
305 BTagging_JER_WP70_R05N5
306 BTagging_JER_WP90_R05N5
307 BTagging_JER_WP50_R05N6
308 BTagging_JER_WP70_R05N6
309 BTagging_JER_WP90_R05N6
310 BTagging_JER_WP50_R07N2
311 BTagging_JER_WP70_R07N2
312 BTagging_JER_WP90_R07N2
313 BTagging_JER_WP50_R07N3
314 BTagging_JER_WP70_R07N3
315 BTagging_JER_WP90_R07N3
316 BTagging_JER_WP50_R07N4
317 BTagging_JER_WP70_R07N4
318 BTagging_JER_WP90_R07N4
319 BTagging_JER_WP50_R07N5
320 BTagging_JER_WP70_R07N5
321 BTagging_JER_WP90_R07N5
322 BTagging_JER_WP50_R07N6
323 BTagging_JER_WP70_R07N6
324 BTagging_JER_WP90_R07N6
325 BTagging_JER_WP50_R10N2
326 BTagging_JER_WP70_R10N2
327 BTagging_JER_WP90_R10N2
328 BTagging_JER_WP50_R10N3
329 BTagging_JER_WP70_R10N3
330 BTagging_JER_WP90_R10N3
331 BTagging_JER_WP50_R10N4
332 BTagging_JER_WP70_R10N4
333 BTagging_JER_WP90_R10N4
334 BTagging_JER_WP50_R10N5
335 BTagging_JER_WP70_R10N5
336 BTagging_JER_WP90_R10N5
337 BTagging_JER_WP50_R10N6
338 BTagging_JER_WP70_R10N6
339 BTagging_JER_WP90_R10N6
340 BTagging_JER_WP50_R12N2
341 BTagging_JER_WP70_R12N2
342 BTagging_JER_WP90_R12N2
343 BTagging_JER_WP50_R12N3
344 BTagging_JER_WP70_R12N3
345 BTagging_JER_WP90_R12N3
346 BTagging_JER_WP50_R12N4
347 BTagging_JER_WP70_R12N4
348 BTagging_JER_WP90_R12N4
349 BTagging_JER_WP50_R12N5
350 BTagging_JER_WP70_R12N5
351 BTagging_JER_WP90_R12N5
352 BTagging_JER_WP50_R12N6
353 BTagging_JER_WP70_R12N6
354 BTagging_JER_WP90_R12N6
355 BTagging_JER_WP50_R15N2
356 BTagging_JER_WP70_R15N2
357 BTagging_JER_WP90_R15N2
358 BTagging_JER_WP50_R15N3
359 BTagging_JER_WP70_R15N3
360 BTagging_JER_WP90_R15N3
361 BTagging_JER_WP50_R15N4
362 BTagging_JER_WP70_R15N4
363 BTagging_JER_WP90_R15N4
364 BTagging_JER_WP50_R15N5
365 BTagging_JER_WP70_R15N5
366 BTagging_JER_WP90_R15N5
367 BTagging_JER_WP50_R15N6
368 BTagging_JER_WP70_R15N6
369 BTagging_JER_WP90_R15N6
370 BTagging_JER_WP50_R05_inclusive
371 BTagging_JER_WP70_R05_inclusive
372 BTagging_JER_WP90_R05_inclusive
373 BTagging_JER_WP50_R07_inclusive
374 BTagging_JER_WP70_R07_inclusive
375 BTagging_JER_WP90_R07_inclusive
376 BTagging_JER_WP50_R10_inclusive
377 BTagging_JER_WP70_R10_inclusive
378 BTagging_JER_WP90_R10_inclusive
379 BTagging_JER_WP50_R12_inclusive
380 BTagging_JER_WP70_R12_inclusive
381 BTagging_JER_WP90_R12_inclusive
382 BTagging_JER_WP50_R15_inclusive
383 BTagging_JER_WP70_R15_inclusive
384 BTagging_JER_WP90_R15_inclusive
385
386 TauTagging_R05N2
387 TauTagging_R05N3
388 TauTagging_R05N4
389 TauTagging_R05N5
390 TauTagging_R05N6
391 TauTagging_R07N2
392 TauTagging_R07N3
393 TauTagging_R07N4
394 TauTagging_R07N5
395 TauTagging_R07N6
396 TauTagging_R10N2
397 TauTagging_R10N3
398 TauTagging_R10N4
399 TauTagging_R10N5
400 TauTagging_R10N6
401 TauTagging_R12N2
402 TauTagging_R12N3
403 TauTagging_R12N4
404 TauTagging_R12N5
405 TauTagging_R12N6
406 TauTagging_R15N2
407 TauTagging_R15N3
408 TauTagging_R15N4
409 TauTagging_R15N5
410 TauTagging_R15N6
411 TauTagging_R05_inclusive
412 TauTagging_R07_inclusive
413 TauTagging_R10_inclusive
414 TauTagging_R12_inclusive
415 TauTagging_R15_inclusive
416
417 TauTagging_JER_R05N2
418 TauTagging_JER_R05N3
419 TauTagging_JER_R05N4
420 TauTagging_JER_R05N5
421 TauTagging_JER_R05N6
422 TauTagging_JER_R07N2
423 TauTagging_JER_R07N3
424 TauTagging_JER_R07N4
425 TauTagging_JER_R07N5
426 TauTagging_JER_R07N6
427 TauTagging_JER_R10N2
428 TauTagging_JER_R10N3
429 TauTagging_JER_R10N4
430 TauTagging_JER_R10N5
431 TauTagging_JER_R10N6
432 TauTagging_JER_R12N2
433 TauTagging_JER_R12N3
434 TauTagging_JER_R12N4
435 TauTagging_JER_R12N5
436 TauTagging_JER_R12N6
437 TauTagging_JER_R15N2
438 TauTagging_JER_R15N3
439 TauTagging_JER_R15N4
440 TauTagging_JER_R15N5
441 TauTagging_JER_R15N6
442 TauTagging_JER_R05_inclusive
443 TauTagging_JER_R07_inclusive
444 TauTagging_JER_R10_inclusive
445 TauTagging_JER_R12_inclusive
446 TauTagging_JER_R15_inclusive
447
448
449 ScalarHT
450
451
452 TreeWriter
453}
454
455#################################
456# Propagate particles in cylinder
457#################################
458
459module ParticlePropagator ParticlePropagator {
460 set InputArray Delphes/stableParticles
461
462 set OutputArray stableParticles
463 set ChargedHadronOutputArray chargedHadrons
464 set ElectronOutputArray electrons
465 set MuonOutputArray muons
466
467 # radius of the magnetic field coverage in the calorimeter, in m
468 set Radius 1.5
469 # half-length of the magnetic field coverage in the calorimeter, in m
470 set HalfLength 2.31
471
472 # magnetic field, in T
473 set Bz 4.0
474}
475
476####################################
477# Charged hadron tracking efficiency
478####################################
479
480module Efficiency ChargedHadronTrackingEfficiency {
481 set InputArray ParticlePropagator/chargedHadrons
482 set OutputArray chargedHadrons
483 # Current full simulation with CLICdet provides for pions:
484
485 set EfficiencyFormula {
486 (abs(eta) > 2.54) * (0.000) +
487 (energy >= 80) * (abs(eta) < 2.54) * (1.000) +
488 (energy < 80 && energy >= 3) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
489 (energy < 80 && energy >= 3) * (abs(eta) <= 2.34) * (1.000) +
490 (energy < 3) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) +
491 (energy < 3) * (abs(eta) <= 0.55 ) * (1.000)
492 }
493}
494
495##############################
496# Electron tracking efficiency
497##############################
498
499module Efficiency ElectronTrackingEfficiency {
500 set InputArray ParticlePropagator/electrons
501 set OutputArray electrons
502
503
504 # Current full simulation with CLICdet provides for electrons:
505 set EfficiencyFormula {
506 (pt <= 1) * (0.000) +
507 (abs(eta) > 2.54) * (0.000) +
508 (energy >= 80) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
509 (energy >= 80) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
510 (energy >= 80) * (abs(eta) <= 2.34 ) * (1.000) +
511 (energy < 80 && energy >= 5) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
512 (energy < 80 && energy >= 5) * (abs(eta) <= 2.17) * (1.000) +
513 (energy < 5) * (abs(eta) > 2.34 ) * (0.000) +
514 (energy < 5) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
515 (energy < 5) * (abs(eta) <= 0.76) * (0.999)
516 }
517}
518
519##########################
520# Muon tracking efficiency
521##########################
522
523module Efficiency MuonTrackingEfficiency {
524 set InputArray ParticlePropagator/muons
525 set OutputArray muons
526
527 # Current full simulation with CLICdet provides for muons:
528 set EfficiencyFormula {
529 (pt < 1) * (0.000) +
530 (abs(eta) > 2.54) * (0.000) +
531 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 80) * (0.994) +
532 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 5 && energy < 80) * (0.996) +
533 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy < 5 ) * (0.996) +
534 (abs(eta) <= 2.44 ) * (energy >= 5 ) * (1.000) +
535 (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy < 5 && pt >=1 ) * (0.999) +
536 (abs(eta) <= 2.25 ) * (energy >= 1) * (1.000)
537
538
539 }
540}
541
542########################################
543# Momentum resolution for charged tracks
544########################################
545
546module MomentumSmearing ChargedHadronMomentumSmearing {
547 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
548 set OutputArray chargedHadrons
549
550
551 # Resolution given in dpT/pT.
552 # CLICdet internal studies
553 set ResolutionFormula {
554 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
555 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
556 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
557 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
558 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
559 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
560
561 }
562}
563
564###################################
565# Momentum resolution for electrons
566###################################
567
568module MomentumSmearing ElectronMomentumSmearing {
569 set InputArray ElectronTrackingEfficiency/electrons
570 set OutputArray electrons
571
572 # Resolution given in dpT/pT.
573 # CLICdet internal studies
574 set ResolutionFormula {
575 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
576 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
577 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
578 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
579 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
580 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
581
582 }
583}
584
585###############################
586# Momentum resolution for muons
587###############################
588
589module MomentumSmearing MuonMomentumSmearing {
590 set InputArray MuonTrackingEfficiency/muons
591 set OutputArray muons
592
593 # Resolution given in dpT/pT.
594
595 # CLICdet internal studies
596 set ResolutionFormula {
597
598 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
599 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
600 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt(9.81626e-06^2 * pt^2 + 0.00379895^2 ) +
601 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt(1.1959e-05^2 * pt^2 + 0.00242417^2 ) +
602 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
603 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
604
605 }
606}
607
608##############
609# Track merger
610##############
611
612module Merger TrackMerger {
613 # add InputArray InputArray
614 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
615 add InputArray ElectronMomentumSmearing/electrons
616 add InputArray MuonMomentumSmearing/muons
617 set OutputArray tracks
618}
619
620#############
621# ECAL
622#############
623
624module SimpleCalorimeter ECal {
625 set ParticleInputArray ParticlePropagator/stableParticles
626 set TrackInputArray TrackMerger/tracks
627
628 set TowerOutputArray ecalTowers
629 set EFlowTrackOutputArray eflowTracks
630 set EFlowTowerOutputArray eflowPhotons
631
632 set IsEcal true
633
634 set EnergyMin 0.5
635 set EnergySignificanceMin 1.0
636
637 set SmearTowerCenter true
638
639 set pi [expr {acos(-1)}]
640
641 # lists of the edges of each tower in eta and phi
642 # each list starts with the lower edge of the first tower
643 # the list ends with the higher edged of the last tower
644
645 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
646 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
647 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
648 #ECAL cell sizes always 5x5 mm^2
649
650 #barrel:
651 #dphi = 0.2 degree towers up to eta <=1.2
652 set PhiBins {}
653 for {set i -900} {$i <= 900} {incr i} {
654 add PhiBins [expr {$i * $pi/900.0 }]
655 }
656 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
657 for {set i -400} {$i <=400} {incr i} {
658 set eta [expr {$i * 0.003}]
659 add EtaPhiBins $eta $PhiBins
660 }
661
662 #endcaps:
663 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
664 set PhiBins {}
665 for {set i -225} {$i <= 225} {incr i} {
666 add PhiBins [expr {$i * $pi/225.}]
667 }
668 #deta=0.02 units for 1.2 < |eta| <=2.5
669 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
670 for {set i 1} {$i <=66} {incr i} {
671 set eta [expr {-2.52 + $i * 0.02}]
672 add EtaPhiBins $eta $PhiBins
673 }
674 #same for 1.2 to 2.5
675 for {set i 1} {$i <=66} {incr i} {
676 set eta [expr {1.18 + $i*0.02}]
677 add EtaPhiBins $eta $PhiBins
678 }
679
680 #plug:
681 #dphi = 1 degree for 2.5 < eta <=3
682 set PhiBins {}
683 for {set i -180} {$i <= 180} {incr i} {
684 add PhiBins [expr {$i * $pi/180.}]
685 }
686 # deta = 0.02 for 2.5 < |eta| <=3
687 # from -3 to -2.5, there will be 25 segments
688 for {set i 1} {$i <= 26} {incr i} {
689 set eta [expr {-3.02 + $i * 0.02}]
690 add EtaPhiBins $eta $PhiBins
691 }
692 #same for 2.5 to 3
693 for {set i 1} {$i <= 26} {incr i} {
694 set eta [expr {2.48 + $i*0.02}]
695 add EtaPhiBins $eta $PhiBins
696 }
697
698
699
700 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
701
702 add EnergyFraction {0} {0.0}
703 # energy fractions for e, gamma and pi0
704 add EnergyFraction {11} {1.0}
705 add EnergyFraction {22} {1.0}
706 add EnergyFraction {111} {1.0}
707 # energy fractions for muon, neutrinos and neutralinos
708 add EnergyFraction {12} {0.0}
709 add EnergyFraction {13} {0.0}
710 add EnergyFraction {14} {0.0}
711 add EnergyFraction {16} {0.0}
712 add EnergyFraction {1000022} {0.0}
713 add EnergyFraction {1000023} {0.0}
714 add EnergyFraction {1000025} {0.0}
715 add EnergyFraction {1000035} {0.0}
716 add EnergyFraction {1000045} {0.0}
717 # energy fractions for K0short and Lambda
718 add EnergyFraction {310} {0.3}
719 add EnergyFraction {3122} {0.3}
720
721 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
722 set ResolutionFormula {
723 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
724 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
725 (abs(eta) <= 3 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
726}
727
728#############
729# HCAL
730#############
731
732module SimpleCalorimeter HCal {
733 set ParticleInputArray ParticlePropagator/stableParticles
734 set TrackInputArray ECal/eflowTracks
735
736 set TowerOutputArray hcalTowers
737 set EFlowTrackOutputArray eflowTracks
738 set EFlowTowerOutputArray eflowNeutralHadrons
739
740 set IsEcal false
741
742 set EnergyMin 1.0
743 set EnergySignificanceMin 1.0
744
745 set SmearTowerCenter true
746
747 set pi [expr {acos(-1)}]
748
749 # lists of the edges of each tower in eta and phi
750 # each list starts with the lower edge of the first tower
751 # the list ends with the higher edged of the last tower
752
753
754 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
755 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
756 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
757 #HCAL cell sizes always 30x30 mm^2
758
759 #barrel and ring:
760 #dphi = 1 degree up to |eta| <=0.9
761 set PhiBins {}
762 for {set i -180} {$i <=180} {incr i} {
763 add PhiBins [expr {$i * $pi/180.0}]
764 }
765 #deta= 0.02 towers up to |eta| <=0.9
766 for {set i -45} {$i <=45} {incr i} {
767 set eta [expr {$i * 0.02}]
768 add EtaPhiBins $eta $PhiBins
769 }
770
771 #endcaps:
772 # dphi = 6 degree
773 set PhiBins {}
774 for {set i -30} {$i <=30} {incr i} {
775 add PhiBins [expr {$i * $pi/30.0}]
776 }
777 # deta =0.1 for 0.9 < |eta| <=3.0
778 #for -3.0 to -0.9, 21 segments
779 for {set i 1} {$i <=22} {incr i} {
780 set eta [expr {-3.1 + $i * 0.1}]
781 add EtaPhiBins $eta $PhiBins
782 }
783 #same for 0.9 to 3.0
784 for {set i 1} {$i <=22} {incr i} {
785 set eta [expr {0.8 + $i * 0.1 }]
786 add EtaPhiBins $eta $PhiBins
787 }
788
789 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
790 add EnergyFraction {0} {1.0}
791 # energy fractions for e, gamma and pi0
792 add EnergyFraction {11} {0.0}
793 add EnergyFraction {22} {0.0}
794 add EnergyFraction {111} {0.0}
795 # energy fractions for muon, neutrinos and neutralinos
796 add EnergyFraction {12} {0.0}
797 add EnergyFraction {13} {0.0}
798 add EnergyFraction {14} {0.0}
799 add EnergyFraction {16} {0.0}
800 add EnergyFraction {1000022} {0.0}
801 add EnergyFraction {1000023} {0.0}
802 add EnergyFraction {1000025} {0.0}
803 add EnergyFraction {1000035} {0.0}
804 add EnergyFraction {1000045} {0.0}
805 # energy fractions for K0short and Lambda
806 add EnergyFraction {310} {0.7}
807 add EnergyFraction {3122} {0.7}
808
809 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
810 #CLICdet internal studies
811 set ResolutionFormula {
812 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
813 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
814 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
815 (abs(eta)<=3 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
816 }
817
818}
819
820#################
821# Electron filter
822#################
823
824module PdgCodeFilter ElectronFilter {
825 set InputArray HCal/eflowTracks
826 set OutputArray electrons
827 set Invert true
828 add PdgCode {11}
829 add PdgCode {-11}
830}
831
832######################
833# ChargedHadronFilter
834######################
835
836module PdgCodeFilter ChargedHadronFilter {
837 set InputArray HCal/eflowTracks
838 set OutputArray chargedHadrons
839
840 add PdgCode {11}
841 add PdgCode {-11}
842 add PdgCode {13}
843 add PdgCode {-13}
844}
845
846
847
848###################################################
849# Tower Merger (in case not using e-flow algorithm)
850###################################################
851
852module Merger Calorimeter {
853 # add InputArray InputArray
854 add InputArray ECal/ecalTowers
855 add InputArray HCal/hcalTowers
856 set OutputArray towers
857}
858
859
860####################
861# Energy flow merger
862####################
863
864module Merger EFlowMerger {
865 # add InputArray InputArray
866 add InputArray HCal/eflowTracks
867 add InputArray ECal/eflowPhotons
868 add InputArray HCal/eflowNeutralHadrons
869 set OutputArray eflow
870}
871
872
873###################
874# Photon efficiency
875###################
876
877module Efficiency PhotonEfficiency {
878 set InputArray ECal/eflowPhotons
879 set OutputArray photons
880
881 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
882
883 # efficiency formula for photons
884 # current full simulation of CLICdet yields:
885 set EfficiencyFormula {
886 (energy < 2.0 ) * (0.000) +
887 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
888 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) }
889
890}
891
892
893##################
894# Photon isolation
895##################
896
897module Isolation PhotonIsolation {
898 set CandidateInputArray PhotonEfficiency/photons
899 set IsolationInputArray EFlowMerger/eflow
900
901 set OutputArray photons
902
903 set DeltaRMax 0.1
904
905 set PTMin 0.5
906
907 set PTRatioMax 0.2
908}
909
910
911#####################
912# Electron efficiency
913#####################
914
915module Efficiency ElectronEfficiency {
916 set InputArray ElectronFilter/electrons
917 set OutputArray electrons
918
919 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
920
921 set EfficiencyFormula {
922 (energy < 3.0) * ( 0.00 ) +
923 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95) * (0.58 ) +
924 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
925 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
926 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
927 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
928 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
929 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95) * ( 0.6 ) +
930 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
931 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
932 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
933 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
934 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
935 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95) * ( 0.6 ) +
936 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
937 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
938 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
939 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
940 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
941 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95) * (0.64 ) +
942 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
943 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
944 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
945 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
946 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
947 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95) * (0.64 ) +
948 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
949 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
950 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
951 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
952 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
953 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95) * (0.67 ) +
954 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
955 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
956 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
957 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
958 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
959 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95) * (0.68 ) +
960 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
961 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
962 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
963 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
964 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
965 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95) * ( 0.7 ) +
966 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
967 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
968 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
969 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
970 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
971 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95) * (0.68 ) +
972 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
973 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
974 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
975 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
976 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
977 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95) * ( 0.68 ) +
978 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
979 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
980 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
981 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
982 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
983 ( energy >=400 ) * (abs(eta) > 1.95) * (0.68 ) +
984 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
985 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
986 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
987 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
988 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
989 }
990}
991
992####################
993# Electron isolation
994####################
995
996module Isolation ElectronIsolation {
997 set CandidateInputArray ElectronEfficiency/electrons
998 set IsolationInputArray EFlowMerger/eflow
999
1000 set OutputArray electrons
1001
1002 set DeltaRMax 0.1
1003
1004 set PTMin 0.5
1005
1006 set PTRatioMax 0.2
1007}
1008
1009#################
1010# Muon efficiency
1011#################
1012
1013module Efficiency MuonEfficiency {
1014 set InputArray MuonMomentumSmearing/muons
1015 set OutputArray muons
1016
1017 # set EfficiencyFormula {efficiency as a function of eta and pt}
1018
1019 # efficiency formula for muons
1020 # current full simulation of CLICdet yields:
1021
1022
1023
1024 set EfficiencyFormula {
1025 (energy < 2.0 ) * (0.00) +
1026 (energy>=2.0) * (0.999)
1027 }
1028}
1029
1030################
1031# Muon isolation
1032################
1033
1034module Isolation MuonIsolation {
1035 set CandidateInputArray MuonEfficiency/muons
1036 set IsolationInputArray EFlowMerger/eflow
1037
1038 set OutputArray muons
1039
1040 set DeltaRMax 0.1
1041
1042 set PTMin 0.5
1043
1044 set PTRatioMax 0.2
1045}
1046
1047
1048
1049###################
1050# Missing ET merger
1051###################
1052
1053module Merger MissingET {
1054 # add InputArray InputArray
1055 add InputArray EFlowMerger/eflow
1056 set MomentumOutputArray momentum
1057}
1058
1059
1060##################
1061# Scalar HT merger
1062##################
1063
1064module Merger ScalarHT {
1065 # add InputArray InputArray
1066 add InputArray EFlowMerger/eflow
1067 set EnergyOutputArray energy
1068}
1069######################
1070# EFlowFilter (UniqueObjectFinder)
1071######################
1072module UniqueObjectFinder EFlowFilter {
1073 add InputArray PhotonIsolation/photons photons
1074 add InputArray ElectronIsolation/electrons electrons
1075 add InputArray MuonIsolation/muons muons
1076 add InputArray EFlowMerger/eflow eflow
1077}
1078
1079#################
1080# Neutrino Filter
1081#################
1082
1083module PdgCodeFilter NeutrinoFilter {
1084
1085 set InputArray Delphes/stableParticles
1086 set OutputArray filteredParticles
1087
1088 set PTMin 0.0
1089
1090 add PdgCode {12}
1091 add PdgCode {14}
1092 add PdgCode {16}
1093 add PdgCode {-12}
1094 add PdgCode {-14}
1095 add PdgCode {-16}
1096
1097}
1098
1099
1100#####################
1101# MC truth jet finder
1102#####################
1103
1104module FastJetFinder GenJetFinder {
1105 set InputArray NeutrinoFilter/filteredParticles
1106
1107 set OutputArray jets
1108
1109 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1110 set JetAlgorithm 9
1111 set ParameterR 0.5
1112
1113 set JetPTMin 20.0
1114}
1115
1116#########################
1117# Gen Missing ET merger
1118########################
1119
1120module Merger GenMissingET {
1121 # add InputArray InputArray
1122 add InputArray NeutrinoFilter/filteredParticles
1123 set MomentumOutputArray momentum
1124}
1125
1126
1127
1128############
1129# Jet finder
1130############
1131
1132module FastJetFinder FastJetFinderKt {
1133 # set InputArray Calorimeter/towers
1134 set InputArray EFlowMerger/eflow
1135
1136 set OutputArray KTjets
1137
1138 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1139 set JetAlgorithm 4
1140 set ParameterR 0.5
1141
1142 set JetPTMin 20.0
1143}
1144
1145################
1146# Jet finder VLC
1147################
1148
1149
1150source CLIC/CLICdet_JetReco.tcl
1151
1152
1153#########################################
1154# Jet Momentum Smearing to mimick overlay
1155#########################################
1156
1157
1158source CLIC/CLICdet_JetSmearing.tcl
1159
1160
1161
1162########################
1163# Jet Flavor Association
1164########################
1165
1166source CLIC/CLICdet_JetFlavorAssociation.tcl
1167source CLIC/CLICdet_JetFlavorAssociation_JER.tcl
1168
1169###########
1170# b-tagging
1171###########
1172# based on CLICdp-Note-2014-002
1173
1174source CLIC/CLICdet_BTagging.tcl
1175source CLIC/CLICdet_BTagging_JER.tcl
1176
1177
1178#############
1179# tau-tagging
1180#############
1181# based on LCD-2010-009
1182
1183source CLIC/CLICdet_TauTagging.tcl
1184source CLIC/CLICdet_TauTagging_JER.tcl
1185
1186
1187
1188##################
1189# ROOT tree writer
1190##################
1191
1192module TreeWriter TreeWriter {
1193 # add Branch InputArray BranchName BranchClass
1194 add Branch Delphes/allParticles Particle GenParticle
1195
1196 add Branch GenJetFinder/jets GenJet Jet
1197
1198 add Branch FastJetFinderKt/KTjets KTjet Jet
1199 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1200 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1201 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1202 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1203 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1204
1205 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1206 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1207 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1208 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1209 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1210
1211 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1212 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1213 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1214 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1215 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1216
1217 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1218 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1219 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1220 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1221 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1222
1223 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1224 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1225 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1226 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1227 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1228
1229 add Branch FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive VLCjetR05_inclusive Jet
1230 add Branch FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive VLCjetR07_inclusive Jet
1231 add Branch FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive VLCjetR10_inclusive Jet
1232 add Branch FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive VLCjetR12_inclusive Jet
1233 add Branch FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive VLCjetR15_inclusive Jet
1234
1235
1236
1237
1238 ###with jet energy smearing
1239 add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 JER_VLCjetR05N2 Jet
1240 add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 JER_VLCjetR05N3 Jet
1241 add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 JER_VLCjetR05N4 Jet
1242 add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 JER_VLCjetR05N5 Jet
1243 add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 JER_VLCjetR05N6 Jet
1244
1245 add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 JER_VLCjetR07N2 Jet
1246 add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 JER_VLCjetR07N3 Jet
1247 add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 JER_VLCjetR07N4 Jet
1248 add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 JER_VLCjetR07N5 Jet
1249 add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 JER_VLCjetR07N6 Jet
1250
1251 add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 JER_VLCjetR10N2 Jet
1252 add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 JER_VLCjetR10N3 Jet
1253 add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 JER_VLCjetR10N4 Jet
1254 add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 JER_VLCjetR10N5 Jet
1255 add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 JER_VLCjetR10N6 Jet
1256
1257 add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 JER_VLCjetR12N2 Jet
1258 add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 JER_VLCjetR12N3 Jet
1259 add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 JER_VLCjetR12N4 Jet
1260 add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 JER_VLCjetR12N5 Jet
1261 add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 JER_VLCjetR12N6 Jet
1262
1263 add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 JER_VLCjetR15N2 Jet
1264 add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 JER_VLCjetR15N3 Jet
1265 add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 JER_VLCjetR15N4 Jet
1266 add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 JER_VLCjetR15N5 Jet
1267 add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 JER_VLCjetR15N6 Jet
1268
1269 add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive JER_VLCjetR05_inclusive Jet
1270 add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive JER_VLCjetR07_inclusive Jet
1271 add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive JER_VLCjetR10_inclusive Jet
1272 add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive JER_VLCjetR12_inclusive Jet
1273 add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive JER_VLCjetR15_inclusive Jet
1274
1275 ####
1276
1277 add Branch GenMissingET/momentum GenMissingET MissingET
1278
1279 add Branch TrackMerger/tracks Track Track
1280 add Branch Calorimeter/towers Tower Tower
1281
1282 add Branch HCal/eflowTracks EFlowTrack Track
1283 add Branch ECal/eflowPhotons EFlowPhoton Tower
1284 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1285
1286 add Branch EFlowFilter/photons Photon Photon
1287 add Branch EFlowFilter/electrons Electron Electron
1288 add Branch EFlowFilter/muons Muon Muon
1289
1290 add Branch MissingET/momentum MissingET MissingET
1291 add Branch ScalarHT/energy ScalarHT ScalarHT
1292}
1293
Note: See TracBrowser for help on using the repository browser.