Fork me on GitHub

Ticket #406: delphes_card_SMEAR.dat

File delphes_card_SMEAR.dat, 18.8 KB (added by Michele Selvaggi, 10 years ago)
Line 
1set RandomSeed 123
2
3#######################################
4# Order of execution of various modules
5#######################################
6
7set ExecutionPath {
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronEnergySmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 TrackVertexSmearing
20
21 Calorimeter
22 EFlowMerger
23
24 PhotonEfficiency
25 PhotonIsolation
26
27 ElectronEfficiency
28 ElectronIsolation
29
30 MuonEfficiency
31 MuonIsolation
32
33 MissingET
34
35 GenJetFinder
36 FastJetFinder
37
38 JetEnergyScale
39
40 BTagging
41 TauTagging
42
43 UniqueObjectFinder
44
45 ScalarHT
46
47 TreeWriter
48}
49
50#################################
51# Propagate particles in cylinder
52#################################
53
54module ParticlePropagator ParticlePropagator {
55 set InputArray Delphes/stableParticles
56
57 set OutputArray stableParticles
58 set ChargedHadronOutputArray chargedHadrons
59 set ElectronOutputArray electrons
60 set MuonOutputArray muons
61
62 # radius of the magnetic field coverage, in m
63 set Radius 1.15
64 # half-length of the magnetic field coverage, in m
65 set HalfLength 3.51
66
67 # magnetic field
68 set Bz 2.0
69}
70
71
72####################################
73# Charged hadron tracking efficiency
74####################################
75
76module Efficiency ChargedHadronTrackingEfficiency {
77 set InputArray ParticlePropagator/chargedHadrons
78 set OutputArray chargedHadrons
79
80 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
81
82 # tracking efficiency formula for charged hadrons
83 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
84 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
85 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
86 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
87 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
88 (abs(eta) > 2.5) * (0.00)}
89}
90
91##############################
92# Electron tracking efficiency
93##############################
94
95module Efficiency ElectronTrackingEfficiency {
96 set InputArray ParticlePropagator/electrons
97 set OutputArray electrons
98
99 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
100
101 # tracking efficiency formula for electrons
102 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
103 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
104 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
105 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
106 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
107 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
108 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
109 (abs(eta) > 2.5) * (0.00)}
110}
111
112##########################
113# Muon tracking efficiency
114##########################
115
116module Efficiency MuonTrackingEfficiency {
117 set InputArray ParticlePropagator/muons
118 set OutputArray muons
119
120 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
121
122 # tracking efficiency formula for muons
123 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
124 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
125 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
126 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
128 (abs(eta) > 2.5) * (0.00)}
129}
130
131########################################
132# Momentum resolution for charged tracks
133########################################
134
135module MomentumSmearing ChargedHadronMomentumSmearing {
136 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
137 set OutputArray chargedHadrons
138
139 # set ResolutionFormula {resolution formula as a function of eta and pt}
140
141 # resolution formula for charged hadrons
142 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
143 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
144 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
145 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
148 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
150}
151
152#################################
153# Energy resolution for electrons
154#################################
155
156module EnergySmearing ElectronEnergySmearing {
157 set InputArray ElectronTrackingEfficiency/electrons
158 set OutputArray electrons
159
160 # set ResolutionFormula {resolution formula as a function of eta and energy}
161
162 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
163 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
164 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
165 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
166
167}
168
169###############################
170# Momentum resolution for muons
171###############################
172
173module MomentumSmearing MuonMomentumSmearing {
174 set InputArray MuonTrackingEfficiency/muons
175 set OutputArray muons
176
177 # set ResolutionFormula {resolution formula as a function of eta and pt}
178
179 # resolution formula for muons
180 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
181 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
182 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
183 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
188}
189
190##############
191# Track merger
192##############
193
194module Merger TrackMerger {
195# add InputArray InputArray
196 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
197 add InputArray ElectronEnergySmearing/electrons
198 add InputArray MuonMomentumSmearing/muons
199 set OutputArray tracks
200}
201
202################################
203# Track vertex smearing
204################################
205
206module TrackVertexSmearing TrackVertexSmearing {
207 set InputArray TrackMerger/tracks
208 set OutputArray tracks
209
210# absolute impact parameter smearing formula (in mm) as a function of pt and eta
211# FNesti 24.2.2015 - built from slide 29 of
212# https://indico.cern.ch/event/300851/session/7/contribution/42/material/slides/0.pdf
213# powerlaw increase with low pt, exponential increase with eta beyond 1.5
214#
215set ResolutionFormula {
216 (100)
217 }
218}
219
220
221#############
222# Calorimeter
223#############
224
225module Calorimeter Calorimeter {
226 set ParticleInputArray ParticlePropagator/stableParticles
227 set TrackInputArray TrackVertexSmearing/tracks
228
229 set TowerOutputArray towers
230 set PhotonOutputArray photons
231
232 set EFlowTrackOutputArray eflowTracks
233 set EFlowPhotonOutputArray eflowPhotons
234 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
235
236 set pi [expr {acos(-1)}]
237
238 # lists of the edges of each tower in eta and phi
239 # each list starts with the lower edge of the first tower
240 # the list ends with the higher edged of the last tower
241
242 # 10 degrees towers
243 set PhiBins {}
244 for {set i -18} {$i <= 18} {incr i} {
245 add PhiBins [expr {$i * $pi/18.0}]
246 }
247 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
248 add EtaPhiBins $eta $PhiBins
249 }
250
251 # 20 degrees towers
252 set PhiBins {}
253 for {set i -9} {$i <= 9} {incr i} {
254 add PhiBins [expr {$i * $pi/9.0}]
255 }
256 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
257 add EtaPhiBins $eta $PhiBins
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 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
281 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
282 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
283 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
284 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
285 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
286
287 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
288 # http://arxiv.org/pdf/hep-ex/0004009v1
289 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
290 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
291 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
292 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
293}
294
295####################
296# Energy flow merger
297####################
298
299module Merger EFlowMerger {
300# add InputArray InputArray
301 add InputArray Calorimeter/eflowTracks
302 add InputArray Calorimeter/eflowPhotons
303 add InputArray Calorimeter/eflowNeutralHadrons
304 set OutputArray eflow
305}
306
307###################
308# Photon efficiency
309###################
310
311module Efficiency PhotonEfficiency {
312 set InputArray Calorimeter/photons
313 set OutputArray photons
314
315 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
316
317 # efficiency formula for photons
318 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
319 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
320 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
321 (abs(eta) > 2.5) * (0.00)}
322}
323
324##################
325# Photon isolation
326##################
327
328module Isolation PhotonIsolation {
329 set CandidateInputArray PhotonEfficiency/photons
330 set IsolationInputArray EFlowMerger/eflow
331
332 set OutputArray photons
333
334 set DeltaRMax 0.5
335
336 set PTMin 0.5
337
338 set PTRatioMax 0.1
339}
340
341#####################
342# Electron efficiency
343#####################
344
345module Efficiency ElectronEfficiency {
346 set InputArray ElectronEnergySmearing/electrons
347 set OutputArray electrons
348
349 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
350
351 # efficiency formula for electrons
352 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
353 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
354 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
355 (abs(eta) > 2.5) * (0.00)}
356}
357
358####################
359# Electron isolation
360####################
361
362module Isolation ElectronIsolation {
363 set CandidateInputArray ElectronEfficiency/electrons
364 set IsolationInputArray EFlowMerger/eflow
365
366 set OutputArray electrons
367
368 set DeltaRMax 0.3
369
370 set PTMin 0.4
371
372 set PTRatioMax 0.1
373}
374
375#################
376# Muon efficiency
377#################
378
379module Efficiency MuonEfficiency {
380 set InputArray MuonMomentumSmearing/muons
381 set OutputArray muons
382
383 # set EfficiencyFormula {efficiency as a function of eta and pt}
384
385 # efficiency formula for muons
386 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
387 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
388 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
389 (abs(eta) > 2.7) * (0.00)}
390}
391
392################
393# Muon isolation
394################
395
396module Isolation MuonIsolation {
397 set CandidateInputArray MuonEfficiency/muons
398 set IsolationInputArray EFlowMerger/eflow
399
400 set OutputArray muons
401
402 set DeltaRMax 0.3
403
404 set PTMin 1
405
406 set PTRatioMax 0.07
407}
408
409###################
410# Missing ET merger
411###################
412
413module Merger MissingET {
414# add InputArray InputArray
415 add InputArray EFlowMerger/eflow
416 set MomentumOutputArray momentum
417}
418
419##################
420# Scalar HT merger
421##################
422
423module Merger ScalarHT {
424# add InputArray InputArray
425 add InputArray UniqueObjectFinder/jets
426 add InputArray UniqueObjectFinder/electrons
427 add InputArray UniqueObjectFinder/photons
428 add InputArray UniqueObjectFinder/muons
429 set EnergyOutputArray energy
430}
431
432#####################
433# MC truth jet finder
434#####################
435
436module FastJetFinder GenJetFinder {
437 set InputArray Delphes/stableParticles
438
439 set OutputArray jets
440
441 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
442 set JetAlgorithm 6
443 set ParameterR 0.4
444
445 set JetPTMin 20.0
446}
447
448############
449# Jet finder
450############
451
452module FastJetFinder FastJetFinder {
453 set InputArray Calorimeter/towers
454
455 set OutputArray jets
456
457 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
458 set JetAlgorithm 6
459 set ParameterR 0.4
460
461 set JetPTMin 20.0
462}
463
464##################
465# Jet Energy Scale
466##################
467
468module EnergyScale JetEnergyScale {
469 set InputArray FastJetFinder/jets
470 set OutputArray jets
471
472 # scale formula for jets
473 set ScaleFormula {1.00}
474}
475
476###########
477# b-tagging
478###########
479
480module BTagging BTagging {
481 set PartonInputArray Delphes/partons
482 set JetInputArray JetEnergyScale/jets
483
484 set BitNumber 0
485
486 set DeltaR 0.5
487
488 set PartonPTMin 1.0
489
490 set PartonEtaMax 2.5
491
492 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
493 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
494 # gluon's PDG code has the lowest priority
495
496 # default efficiency formula (misidentification rate)
497 add EfficiencyFormula {0} {0.001}
498
499 # efficiency formula for c-jets (misidentification rate)
500 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
501 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
502 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
503 (abs(eta) > 2.5) * (0.000)}
504
505 # efficiency formula for b-jets
506 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
507 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
508 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
509 (abs(eta) > 2.5) * (0.000)}
510}
511
512module TauTagging TauTagging {
513 set ParticleInputArray Delphes/allParticles
514 set PartonInputArray Delphes/partons
515 set JetInputArray JetEnergyScale/jets
516
517 set DeltaR 0.5
518
519 set TauPTMin 1.0
520
521 set TauEtaMax 2.5
522
523 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
524
525 # default efficiency formula (misidentification rate)
526 add EfficiencyFormula {0} {0.001}
527 # efficiency formula for tau-jets
528 add EfficiencyFormula {15} {0.4}
529}
530
531#####################################################
532# Find uniquely identified photons/electrons/tau/jets
533#####################################################
534
535module UniqueObjectFinder UniqueObjectFinder {
536# earlier arrays take precedence over later ones
537# add InputArray InputArray OutputArray
538 add InputArray PhotonIsolation/photons photons
539 add InputArray ElectronIsolation/electrons electrons
540 add InputArray MuonIsolation/muons muons
541 add InputArray JetEnergyScale/jets jets
542}
543
544##################
545# ROOT tree writer
546##################
547
548# tracks, towers and eflow objects are not stored by default in the output.
549# if needed (for jet constituent or other studies), uncomment the relevant
550# "add Branch ..." lines.
551
552module TreeWriter TreeWriter {
553# add Branch InputArray BranchName BranchClass
554 add Branch Delphes/allParticles Particle GenParticle
555
556 add Branch TrackMerger/tracks Track Track
557 add Branch TrackVertexSmearing/tracks TrackSmeared Track
558
559 add Branch Calorimeter/towers Tower Tower
560
561 add Branch Calorimeter/eflowTracks EFlowTrack Track
562 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
563 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
564
565 add Branch GenJetFinder/jets GenJet Jet
566 add Branch UniqueObjectFinder/jets Jet Jet
567 add Branch UniqueObjectFinder/electrons Electron Electron
568 add Branch UniqueObjectFinder/photons Photon Photon
569 add Branch UniqueObjectFinder/muons Muon Muon
570 add Branch MissingET/momentum MissingET MissingET
571 add Branch ScalarHT/energy ScalarHT ScalarHT
572}
573