Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage3_fcal.tcl@ 25291dd

Last change on this file since 25291dd was 8dd735c, checked in by Aleksander Filip Zarnecki <zarnecki@…>, 4 years ago

File header updated

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