Fork me on GitHub

source: git/cards/delphes_card_LHeC.tcl@ a6e8029

Last change on this file since a6e8029 was c856b3a, checked in by Michele Selvaggi <michele.selvaggi@…>, 4 years ago

added LHeC and FCCeh cards
  • Property mode set to 100644
File size: 23.1 KB
Line 
1#######################################
2# Order of execution of various modules
3# 15.9.2014 try a fine granulated calorimeter for LHeC
4# 0.025 in eta and phi
5# 19.12.16 add paramterisations according to Peter Kostka's design
6# muons not considered here - take CMS/FCC values
7# without CTagging, with electrons, muons and b-tagging issues considered
8# 30.1.2018 hammad add Muon Filter module
9#######################################
10
11set ExecutionPath {
12 ParticlePropagator
13
14 ChargedHadronTrackingEfficiency
15 ElectronTrackingEfficiency
16 MuonTrackingEfficiency
17
18 ChargedHadronMomentumSmearing
19 ElectronMomentumSmearing
20 MuonMomentumSmearing
21
22 TrackMerger
23 Calorimeter
24 EFlowMerger
25
26 PhotonEfficiency
27 PhotonIsolation
28
29 ElectronFilter
30 MuonFilter
31 ElectronEfficiency
32 ElectronIsolation
33
34 MuonEfficiency
35 MuonIsolation
36
37 MissingET
38
39 GenJetFinder
40 FastJetFinder
41
42 JetEnergyScale
43
44 GenJetFlavorAssociation
45 JetFlavorAssociation
46
47 GenBTagging
48 BTagging
49
50 GenTauTagging
51 TauTagging
52
53 UniqueObjectFinder
54
55 ScalarHT
56
57 TreeWriter
58}
59
60# exclude since it adds also an entry to Jet->BTag
61# TrackCountingBTagging
62
63#################################
64# Propagate particles in cylinder
65#################################
66
67module ParticlePropagator ParticlePropagator {
68 set InputArray Delphes/stableParticles
69
70 set OutputArray stableParticles
71 set ChargedHadronOutputArray chargedHadrons
72 set ElectronOutputArray electrons
73 set MuonOutputArray muons
74
75# LHeC
76# radius of the magnetic field coverage, in m
77 set Radius 1.14
78 # half-length of the magnetic field coverage, in m
79 set HalfLength 2.9
80
81 # magnetic field
82 set Bz 3.8
83
84}
85
86####################################
87# Charged hadron tracking efficiency
88####################################
89
90module Efficiency ChargedHadronTrackingEfficiency {
91 set InputArray ParticlePropagator/chargedHadrons
92 set OutputArray chargedHadrons
93
94# add EfficiencyFormula {efficiency formula as a function of eta and pt}
95
96# tracking efficiency formula for charged hadrons
97# set to 1 for full range (-0.2 eta already subtracted from edges)
98 set EfficiencyFormula { ( eta <= 4.9 && eta >= -4.3) * (1.0) + \
99 (eta > 4.9 || eta < -4.3 ) * (0.00)}
100}
101##############################
102# Electron tracking efficiency
103##############################
104
105module Efficiency ElectronTrackingEfficiency {
106 set InputArray ParticlePropagator/electrons
107 set OutputArray electrons
108
109 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
110
111 # tracking efficiency formula for electrons
112# set to 1
113 set EfficiencyFormula { ( (eta <=4.9 && eta >= -4.3) ) * (1.0) + \
114 ((eta > 4.9 || eta < -4.3)) * (0.00)}
115 }
116##########################
117# Muon tracking efficiency
118##########################
119
120module Efficiency MuonTrackingEfficiency {
121 set InputArray ParticlePropagator/muons
122 set OutputArray muons
123
124 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
125
126 # tracking efficiency formula for muons
127# set to 1
128 set EfficiencyFormula { (abs(eta) <= 4) * (1.0) + \
129 (abs(eta) > 4) * (0.00)}
130}
131
132
133########################################
134# Momentum resolution for charged tracks
135########################################
136
137module MomentumSmearing ChargedHadronMomentumSmearing {
138 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
139 set OutputArray chargedHadrons
140
141 # set ResolutionFormula {resolution formula as a function of eta and pt}
142
143 # resolution formula for charged hadrons
144# taken from FCC but adjust to ep
145#the abs(eta) <= 2.0 should be the "central best" and abs(eta) >
146#2.0 and <= 4.9 the fwd / <= 4.3 the bwd (LHeC)
147# and abs(eta) > 2.0 and <= 5.2 fwd / <= 5.0 the bwd (FCC)
148
149 set ResolutionFormula { (abs(eta) <= 2.0) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
150 (abs(eta) <= 2.0) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
151 (abs(eta) <= 2.0) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
152 (abs(eta) <= 2.0) * (pt > 2.0e2) * (0.05) + \
153 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
154 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
155 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
156 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 2.0e2) * (0.05)}
157}
158
159#################################
160# Energy resolution for electrons
161#################################
162
163module MomentumSmearing ElectronMomentumSmearing {
164 set InputArray ElectronTrackingEfficiency/electrons
165 set OutputArray electrons
166
167 # set ResolutionFormula {resolution formula as a function of eta and pt}
168
169 # same as resolution formula for charged hadrons
170# taken from FCC but adjust to ep
171#the abs(eta) <= 2.0 should be the "central best" and abs(eta) >
172#2.0 and <= 4.9 the fwd / <= 4.3 the bwd (LHeC)
173# and abs(eta) > 2.0 and <= 5.2 fwd / <= 5.0 the bwd (FCC)
174
175 set ResolutionFormula { (abs(eta) <= 2.0) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
176 (abs(eta) <= 2.0) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
177 (abs(eta) <= 2.0) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
178 (abs(eta) <= 2.0) * (pt > 2.0e2) * (0.05) + \
179 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
180 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
181 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
182 (abs(eta) > 2.0 && ( eta <= 4.9 && eta >= -4.3 ) ) * (pt > 2.0e2) * (0.05)}
183}
184
185
186###############################
187# Momentum resolution for muons
188###############################
189
190module MomentumSmearing MuonMomentumSmearing {
191 set InputArray MuonTrackingEfficiency/muons
192 set OutputArray muons
193
194 # set ResolutionFormula {resolution formula as a function of eta and pt}
195
196 # resolution formula for muons
197# taken from FCC
198 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
199 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
200 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
201 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
202 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
203 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
204 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
205 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
206 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
207 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
208 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
209 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
210}
211
212
213
214##############
215# Track merger
216##############
217
218module Merger TrackMerger {
219# add InputArray InputArray
220 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
221 add InputArray ElectronMomentumSmearing/electrons
222 add InputArray MuonMomentumSmearing/muons
223 set OutputArray tracks
224}
225
226##################
227# LHeC Calorimeter
228##################
229
230module Calorimeter Calorimeter {
231 set ParticleInputArray ParticlePropagator/stableParticles
232 set TrackInputArray TrackMerger/tracks
233
234 set TowerOutputArray towers
235 set PhotonOutputArray photons
236
237 set EFlowTrackOutputArray eflowTracks
238 set EFlowPhotonOutputArray eflowPhotons
239 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
240
241 set pi [expr {acos(-1)}]
242
243 # lists of the edges of each tower in eta and phi
244 # each list starts with the lower edge of the first tower
245 # the list ends with the higher edged of the last tower
246
247 # use a delta_phi=0.025 for 2*phi range == 252 cells
248 # use a delta_eta-0.025 for eta = -5 to 5 == 400 cells
249 # Uta 22.12.16 use eta range from Peter Kostka -4.7 - 5.2
250 set PhiBins {}
251 for {set i -126} {$i <= 126} {incr i} {
252 add PhiBins [expr {$i * $pi/126.0}]
253 }
254 foreach eta { -4.700 -4.675 -4.650 -4.625 -4.600 -4.575 -4.550 -4.525 -4.500 -4.475 -4.450 -4.425 -4.400 -4.375 -4.350 -4.325 -4.300 -4.275 -4.250 -4.225 -4.200 -4.175 -4.150 -4.125 -4.100 -4.075 -4.050 -4.025 -4.000 -3.975 -3.950 -3.925 -3.900 -3.875 -3.850 -3.825 -3.800 -3.775 -3.750 -3.725 -3.700 -3.675 -3.650 -3.625 -3.600 -3.575 -3.550 -3.525 -3.500 -3.475 -3.450 -3.425 -3.400 -3.375 -3.350 -3.325 -3.300 -3.275 -3.250 -3.225 -3.200 -3.175 -3.150 -3.125 -3.100 -3.075 -3.050 -3.025 -3.000 -2.975 -2.950 -2.925 -2.900 -2.875 -2.850 -2.825 -2.800 -2.775 -2.750 -2.725 -2.700 -2.675 -2.650 -2.625 -2.600 -2.575 -2.550 -2.525 -2.500 -2.475 -2.450 -2.425 -2.400 -2.375 -2.350 -2.325 -2.300 -2.275 -2.250 -2.225 -2.200 -2.175 -2.150 -2.125 -2.100 -2.075 -2.050 -2.025 -2.000 -1.975 -1.950 -1.925 -1.900 -1.875 -1.850 -1.825 -1.800 -1.775 -1.750 -1.725 -1.700 -1.675 -1.650 -1.625 -1.600 -1.575 -1.550 -1.525 -1.500 -1.475 -1.450 -1.425 -1.400 -1.375 -1.350 -1.325 -1.300 -1.275 -1.250 -1.225 -1.200 -1.175 -1.150 -1.125 -1.100 -1.075 -1.050 -1.025 -1.000 -0.975 -0.950 -0.925 -0.900 -0.875 -0.850 -0.825 -0.800 -0.775 -0.750 -0.725 -0.700 -0.675 -0.650 -0.625 -0.600 -0.575 -0.550 -0.525 -0.500 -0.475 -0.450 -0.425 -0.400 -0.375 -0.350 -0.325 -0.300 -0.275 -0.250 -0.225 -0.200 -0.175 -0.150 -0.125 -0.100 -0.075 -0.050 -0.025 0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 0.225 0.250 0.275 0.300 0.325 0.350 0.375 0.400 0.425 0.450 0.475 0.500 0.525 0.550 0.575 0.600 0.625 0.650 0.675 0.700 0.725 0.750 0.775 0.800 0.825 0.850 0.875 0.900 0.925 0.950 0.975 1.000 1.025 1.050 1.075 1.100 1.125 1.150 1.175 1.200 1.225 1.250 1.275 1.300 1.325 1.350 1.375 1.400 1.425 1.450 1.475 1.500 1.525 1.550 1.575 1.600 1.625 1.650 1.675 1.700 1.725 1.750 1.775 1.800 1.825 1.850 1.875 1.900 1.925 1.950 1.975 2.000 2.025 2.050 2.075 2.100 2.125 2.150 2.175 2.200 2.225 2.250 2.275 2.300 2.325 2.350 2.375 2.400 2.425 2.450 2.475 2.500 2.525 2.550 2.575 2.600 2.625 2.650 2.675 2.700 2.725 2.750 2.775 2.800 2.825 2.850 2.875 2.900 2.925 2.950 2.975 3.000 3.025 3.050 3.075 3.100 3.125 3.150 3.175 3.200 3.225 3.250 3.275 3.300 3.325 3.350 3.375 3.400 3.425 3.450 3.475 3.500 3.525 3.550 3.575 3.600 3.625 3.650 3.675 3.700 3.725 3.750 3.775 3.800 3.825 3.850 3.875 3.900 3.925 3.950 3.975 4.000 4.025 4.050 4.075 4.100 4.125 4.150 4.175 4.200 4.225 4.250 4.275 4.300 4.325 4.350 4.375 4.400 4.425 4.450 4.475 4.500 4.525 4.550 4.575 4.600 4.625 4.650 4.675 4.700 4.725 4.750 4.775 4.800 4.825 4.850 4.875 4.900 4.925 4.950 4.975 5.00 5.025 5.050 5.075 5.100 5.125 5.150 5.175 5.200
255 } {
256 add EtaPhiBins $eta $PhiBins
257 }
258
259
260 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
261 add EnergyFraction {0} {0.0 1.0}
262 # energy fractions for e, gamma and pi0
263 add EnergyFraction {11} {1.0 0.0}
264 add EnergyFraction {22} {1.0 0.0}
265 add EnergyFraction {111} {1.0 0.0}
266 # energy fractions for muon, neutrinos and neutralinos
267 add EnergyFraction {12} {0.0 0.0}
268 add EnergyFraction {13} {0.0 0.0}
269 add EnergyFraction {14} {0.0 0.0}
270 add EnergyFraction {16} {0.0 0.0}
271 add EnergyFraction {1000022} {0.0 0.0}
272 add EnergyFraction {1000023} {0.0 0.0}
273 add EnergyFraction {1000025} {0.0 0.0}
274 add EnergyFraction {1000035} {0.0 0.0}
275 add EnergyFraction {1000045} {0.0 0.0}
276 # energy fractions for K0short and Lambda
277 add EnergyFraction {310} {0.3 0.7}
278 add EnergyFraction {3122} {0.3 0.7}
279
280## no minimum p range given (Uta 22.12.16)
281## sigma_E = E* (a/sqrt(E) + b) = sqrt ( E*a^2 E +E^2*b^2 )
282## a = sampling term, b = constant term
283# ECAL eta : 5.0 - 2.7 - -2.1 - -4.4 (subtract -0.1 in eta as safety distance to edges)
284 set ECalResolutionFormula { ( eta <= 5.0 && eta > 2.7 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) + \
285 ( eta <=2.7 && eta >= -2.1 ) * sqrt(energy*0.09^2 + energy^2*0.02^2) + \
286 ( eta < -2.1 && eta >= -4.4 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) }
287
288#### HCAL eta : 5.1 - 2.1 - -1.7 - -4.6 (subtract -0.1 in eta as safety distance to edges)
289#set HCalResolutionFormula {resolution formula as a function of eta and energy}
290
291 set HCalResolutionFormula { ( eta <= 5.1 && eta > 2.1 ) * sqrt(energy*0.4^2 + energy^2*0.02^2) + \
292 ( eta <=2.1 && eta >= -1.7 ) * sqrt(energy*0.4^2 + energy^2*0.02^2) + \
293 ( eta < -1.7 && eta >= -4.7 ) * sqrt(energy*0.4^2 + energy^2*0.04^2) }
294
295}
296####################
297# Energy flow merger
298####################
299
300module Merger EFlowMerger {
301# add InputArray InputArray
302# exchange UK: 25.7.16 change back to all eflowtracks!
303# add InputArray ImpactParameterSmearing/tracks
304## eflowTracks kills electrons for Delphes 3.3.3 and ep??
305 add InputArray Calorimeter/eflowTracks
306 add InputArray Calorimeter/eflowPhotons
307 add InputArray Calorimeter/eflowNeutralHadrons
308 set OutputArray eflow
309}
310
311###################
312# Photon efficiency
313###################
314
315module Efficiency PhotonEfficiency {
316# set InputArray Calorimeter/eflowPhotons
317# 21.7.16 use photons
318 set InputArray Calorimeter/photons
319 set OutputArray photons
320
321 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
322
323# efficiency formula for photons - same as for electrons
324# expand ad set flat
325## OLD EFFICIENCY
326## set EfficiencyFormula { ( (eta <=4.9 && eta >= -4.3) ) * (1.0) + \
327## ((eta > 4.9 || eta < -4.3)) * (0.00)}
328#########################################################
329# OF, 20.03.24: NEW EFFICIENCY TO AVOID SOFT PHOTONS ##
330#########################################################
331 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( (eta <=4.9 && eta >= -4.3) ) * (1.0) + \
332 ((eta > 4.9 || eta < -4.3)) * (0.00)}
333
334}
335
336##################
337# Photon isolation
338##################
339
340module Isolation PhotonIsolation {
341 set CandidateInputArray PhotonEfficiency/photons
342 set IsolationInputArray EFlowMerger/eflow
343
344 set OutputArray photons
345
346 set DeltaRMax 0.4
347
348 set PTMin 0.5
349
350 set PTRatioMax 0.1
351}
352
353#################
354# Electron filter
355#################
356
357module PdgCodeFilter ElectronFilter {
358 set InputArray Calorimeter/eflowTracks
359 set OutputArray electrons
360 set Invert true
361 add PdgCode {11}
362 add PdgCode {-11}
363}
364
365#################
366# Muon filter
367#################
368
369module PdgCodeFilter MuonFilter {
370 set InputArray Calorimeter/eflowTracks
371 set OutputArray muons
372 set Invert true
373 add PdgCode {13}
374 add PdgCode {-13}
375}
376
377#####################
378# Electron efficiency
379#####################
380
381module Efficiency ElectronEfficiency {
382 set InputArray ElectronFilter/electrons
383 set OutputArray electrons
384
385 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
386
387 # efficiency formula for electrons
388 # 1 degree = eta<4.7
389## OLD EFFICIENCY
390## set EfficiencyFormula { ( (eta <=4.9 && eta >= -4.3) ) * (1.0) + \
391## ((eta > 4.9 || eta < -4.3)) * (0.00)}
392#########################################################
393# OF, 20.03.24: NEW EFFICIENCY TO AVOID SOFT ELECTRONS ##
394#########################################################
395 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( (eta <=4.9 && eta >= -4.3) ) * (1.0) + \
396 ((eta > 4.9 || eta < -4.3)) * (0.00)}
397
398}
399
400####################
401# Electron isolation
402####################
403
404module Isolation ElectronIsolation {
405 set CandidateInputArray ElectronEfficiency/electrons
406 set IsolationInputArray EFlowMerger/eflow
407
408 set OutputArray electrons
409
410 set DeltaRMax 0.4
411
412 set PTMin 0.5
413
414 set PTRatioMax 0.1
415}
416
417#################
418# Muon efficiency
419#################
420
421module Efficiency MuonEfficiency {
422 set InputArray MuonMomentumSmearing/muons
423 set InputArray MuonFilter/muons
424 set OutputArray muons
425
426 # set EfficiencyFormula {efficiency as a function of eta and pt}
427
428 # efficiency formula for muons
429# set to 1 for |eta|<4
430 set EfficiencyFormula { (abs(eta) <= 4) * (1.0) + \
431 (abs(eta) > 4) * (0.00)}
432}
433
434################
435# Muon isolation
436################
437
438module Isolation MuonIsolation {
439 set CandidateInputArray MuonEfficiency/muons
440 set IsolationInputArray EFlowMerger/eflow
441
442 set OutputArray muons
443
444 set DeltaRMax 0.4
445
446 set PTMin 0.5
447
448 set PTRatioMax 0.1
449}
450
451###################
452# Missing ET merger
453###################
454
455module Merger MissingET {
456# add InputArray InputArray
457 add InputArray EFlowMerger/eflow
458 set MomentumOutputArray momentum
459}
460
461##################
462# Scalar HT merger
463##################
464
465module Merger ScalarHT {
466# add InputArray InputArray
467 add InputArray UniqueObjectFinder/jets
468 add InputArray UniqueObjectFinder/electrons
469 add InputArray UniqueObjectFinder/photons
470 add InputArray UniqueObjectFinder/muons
471 set EnergyOutputArray energy
472}
473
474#####################
475# MC truth jet finder
476#####################
477
478module FastJetFinder GenJetFinder {
479 set InputArray Delphes/stableParticles
480
481 set OutputArray jets
482
483 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
484 set JetAlgorithm 6
485 set ParameterR 0.7
486
487 set JetPTMin 3.0
488}
489
490############
491# Jet finder
492############
493
494module FastJetFinder FastJetFinder {
495# ATLAS uses towers
496# set InputArray Calorimeter/towers
497# CMS uses eflow?
498 set InputArray EFlowMerger/eflow
499
500 set OutputArray jets
501
502 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
503 set JetAlgorithm 6
504 set ParameterR 0.4
505
506 set JetPTMin 1.0
507}
508
509
510##################
511# Jet Energy Scale
512##################
513
514module EnergyScale JetEnergyScale {
515 set InputArray FastJetFinder/jets
516 set OutputArray jets
517
518 # scale formula for jets
519 set ScaleFormula {1.00}
520}
521
522########################
523# GenJet Flavor Association
524########################
525
526module JetFlavorAssociation GenJetFlavorAssociation {
527 set PartonInputArray Delphes/partons
528 set ParticleInputArray Delphes/allParticles
529 set ParticleLHEFInputArray Delphes/allParticlesLHEF
530 set JetInputArray GenJetFinder/jets
531
532 set DeltaR 0.4
533 set PartonPTMin 0.5
534 set PartonEtaMax 3.0
535
536}
537
538########################
539# Jet Flavor Association
540########################
541
542module JetFlavorAssociation JetFlavorAssociation {
543
544 set PartonInputArray Delphes/partons
545 set ParticleInputArray Delphes/allParticles
546 set ParticleLHEFInputArray Delphes/allParticlesLHEF
547 set JetInputArray JetEnergyScale/jets
548
549 set DeltaR 0.4
550 set PartonPTMin 0.5
551 set PartonEtaMax 3.0
552
553}
554
555###########
556# b-tagging
557###########
558module BTagging GenBTagging {
559 set PartonInputArray Delphes/partons
560 set JetInputArray GenJetFinder/jets
561
562 set BitNumber 0
563
564 set DeltaR 0.4
565
566 set PartonPTMin 0.5
567
568 set PartonEtaMax 6.0
569
570 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
571 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
572 # gluon's PDG code has the lowest priority
573
574 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
575 # default efficiency formula (misidentification rate)
576 add EfficiencyFormula {0} {0.001}
577
578 # efficiency formula for c-jets (misidentification rate)
579 add EfficiencyFormula {4} {0.05}
580
581 # efficiency formula for b-jets
582 add EfficiencyFormula {5} {0.75}
583}
584
585
586####################
587# standard B-tagging
588## ATTENTION : idealised values - use for cross check only
589####################
590module BTagging BTagging {
591 set PartonInputArray Delphes/partons
592 set JetInputArray JetEnergyScale/jets
593
594 set BitNumber 0
595
596 set DeltaR 0.4
597
598 set PartonPTMin 0.5
599
600 set PartonEtaMax 6.0
601
602 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
603 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
604 # gluon's PDG code has the lowest priority
605
606 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
607 # default efficiency formula (misidentification rate)
608 add EfficiencyFormula {0} {0.001}
609
610
611 # efficiency formula for c-jets (misidentification rate)
612 add EfficiencyFormula {4} {0.05}
613
614 # efficiency formula for b-jets
615 add EfficiencyFormula {5} {0.75}
616
617}
618
619
620
621module TauTagging TauTagging {
622 set ParticleInputArray Delphes/allParticles
623 set PartonInputArray Delphes/partons
624 set JetInputArray JetEnergyScale/jets
625
626 set DeltaR 0.4
627
628 set TauPTMin 0.5
629
630 set TauEtaMax 6.0
631
632 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
633
634 # default efficiency formula (misidentification rate)
635 add EfficiencyFormula {0} {0.001}
636 # efficiency formula for tau-jets
637 add EfficiencyFormula {15} {0.4}
638}
639
640module TauTagging GenTauTagging {
641 set ParticleInputArray Delphes/allParticles
642 set PartonInputArray Delphes/partons
643 set JetInputArray GenJetFinder/jets
644
645 set DeltaR 0.4
646
647 set TauPTMin 0.5
648
649 set TauEtaMax 6.0
650
651 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
652
653 # default efficiency formula (misidentification rate)
654 add EfficiencyFormula {0} {0.001}
655 # efficiency formula for tau-jets
656 add EfficiencyFormula {15} {0.4}
657}
658
659#####################################################
660# Find uniquely identified photons/electrons/tau/jets
661#####################################################
662
663module UniqueObjectFinder UniqueObjectFinder {
664# earlier arrays take precedence over later ones
665# add InputArray InputArray OutputArray
666 add InputArray PhotonIsolation/photons photons
667 add InputArray ElectronIsolation/electrons electrons
668 add InputArray MuonIsolation/muons muons
669 add InputArray JetEnergyScale/jets jets
670}
671
672##################
673# ROOT tree writer
674##################
675
676# tracks, towers and eflow objects are not stored by default in the output.
677# if needed (for jet constituent or other studies), uncomment the relevant
678# "add Branch ..." lines.
679
680module TreeWriter TreeWriter {
681# add Branch InputArray BranchName BranchClass
682 add Branch Delphes/allParticles Particle GenParticle
683
684# add Branch TrackMerger/tracks Track Track
685# add Branch Calorimeter/towers Tower Tower
686
687# add Branch Calorimeter/eflowTracks EFlowTrack Track
688# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
689# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
690
691 add Branch GenJetFinder/jets GenJet Jet
692 add Branch UniqueObjectFinder/jets Jet Jet
693 add Branch UniqueObjectFinder/electrons Electron Electron
694 add Branch UniqueObjectFinder/photons Photon Photon
695 add Branch UniqueObjectFinder/muons Muon Muon
696 add Branch MissingET/momentum MissingET MissingET
697 add Branch ScalarHT/energy ScalarHT ScalarHT
698}
Note: See TracBrowser for help on using the repository browser.