Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ 283bed1

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

fixed electron, photon and muon efficiencies for low energy

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