Fork me on GitHub

Ticket #1398: cms_fake.tcl

File cms_fake.tcl, 23.4 KB (added by padhan_rojalin, 5 years ago)

CMS card that I am using.

Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23 EFlowFilter
24
25 PhotonEfficiency
26 PhotonIsolation
27
28 ElectronFilter
29 ElectronEfficiency
30 ElectronIsolation
31
32 ChargedHadronFilter
33
34 MuonEfficiency
35 MuonIsolation
36
37 MissingET
38
39 NeutrinoFilter
40 GenJetFinder
41 GenMissingET
42
43 FastJetFinder
44 JetFakeParticle
45 ElectronMerger
46 MuonMerger
47 PhotonMerger
48 FatJetFinder
49
50 JetEnergyScale
51
52 JetFlavorAssociation
53
54 BTagging
55 TauTagging
56
57 UniqueObjectFinder
58
59 ScalarHT
60
61 TreeWriter
62}
63
64#################################
65# Propagate particles in cylinder
66#################################
67
68module ParticlePropagator ParticlePropagator {
69 set InputArray Delphes/stableParticles
70
71 set OutputArray stableParticles
72 set ChargedHadronOutputArray chargedHadrons
73 set ElectronOutputArray electrons
74 set MuonOutputArray muons
75
76 # radius of the magnetic field coverage, in m
77 set Radius 1.29
78 # half-length of the magnetic field coverage, in m
79 set HalfLength 3.00
80
81 # magnetic field
82 set Bz 3.8
83}
84
85####################################
86# Charged hadron tracking efficiency
87####################################
88
89module Efficiency ChargedHadronTrackingEfficiency {
90 set InputArray ParticlePropagator/chargedHadrons
91 set OutputArray chargedHadrons
92
93 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
94
95 # tracking efficiency formula for charged hadrons
96 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
97 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
98 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
99 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
100 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
101 (abs(eta) > 2.5) * (0.00)}
102}
103
104##############################
105# Electron tracking efficiency
106##############################
107
108module Efficiency ElectronTrackingEfficiency {
109 set InputArray ParticlePropagator/electrons
110 set OutputArray electrons
111
112 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
113
114 # tracking efficiency formula for electrons
115 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
116 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
117 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
118 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
119 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
120 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.95) +
122 (abs(eta) > 2.5) * (0.0)}
123}
124
125##########################
126# Muon tracking efficiency
127##########################
128
129module Efficiency MuonTrackingEfficiency {
130 set InputArray ParticlePropagator/muons
131 set OutputArray muons
132
133 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
134
135 # tracking efficiency formula for muons
136 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
137 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
138 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
139 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
140
141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
144 (abs(eta) > 2.5) * (0.90)}
145}
146
147########################################
148# Momentum resolution for charged tracks
149########################################
150
151module MomentumSmearing ChargedHadronMomentumSmearing {
152 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
153 set OutputArray chargedHadrons
154
155 # set ResolutionFormula {resolution formula as a function of eta and pt}
156
157 # resolution formula for charged hadrons
158 # based on arXiv:1405.6569
159 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
160 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
161 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
162}
163
164###################################
165# Momentum resolution for electrons
166###################################
167
168module MomentumSmearing ElectronMomentumSmearing {
169 set InputArray ElectronTrackingEfficiency/electrons
170 set OutputArray electrons
171
172 # set ResolutionFormula {resolution formula as a function of eta and energy}
173
174 # resolution formula for electrons
175 # based on arXiv:1405.6569
176 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
177 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
178 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
179}
180
181###############################
182# Momentum resolution for muons
183###############################
184
185module MomentumSmearing MuonMomentumSmearing {
186 set InputArray MuonTrackingEfficiency/muons
187 set OutputArray muons
188
189 # set ResolutionFormula {resolution formula as a function of eta and pt}
190
191 # resolution formula for muons
192 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
193 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
194 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
195}
196
197##############
198# Track merger
199##############
200
201module Merger TrackMerger {
202# add InputArray InputArray
203 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
204 add InputArray ElectronMomentumSmearing/electrons
205 add InputArray MuonMomentumSmearing/muons
206 set OutputArray tracks
207}
208
209
210
211#############
212# ECAL
213#############
214
215module SimpleCalorimeter ECal {
216 set ParticleInputArray ParticlePropagator/stableParticles
217 set TrackInputArray TrackMerger/tracks
218
219 set TowerOutputArray ecalTowers
220 set EFlowTrackOutputArray eflowTracks
221 set EFlowTowerOutputArray eflowPhotons
222
223 set IsEcal true
224
225 set EnergyMin 0.5
226 set EnergySignificanceMin 2.0
227
228 set SmearTowerCenter true
229
230 set pi [expr {acos(-1)}]
231
232 # lists of the edges of each tower in eta and phi
233 # each list starts with the lower edge of the first tower
234 # the list ends with the higher edged of the last tower
235
236 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
237
238 set PhiBins {}
239 for {set i -180} {$i <= 180} {incr i} {
240 add PhiBins [expr {$i * $pi/180.0}]
241 }
242
243 # 0.02 unit in eta up to eta = 1.5 (barrel)
244 for {set i -85} {$i <= 86} {incr i} {
245 set eta [expr {$i * 0.0174}]
246 add EtaPhiBins $eta $PhiBins
247 }
248
249 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
250
251 set PhiBins {}
252 for {set i -180} {$i <= 180} {incr i} {
253 add PhiBins [expr {$i * $pi/180.0}]
254 }
255
256 # 0.02 unit in eta up to eta = 3
257 for {set i 1} {$i <= 84} {incr i} {
258 set eta [expr { -2.958 + $i * 0.0174}]
259 add EtaPhiBins $eta $PhiBins
260 }
261
262 for {set i 1} {$i <= 84} {incr i} {
263 set eta [expr { 1.4964 + $i * 0.0174}]
264 add EtaPhiBins $eta $PhiBins
265 }
266
267 # take present CMS granularity for HF
268
269 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
270 set PhiBins {}
271 for {set i -18} {$i <= 18} {incr i} {
272 add PhiBins [expr {$i * $pi/18.0}]
273 }
274
275 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
276 add EtaPhiBins $eta $PhiBins
277 }
278
279
280 add EnergyFraction {0} {0.0}
281 # energy fractions for e, gamma and pi0
282 add EnergyFraction {11} {1.0}
283 add EnergyFraction {22} {1.0}
284 add EnergyFraction {111} {1.0}
285 # energy fractions for muon, neutrinos and neutralinos
286 add EnergyFraction {12} {0.0}
287 add EnergyFraction {13} {0.0}
288 add EnergyFraction {14} {0.0}
289 add EnergyFraction {16} {0.0}
290 add EnergyFraction {1000022} {0.0}
291 add EnergyFraction {1000023} {0.0}
292 add EnergyFraction {1000025} {0.0}
293 add EnergyFraction {1000035} {0.0}
294 add EnergyFraction {1000045} {0.0}
295 # energy fractions for K0short and Lambda
296 add EnergyFraction {310} {0.3}
297 add EnergyFraction {3122} {0.3}
298
299 # set ResolutionFormula {resolution formula as a function of eta and energy}
300
301 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
302
303 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
304 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
305 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
306 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
307 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
308
309}
310
311
312#############
313# HCAL
314#############
315
316module SimpleCalorimeter HCal {
317 set ParticleInputArray ParticlePropagator/stableParticles
318 set TrackInputArray ECal/eflowTracks
319
320 set TowerOutputArray hcalTowers
321 set EFlowTrackOutputArray eflowTracks
322 set EFlowTowerOutputArray eflowNeutralHadrons
323
324 set IsEcal false
325
326 set EnergyMin 1.0
327 set EnergySignificanceMin 1.0
328
329 set SmearTowerCenter true
330
331 set pi [expr {acos(-1)}]
332
333 # lists of the edges of each tower in eta and phi
334 # each list starts with the lower edge of the first tower
335 # the list ends with the higher edged of the last tower
336
337 # 5 degrees towers
338 set PhiBins {}
339 for {set i -36} {$i <= 36} {incr i} {
340 add PhiBins [expr {$i * $pi/36.0}]
341 }
342 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
343 add EtaPhiBins $eta $PhiBins
344 }
345
346 # 10 degrees towers
347 set PhiBins {}
348 for {set i -18} {$i <= 18} {incr i} {
349 add PhiBins [expr {$i * $pi/18.0}]
350 }
351 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
352 add EtaPhiBins $eta $PhiBins
353 }
354
355 # 20 degrees towers
356 set PhiBins {}
357 for {set i -9} {$i <= 9} {incr i} {
358 add PhiBins [expr {$i * $pi/9.0}]
359 }
360 foreach eta {-5 -4.7 -4.525 4.7 5} {
361 add EtaPhiBins $eta $PhiBins
362 }
363
364 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
365 add EnergyFraction {0} {1.0}
366 # energy fractions for e, gamma and pi0
367 add EnergyFraction {11} {0.0}
368 add EnergyFraction {22} {0.0}
369 add EnergyFraction {111} {0.0}
370 # energy fractions for muon, neutrinos and neutralinos
371 add EnergyFraction {12} {0.0}
372 add EnergyFraction {13} {0.0}
373 add EnergyFraction {14} {0.0}
374 add EnergyFraction {16} {0.0}
375 add EnergyFraction {1000022} {0.0}
376 add EnergyFraction {1000023} {0.0}
377 add EnergyFraction {1000025} {0.0}
378 add EnergyFraction {1000035} {0.0}
379 add EnergyFraction {1000045} {0.0}
380 # energy fractions for K0short and Lambda
381 add EnergyFraction {310} {0.7}
382 add EnergyFraction {3122} {0.7}
383
384 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
385 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
386 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
387
388}
389
390
391#################
392# Electron filter
393#################
394
395module PdgCodeFilter ElectronFilter {
396 set InputArray HCal/eflowTracks
397 set OutputArray electrons
398 set Invert true
399 add PdgCode {11}
400 add PdgCode {-11}
401}
402
403######################
404# ChargedHadronFilter
405######################
406
407module PdgCodeFilter ChargedHadronFilter {
408 set InputArray HCal/eflowTracks
409 set OutputArray chargedHadrons
410
411 add PdgCode {11}
412 add PdgCode {-11}
413 add PdgCode {13}
414 add PdgCode {-13}
415}
416
417
418###################################################
419# Tower Merger (in case not using e-flow algorithm)
420###################################################
421
422module Merger Calorimeter {
423# add InputArray InputArray
424 add InputArray ECal/ecalTowers
425 add InputArray HCal/hcalTowers
426 set OutputArray towers
427}
428
429
430
431####################
432# Energy flow merger
433####################
434
435module Merger EFlowMerger {
436# add InputArray InputArray
437 add InputArray HCal/eflowTracks
438 add InputArray ECal/eflowPhotons
439 add InputArray HCal/eflowNeutralHadrons
440 set OutputArray eflow
441}
442
443######################
444# EFlowFilter
445######################
446
447module PdgCodeFilter EFlowFilter {
448 set InputArray EFlowMerger/eflow
449 set OutputArray eflow
450
451 add PdgCode {11}
452 add PdgCode {-11}
453 add PdgCode {13}
454 add PdgCode {-13}
455}
456
457
458###################
459# Photon efficiency
460###################
461
462module Efficiency PhotonEfficiency {
463 set InputArray ECal/eflowPhotons
464 set OutputArray photons
465
466 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
467
468 # efficiency formula for photons
469 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
470 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
471 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
472 (abs(eta) > 2.5) * (0.00)}
473}
474
475##################
476# Photon isolation
477##################
478
479module Isolation PhotonIsolation {
480 set CandidateInputArray PhotonEfficiency/photons
481 set IsolationInputArray EFlowFilter/eflow
482
483 set OutputArray photons
484
485 set DeltaRMax 0.5
486
487 set PTMin 0.5
488
489 set PTRatioMax 0.12
490}
491
492
493#####################
494# Electron efficiency
495#####################
496
497module Efficiency ElectronEfficiency {
498 set InputArray ElectronFilter/electrons
499 set OutputArray electrons
500
501 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
502
503 # efficiency formula for electrons
504 set EfficiencyFormula { (pt <= 10.0) * (0.0) +
505 (abs(eta) <= 1.5) * (pt > 10.0) * (0.97) +
506 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.97) +
507 (abs(eta) > 2.5) * (0.0)}
508}
509
510####################
511# Electron isolation
512####################
513
514module Isolation ElectronIsolation {
515 set CandidateInputArray ElectronEfficiency/electrons
516 set IsolationInputArray EFlowFilter/eflow
517
518 set OutputArray electrons
519
520 set DeltaRMax 0.4
521
522 set PTMin 0.1
523
524 set PTRatioMax 0.2
525}
526
527#################
528# Muon efficiency
529#################
530
531module Efficiency MuonEfficiency {
532 set InputArray MuonMomentumSmearing/muons
533 set OutputArray muons
534
535 # set EfficiencyFormula {efficiency as a function of eta and pt}
536
537 # efficiency formula for muons
538 set EfficiencyFormula { (pt <= 10.0) * (0.0) +
539 (abs(eta) <= 1.5) * (pt > 10.0) * (0.97) +
540 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.97) +
541 (abs(eta) > 2.4) * (0.0)}
542}
543
544################
545# Muon isolation
546################
547
548module Isolation MuonIsolation {
549 set CandidateInputArray MuonEfficiency/muons
550 set IsolationInputArray EFlowFilter/eflow
551
552 set OutputArray muons
553
554 set DeltaRMax 0.4
555
556 set PTMin 0.1
557
558 set PTRatioMax 0.25
559}
560
561###################
562# Missing ET merger
563###################
564
565module Merger MissingET {
566# add InputArray InputArray
567 add InputArray EFlowMerger/eflow
568 set MomentumOutputArray momentum
569}
570
571##################
572# Scalar HT merger
573##################
574
575module Merger ScalarHT {
576# add InputArray InputArray
577 add InputArray UniqueObjectFinder/jets
578 add InputArray UniqueObjectFinder/electrons
579 add InputArray UniqueObjectFinder/photons
580 add InputArray UniqueObjectFinder/muons
581 set EnergyOutputArray energy
582}
583
584
585#####################
586# Neutrino Filter
587#####################
588
589module PdgCodeFilter NeutrinoFilter {
590
591 set InputArray Delphes/stableParticles
592 set OutputArray filteredParticles
593
594 set PTMin 0.0
595
596 add PdgCode {12}
597 add PdgCode {14}
598 add PdgCode {16}
599 add PdgCode {-12}
600 add PdgCode {-14}
601 add PdgCode {-16}
602
603}
604
605
606#####################
607# MC truth jet finder
608#####################
609
610module FastJetFinder GenJetFinder {
611 set InputArray NeutrinoFilter/filteredParticles
612
613 set OutputArray jets
614
615 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
616 set JetAlgorithm 6
617 set ParameterR 0.5
618
619 set JetPTMin 20.0
620}
621
622#########################
623# Gen Missing ET merger
624########################
625
626module Merger GenMissingET {
627# add InputArray InputArray
628 add InputArray NeutrinoFilter/filteredParticles
629 set MomentumOutputArray momentum
630}
631
632
633
634############
635# Jet finder
636############
637
638module FastJetFinder FastJetFinder {
639# set InputArray Calorimeter/towers
640 set InputArray EFlowMerger/eflow
641
642 set OutputArray jets
643
644 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
645 set JetAlgorithm 6
646 set ParameterR 0.5
647
648 set JetPTMin 20.0
649}
650
651
652
653###############
654# Jet Fakes
655###############
656
657
658module JetFakeParticle JetFakeParticle {
659set InputArray FastJetFinder/jets
660set PhotonOutputArray photons
661set ElectronOutputArray electrons
662set MuonOutputArray muons
663set OutputArray jets
664# {PID out} {formula}
665# make sure output particles are only electrons, muons or photons
666
667 add EfficiencyFormula {11} { (abs(eta) <= 2.5) * (pt > 20.0) * (0.001)}
668 add EfficiencyFormula {13} { (abs(eta) <= 2.5) * (pt > 20.0) * (0.001)}
669 add EfficiencyFormula {22} { (abs(eta) <= 2.5) * (pt > 20.0) * (0.001)}
670 }
671
672
673##################
674# Electron merger
675##################
676
677module Merger ElectronMerger {
678# add InputArray InputArray
679 add InputArray ElectronIsolation/electrons
680 add InputArray JetFakeParticle/electrons
681 set OutputArray electrons
682 }
683
684##################
685# Muon merger
686##################
687
688 module Merger MuonMerger {
689 # add InputArray InputArray
690 add InputArray MuonIsolation/muons
691 add InputArray JetFakeParticle/muons
692 set OutputArray muons
693}
694
695##################
696# Photon merger
697 ##################
698 module Merger PhotonMerger {
699 # add InputArray InputArray
700 add InputArray PhotonIsolation/photons
701 add InputArray JetFakeParticle/photons
702 set OutputArray photons
703 }
704
705##################
706# Fat Jet finder
707##################
708
709module FastJetFinder FatJetFinder {
710 set InputArray EFlowMerger/eflow
711
712 set OutputArray jets
713
714 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
715 set JetAlgorithm 6
716 set ParameterR 0.8
717
718 set ComputeNsubjettiness 1
719 set Beta 1.0
720 set AxisMode 4
721
722 set ComputeTrimming 1
723 set RTrim 0.2
724 set PtFracTrim 0.05
725
726 set ComputePruning 1
727 set ZcutPrun 0.1
728 set RcutPrun 0.5
729 set RPrun 0.8
730
731 set ComputeSoftDrop 1
732 set BetaSoftDrop 0.0
733 set SymmetryCutSoftDrop 0.1
734 set R0SoftDrop 0.8
735
736 set JetPTMin 200.0
737}
738
739
740
741
742##################
743# Jet Energy Scale
744##################
745
746module EnergyScale JetEnergyScale {
747 set InputArray FastJetFinder/jets
748 set InputArray JetFakeParticle/jets
749 set OutputArray jets
750
751 # scale formula for jets
752 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
753}
754
755########################
756# Jet Flavor Association
757########################
758
759module JetFlavorAssociation JetFlavorAssociation {
760
761 set PartonInputArray Delphes/partons
762 set ParticleInputArray Delphes/allParticles
763 set ParticleLHEFInputArray Delphes/allParticlesLHEF
764 set JetInputArray JetEnergyScale/jets
765
766 set DeltaR 0.5
767 set PartonPTMin 1.0
768 set PartonEtaMax 2.5
769
770}
771
772###########
773# b-tagging
774###########
775
776module BTagging BTagging {
777 set JetInputArray JetEnergyScale/jets
778
779 set BitNumber 0
780
781 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
782 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
783 # gluon's PDG code has the lowest priority
784
785 # based on arXiv:1211.4462
786
787 # default efficiency formula (misidentification rate)
788 add EfficiencyFormula {0} {0.01+0.000038*pt}
789
790 # efficiency formula for c-jets (misidentification rate)
791 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
792
793 # efficiency formula for b-jets
794 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
795}
796
797#############
798# tau-tagging
799#############
800
801module TauTagging TauTagging {
802 set ParticleInputArray Delphes/allParticles
803 set PartonInputArray Delphes/partons
804 set JetInputArray JetEnergyScale/jets
805
806 set DeltaR 0.5
807
808 set TauPTMin 1.0
809
810 set TauEtaMax 2.5
811
812 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
813
814 # default efficiency formula (misidentification rate)
815 add EfficiencyFormula {0} {0.01}
816 # efficiency formula for tau-jets
817 add EfficiencyFormula {15} {0.6}
818}
819
820#####################################################
821# Find uniquely identified photons/electrons/tau/jets
822#####################################################
823
824module UniqueObjectFinder UniqueObjectFinder {
825# earlier arrays take precedence over later ones
826# add InputArray InputArray OutputArray
827 add InputArray PhotonMerger/photons photons
828 add InputArray ElectronMerger/electrons electrons
829 add InputArray MuonMerger/muons muons
830 add InputArray JetFakeParticle/jets jets
831 # add InputArray JetEnergyScale/jets jets
832}
833
834##################
835# ROOT tree writer
836##################
837
838# tracks, towers and eflow objects are not stored by default in the output.
839# if needed (for jet constituent or other studies), uncomment the relevant
840# "add Branch ..." lines.
841
842module TreeWriter TreeWriter {
843# add Branch InputArray BranchName BranchClass
844 add Branch Delphes/allParticles Particle GenParticle
845
846 add Branch TrackMerger/tracks Track Track
847 add Branch Calorimeter/towers Tower Tower
848
849 add Branch HCal/eflowTracks EFlowTrack Track
850 add Branch ECal/eflowPhotons EFlowPhoton Tower
851 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
852
853 add Branch GenJetFinder/jets GenJet Jet
854 add Branch GenMissingET/momentum GenMissingET MissingET
855
856 add Branch UniqueObjectFinder/jets Jet Jet
857 add Branch UniqueObjectFinder/electrons Electron Electron
858 add Branch UniqueObjectFinder/photons Photon Photon
859 add Branch UniqueObjectFinder/muons Muon Muon
860
861 add Branch FatJetFinder/jets FatJet Jet
862
863 add Branch MissingET/momentum MissingET MissingET
864 add Branch ScalarHT/energy ScalarHT ScalarHT
865}