Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage2.tcl@ 8dd735c

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