Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ 3037765

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

added tracking efficiency numbers and momentum resolutions for e, mu, pions
  • Property mode set to 100644
File size: 31.5 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 JetEnergyScale
83
84 JetFlavorAssociation
85 BTagging
86
87 TauTagging
88
89 ScalarHT
90
91
92
93 TreeWriter
94}
95
96#################################
97# Propagate particles in cylinder
98#################################
99
100module ParticlePropagator ParticlePropagator {
101 set InputArray Delphes/stableParticles
102
103 set OutputArray stableParticles
104 set ChargedHadronOutputArray chargedHadrons
105 set ElectronOutputArray electrons
106 set MuonOutputArray muons
107
108 # radius of the magnetic field coverage in the calorimeter, in m
109 set Radius 1.5
110 # half-length of the magnetic field coverage in the calorimeter, in m
111 set HalfLength 2.31
112
113 # magnetic field, in T
114 set Bz 4.0
115}
116
117####################################
118# Charged hadron tracking efficiency
119####################################
120
121module Efficiency ChargedHadronTrackingEfficiency {
122 set InputArray ParticlePropagator/chargedHadrons
123 set OutputArray chargedHadrons
124 # Current full simulation with CLICdet provides for pions:
125
126 set EfficiencyFormula {
127(abs(eta) > 2.54) * (0.000) +
128(energy >= 100) * (abs(eta) < 2.54) * (1.000) +
129(energy < 100 && energy >= 10) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
130(energy < 100 && energy >= 10) * (abs(eta) <= 2.34) * (1.000) +
131(energy < 10 && energy >= 1) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) +
132(energy < 10 && energy >= 1) * (abs(eta) <= 0.55 ) * (1.000)
133 }
134}
135
136##############################
137# Electron tracking efficiency
138##############################
139
140module Efficiency ElectronTrackingEfficiency {
141 set InputArray ParticlePropagator/electrons
142 set OutputArray electrons
143
144
145 # Current full simulation with CLICdet provides for electrons:
146 set EfficiencyFormula {
147 (pt <= 1) * (0.000) +
148 (abs(eta) > 2.54) * (0.000) +
149 (energy > 100) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
150 (energy > 100) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
151 (energy > 100) * (abs(eta) <= 2.34 ) * (1.000) +
152 (energy <= 100 && energy > 10) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
153 (energy <= 100 && energy > 10) * (abs(eta) <= 2.17) * (1.000) +
154 (energy <= 10 && energy > 1) * (abs(eta) > 2.34 ) * (0.000) +
155 (energy <= 10 && energy > 1) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
156 (energy <= 10 && energy > 1) * (abs(eta) <= 0.76) * (0.999)
157 }
158}
159
160##########################
161# Muon tracking efficiency
162##########################
163
164module Efficiency MuonTrackingEfficiency {
165 set InputArray ParticlePropagator/muons
166 set OutputArray muons
167
168 # Current full simulation with CLICdet provides for muons:
169 set EfficiencyFormula {
170 (pt < 1) * (0.000)+
171 (abs(eta) > 2.54) * (0.0000) +
172 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 100) * (0.994) +
173 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 10 && energy < 100) * (0.997) +
174 (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 1 && energy < 10) * (0.996) +
175 (abs(eta) <= 2.44 ) * (energy >= 10) * (1.000) +
176 (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy >= 1 && energy < 10) * (0.999) +
177 (abs(eta) <= 2.25 ) * (energy >= 1) * (1.000)
178 }
179}
180
181########################################
182# Momentum resolution for charged tracks
183########################################
184
185module MomentumSmearing ChargedHadronMomentumSmearing {
186 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
187 set OutputArray chargedHadrons
188
189
190 # Resolution given in dpT/pT.
191 set ResolutionFormula {
192 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
193 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
194 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
195 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
196 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
197
198 }
199}
200
201###################################
202# Momentum resolution for electrons
203###################################
204
205module MomentumSmearing ElectronMomentumSmearing {
206 set InputArray ElectronTrackingEfficiency/electrons
207 set OutputArray electrons
208
209 # Resolution given in dpT/pT.
210 set ResolutionFormula {
211 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
212 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
213 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
214 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
215 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
216
217 }
218}
219
220###############################
221# Momentum resolution for muons
222###############################
223
224module MomentumSmearing MuonMomentumSmearing {
225 set InputArray MuonTrackingEfficiency/muons
226 set OutputArray muons
227
228 # Resolution given in dpT/pT.
229 set ResolutionFormula {
230
231 (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt(4.57439e-05^2 * pt^2* + 0.0149328^2 ) +
232 (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt(9.81626e-06^2 * pt^2* + 0.00379895^2 ) +
233 (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt(1.1959e-05^2 * pt^2* + 0.00242417^2 ) +
234 (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
235 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
236
237 }
238}
239
240##############
241# Track merger
242##############
243
244module Merger TrackMerger {
245 # add InputArray InputArray
246 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
247 add InputArray ElectronMomentumSmearing/electrons
248 add InputArray MuonMomentumSmearing/muons
249 set OutputArray tracks
250}
251
252#############
253# ECAL
254#############
255
256module SimpleCalorimeter ECal {
257 set ParticleInputArray ParticlePropagator/stableParticles
258 set TrackInputArray TrackMerger/tracks
259
260 set TowerOutputArray ecalTowers
261 set EFlowTrackOutputArray eflowTracks
262 set EFlowTowerOutputArray eflowPhotons
263
264 set IsEcal true
265
266 set EnergyMin 0.5
267 set EnergySignificanceMin 1.0
268
269 set SmearTowerCenter true
270
271 set pi [expr {acos(-1)}]
272
273 # lists of the edges of each tower in eta and phi
274 # each list starts with the lower edge of the first tower
275 # the list ends with the higher edged of the last tower
276
277 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
278 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
279 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
280 #ECAL cell sizes always 5x5 mm^2
281
282 #barrel:
283 #dphi = 0.2 degree towers up to eta <=1.2
284 set PhiBins {}
285 for {set i -900} {$i <= 900} {incr i} {
286 add PhiBins [expr {$i * $pi/900.0 }]
287 }
288 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
289 for {set i -400} {$i <=400} {incr i} {
290 set eta [expr {$i * 0.003}]
291 add EtaPhiBins $eta $PhiBins
292 }
293
294 #endcaps:
295 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
296 set PhiBins {}
297 for {set i -225} {$i <= 225} {incr i} {
298 add PhiBins [expr {$i * $pi/225.}]
299 }
300 #deta=0.02 units for 1.2 < |eta| <=2.5
301 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
302 for {set i 1} {$i <=66} {incr i} {
303 set eta [expr {-2.52 + $i * 0.02}]
304 add EtaPhiBins $eta $PhiBins
305 }
306 #same for 1.2 to 2.5
307 for {set i 1} {$i <=66} {incr i} {
308 set eta [expr {1.18 + $i*0.02}]
309 add EtaPhiBins $eta $PhiBins
310 }
311
312 #plug:
313 #dphi = 1 degree for 2.5 < eta <=3
314 set PhiBins {}
315 for {set i -180} {$i <= 180} {incr i} {
316 add PhiBins [expr {$i * $pi/180.}]
317 }
318 # deta = 0.02 for 2.5 < |eta| <=3
319 # from -3 to -2.5, there will be 25 segments
320 for {set i 1} {$i <= 26} {incr i} {
321 set eta [expr {-3.02 + $i * 0.02}]
322 add EtaPhiBins $eta $PhiBins
323 }
324 #same for 2.5 to 3
325 for {set i 1} {$i <= 26} {incr i} {
326 set eta [expr {2.48 + $i*0.02}]
327 add EtaPhiBins $eta $PhiBins
328 }
329
330
331
332 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
333
334 add EnergyFraction {0} {0.0}
335 # energy fractions for e, gamma and pi0
336 add EnergyFraction {11} {1.0}
337 add EnergyFraction {22} {1.0}
338 add EnergyFraction {111} {1.0}
339 # energy fractions for muon, neutrinos and neutralinos
340 add EnergyFraction {12} {0.0}
341 add EnergyFraction {13} {0.0}
342 add EnergyFraction {14} {0.0}
343 add EnergyFraction {16} {0.0}
344 add EnergyFraction {1000022} {0.0}
345 add EnergyFraction {1000023} {0.0}
346 add EnergyFraction {1000025} {0.0}
347 add EnergyFraction {1000035} {0.0}
348 add EnergyFraction {1000045} {0.0}
349 # energy fractions for K0short and Lambda
350 add EnergyFraction {310} {0.3}
351 add EnergyFraction {3122} {0.3}
352
353 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
354
355 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.15^2) }
356
357}
358
359#############
360# HCAL
361#############
362
363module SimpleCalorimeter HCal {
364 set ParticleInputArray ParticlePropagator/stableParticles
365 set TrackInputArray ECal/eflowTracks
366
367 set TowerOutputArray hcalTowers
368 set EFlowTrackOutputArray eflowTracks
369 set EFlowTowerOutputArray eflowNeutralHadrons
370
371 set IsEcal false
372
373 set EnergyMin 1.0
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
385 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
386 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
387 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
388 #HCAL cell sizes always 30x30 mm^2
389
390 #barrel and ring:
391 #dphi = 1 degree up to |eta| <=0.9
392 set PhiBins {}
393 for {set i -180} {$i <=180} {incr i} {
394 add PhiBins [expr {$i * $pi/180.0}]
395 }
396 #deta= 0.02 towers up to |eta| <=0.9
397 for {set i -45} {$i <=45} {incr i} {
398 set eta [expr {$i * 0.02}]
399 add EtaPhiBins $eta $PhiBins
400 }
401
402 #endcaps:
403 # dphi = 6 degree
404 set PhiBins {}
405 for {set i -30} {$i <=30} {incr i} {
406 add PhiBins [expr {$i * $pi/30.0}]
407 }
408 # deta =0.1 for 0.9 < |eta| <=3.5
409 #for -3.5 to -0.9, 26 segments
410 for {set i 1} {$i <=27} {incr i} {
411 set eta [expr {-3.6 + $i * 0.1}]
412 add EtaPhiBins $eta $PhiBins
413 }
414 #same for 0.9 to 3.5
415 for {set i 1} {$i <=27} {incr i} {
416 set eta [expr {0.8 + $i * 0.1 }]
417 add EtaPhiBins $eta $PhiBins
418 }
419
420 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
421 add EnergyFraction {0} {1.0}
422 # energy fractions for e, gamma and pi0
423 add EnergyFraction {11} {0.0}
424 add EnergyFraction {22} {0.0}
425 add EnergyFraction {111} {0.0}
426 # energy fractions for muon, neutrinos and neutralinos
427 add EnergyFraction {12} {0.0}
428 add EnergyFraction {13} {0.0}
429 add EnergyFraction {14} {0.0}
430 add EnergyFraction {16} {0.0}
431 add EnergyFraction {1000022} {0.0}
432 add EnergyFraction {1000023} {0.0}
433 add EnergyFraction {1000025} {0.0}
434 add EnergyFraction {1000035} {0.0}
435 add EnergyFraction {1000045} {0.0}
436 # energy fractions for K0short and Lambda
437 add EnergyFraction {310} {0.7}
438 add EnergyFraction {3122} {0.7}
439
440 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
441
442 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.015^2 + energy*0.50^2)}
443
444}
445
446#################
447# Electron filter
448#################
449
450module PdgCodeFilter ElectronFilter {
451 set InputArray HCal/eflowTracks
452 set OutputArray electrons
453 set Invert true
454 add PdgCode {11}
455 add PdgCode {-11}
456}
457
458######################
459# ChargedHadronFilter
460######################
461
462module PdgCodeFilter ChargedHadronFilter {
463 set InputArray HCal/eflowTracks
464 set OutputArray chargedHadrons
465
466 add PdgCode {11}
467 add PdgCode {-11}
468 add PdgCode {13}
469 add PdgCode {-13}
470}
471
472
473
474###################################################
475# Tower Merger (in case not using e-flow algorithm)
476###################################################
477
478module Merger Calorimeter {
479 # add InputArray InputArray
480 add InputArray ECal/ecalTowers
481 add InputArray HCal/hcalTowers
482 set OutputArray towers
483}
484
485
486####################
487# Energy flow merger
488####################
489
490module Merger EFlowMerger {
491 # add InputArray InputArray
492 add InputArray HCal/eflowTracks
493 add InputArray ECal/eflowPhotons
494 add InputArray HCal/eflowNeutralHadrons
495 set OutputArray eflow
496}
497
498##here
499
500###################
501# Photon efficiency
502###################
503
504module Efficiency PhotonEfficiency {
505 set InputArray ECal/eflowPhotons
506 set OutputArray photons
507
508 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
509
510 # efficiency formula for photons
511 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
512 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
513 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
514 (abs(eta) > 2.5) * (0.00)}
515}
516
517##################
518# Photon isolation
519##################
520
521module Isolation PhotonIsolation {
522 set CandidateInputArray PhotonEfficiency/photons
523 set IsolationInputArray EFlowMerger/eflow
524
525 set OutputArray photons
526
527 set DeltaRMax 0.5
528
529 set PTMin 0.5
530
531 set PTRatioMax 0.12
532}
533
534#####################
535# Electron efficiency
536#####################
537
538module Efficiency ElectronEfficiency {
539 set InputArray ElectronFilter/electrons
540 set OutputArray electrons
541
542 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
543
544 # efficiency formula for electrons
545 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
546 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
547 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
548 (abs(eta) > 2.5) * (0.00)}
549}
550
551####################
552# Electron isolation
553####################
554
555module Isolation ElectronIsolation {
556 set CandidateInputArray ElectronEfficiency/electrons
557 set IsolationInputArray EFlowMerger/eflow
558
559 set OutputArray electrons
560
561 set DeltaRMax 0.5
562
563 set PTMin 0.5
564
565 set PTRatioMax 0.12
566}
567
568#################
569# Muon efficiency
570#################
571
572module Efficiency MuonEfficiency {
573 set InputArray MuonMomentumSmearing/muons
574 set OutputArray muons
575
576 # set EfficiencyFormula {efficiency as a function of eta and pt}
577
578 # efficiency formula for muons
579 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
580 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
581 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
582 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
583 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
584 (abs(eta) > 2.4) * (0.00)}
585}
586
587################
588# Muon isolation
589################
590
591module Isolation MuonIsolation {
592 set CandidateInputArray MuonEfficiency/muons
593 set IsolationInputArray EFlowMerger/eflow
594
595 set OutputArray muons
596
597 set DeltaRMax 0.5
598
599 set PTMin 0.5
600
601 set PTRatioMax 0.25
602}
603
604
605
606###################
607# Missing ET merger
608###################
609
610module Merger MissingET {
611 # add InputArray InputArray
612 add InputArray EFlowMerger/eflow
613 set MomentumOutputArray momentum
614}
615
616
617##################
618# Scalar HT merger
619##################
620
621module Merger ScalarHT {
622 # add InputArray InputArray
623 add InputArray EFlowMerger/eflow
624 set EnergyOutputArray energy
625}
626######################
627# EFlowFilter (UniqueObjectFinder)
628######################
629module UniqueObjectFinder EFlowFilter {
630 add InputArray PhotonIsolation/photons photons
631 add InputArray ElectronIsolation/electrons electrons
632 add InputArray MuonIsolation/muons muons
633 add InputArray EFlowMerger/eflow eflow
634}
635
636
637#
638#module PdgCodeFilter EFlowFilter {
639# set InputArray EFlowMerger/eflow
640# set OutputArray eflow
641#
642# add PdgCode {11}
643# add PdgCode {-11}
644# add PdgCode {13}
645# add PdgCode {-13}
646#}
647
648#################
649# Neutrino Filter
650#################
651
652module PdgCodeFilter NeutrinoFilter {
653
654 set InputArray Delphes/stableParticles
655 set OutputArray filteredParticles
656
657 set PTMin 0.0
658
659 add PdgCode {12}
660 add PdgCode {14}
661 add PdgCode {16}
662 add PdgCode {-12}
663 add PdgCode {-14}
664 add PdgCode {-16}
665
666}
667
668
669#####################
670# MC truth jet finder
671#####################
672
673module FastJetFinder GenJetFinder {
674 set InputArray NeutrinoFilter/filteredParticles
675
676 set OutputArray jets
677
678 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
679 set JetAlgorithm 9
680 set ParameterR 0.5
681
682 set JetPTMin 20.0
683}
684
685#########################
686# Gen Missing ET merger
687########################
688
689module Merger GenMissingET {
690 # add InputArray InputArray
691 add InputArray NeutrinoFilter/filteredParticles
692 set MomentumOutputArray momentum
693}
694
695
696
697############
698# Jet finder
699############
700
701module FastJetFinder FastJetFinderKt {
702 # set InputArray Calorimeter/towers
703 set InputArray EFlowMerger/eflow
704
705 set OutputArray KTjets
706
707 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
708 set JetAlgorithm 4
709 set ParameterR 0.5
710
711 set JetPTMin 20.0
712}
713
714############
715# Jet finder VLC
716############
717#R05 N2
718module FastJetFinder FastJetFinderVLC_R05_N2 {
719 # set InputArray Calorimeter/towers
720 set InputArray EFlowFilter/eflow
721
722 set OutputArray VLCjetsR05N2
723
724 # 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
725 set NJets 2
726 set ExclusiveClustering true
727 set JetAlgorithm 9
728 set ParameterR 0.5
729 set Beta 1.0
730 set Gamma 1.0
731
732 set JetPTMin 20.0
733}
734#R05 N3
735module FastJetFinder FastJetFinderVLC_R05_N3 {
736 set InputArray EFlowFilter/eflow
737 set OutputArray VLCjetsR05N3
738 set NJets 3
739 set ExclusiveClustering true
740 set JetAlgorithm 9
741 set ParameterR 0.5
742 set Beta 1.0
743 set Gamma 1.0
744 set JetPTMin 20.0
745}
746#R05 N4
747module FastJetFinder FastJetFinderVLC_R05_N4 {
748 set InputArray EFlowFilter/eflow
749 set OutputArray VLCjetsR05N4
750 set NJets 4
751 set ExclusiveClustering true
752 set JetAlgorithm 9
753 set ParameterR 0.5
754 set Beta 1.0
755 set Gamma 1.0
756 set JetPTMin 20.0
757}
758#R05 N5
759module FastJetFinder FastJetFinderVLC_R05_N5 {
760 set InputArray EFlowFilter/eflow
761 set OutputArray VLCjetsR05N5
762 set NJets 5
763 set ExclusiveClustering true
764 set JetAlgorithm 9
765 set ParameterR 0.5
766 set Beta 1.0
767 set Gamma 1.0
768 set JetPTMin 20.0
769}
770#R05 N6
771module FastJetFinder FastJetFinderVLC_R05_N6 {
772 set InputArray EFlowFilter/eflow
773 set OutputArray VLCjetsR05N6
774 set NJets 6
775 set ExclusiveClustering true
776 set JetAlgorithm 9
777 set ParameterR 0.5
778 set Beta 1.0
779 set Gamma 1.0
780 set JetPTMin 20.0
781}
782#R07 N2
783module FastJetFinder FastJetFinderVLC_R07_N2 {
784 set InputArray EFlowFilter/eflow
785 set OutputArray VLCjetsR07N2
786 set NJets 2
787 set ExclusiveClustering true
788 set JetAlgorithm 9
789 set ParameterR 0.7
790 set Beta 1.0
791 set Gamma 1.0
792 set JetPTMin 20.0
793}
794#R07 N3
795module FastJetFinder FastJetFinderVLC_R07_N3 {
796 set InputArray EFlowFilter/eflow
797 set OutputArray VLCjetsR07N3
798 set NJets 3
799 set ExclusiveClustering true
800 set JetAlgorithm 9
801 set ParameterR 0.7
802 set Beta 1.0
803 set Gamma 1.0
804 set JetPTMin 20.0
805}
806#R07 N4
807module FastJetFinder FastJetFinderVLC_R07_N4 {
808 set InputArray EFlowFilter/eflow
809 set OutputArray VLCjetsR07N4
810 set NJets 4
811 set ExclusiveClustering true
812 set JetAlgorithm 9
813 set ParameterR 0.7
814 set Beta 1.0
815 set Gamma 1.0
816 set JetPTMin 20.0
817}
818#R07 N5
819module FastJetFinder FastJetFinderVLC_R07_N5 {
820 set InputArray EFlowFilter/eflow
821 set OutputArray VLCjetsR07N5
822 set NJets 5
823 set ExclusiveClustering true
824 set JetAlgorithm 9
825 set ParameterR 0.7
826 set Beta 1.0
827 set Gamma 1.0
828 set JetPTMin 20.0
829}
830#R07 N6
831module FastJetFinder FastJetFinderVLC_R07_N6 {
832 set InputArray EFlowFilter/eflow
833 set OutputArray VLCjetsR07N6
834 set NJets 6
835 set ExclusiveClustering true
836 set JetAlgorithm 9
837 set ParameterR 0.7
838 set Beta 1.0
839 set Gamma 1.0
840 set JetPTMin 20.0
841}
842
843#R10N2
844module FastJetFinder FastJetFinderVLC_R10_N2 {
845 set InputArray EFlowFilter/eflow
846 set OutputArray VLCjetsR10N2
847 set NJets 2
848 set ExclusiveClustering true
849 set JetAlgorithm 9
850 set ParameterR 1.0
851 set Beta 1.0
852 set Gamma 1.0
853 set JetPTMin 20.0
854}
855#R10 N3
856module FastJetFinder FastJetFinderVLC_R10_N3 {
857 set InputArray EFlowFilter/eflow
858 set OutputArray VLCjetsR10N3
859 set NJets 3
860 set ExclusiveClustering true
861 set JetAlgorithm 9
862 set ParameterR 1.0
863 set Beta 1.0
864 set Gamma 1.0
865 set JetPTMin 20.0
866}
867#R10 N4
868module FastJetFinder FastJetFinderVLC_R10_N4 {
869 set InputArray EFlowFilter/eflow
870 set OutputArray VLCjetsR10N4
871 set NJets 4
872 set ExclusiveClustering true
873 set JetAlgorithm 9
874 set ParameterR 1.0
875 set Beta 1.0
876 set Gamma 1.0
877 set JetPTMin 20.0
878}
879#R10 N5
880module FastJetFinder FastJetFinderVLC_R10_N5 {
881 set InputArray EFlowFilter/eflow
882 set OutputArray VLCjetsR10N5
883 set NJets 5
884 set ExclusiveClustering true
885 set JetAlgorithm 9
886 set ParameterR 1.0
887 set Beta 1.0
888 set Gamma 1.0
889 set JetPTMin 20.0
890}
891#R10 N6
892module FastJetFinder FastJetFinderVLC_R10_N6 {
893 set InputArray EFlowFilter/eflow
894 set OutputArray VLCjetsR10N6
895 set NJets 6
896 set ExclusiveClustering true
897 set JetAlgorithm 9
898 set ParameterR 1.0
899 set Beta 1.0
900 set Gamma 1.0
901 set JetPTMin 20.0
902}
903
904#R12 N2
905module FastJetFinder FastJetFinderVLC_R12_N2 {
906 set InputArray EFlowFilter/eflow
907 set OutputArray VLCjetsR12N2
908 set NJets 2
909 set ExclusiveClustering true
910 set JetAlgorithm 9
911 set ParameterR 1.2
912 set Beta 1.0
913 set Gamma 1.0
914 set JetPTMin 20.0
915}
916#R12 N3
917module FastJetFinder FastJetFinderVLC_R12_N3 {
918 set InputArray EFlowFilter/eflow
919 set OutputArray VLCjetsR12N3
920 set NJets 3
921 set ExclusiveClustering true
922 set JetAlgorithm 9
923 set ParameterR 1.2
924 set Beta 1.0
925 set Gamma 1.0
926 set JetPTMin 20.0
927}
928#R12 N4
929module FastJetFinder FastJetFinderVLC_R12_N4 {
930 set InputArray EFlowFilter/eflow
931 set OutputArray VLCjetsR12N4
932 set NJets 4
933 set ExclusiveClustering true
934 set JetAlgorithm 9
935 set ParameterR 1.2
936 set Beta 1.0
937 set Gamma 1.0
938 set JetPTMin 20.0
939}
940#R12 N5
941module FastJetFinder FastJetFinderVLC_R12_N5 {
942 set InputArray EFlowFilter/eflow
943 set OutputArray VLCjetsR12N5
944 set NJets 5
945 set ExclusiveClustering true
946 set JetAlgorithm 9
947 set ParameterR 1.2
948 set Beta 1.0
949 set Gamma 1.0
950 set JetPTMin 20.0
951}
952#R12 N6
953module FastJetFinder FastJetFinderVLC_R12_N6 {
954 set InputArray EFlowFilter/eflow
955 set OutputArray VLCjetsR12N6
956 set NJets 6
957 set ExclusiveClustering true
958 set JetAlgorithm 9
959 set ParameterR 1.2
960 set Beta 1.0
961 set Gamma 1.0
962 set JetPTMin 20.0
963}
964
965
966#R15 N2
967module FastJetFinder FastJetFinderVLC_R15_N2 {
968 set InputArray EFlowFilter/eflow
969 set OutputArray VLCjetsR15N2
970 set NJets 2
971 set ExclusiveClustering true
972 set JetAlgorithm 9
973 set ParameterR 1.5
974 set Beta 1.0
975 set Gamma 1.0
976 set JetPTMin 20.0
977}
978#R15 N3
979module FastJetFinder FastJetFinderVLC_R15_N3 {
980 set InputArray EFlowFilter/eflow
981 set OutputArray VLCjetsR15N3
982 set NJets 3
983 set ExclusiveClustering true
984 set JetAlgorithm 9
985 set ParameterR 1.5
986 set Beta 1.0
987 set Gamma 1.0
988 set JetPTMin 20.0
989}
990#R15 N4
991module FastJetFinder FastJetFinderVLC_R15_N4 {
992 set InputArray EFlowFilter/eflow
993 set OutputArray VLCjetsR15N4
994 set NJets 4
995 set ExclusiveClustering true
996 set JetAlgorithm 9
997 set ParameterR 1.5
998 set Beta 1.0
999 set Gamma 1.0
1000 set JetPTMin 20.0
1001}
1002#R15 N5
1003module FastJetFinder FastJetFinderVLC_R15_N5 {
1004 set InputArray EFlowFilter/eflow
1005 set OutputArray VLCjetsR15N5
1006 set NJets 5
1007 set ExclusiveClustering true
1008 set JetAlgorithm 9
1009 set ParameterR 1.5
1010 set Beta 1.0
1011 set Gamma 1.0
1012 set JetPTMin 20.0
1013}
1014#R15 N6
1015module FastJetFinder FastJetFinderVLC_R15_N6 {
1016 set InputArray EFlowFilter/eflow
1017 set OutputArray VLCjetsR15N6
1018 set NJets 6
1019 set ExclusiveClustering true
1020 set JetAlgorithm 9
1021 set ParameterR 1.5
1022 set Beta 1.0
1023 set Gamma 1.0
1024 set JetPTMin 20.0
1025}
1026##################
1027# Jet Energy Scale
1028##################
1029
1030module EnergyScale JetEnergyScale {
1031 set InputArray FastJetFinderVLC_R10_N4/VLCjetsR10N4
1032 set OutputArray jets
1033
1034 # scale formula for jets
1035 set ScaleFormula {1.00}
1036}
1037
1038
1039########################
1040# Jet Flavor Association
1041########################
1042
1043module JetFlavorAssociation JetFlavorAssociation {
1044
1045 set PartonInputArray Delphes/partons
1046 set ParticleInputArray Delphes/allParticles
1047 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1048 set JetInputArray JetEnergyScale/jets
1049
1050 set DeltaR 0.5
1051 set PartonPTMin 1.0
1052 set PartonEtaMax 2.5
1053
1054}
1055
1056
1057###########
1058# b-tagging
1059###########
1060
1061module BTagging BTagging {
1062 set JetInputArray JetEnergyScale/jets
1063
1064 set BitNumber 0
1065
1066 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1067 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
1068 # gluon's PDG code has the lowest priority
1069
1070 # based on arXiv:1211.4462 (CMS)
1071
1072 # default efficiency formula (misidentification rate)
1073 add EfficiencyFormula {0} {0.01+0.000038*pt}
1074
1075 # efficiency formula for c-jets (misidentification rate)
1076 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
1077
1078 # efficiency formula for b-jets
1079 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
1080}
1081
1082#############
1083# tau-tagging
1084#############
1085
1086
1087module TauTagging TauTagging {
1088 set ParticleInputArray Delphes/allParticles
1089 set PartonInputArray Delphes/partons
1090 set JetInputArray JetEnergyScale/jets
1091
1092 set DeltaR 0.5
1093
1094 set TauPTMin 1.0
1095
1096 set TauEtaMax 4.0
1097
1098 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1099
1100 # default efficiency formula (misidentification rate)
1101 add EfficiencyFormula {0} {0.001}
1102 # efficiency formula for tau-jets
1103 add EfficiencyFormula {15} {0.4}
1104}
1105
1106#####################################################
1107# Find uniquely identified photons/electrons/tau/jets
1108#####################################################
1109
1110#module UniqueObjectFinder UniqueObjectFinder {
1111# # earlier arrays take precedence over later ones
1112# # add InputArray InputArray OutputArray
1113# add InputArray PhotonIsolation/photons photons
1114# add InputArray ElectronIsolation/electrons electrons
1115# add InputArray MuonIsolation/muons muons
1116# add InputArray JetEnergyScale/jets jets
1117####not necessary any more, since these objects have already been removed if necessary??
1118#}
1119#
1120
1121##################
1122# ROOT tree writer
1123##################
1124
1125module TreeWriter TreeWriter {
1126 # add Branch InputArray BranchName BranchClass
1127 add Branch Delphes/allParticles Particle GenParticle
1128
1129 add Branch GenJetFinder/jets GenJet Jet
1130
1131 add Branch FastJetFinderKt/KTjets KTjet Jet
1132 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1133 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1134 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1135 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1136 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1137
1138 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1139 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1140 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1141 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1142 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1143
1144 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1145 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1146 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1147 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1148 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1149
1150 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1151 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1152 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1153 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1154 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1155
1156 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1157 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1158 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1159 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1160 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1161
1162
1163 add Branch GenMissingET/momentum GenMissingET MissingET
1164
1165 add Branch TrackMerger/tracks Track Track
1166 add Branch Calorimeter/towers Tower Tower
1167
1168 add Branch HCal/eflowTracks EFlowTrack Track
1169 add Branch ECal/eflowPhotons EFlowPhoton Tower
1170 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1171
1172# add Branch UniqueObjectFinder/photons Photon Photon
1173# add Branch UniqueObjectFinder/electrons Electron Electron
1174# add Branch UniqueObjectFinder/muons Muon Muon
1175# add Branch UniqueObjectFinder/jets Jet Jet
1176#
1177
1178# add Branch PhotonIsolation/photons Photon Photon
1179# add Branch ElectronIsolation/electrons Electron Electron
1180# add Branch MuonIsolation/muons Muon Muon
1181
1182 add Branch EFlowFilter/photons Photon Photon
1183 add Branch EFlowFilter/electrons Electron Electron
1184 add Branch EFlowFilter/muons Muon Muon
1185
1186
1187
1188 add Branch MissingET/momentum MissingET MissingET
1189 add Branch ScalarHT/energy ScalarHT ScalarHT
1190}
1191
Note: See TracBrowser for help on using the repository browser.