Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ d9b8093

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

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