Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ f521c11

ImprovedOutputFile Timing dual_readout llp
Last change on this file since f521c11 was f521c11, checked in by Ulrike Schnoor <schnooru@…>, 7 years ago

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