Fork me on GitHub

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

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

isolation criteria adjusted for CLICdet
  • Property mode set to 100644
File size: 32.0 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 low energy stage of
9# CLIC: 380 (350) GeV
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
99
100 JetFlavorAssociation_R05N2
101 JetFlavorAssociation_R05N3
102 JetFlavorAssociation_R05N4
103 JetFlavorAssociation_R05N5
104 JetFlavorAssociation_R05N6
105
106 JetFlavorAssociation_R07N2
107 JetFlavorAssociation_R07N3
108 JetFlavorAssociation_R07N4
109 JetFlavorAssociation_R07N5
110 JetFlavorAssociation_R07N6
111
112 JetFlavorAssociation_R10N2
113 JetFlavorAssociation_R10N3
114 JetFlavorAssociation_R10N4
115 JetFlavorAssociation_R10N5
116 JetFlavorAssociation_R10N6
117
118 JetFlavorAssociation_R12N2
119 JetFlavorAssociation_R12N3
120 JetFlavorAssociation_R12N4
121 JetFlavorAssociation_R12N5
122 JetFlavorAssociation_R12N6
123
124 JetFlavorAssociation_R15N2
125 JetFlavorAssociation_R15N3
126 JetFlavorAssociation_R15N4
127 JetFlavorAssociation_R15N5
128 JetFlavorAssociation_R15N6
129
130 JetFlavorAssociation_R05_inclusive
131 JetFlavorAssociation_R07_inclusive
132 JetFlavorAssociation_R10_inclusive
133 JetFlavorAssociation_R12_inclusive
134 JetFlavorAssociation_R15_inclusive
135
136
137 BTaggingWP50_R05N2
138 BTaggingWP70_R05N2
139 BTaggingWP90_R05N2
140 BTaggingWP50_R05N3
141 BTaggingWP70_R05N3
142 BTaggingWP90_R05N3
143 BTaggingWP50_R05N4
144 BTaggingWP70_R05N4
145 BTaggingWP90_R05N4
146 BTaggingWP50_R05N5
147 BTaggingWP70_R05N5
148 BTaggingWP90_R05N5
149 BTaggingWP50_R05N6
150 BTaggingWP70_R05N6
151 BTaggingWP90_R05N6
152 BTaggingWP50_R07N2
153 BTaggingWP70_R07N2
154 BTaggingWP90_R07N2
155 BTaggingWP50_R07N3
156 BTaggingWP70_R07N3
157 BTaggingWP90_R07N3
158 BTaggingWP50_R07N4
159 BTaggingWP70_R07N4
160 BTaggingWP90_R07N4
161 BTaggingWP50_R07N5
162 BTaggingWP70_R07N5
163 BTaggingWP90_R07N5
164 BTaggingWP50_R07N6
165 BTaggingWP70_R07N6
166 BTaggingWP90_R07N6
167 BTaggingWP50_R10N2
168 BTaggingWP70_R10N2
169 BTaggingWP90_R10N2
170 BTaggingWP50_R10N3
171 BTaggingWP70_R10N3
172 BTaggingWP90_R10N3
173 BTaggingWP50_R10N4
174 BTaggingWP70_R10N4
175 BTaggingWP90_R10N4
176 BTaggingWP50_R10N5
177 BTaggingWP70_R10N5
178 BTaggingWP90_R10N5
179 BTaggingWP50_R10N6
180 BTaggingWP70_R10N6
181 BTaggingWP90_R10N6
182 BTaggingWP50_R12N2
183 BTaggingWP70_R12N2
184 BTaggingWP90_R12N2
185 BTaggingWP50_R12N3
186 BTaggingWP70_R12N3
187 BTaggingWP90_R12N3
188 BTaggingWP50_R12N4
189 BTaggingWP70_R12N4
190 BTaggingWP90_R12N4
191 BTaggingWP50_R12N5
192 BTaggingWP70_R12N5
193 BTaggingWP90_R12N5
194 BTaggingWP50_R12N6
195 BTaggingWP70_R12N6
196 BTaggingWP90_R12N6
197 BTaggingWP50_R15N2
198 BTaggingWP70_R15N2
199 BTaggingWP90_R15N2
200 BTaggingWP50_R15N3
201 BTaggingWP70_R15N3
202 BTaggingWP90_R15N3
203 BTaggingWP50_R15N4
204 BTaggingWP70_R15N4
205 BTaggingWP90_R15N4
206 BTaggingWP50_R15N5
207 BTaggingWP70_R15N5
208 BTaggingWP90_R15N5
209 BTaggingWP50_R15N6
210 BTaggingWP70_R15N6
211 BTaggingWP90_R15N6
212 BTaggingWP50_R05_inclusive
213 BTaggingWP70_R05_inclusive
214 BTaggingWP90_R05_inclusive
215 BTaggingWP50_R07_inclusive
216 BTaggingWP70_R07_inclusive
217 BTaggingWP90_R07_inclusive
218 BTaggingWP50_R10_inclusive
219 BTaggingWP70_R10_inclusive
220 BTaggingWP90_R10_inclusive
221 BTaggingWP50_R12_inclusive
222 BTaggingWP70_R12_inclusive
223 BTaggingWP90_R12_inclusive
224 BTaggingWP50_R15_inclusive
225 BTaggingWP70_R15_inclusive
226 BTaggingWP90_R15_inclusive
227
228
229 TauTagging_R05N2
230 TauTagging_R05N3
231 TauTagging_R05N4
232 TauTagging_R05N5
233 TauTagging_R05N6
234 TauTagging_R07N2
235 TauTagging_R07N3
236 TauTagging_R07N4
237 TauTagging_R07N5
238 TauTagging_R07N6
239 TauTagging_R10N2
240 TauTagging_R10N3
241 TauTagging_R10N4
242 TauTagging_R10N5
243 TauTagging_R10N6
244 TauTagging_R12N2
245 TauTagging_R12N3
246 TauTagging_R12N4
247 TauTagging_R12N5
248 TauTagging_R12N6
249 TauTagging_R15N2
250 TauTagging_R15N3
251 TauTagging_R15N4
252 TauTagging_R15N5
253 TauTagging_R15N6
254
255 TauTagging_R05_inclusive
256 TauTagging_R07_inclusive
257 TauTagging_R10_inclusive
258 TauTagging_R12_inclusive
259 TauTagging_R15_inclusive
260
261
262
263 ScalarHT
264
265
266 TreeWriter
267}
268
269#################################
270# Propagate particles in cylinder
271#################################
272
273module ParticlePropagator ParticlePropagator {
274 set InputArray Delphes/stableParticles
275
276 set OutputArray stableParticles
277 set ChargedHadronOutputArray chargedHadrons
278 set ElectronOutputArray electrons
279 set MuonOutputArray muons
280
281 # radius of the magnetic field coverage in the calorimeter, in m
282 set Radius 1.5
283 # half-length of the magnetic field coverage in the calorimeter, in m
284 set HalfLength 2.31
285
286 # magnetic field, in T
287 set Bz 4.0
288}
289
290####################################
291# Charged hadron tracking efficiency
292####################################
293
294module Efficiency ChargedHadronTrackingEfficiency {
295 set InputArray ParticlePropagator/chargedHadrons
296 set OutputArray chargedHadrons
297 # Current full simulation with CLICdet provides for pions:
298
299 set EfficiencyFormula {
300 (abs(eta) > 2.54) * (0.000) +
301 (energy >= 80) * (abs(eta) < 2.54) * (1.000) +
302 (energy < 80 && energy >= 3) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
303 (energy < 80 && energy >= 3) * (abs(eta) <= 2.34) * (1.000) +
304 (energy < 3) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) +
305 (energy < 3) * (abs(eta) <= 0.55 ) * (1.000)
306 }
307}
308
309##############################
310# Electron tracking efficiency
311##############################
312
313module Efficiency ElectronTrackingEfficiency {
314 set InputArray ParticlePropagator/electrons
315 set OutputArray electrons
316
317
318 # Current full simulation with CLICdet provides for electrons:
319 set EfficiencyFormula {
320 (pt <= 1) * (0.000) +
321 (abs(eta) > 2.54) * (0.000) +
322 (energy >= 80) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
323 (energy >= 80) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
324 (energy >= 80) * (abs(eta) <= 2.34 ) * (1.000) +
325 (energy < 80 && energy >= 5) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
326 (energy < 80 && energy >= 5) * (abs(eta) <= 2.17) * (1.000) +
327 (energy < 5) * (abs(eta) > 2.34 ) * (0.000) +
328 (energy < 5) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
329 (energy < 5) * (abs(eta) <= 0.76) * (0.999)
330 }
331}
332
333##########################
334# Muon tracking efficiency
335##########################
336
337module Efficiency MuonTrackingEfficiency {
338 set InputArray ParticlePropagator/muons
339 set OutputArray muons
340
341 # Current full simulation with CLICdet provides for muons:
342 set EfficiencyFormula {
343 (pt < 1) * (0.000) +
344 (abs(eta) > 2.54) * (0.000) +
345 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 80) * (0.994) +
346 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 5 && energy < 80) * (0.996) +
347 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy < 5 ) * (0.996) +
348 (abs(eta) <= 2.44 ) * (energy >= 5 ) * (1.000) +
349 (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy < 5 && pt >=1 ) * (0.999) +
350 (abs(eta) <= 2.25 ) * (energy >= 1) * (1.000)
351
352
353 }
354}
355
356########################################
357# Momentum resolution for charged tracks
358########################################
359
360module MomentumSmearing ChargedHadronMomentumSmearing {
361 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
362 set OutputArray chargedHadrons
363
364
365 # Resolution given in dpT/pT.
366 # CLICdet internal studies
367 set ResolutionFormula {
368 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
369 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
370 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
371 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
372 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
373 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
374
375 }
376}
377
378###################################
379# Momentum resolution for electrons
380###################################
381
382module MomentumSmearing ElectronMomentumSmearing {
383 set InputArray ElectronTrackingEfficiency/electrons
384 set OutputArray electrons
385
386 # Resolution given in dpT/pT.
387 # CLICdet internal studies
388 set ResolutionFormula {
389 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
390 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
391 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
392 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
393 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
394 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
395
396 }
397}
398
399###############################
400# Momentum resolution for muons
401###############################
402
403module MomentumSmearing MuonMomentumSmearing {
404 set InputArray MuonTrackingEfficiency/muons
405 set OutputArray muons
406
407 # Resolution given in dpT/pT.
408
409 # CLICdet internal studies
410 set ResolutionFormula {
411
412 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
413 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
414 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt(9.81626e-06^2 * pt^2 + 0.00379895^2 ) +
415 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt(1.1959e-05^2 * pt^2 + 0.00242417^2 ) +
416 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
417 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
418
419 }
420}
421
422##############
423# Track merger
424##############
425
426module Merger TrackMerger {
427 # add InputArray InputArray
428 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
429 add InputArray ElectronMomentumSmearing/electrons
430 add InputArray MuonMomentumSmearing/muons
431 set OutputArray tracks
432}
433
434#############
435# ECAL
436#############
437
438module SimpleCalorimeter ECal {
439 set ParticleInputArray ParticlePropagator/stableParticles
440 set TrackInputArray TrackMerger/tracks
441
442 set TowerOutputArray ecalTowers
443 set EFlowTrackOutputArray eflowTracks
444 set EFlowTowerOutputArray eflowPhotons
445
446 set IsEcal true
447
448 set EnergyMin 0.5
449 set EnergySignificanceMin 1.0
450
451 set SmearTowerCenter true
452
453 set pi [expr {acos(-1)}]
454
455 # lists of the edges of each tower in eta and phi
456 # each list starts with the lower edge of the first tower
457 # the list ends with the higher edged of the last tower
458
459 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
460 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
461 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
462 #ECAL cell sizes always 5x5 mm^2
463
464 #barrel:
465 #dphi = 0.2 degree towers up to eta <=1.2
466 set PhiBins {}
467 for {set i -900} {$i <= 900} {incr i} {
468 add PhiBins [expr {$i * $pi/900.0 }]
469 }
470 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
471 for {set i -400} {$i <=400} {incr i} {
472 set eta [expr {$i * 0.003}]
473 add EtaPhiBins $eta $PhiBins
474 }
475
476 #endcaps:
477 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
478 set PhiBins {}
479 for {set i -225} {$i <= 225} {incr i} {
480 add PhiBins [expr {$i * $pi/225.}]
481 }
482 #deta=0.02 units for 1.2 < |eta| <=2.5
483 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
484 for {set i 1} {$i <=66} {incr i} {
485 set eta [expr {-2.52 + $i * 0.02}]
486 add EtaPhiBins $eta $PhiBins
487 }
488 #same for 1.2 to 2.5
489 for {set i 1} {$i <=66} {incr i} {
490 set eta [expr {1.18 + $i*0.02}]
491 add EtaPhiBins $eta $PhiBins
492 }
493
494 #plug:
495 #dphi = 1 degree for 2.5 < eta <=3
496 set PhiBins {}
497 for {set i -180} {$i <= 180} {incr i} {
498 add PhiBins [expr {$i * $pi/180.}]
499 }
500 # deta = 0.02 for 2.5 < |eta| <=3
501 # from -3 to -2.5, there will be 25 segments
502 for {set i 1} {$i <= 26} {incr i} {
503 set eta [expr {-3.02 + $i * 0.02}]
504 add EtaPhiBins $eta $PhiBins
505 }
506 #same for 2.5 to 3
507 for {set i 1} {$i <= 26} {incr i} {
508 set eta [expr {2.48 + $i*0.02}]
509 add EtaPhiBins $eta $PhiBins
510 }
511
512
513
514 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
515
516 add EnergyFraction {0} {0.0}
517 # energy fractions for e, gamma and pi0
518 add EnergyFraction {11} {1.0}
519 add EnergyFraction {22} {1.0}
520 add EnergyFraction {111} {1.0}
521 # energy fractions for muon, neutrinos and neutralinos
522 add EnergyFraction {12} {0.0}
523 add EnergyFraction {13} {0.0}
524 add EnergyFraction {14} {0.0}
525 add EnergyFraction {16} {0.0}
526 add EnergyFraction {1000022} {0.0}
527 add EnergyFraction {1000023} {0.0}
528 add EnergyFraction {1000025} {0.0}
529 add EnergyFraction {1000035} {0.0}
530 add EnergyFraction {1000045} {0.0}
531 # energy fractions for K0short and Lambda
532 add EnergyFraction {310} {0.3}
533 add EnergyFraction {3122} {0.3}
534
535 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
536 set ResolutionFormula {
537 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
538 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
539 (abs(eta) <= 3 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
540}
541
542#############
543# HCAL
544#############
545
546module SimpleCalorimeter HCal {
547 set ParticleInputArray ParticlePropagator/stableParticles
548 set TrackInputArray ECal/eflowTracks
549
550 set TowerOutputArray hcalTowers
551 set EFlowTrackOutputArray eflowTracks
552 set EFlowTowerOutputArray eflowNeutralHadrons
553
554 set IsEcal false
555
556 set EnergyMin 1.0
557 set EnergySignificanceMin 1.0
558
559 set SmearTowerCenter true
560
561 set pi [expr {acos(-1)}]
562
563 # lists of the edges of each tower in eta and phi
564 # each list starts with the lower edge of the first tower
565 # the list ends with the higher edged of the last tower
566
567
568 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
569 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
570 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
571 #HCAL cell sizes always 30x30 mm^2
572
573 #barrel and ring:
574 #dphi = 1 degree up to |eta| <=0.9
575 set PhiBins {}
576 for {set i -180} {$i <=180} {incr i} {
577 add PhiBins [expr {$i * $pi/180.0}]
578 }
579 #deta= 0.02 towers up to |eta| <=0.9
580 for {set i -45} {$i <=45} {incr i} {
581 set eta [expr {$i * 0.02}]
582 add EtaPhiBins $eta $PhiBins
583 }
584
585 #endcaps:
586 # dphi = 6 degree
587 set PhiBins {}
588 for {set i -30} {$i <=30} {incr i} {
589 add PhiBins [expr {$i * $pi/30.0}]
590 }
591 # deta =0.1 for 0.9 < |eta| <=3.0
592 #for -3.0 to -0.9, 21 segments
593 for {set i 1} {$i <=22} {incr i} {
594 set eta [expr {-3.1 + $i * 0.1}]
595 add EtaPhiBins $eta $PhiBins
596 }
597 #same for 0.9 to 3.0
598 for {set i 1} {$i <=22} {incr i} {
599 set eta [expr {0.8 + $i * 0.1 }]
600 add EtaPhiBins $eta $PhiBins
601 }
602
603 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
604 add EnergyFraction {0} {1.0}
605 # energy fractions for e, gamma and pi0
606 add EnergyFraction {11} {0.0}
607 add EnergyFraction {22} {0.0}
608 add EnergyFraction {111} {0.0}
609 # energy fractions for muon, neutrinos and neutralinos
610 add EnergyFraction {12} {0.0}
611 add EnergyFraction {13} {0.0}
612 add EnergyFraction {14} {0.0}
613 add EnergyFraction {16} {0.0}
614 add EnergyFraction {1000022} {0.0}
615 add EnergyFraction {1000023} {0.0}
616 add EnergyFraction {1000025} {0.0}
617 add EnergyFraction {1000035} {0.0}
618 add EnergyFraction {1000045} {0.0}
619 # energy fractions for K0short and Lambda
620 add EnergyFraction {310} {0.7}
621 add EnergyFraction {3122} {0.7}
622
623 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
624 #CLICdet internal studies
625 set ResolutionFormula {
626 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
627 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
628 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
629 (abs(eta)<=3 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
630 }
631
632}
633
634#################
635# Electron filter
636#################
637
638module PdgCodeFilter ElectronFilter {
639 set InputArray HCal/eflowTracks
640 set OutputArray electrons
641 set Invert true
642 add PdgCode {11}
643 add PdgCode {-11}
644}
645
646######################
647# ChargedHadronFilter
648######################
649
650module PdgCodeFilter ChargedHadronFilter {
651 set InputArray HCal/eflowTracks
652 set OutputArray chargedHadrons
653
654 add PdgCode {11}
655 add PdgCode {-11}
656 add PdgCode {13}
657 add PdgCode {-13}
658}
659
660
661
662###################################################
663# Tower Merger (in case not using e-flow algorithm)
664###################################################
665
666module Merger Calorimeter {
667 # add InputArray InputArray
668 add InputArray ECal/ecalTowers
669 add InputArray HCal/hcalTowers
670 set OutputArray towers
671}
672
673
674####################
675# Energy flow merger
676####################
677
678module Merger EFlowMerger {
679 # add InputArray InputArray
680 add InputArray HCal/eflowTracks
681 add InputArray ECal/eflowPhotons
682 add InputArray HCal/eflowNeutralHadrons
683 set OutputArray eflow
684}
685
686
687###################
688# Photon efficiency
689###################
690
691module Efficiency PhotonEfficiency {
692 set InputArray ECal/eflowPhotons
693 set OutputArray photons
694
695 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
696
697 # efficiency formula for photons
698 # current full simulation of CLICdet yields:
699 set EfficiencyFormula {
700 (energy < 2.0 ) * (0.000) +
701 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
702 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) }
703
704}
705
706##################
707# Photon isolation
708##################
709
710module Isolation PhotonIsolation {
711 set CandidateInputArray PhotonEfficiency/photons
712 set IsolationInputArray EFlowMerger/eflow
713
714 set OutputArray photons
715
716 set DeltaRMax 0.1
717
718 set PTMin 0.5
719
720 set PTRatioMax 0.2
721}
722
723#####################
724# Electron efficiency
725#####################
726
727module Efficiency ElectronEfficiency {
728 set InputArray ElectronFilter/electrons
729 set OutputArray electrons
730
731 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
732
733 set EfficiencyFormula {
734 (energy < 3.0) * ( 0.00 ) +
735 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95) * (0.58 ) +
736 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
737 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
738 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
739 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
740 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
741 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95) * ( 0.6 ) +
742 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
743 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
744 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
745 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
746 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
747 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95) * ( 0.6 ) +
748 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
749 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
750 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
751 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
752 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
753 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95) * (0.64 ) +
754 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
755 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
756 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
757 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
758 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
759 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95) * (0.64 ) +
760 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
761 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
762 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
763 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
764 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
765 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95) * (0.67 ) +
766 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
767 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
768 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
769 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
770 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
771 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95) * (0.68 ) +
772 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
773 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
774 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
775 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
776 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
777 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95) * ( 0.7 ) +
778 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
779 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
780 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
781 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
782 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
783 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95) * (0.68 ) +
784 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
785 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
786 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
787 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
788 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
789 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95) * ( 0.68 ) +
790 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
791 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
792 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
793 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
794 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
795 ( energy >=400 ) * (abs(eta) > 1.95) * (0.68 ) +
796 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
797 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
798 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
799 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
800 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
801 }
802}
803
804####################
805# Electron isolation
806####################
807
808module Isolation ElectronIsolation {
809 set CandidateInputArray ElectronEfficiency/electrons
810 set IsolationInputArray EFlowMerger/eflow
811
812 set OutputArray electrons
813
814 set DeltaRMax 0.1
815
816 set PTMin 0.5
817
818 set PTRatioMax 0.2
819}
820
821#################
822# Muon efficiency
823#################
824
825module Efficiency MuonEfficiency {
826 set InputArray MuonMomentumSmearing/muons
827 set OutputArray muons
828
829 # set EfficiencyFormula {efficiency as a function of eta and pt}
830
831 # efficiency formula for muons
832 # current full simulation of CLICdet yields:
833
834
835
836 set EfficiencyFormula {
837 (energy < 2.0 ) * (0.00) +
838 (energy>=2.0) * (0.999)
839
840 }
841}
842
843################
844# Muon isolation
845################
846
847module Isolation MuonIsolation {
848 set CandidateInputArray MuonEfficiency/muons
849 set IsolationInputArray EFlowMerger/eflow
850
851 set OutputArray muons
852
853 set DeltaRMax 0.1
854
855 set PTMin 0.5
856
857 set PTRatioMax 0.2
858}
859
860
861
862###################
863# Missing ET merger
864###################
865
866module Merger MissingET {
867 # add InputArray InputArray
868 add InputArray EFlowMerger/eflow
869 set MomentumOutputArray momentum
870}
871
872
873##################
874# Scalar HT merger
875##################
876
877module Merger ScalarHT {
878 # add InputArray InputArray
879 add InputArray EFlowMerger/eflow
880 set EnergyOutputArray energy
881}
882######################
883# EFlowFilter (UniqueObjectFinder)
884######################
885module UniqueObjectFinder EFlowFilter {
886 add InputArray PhotonIsolation/photons photons
887 add InputArray ElectronIsolation/electrons electrons
888 add InputArray MuonIsolation/muons muons
889 add InputArray EFlowMerger/eflow eflow
890}
891
892#################
893# Neutrino Filter
894#################
895
896module PdgCodeFilter NeutrinoFilter {
897
898 set InputArray Delphes/stableParticles
899 set OutputArray filteredParticles
900
901 set PTMin 0.0
902
903 add PdgCode {12}
904 add PdgCode {14}
905 add PdgCode {16}
906 add PdgCode {-12}
907 add PdgCode {-14}
908 add PdgCode {-16}
909
910}
911
912
913#####################
914# MC truth jet finder
915#####################
916
917module FastJetFinder GenJetFinder {
918 set InputArray NeutrinoFilter/filteredParticles
919
920 set OutputArray jets
921
922 # 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
923 set JetAlgorithm 9
924 set ParameterR 0.5
925
926 set JetPTMin 20.0
927}
928
929#########################
930# Gen Missing ET merger
931########################
932
933module Merger GenMissingET {
934 # add InputArray InputArray
935 add InputArray NeutrinoFilter/filteredParticles
936 set MomentumOutputArray momentum
937}
938
939
940
941############
942# Jet finder
943############
944
945module FastJetFinder FastJetFinderKt {
946 # set InputArray Calorimeter/towers
947 set InputArray EFlowMerger/eflow
948
949 set OutputArray KTjets
950
951 # 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
952 set JetAlgorithm 4
953 set ParameterR 0.5
954
955 set JetPTMin 20.0
956}
957
958################
959# Jet finder VLC
960################
961
962
963source CLIC/CLICdet_JetReco.tcl
964
965
966########################
967# Jet Flavor Association
968########################
969
970source CLIC/CLICdet_JetFlavorAssociation.tcl
971
972###########
973# b-tagging
974###########
975# based on CLICdp-Note-2014-002
976
977source CLIC/CLICdet_BTagging.tcl
978
979
980#############
981# tau-tagging
982#############
983# based on LCD-2010-009
984
985source CLIC/CLICdet_TauTagging.tcl
986
987
988
989##################
990# ROOT tree writer
991##################
992
993module TreeWriter TreeWriter {
994 # add Branch InputArray BranchName BranchClass
995 add Branch Delphes/allParticles Particle GenParticle
996
997 add Branch GenJetFinder/jets GenJet Jet
998
999 add Branch FastJetFinderKt/KTjets KTjet Jet
1000 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1001 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1002 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1003 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1004 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1005
1006 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1007 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1008 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1009 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1010 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1011
1012 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1013 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1014 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1015 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1016 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1017
1018 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1019 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1020 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1021 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1022 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1023
1024 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1025 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1026 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1027 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1028 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1029
1030 add Branch FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive VLCjetR05_inclusive Jet
1031 add Branch FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive VLCjetR07_inclusive Jet
1032 add Branch FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive VLCjetR10_inclusive Jet
1033 add Branch FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive VLCjetR12_inclusive Jet
1034 add Branch FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive VLCjetR15_inclusive Jet
1035
1036
1037
1038 ####
1039
1040 add Branch GenMissingET/momentum GenMissingET MissingET
1041
1042 add Branch TrackMerger/tracks Track Track
1043 add Branch Calorimeter/towers Tower Tower
1044
1045 add Branch HCal/eflowTracks EFlowTrack Track
1046 add Branch ECal/eflowPhotons EFlowPhoton Tower
1047 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1048
1049 add Branch EFlowFilter/photons Photon Photon
1050 add Branch EFlowFilter/electrons Electron Electron
1051 add Branch EFlowFilter/muons Muon Muon
1052
1053 add Branch MissingET/momentum MissingET MissingET
1054 add Branch ScalarHT/energy ScalarHT ScalarHT
1055}
1056
Note: See TracBrowser for help on using the repository browser.