Fork me on GitHub

Ticket #744: info_v2.txt

File info_v2.txt, 20.5 KB (added by Michele Selvaggi, 9 years ago)
Line 
1# FCC proposed detector. Sep. 2015
2# taus reconstruction is removed to deal with input slimmed records
3# based on Delphes 3.3 (Oct. 2015)
4# Created by Heather Gray and Filip Moortgat
5# Modified for HepSim by S.Chekanov (ANL)
6#######################################
7# Order of execution of various modules
8#######################################
9
10set ExecutionPath {
11ParticlePropagator
12StatusPid
13
14ChargedHadronTrackingEfficiency
15
16ElectronEfficiency
17MuonEfficiency
18
19ChargedHadronMomentumSmearing
20ElectronEnergySmearing
21MuonMomentumSmearing
22
23TrackMerger
24ImpactParameterSmearing
25
26Ecal
27Hcal
28
29TowerMerger
30EFlowMerger
31
32PhotonEfficiency
33MuonIsolation
34ElectronIsolation
35PhotonIsolation
36
37MissingET
38
39GenJetFinder
40FastJetFinder
41
42JetEnergyScale
43
44JetFlavorAssociation
45
46BTagging
47CTagging
48
49UniqueObjectFinder
50
51ScalarHT
52
53TreeWriter
54}
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61set InputArray Delphes/stableParticles
62
63set OutputArray stableParticles
64set ChargedHadronOutputArray chargedHadrons
65set ElectronOutputArray electrons
66set MuonOutputArray muons
67
68# radius of the magnetic field coverage, in m
69set Radius 6.00
70# half-length of the magnetic field coverage, in m
71set HalfLength 11.5
72
73# magnetic field
74set Bz 6.0
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82set InputArray ParticlePropagator/chargedHadrons
83set OutputArray chargedHadrons
84
85# add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87# tracking efficiency formula for charged hadrons
88set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
89(abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
90(abs(eta) <= 1.5) * (pt > 1) * (0.95) + \
91(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
92(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
93(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
94(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
95(abs(eta) > 6.0) * (0.00)}
96}
97
98##############################
99# Electron tracking efficiency
100##############################
101
102module Efficiency ElectronEfficiency {
103set InputArray ParticlePropagator/electrons
104set OutputArray electrons
105
106# set EfficiencyFormula {efficiency formula as a function of eta and pt}
107
108# tracking efficiency formula for electrons
109# tracking efficiency formula for electrons
110set EfficiencyFormula { (pt <= 10) * (0.00) + \
111(abs(eta) <= 1.5) * (pt > 10 && pt <= 50) * (0.80) + \
112(abs(eta) <= 1.5) * (pt > 50) * (0.90) + \
113(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 10 && pt <= 50) * (0.80) + \
114(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 50) * (0.90) + \
115(abs(eta) > 4 && abs(eta) <= 6) * (pt > 10 && pt <= 50) * (0.70) + \
116(abs(eta) > 4 && abs(eta) <= 6) * (pt > 50) * (0.80) + \
117(abs(eta) > 6) * (0.00)}
118}
119
120##########################
121# Muon tracking efficiency
122##########################
123
124module Efficiency MuonEfficiency {
125set InputArray ParticlePropagator/muons
126set OutputArray muons
127
128# set EfficiencyFormula {efficiency formula as a function of eta and pt}
129
130# tracking efficiency formula for muons
131set EfficiencyFormula { (pt <= 5) * (0.00) + \
132(abs(eta) <= 6) * (pt > 5 ) * (0.99) + \
133(abs(eta) > 6) * (0.00)}
134}
135
136########################################
137# Momentum resolution for charged tracks
138########################################
139
140module MomentumSmearing ChargedHadronMomentumSmearing {
141set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
142set OutputArray chargedHadrons
143
144# set ResolutionFormula {resolution formula as a function of eta and pt}
145
146# resolution formula for charged hadrons
147set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1.0) * (0.02) + \
148 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
149 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
150 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
151 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1.0) * (0.03) + \
152 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
153 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
154 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 2.0e2) * (0.05) + \
155 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1.0) * (0.06) + \
156 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0 && pt <= 1.0e1) * (0.04) + \
157 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0e1 && pt <= 2.0e2) * (0.08) + \
158 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 2.0e2) * (0.08) + \
159 (abs(eta) > 6) * (0.00)}
160}
161#################################
162# Energy resolution for electrons
163#################################
164
165module EnergySmearing ElectronEnergySmearing {
166set InputArray ElectronEfficiency/electrons
167set OutputArray electrons
168
169# set ResolutionFormula {resolution formula as a function of eta and energy}
170
171# resolution formula for electrons
172set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) + \
173(abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) + \
174(abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
175}
176
177###############################
178# Momentum resolution for muons
179###############################
180
181module MomentumSmearing MuonMomentumSmearing {
182set InputArray MuonEfficiency/muons
183set OutputArray muons
184
185# set ResolutionFormula {resolution formula as a function of eta and pt}
186
187# resolution formula for muons
188set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
189(abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
190(abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
191(abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
192(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
193(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
194(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
195(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
196(abs(eta) > 1.5 && abs(eta) <= 6.0) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
197(abs(eta) > 1.5 && abs(eta) <= 6.0) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
198(abs(eta) > 1.5 && abs(eta) <= 6.0) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
199(abs(eta) > 1.5 && abs(eta) <= 6.0) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
200}
201
202####################
203# Electron isolation
204####################
205
206module Isolation ElectronIsolation {
207set CandidateInputArray ElectronEfficiency/electrons
208set IsolationInputArray EFlowMerger/eflow
209
210set OutputArray electrons
211
212set DeltaRMax 0.3
213
214set PTMin 50
215
216set PTRatioMax 0.1
217}
218
219################
220# Muon isolation
221################
222
223module Isolation MuonIsolation {
224set CandidateInputArray MuonEfficiency/muons
225set IsolationInputArray EFlowMerger/eflow
226
227set OutputArray muons
228
229set DeltaRMax 0.3
230
231set PTMin 50
232
233set PTRatioMax 0.1
234}
235
236###################
237# Photon efficiency
238###################
239
240module Efficiency PhotonEfficiency {
241set InputArray Ecal/eflowPhotons
242set OutputArray photons
243
244# set EfficiencyFormula {efficiency formula as a function of eta and pt}
245# efficiency formula for photons
246set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
247(abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
248(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.9624) + \
249(abs(eta) > 2.5) * (0.00)}
250}
251
252##################
253# Photon isolation
254##################
255
256module Isolation PhotonIsolation {
257set CandidateInputArray Ecal/eflowPhotons
258set IsolationInputArray EFlowMerger/eflow
259
260set OutputArray photons
261
262set DeltaRMax 0.3
263
264set PTMin 50.0
265
266set PTRatioMax 0.1
267}
268
269##############
270# Track merger
271##############
272
273module Merger TrackMerger {
274# add InputArray InputArray
275add InputArray ChargedHadronMomentumSmearing/chargedHadrons
276add InputArray ElectronEnergySmearing/electrons
277add InputArray MuonMomentumSmearing/muons
278set OutputArray tracks
279}
280
281################################
282# Track impact parameter smearing
283################################
284
285module ImpactParameterSmearing ImpactParameterSmearing {
286set InputArray TrackMerger/tracks
287set OutputArray tracks
288
289
290# absolute impact parameter smearing formula (in mm) as a function of pt and eta
291set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \
292(pt > 5.0) * (0.005)}
293
294}
295
296#############
297# ECAL
298#############
299
300module SimpleCalorimeter Ecal {
301 set ParticleInputArray ParticlePropagator/stableParticles
302 set TrackInputArray ImpactParameterSmearing/tracks
303
304 set TowerOutputArray ecalTowers
305 set EFlowTowerOutputArray eflowPhotons
306
307 set pi [expr {acos(-1)}]
308
309 # lists of the edges of each tower in eta and phi
310 # each list starts with the lower edge of the first tower
311 # the list ends with the higher edged of the last tower
312 # ECAL is x4 better than HCAL
313 # 0.02453125 rad (1.40625 deg) in Phi (256 modules)
314 # Eta: 0.025
315 set PhiBins {}
316 for {set i -128} {$i <= 128} {incr i} {
317 add PhiBins [expr {$i * $pi/128.0}]
318 }
319 foreach eta {-1.6 -1.575 -1.55 -1.525 -1.5 -1.475 -1.45 -1.425 -1.4 -1.375 -1.35 -1.325 -1.3 -1.275 -1.25 -1.225 -1.2 -1.175 -1.15 -1.125 -1.1 -1.075 -1.05 -1.025 -1.0 -0.975 -0.95 -0.925 -0.9 -0.875 -0.85 -0.825 -0.8 -0.775 -0.75 -0.725 -0.7 -0.675 -0.65 -0.625 -0.6 -0.575 -0.55 -0.525 -0.5 -0.475 -0.45 -0.425 -0.4 -0.375 -0.35 -0.325 -0.3 -0.275 -0.25 -0.225 -0.2 -0.175 -0.15 -0.125 -0.1 -0.075 -0.05 -0.025 0.0 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 0.3 0.325 0.35 0.375 0.4 0.425 0.45 0.475 0.5 0.525 0.55 0.575 0.6 0.625 0.65 0.675 0.7 0.725 0.75 0.775 0.8 0.825 0.85 0.875 0.9 0.925 0.95 0.975 1.0 1.025 1.05 1.075 1.1 1.125 1.15 1.175 1.2 1.225 1.25 1.275 1.3 1.325 1.35 1.375 1.4 1.425 1.45 1.475 1.5 1.525 1.55 1.575 1.6} {
320 add EtaPhiBins $eta $PhiBins
321 }
322
323 # 0.05 eta x 0.05 deg phi
324 set PhiBins {}
325 for {set i -64} {$i <= 64} {incr i} {
326 add PhiBins [expr {$i * $pi/64.0}]
327 }
328 foreach eta {-4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 -4.0 -3.95 -3.9 -3.85 -3.8 -3.75 -3.7 -3.65 -3.6 -3.55 -3.5 -3.45 -3.4 -3.35 -3.3 -3.25 -3.2 -3.15 -3.1 -3.05 -3.0 -2.95 -2.9 -2.85 -2.8 -2.75 -2.7 -2.65 -2.6 -2.55 -2.5 -2.45 -2.4 -2.35 -2.3 -2.25 -2.2 -2.15 -2.1 -2.05 -2.0 -1.95 -1.9 -1.85 -1.8 -1.75 -1.7 1.7 1.75 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 3.0 3.05 3.1 3.15 3.2 3.25 3.3 3.35 3.4 3.45 3.5 3.55 3.6 3.65 3.7 3.75 3.8 3.85 3.9 3.95 4.0 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5} {
329 add EtaPhiBins $eta $PhiBins
330 }
331
332 # 0.2 degrees towers, 20 degree
333 set PhiBins {}
334 for {set i -9} {$i <= 9} {incr i} {
335 add PhiBins [expr {$i * $pi/9.0}]
336 }
337 foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
338 add EtaPhiBins $eta $PhiBins
339 }
340 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
341
342 add EnergyFraction {0} {0.4}
343 # energy fractions for e, gamma and pi0
344 add EnergyFraction {11} {1.0}
345 add EnergyFraction {22} {1.0}
346 add EnergyFraction {111} {1.0}
347 # energy fractions for muon, neutrinos and neutralinos
348 add EnergyFraction {12} {0.0}
349 add EnergyFraction {13} {0.0}
350 add EnergyFraction {14} {0.0}
351 add EnergyFraction {16} {0.0}
352 add EnergyFraction {1000022} {0.0}
353 add EnergyFraction {1000023} {0.0}
354 add EnergyFraction {1000025} {0.0}
355 add EnergyFraction {1000035} {0.0}
356 add EnergyFraction {1000045} {0.0}
357 # energy fractions for K0short and Lambda
358 # add EnergyFraction {310} {0.3}
359 # add EnergyFraction {3122} {0.3}
360
361 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
362 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
363 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
364
365
366}
367
368
369#############
370# HCAL
371#############
372
373module SimpleCalorimeter Hcal {
374 set ParticleInputArray ParticlePropagator/stableParticles
375 set TrackInputArray ImpactParameterSmearing/tracks
376
377 set TowerOutputArray hcalTowers
378 set EFlowTowerOutputArray eflowNeutralHadrons
379
380 set pi [expr {acos(-1)}]
381
382 # lists of the edges of each tower in eta and phi
383 # each list starts with the lower edge of the first tower
384 # the list ends with the higher edged of the last tower
385
386 # 128 modules in phi like in ATLAS
387 # 1.40625 degree in Phi
388 set PhiBins {}
389 for {set i -128} {$i <= 128} {incr i} {
390 add PhiBins [expr {$i * $pi/128.0}]
391 }
392 # Use the 0.025 bin in Eta (better than ATLAS 0.1)
393 foreach eta {-1.6 -1.575 -1.55 -1.525 -1.5 -1.475 -1.45 -1.425 -1.4 -1.375 -1.35 -1.325 -1.3 -1.275 -1.25 -1.225 -1.2 -1.175 -1.15 -1.125 -1.1 -1.075 -1.05 -1.025 -1.0 -0.975 -0.95 -0.925 -0.9 -0.875 -0.85 -0.825 -0.8 -0.775 -0.75 -0.725 -0.7 -0.675 -0.65 -0.625 -0.6 -0.575 -0.55 -0.525 -0.5 -0.475 -0.45 -0.425 -0.4 -0.375 -0.35 -0.325 -0.3 -0.275 -0.25 -0.225 -0.2 -0.175 -0.15 -0.125 -0.1 -0.075 -0.05 -0.025 0.0 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 0.3 0.325 0.35 0.375 0.4 0.425 0.45 0.475 0.5 0.525 0.55 0.575 0.6 0.625 0.65 0.675 0.7 0.725 0.75 0.775 0.8 0.825 0.85 0.875 0.9 0.925 0.95 0.975 1.0 1.025 1.05 1.075 1.1 1.125 1.15 1.175 1.2 1.225 1.25 1.275 1.3 1.325 1.35 1.375 1.4 1.425 1.45 1.475 1.5 1.525 1.55 1.575 1.6 } {
394 add EtaPhiBins $eta $PhiBins
395 }
396
397 # 5.625 degrees x 0.1 in eta towers
398 set PhiBins {}
399 for {set i -32} {$i <= 32} {incr i} {
400 add PhiBins [expr {$i * $pi/32.0}]
401 }
402 foreach eta {-4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 -4.0 -3.95 -3.9 -3.85 -3.8 -3.75 -3.7 -3.65 -3.6 -3.55 -3.5 -3.45 -3.4 -3.35 -3.3 -3.25 -3.2 -3.15 -3.1 -3.05 -3.0 -2.95 -2.9 -2.85 -2.8 -2.75 -2.7 -2.65 -2.6 -2.55 -2.5 -2.45 -2.4 -2.35 -2.3 -2.25 -2.2 -2.15 -2.1 -2.05 -2.0 -1.95 -1.9 -1.85 -1.8 -1.75 -1.7 1.7 1.75 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 3.0 3.05 3.1 3.15 3.2 3.25 3.3 3.35 3.4 3.45 3.5 3.55 3.6 3.65 3.7 3.75 3.8 3.85 3.9 3.95 4.0 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5} {
403 add EtaPhiBins $eta $PhiBins
404 }
405
406 # 20 degrees towers
407 set PhiBins {}
408 for {set i -9} {$i <= 9} {incr i} {
409 add PhiBins [expr {$i * $pi/9.0}]
410 }
411 foreach eta {6 -5.6 -5.3 -5 -4.7 -4.525 4.7 5 5.3 5.6 6} {
412 add EtaPhiBins $eta $PhiBins
413 }
414
415 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
416 add EnergyFraction {0} {0.6}
417 # energy fractions for e, gamma and pi0
418 add EnergyFraction {11} {0.0}
419 add EnergyFraction {22} {0.0}
420 add EnergyFraction {111} {0.0}
421 # energy fractions for muon, neutrinos and neutralinos
422 add EnergyFraction {12} {0.0}
423 add EnergyFraction {13} {0.0}
424 add EnergyFraction {14} {0.0}
425 add EnergyFraction {16} {0.0}
426 add EnergyFraction {1000022} {0.0}
427 add EnergyFraction {1000023} {0.0}
428 add EnergyFraction {1000025} {0.0}
429 add EnergyFraction {1000035} {0.0}
430 add EnergyFraction {1000045} {0.0}
431 # energy fractions for K0short and Lambda
432 # add EnergyFraction {310} {0.7}
433 # add EnergyFraction {3122} {0.7}
434
435 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
436 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
437 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
438}
439
440
441
442
443####################
444# Tower Merger (in case not using e-flow algorithm)
445####################
446
447module Merger TowerMerger {
448# add InputArray InputArray
449add InputArray Ecal/ecalTowers
450add InputArray Hcal/hcalTowers
451set OutputArray towers
452}
453
454####################
455# Energy flow merger
456####################
457
458module Merger EFlowMerger {
459# add InputArray InputArray
460add InputArray ImpactParameterSmearing/tracks
461add InputArray Ecal/eflowPhotons
462add InputArray Hcal/eflowNeutralHadrons
463set OutputArray eflow
464}
465
466
467###################
468# Missing ET merger
469###################
470
471module Merger MissingET {
472# add InputArray InputArray
473add InputArray EFlowMerger/eflow
474set MomentumOutputArray momentum
475}
476
477
478##################
479# Scalar HT merger
480##################
481
482module Merger ScalarHT {
483# add InputArray InputArray
484add InputArray EFlowMerger/eflow
485set EnergyOutputArray energy
486}
487
488#####################
489# MC truth jet finder
490#####################
491
492module FastJetFinder GenJetFinder {
493set InputArray Delphes/stableParticles
494
495set OutputArray jets
496
497# To compute N-subjettiness variables (see JHEP 1103:015 (2011) and JHEP 1202:093 (2012)), use option 7
498# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
499set JetAlgorithm 7
500set ParameterR 0.5
501set ComputeNsubjettiness 1
502set Beta 2.0
503set AxisMode 1
504
505set JetPTMin 50.0
506}
507
508############
509# Jet finder
510############
511
512module FastJetFinder FastJetFinder {
513# set InputArray Calorimeter/towers
514set InputArray EFlowMerger/eflow
515
516set OutputArray jets
517
518# To compute N-subjettiness variables (see JHEP 1103:015 (2011) and JHEP 1202:093 (2012)), use option 7
519# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
520set JetAlgorithm 7
521set ParameterR 0.5
522set ComputeNsubjettiness 1
523set Beta 2.0
524set AxisMode 1
525
526
527set JetPTMin 50.0
528}
529
530##################
531# Jet Energy Scale
532##################
533
534module EnergyScale JetEnergyScale {
535set InputArray FastJetFinder/jets
536set OutputArray jets
537
538# scale formula for jets
539set ScaleFormula {1.00}
540}
541
542
543########################
544# Jet Flavor Association
545########################
546
547module JetFlavorAssociation JetFlavorAssociation {
548
549 set PartonInputArray Delphes/partons
550 set ParticleInputArray Delphes/allParticles
551 set ParticleLHEFInputArray Delphes/allParticlesLHEF
552 set JetInputArray JetEnergyScale/jets
553
554 set DeltaR 0.5
555 set PartonPTMin 1.0
556 set PartonEtaMax 6.0
557
558}
559
560##########################
561# tau-tagging
562##########################
563
564
565module TauTagging TauTagging {
566set ParticleInputArray Delphes/allParticles
567set PartonInputArray Delphes/partons
568set JetInputArray JetEnergyScale/jets
569set DeltaR 0.4
570set TauPTMin 50.0
571set TauEtaMax 6
572# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
573add EfficiencyFormula {0} {0.01}
574add EfficiencyFormula {11} {0.005}
575add EfficiencyFormula {15} {0.6}
576}
577
578
579##########################
580# b-tagging
581##########################
582
583module BTagging BTagging { # need to check that this module exists or write it
584set ParticleInputArray Delphes/allParticles
585set PartonInputArray Delphes/partons
586set JetInputArray JetEnergyScale/jets
587
588set DeltaR 0.4
589set BPTMin 50.0
590set BEtaMax 6.0
591
592# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
593
594# default efficiency formula (misidentification rate)
595add EfficiencyFormula {0} {0.001}
596add EfficiencyFormula {4} {0.04}
597add EfficiencyFormula {5} {0.75}
598}
599
600##########################
601# c-tagging
602##########################
603
604
605module BTagging CTagging {
606set ParticleInputArray Delphes/allParticles
607set PartonInputArray Delphes/partons
608set JetInputArray JetEnergyScale/jets
609
610set DeltaR 0.4
611set BPTMin 30.0
612set BEtaMax 6.0
613
614# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
615add EfficiencyFormula {0} {0.01}
616add EfficiencyFormula {5} {0.10}
617add EfficiencyFormula {4} {0.25}
618}
619
620####################################
621# StatusPidFilter
622# This module removes all generated particles
623# except electrons, muons, taus, and status==3
624####################################
625
626module StatusPidFilter StatusPid {
627# set InputArray Delphes/stableParticles
628 set InputArray Delphes/allParticles
629 set OutputArray filteredParticles
630
631 set PTMin 0.5
632}
633
634
635#####################################################
636# Find uniquely identified photons/electrons/tau/jets
637#####################################################
638
639module UniqueObjectFinder UniqueObjectFinder {
640# earlier arrays take precedence over later ones
641# add InputArray InputArray OutputArray
642add InputArray PhotonIsolation/photons photons
643add InputArray ElectronIsolation/electrons electrons
644add InputArray MuonIsolation/muons muons
645add InputArray FastJetFinder/jets jets
646}
647
648##################
649# ROOT tree writer
650##################
651
652module TreeWriter TreeWriter {
653# add Branch InputArray BranchName BranchClass
654# add Branch Delphes/allParticles Particle GenParticle
655add Branch StatusPid/filteredParticles Particle GenParticle
656add Branch GenJetFinder/jets GenJet Jet
657
658# add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
659# add Branch Hcal/eflowNeutralHadrons NeutralHadron Tower
660
661add Branch Ecal/eflowPhotons Photon Photon
662add Branch ElectronEnergySmearing/electrons Electron Electron
663add Branch MuonMomentumSmearing/muons Muon Muon
664add Branch JetEnergyScale/jets Jet Jet
665add Branch MissingET/momentum MissingET MissingET
666add Branch ScalarHT/energy ScalarHT ScalarHT
667}
668
669
670