Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ 5953e02

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

added also electron efficiencies (preliminary)
  • Property mode set to 100644
File size: 26.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 >= 100) * (abs(eta) < 2.54) * (1.000) +
232(energy < 100 && energy >= 10) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
233(energy < 100 && energy >= 10) * (abs(eta) <= 2.34) * (1.000) +
234(energy < 10 && energy >= 1) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) +
235(energy < 10 && energy >= 1) * (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 > 100) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
253 (energy > 100) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
254 (energy > 100) * (abs(eta) <= 2.34 ) * (1.000) +
255 (energy <= 100 && energy > 10) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
256 (energy <= 100 && energy > 10) * (abs(eta) <= 2.17) * (1.000) +
257 (energy <= 10 && energy > 1) * (abs(eta) > 2.34 ) * (0.000) +
258 (energy <= 10 && energy > 1) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
259 (energy <= 10 && energy > 1) * (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.0000) +
275 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 100) * (0.994) +
276 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 10 && energy < 100) * (0.997) +
277 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 1 && energy < 10) * (0.996) +
278 (abs(eta) <= 2.44 ) * (energy >= 10) * (1.000) +
279 (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy >= 1 && energy < 10) * (0.999) +
280 (abs(eta) <= 2.25 ) * (energy >= 1) * (1.000)
281 }
282}
283
284########################################
285# Momentum resolution for charged tracks
286########################################
287
288module MomentumSmearing ChargedHadronMomentumSmearing {
289 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
290 set OutputArray chargedHadrons
291
292
293 # Resolution given in dpT/pT.
294 set ResolutionFormula {
295 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
296 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
297 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
298 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
299 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
300
301 }
302}
303
304###################################
305# Momentum resolution for electrons
306###################################
307
308module MomentumSmearing ElectronMomentumSmearing {
309 set InputArray ElectronTrackingEfficiency/electrons
310 set OutputArray electrons
311
312 # Resolution given in dpT/pT.
313 set ResolutionFormula {
314 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
315 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
316 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
317 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
318 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
319
320 }
321}
322
323###############################
324# Momentum resolution for muons
325###############################
326
327module MomentumSmearing MuonMomentumSmearing {
328 set InputArray MuonTrackingEfficiency/muons
329 set OutputArray muons
330
331 # Resolution given in dpT/pT.
332 set ResolutionFormula {
333
334 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt(4.57439e-05^2 * pt^2* + 0.0149328^2 ) +
335 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt(9.81626e-06^2 * pt^2* + 0.00379895^2 ) +
336 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt(1.1959e-05^2 * pt^2* + 0.00242417^2 ) +
337 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
338 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
339
340 }
341}
342
343##############
344# Track merger
345##############
346
347module Merger TrackMerger {
348 # add InputArray InputArray
349 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
350 add InputArray ElectronMomentumSmearing/electrons
351 add InputArray MuonMomentumSmearing/muons
352 set OutputArray tracks
353}
354
355#############
356# ECAL
357#############
358
359module SimpleCalorimeter ECal {
360 set ParticleInputArray ParticlePropagator/stableParticles
361 set TrackInputArray TrackMerger/tracks
362
363 set TowerOutputArray ecalTowers
364 set EFlowTrackOutputArray eflowTracks
365 set EFlowTowerOutputArray eflowPhotons
366
367 set IsEcal true
368
369 set EnergyMin 0.5
370 set EnergySignificanceMin 1.0
371
372 set SmearTowerCenter true
373
374 set pi [expr {acos(-1)}]
375
376 # lists of the edges of each tower in eta and phi
377 # each list starts with the lower edge of the first tower
378 # the list ends with the higher edged of the last tower
379
380 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
381 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
382 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
383 #ECAL cell sizes always 5x5 mm^2
384
385 #barrel:
386 #dphi = 0.2 degree towers up to eta <=1.2
387 set PhiBins {}
388 for {set i -900} {$i <= 900} {incr i} {
389 add PhiBins [expr {$i * $pi/900.0 }]
390 }
391 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
392 for {set i -400} {$i <=400} {incr i} {
393 set eta [expr {$i * 0.003}]
394 add EtaPhiBins $eta $PhiBins
395 }
396
397 #endcaps:
398 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
399 set PhiBins {}
400 for {set i -225} {$i <= 225} {incr i} {
401 add PhiBins [expr {$i * $pi/225.}]
402 }
403 #deta=0.02 units for 1.2 < |eta| <=2.5
404 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
405 for {set i 1} {$i <=66} {incr i} {
406 set eta [expr {-2.52 + $i * 0.02}]
407 add EtaPhiBins $eta $PhiBins
408 }
409 #same for 1.2 to 2.5
410 for {set i 1} {$i <=66} {incr i} {
411 set eta [expr {1.18 + $i*0.02}]
412 add EtaPhiBins $eta $PhiBins
413 }
414
415 #plug:
416 #dphi = 1 degree for 2.5 < eta <=3
417 set PhiBins {}
418 for {set i -180} {$i <= 180} {incr i} {
419 add PhiBins [expr {$i * $pi/180.}]
420 }
421 # deta = 0.02 for 2.5 < |eta| <=3
422 # from -3 to -2.5, there will be 25 segments
423 for {set i 1} {$i <= 26} {incr i} {
424 set eta [expr {-3.02 + $i * 0.02}]
425 add EtaPhiBins $eta $PhiBins
426 }
427 #same for 2.5 to 3
428 for {set i 1} {$i <= 26} {incr i} {
429 set eta [expr {2.48 + $i*0.02}]
430 add EtaPhiBins $eta $PhiBins
431 }
432
433
434
435 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
436
437 add EnergyFraction {0} {0.0}
438 # energy fractions for e, gamma and pi0
439 add EnergyFraction {11} {1.0}
440 add EnergyFraction {22} {1.0}
441 add EnergyFraction {111} {1.0}
442 # energy fractions for muon, neutrinos and neutralinos
443 add EnergyFraction {12} {0.0}
444 add EnergyFraction {13} {0.0}
445 add EnergyFraction {14} {0.0}
446 add EnergyFraction {16} {0.0}
447 add EnergyFraction {1000022} {0.0}
448 add EnergyFraction {1000023} {0.0}
449 add EnergyFraction {1000025} {0.0}
450 add EnergyFraction {1000035} {0.0}
451 add EnergyFraction {1000045} {0.0}
452 # energy fractions for K0short and Lambda
453 add EnergyFraction {310} {0.3}
454 add EnergyFraction {3122} {0.3}
455
456 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
457
458 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.15^2) }
459
460}
461
462#############
463# HCAL
464#############
465
466module SimpleCalorimeter HCal {
467 set ParticleInputArray ParticlePropagator/stableParticles
468 set TrackInputArray ECal/eflowTracks
469
470 set TowerOutputArray hcalTowers
471 set EFlowTrackOutputArray eflowTracks
472 set EFlowTowerOutputArray eflowNeutralHadrons
473
474 set IsEcal false
475
476 set EnergyMin 1.0
477 set EnergySignificanceMin 1.0
478
479 set SmearTowerCenter true
480
481 set pi [expr {acos(-1)}]
482
483 # lists of the edges of each tower in eta and phi
484 # each list starts with the lower edge of the first tower
485 # the list ends with the higher edged of the last tower
486
487
488 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
489 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
490 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
491 #HCAL cell sizes always 30x30 mm^2
492
493 #barrel and ring:
494 #dphi = 1 degree up to |eta| <=0.9
495 set PhiBins {}
496 for {set i -180} {$i <=180} {incr i} {
497 add PhiBins [expr {$i * $pi/180.0}]
498 }
499 #deta= 0.02 towers up to |eta| <=0.9
500 for {set i -45} {$i <=45} {incr i} {
501 set eta [expr {$i * 0.02}]
502 add EtaPhiBins $eta $PhiBins
503 }
504
505 #endcaps:
506 # dphi = 6 degree
507 set PhiBins {}
508 for {set i -30} {$i <=30} {incr i} {
509 add PhiBins [expr {$i * $pi/30.0}]
510 }
511 # deta =0.1 for 0.9 < |eta| <=3.5
512 #for -3.5 to -0.9, 26 segments
513 for {set i 1} {$i <=27} {incr i} {
514 set eta [expr {-3.6 + $i * 0.1}]
515 add EtaPhiBins $eta $PhiBins
516 }
517 #same for 0.9 to 3.5
518 for {set i 1} {$i <=27} {incr i} {
519 set eta [expr {0.8 + $i * 0.1 }]
520 add EtaPhiBins $eta $PhiBins
521 }
522
523 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
524 add EnergyFraction {0} {1.0}
525 # energy fractions for e, gamma and pi0
526 add EnergyFraction {11} {0.0}
527 add EnergyFraction {22} {0.0}
528 add EnergyFraction {111} {0.0}
529 # energy fractions for muon, neutrinos and neutralinos
530 add EnergyFraction {12} {0.0}
531 add EnergyFraction {13} {0.0}
532 add EnergyFraction {14} {0.0}
533 add EnergyFraction {16} {0.0}
534 add EnergyFraction {1000022} {0.0}
535 add EnergyFraction {1000023} {0.0}
536 add EnergyFraction {1000025} {0.0}
537 add EnergyFraction {1000035} {0.0}
538 add EnergyFraction {1000045} {0.0}
539 # energy fractions for K0short and Lambda
540 add EnergyFraction {310} {0.7}
541 add EnergyFraction {3122} {0.7}
542
543 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
544
545 # set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.015^2 + energy*0.50^2)}
546 set ResolutionFormula { (abs(eta)< 0.8) * sqrt(1.48^2 + 0.251^2*energy + 0.0537^2*energy^2) +
547 (abs(eta)>= 0.8 && abs(eta)<=3) * sqrt( 1.09^2 + 0.321^2*energy + 0.0517^2*energy^2 )}
548
549}
550
551#################
552# Electron filter
553#################
554
555module PdgCodeFilter ElectronFilter {
556 set InputArray HCal/eflowTracks
557 set OutputArray electrons
558 set Invert true
559 add PdgCode {11}
560 add PdgCode {-11}
561}
562
563######################
564# ChargedHadronFilter
565######################
566
567module PdgCodeFilter ChargedHadronFilter {
568 set InputArray HCal/eflowTracks
569 set OutputArray chargedHadrons
570
571 add PdgCode {11}
572 add PdgCode {-11}
573 add PdgCode {13}
574 add PdgCode {-13}
575}
576
577
578
579###################################################
580# Tower Merger (in case not using e-flow algorithm)
581###################################################
582
583module Merger Calorimeter {
584 # add InputArray InputArray
585 add InputArray ECal/ecalTowers
586 add InputArray HCal/hcalTowers
587 set OutputArray towers
588}
589
590
591####################
592# Energy flow merger
593####################
594
595module Merger EFlowMerger {
596 # add InputArray InputArray
597 add InputArray HCal/eflowTracks
598 add InputArray ECal/eflowPhotons
599 add InputArray HCal/eflowNeutralHadrons
600 set OutputArray eflow
601}
602
603##here
604
605###################
606# Photon efficiency
607###################
608
609module Efficiency PhotonEfficiency {
610 set InputArray ECal/eflowPhotons
611 set OutputArray photons
612
613 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
614
615 # efficiency formula for photons
616 # including converted photons
617 #0.94 for |costheta| < 0.6, i.e. |eta|< 0.7
618 #0.9 for |costheta| > 0.6, i.e. |eta|>0.7 but smaller than 3
619
620
621 set EfficiencyFormula {(abs(eta) < 0.7)*(0.94) + (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) }
622 #old:
623 #set EfficiencyFormula { (pt <= 10.0) * (0.00) +
624 #(abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
625 #(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
626 #(abs(eta) > 2.5) * (0.00)}
627}
628
629##################
630# Photon isolation
631##################
632
633module Isolation PhotonIsolation {
634 set CandidateInputArray PhotonEfficiency/photons
635 set IsolationInputArray EFlowMerger/eflow
636
637 set OutputArray photons
638
639 set DeltaRMax 0.5
640
641 set PTMin 0.5
642
643 set PTRatioMax 0.12
644}
645
646#####################
647# Electron efficiency
648#####################
649
650module Efficiency ElectronEfficiency {
651 set InputArray ElectronFilter/electrons
652 set OutputArray electrons
653
654 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
655 set EfficiencyFormula {(abs( eta)<0.7 )*( 0.89)+
656 (abs(eta) > 0.7 && abs(eta) < 3)*( 0.8)}
657
658 # efficiency formula for electrons
659# set EfficiencyFormula { (pt <= 10.0) * (0.00) +
660# (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
661# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
662# (abs(eta) > 2.5) * (0.00)}
663}
664
665####################
666# Electron isolation
667####################
668
669module Isolation ElectronIsolation {
670 set CandidateInputArray ElectronEfficiency/electrons
671 set IsolationInputArray EFlowMerger/eflow
672
673 set OutputArray electrons
674
675 set DeltaRMax 0.5
676
677 set PTMin 0.5
678
679 set PTRatioMax 0.12
680}
681
682#################
683# Muon efficiency
684#################
685
686module Efficiency MuonEfficiency {
687 set InputArray MuonMomentumSmearing/muons
688 set OutputArray muons
689
690 # set EfficiencyFormula {efficiency as a function of eta and pt}
691
692 # efficiency formula for muons
693 set EfficiencyFormula {(energy < 50 )*(0.98) + (energy>=50) *(0.999)}
694 #set EfficiencyFormula { (pt <= 10.0) * (0.00) +
695 #(abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
696 #(abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
697 #(abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
698 #(abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
699 #(abs(eta) > 2.4) * (0.00)}
700}
701
702################
703# Muon isolation
704################
705
706module Isolation MuonIsolation {
707 set CandidateInputArray MuonEfficiency/muons
708 set IsolationInputArray EFlowMerger/eflow
709
710 set OutputArray muons
711
712 set DeltaRMax 0.5
713
714 set PTMin 0.5
715
716 set PTRatioMax 0.25
717}
718
719
720
721###################
722# Missing ET merger
723###################
724
725module Merger MissingET {
726 # add InputArray InputArray
727 add InputArray EFlowMerger/eflow
728 set MomentumOutputArray momentum
729}
730
731
732##################
733# Scalar HT merger
734##################
735
736module Merger ScalarHT {
737 # add InputArray InputArray
738 add InputArray EFlowMerger/eflow
739 set EnergyOutputArray energy
740}
741######################
742# EFlowFilter (UniqueObjectFinder)
743######################
744module UniqueObjectFinder EFlowFilter {
745 add InputArray PhotonIsolation/photons photons
746 add InputArray ElectronIsolation/electrons electrons
747 add InputArray MuonIsolation/muons muons
748 add InputArray EFlowMerger/eflow eflow
749}
750
751
752#
753#module PdgCodeFilter EFlowFilter {
754# set InputArray EFlowMerger/eflow
755# set OutputArray eflow
756#
757# add PdgCode {11}
758# add PdgCode {-11}
759# add PdgCode {13}
760# add PdgCode {-13}
761#}
762
763#################
764# Neutrino Filter
765#################
766
767module PdgCodeFilter NeutrinoFilter {
768
769 set InputArray Delphes/stableParticles
770 set OutputArray filteredParticles
771
772 set PTMin 0.0
773
774 add PdgCode {12}
775 add PdgCode {14}
776 add PdgCode {16}
777 add PdgCode {-12}
778 add PdgCode {-14}
779 add PdgCode {-16}
780
781}
782
783
784#####################
785# MC truth jet finder
786#####################
787
788module FastJetFinder GenJetFinder {
789 set InputArray NeutrinoFilter/filteredParticles
790
791 set OutputArray jets
792
793 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
794 set JetAlgorithm 9
795 set ParameterR 0.5
796
797 set JetPTMin 20.0
798}
799
800#########################
801# Gen Missing ET merger
802########################
803
804module Merger GenMissingET {
805 # add InputArray InputArray
806 add InputArray NeutrinoFilter/filteredParticles
807 set MomentumOutputArray momentum
808}
809
810
811
812############
813# Jet finder
814############
815
816module FastJetFinder FastJetFinderKt {
817 # set InputArray Calorimeter/towers
818 set InputArray EFlowMerger/eflow
819
820 set OutputArray KTjets
821
822 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
823 set JetAlgorithm 4
824 set ParameterR 0.5
825
826 set JetPTMin 20.0
827}
828
829source CLICdet_JetReco.tcl
830
831
832###################
833## Jet Energy Scale
834###################
835#
836#module EnergyScale JetEnergyScale {
837# set InputArray FastJetFinderVLC_R10_N4/VLCjetsR10N4
838# set OutputArray jets
839#
840# # scale formula for jets
841# set ScaleFormula {1.00}
842#}
843
844
845########################
846# Jet Flavor Association
847########################
848
849source CLICdet_JetFlavorAssociation.tcl
850
851###########
852# b-tagging
853###########
854# based on CLICdp-Note-2014-002
855
856source CLICdet_BTagging.tcl
857
858
859#############
860# tau-tagging
861#############
862
863
864module TauTagging TauTagging_R05N2 {
865 set ParticleInputArray Delphes/allParticles
866 set PartonInputArray Delphes/partons
867 set JetInputArray FastJetFinderVLC_R05_N2/VLCjetsR05N2
868
869 set DeltaR 0.5
870
871 set TauPTMin 1.0
872
873 set TauEtaMax 4.0
874
875 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
876
877 # default efficiency formula (misidentification rate)
878 add EfficiencyFormula {0} {0.001}
879 # efficiency formula for tau-jets
880 add EfficiencyFormula {15} {0.4}
881}
882
883#####################################################
884# Find uniquely identified photons/electrons/tau/jets
885#####################################################
886
887#module UniqueObjectFinder UniqueObjectFinder {
888# # earlier arrays take precedence over later ones
889# # add InputArray InputArray OutputArray
890# add InputArray PhotonIsolation/photons photons
891# add InputArray ElectronIsolation/electrons electrons
892# add InputArray MuonIsolation/muons muons
893# add InputArray JetEnergyScale/jets jets
894####not necessary any more, since these objects have already been removed if necessary??
895#}
896#
897
898##################
899# ROOT tree writer
900##################
901
902module TreeWriter TreeWriter {
903 # add Branch InputArray BranchName BranchClass
904 add Branch Delphes/allParticles Particle GenParticle
905
906 add Branch GenJetFinder/jets GenJet Jet
907
908 add Branch FastJetFinderKt/KTjets KTjet Jet
909 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
910 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
911 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
912 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
913 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
914
915 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
916 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
917 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
918 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
919 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
920
921 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
922 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
923 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
924 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
925 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
926
927 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
928 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
929 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
930 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
931 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
932
933 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
934 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
935 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
936 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
937 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
938
939# add Branch BTaggingWP50/jets BTagged Jet
940# add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 BTagged Jet
941
942 add Branch GenMissingET/momentum GenMissingET MissingET
943
944 add Branch TrackMerger/tracks Track Track
945 add Branch Calorimeter/towers Tower Tower
946
947 add Branch HCal/eflowTracks EFlowTrack Track
948 add Branch ECal/eflowPhotons EFlowPhoton Tower
949 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
950
951# add Branch UniqueObjectFinder/photons Photon Photon
952# add Branch UniqueObjectFinder/electrons Electron Electron
953# add Branch UniqueObjectFinder/muons Muon Muon
954# add Branch UniqueObjectFinder/jets Jet Jet
955#
956
957# add Branch PhotonIsolation/photons Photon Photon
958# add Branch ElectronIsolation/electrons Electron Electron
959# add Branch MuonIsolation/muons Muon Muon
960
961 add Branch EFlowFilter/photons Photon Photon
962 add Branch EFlowFilter/electrons Electron Electron
963 add Branch EFlowFilter/muons Muon Muon
964
965
966
967 add Branch MissingET/momentum MissingET MissingET
968 add Branch ScalarHT/energy ScalarHT ScalarHT
969}
970
Note: See TracBrowser for help on using the repository browser.