Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ 96d6bcf

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 96d6bcf was 96d6bcf, checked in by Michele Selvaggi <michele.selvaggi@…>, 7 years ago

added jet substr. in default reco sequence
  • Property mode set to 100644
File size: 25.8 KB
Line 
1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Heather Gray (CERN)
5# Filip Moortgat (CERN)
6# Michele Selvaggi (CERN)
7#
8# Released on: Nov 14th, 2016
9#
10# Configuration: FCC-hh baseline detector
11#
12#######################################
13# Order of execution of various modules
14#######################################
15
16set ExecutionPath {
17
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
24 ChargedHadronMomentumSmearing
25 ElectronMomentumSmearing
26 MuonMomentumSmearing
27
28 TrackMerger
29
30 ECal
31 HCal
32
33 Calorimeter
34 EFlowMerger
35 EFlowFilter
36
37 PhotonEfficiency
38 PhotonIsolation
39
40 ElectronFilter
41 ElectronIsolation
42
43 ChargedHadronFilter
44
45 MuonIsolation
46
47 NeutrinoFilter
48
49 MissingET
50 GenMissingET
51
52 GenJetFinder
53 FastJetFinder
54 FatJetFinder
55
56 JetEnergyScale
57
58 JetFlavorAssociation
59
60 BTagging
61 CTagging
62 TauTagging
63
64 ScalarHT
65
66 UniqueObjectFinder
67
68 TreeWriter
69}
70
71#################################
72# Propagate particles in cylinder
73#################################
74
75module ParticlePropagator ParticlePropagator {
76 set InputArray Delphes/stableParticles
77
78 set OutputArray stableParticles
79 set ChargedHadronOutputArray chargedHadrons
80 set ElectronOutputArray electrons
81 set MuonOutputArray muons
82
83 # radius of the magnetic field coverage, in m
84 set Radius 1.5
85 # half-length of the magnetic field coverage, in m
86 set HalfLength 5
87
88 # magnetic field
89 set Bz 4.0
90}
91
92####################################
93# Charged hadron tracking efficiency
94####################################
95
96module Efficiency ChargedHadronTrackingEfficiency {
97 set InputArray ParticlePropagator/chargedHadrons
98 set OutputArray chargedHadrons
99
100 # TBC (which eta_max ? which pT min?)
101
102 # tracking efficiency formula for charged hadrons
103
104 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
105(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
106(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
107(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
108(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
109(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
110(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
111(abs(eta) > 6.0) * (0.00)}
112
113}
114
115##############################
116# Electron tracking efficiency
117##############################
118
119module Efficiency ElectronTrackingEfficiency {
120 set InputArray ParticlePropagator/electrons
121 set OutputArray electrons
122
123# TBC (which eta_max ?)
124# putting same as charged hadrons for now...
125
126 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
127 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
128 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
129 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
130 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
131 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
132 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
133 (abs(eta) > 6.0) * (0.00)}
134
135}
136##########################
137# Muon tracking efficiency
138##########################
139
140module Efficiency MuonTrackingEfficiency {
141 set InputArray ParticlePropagator/muons
142 set OutputArray muons
143
144# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
145# what about high pT ?
146 # tracking efficiency formula for muons
147 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
148 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
149 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
150 (abs(eta) > 6.0) * (0.00)}
151
152}
153
154########################################
155# Momentum resolution for charged tracks
156########################################
157
158module MomentumSmearing ChargedHadronMomentumSmearing {
159 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
160 set OutputArray chargedHadrons
161
162 source momentumResolutionVsP.tcl
163}
164
165
166###################################
167# Momentum resolution for electrons
168###################################
169
170module MomentumSmearing ElectronMomentumSmearing {
171 set InputArray ElectronTrackingEfficiency/electrons
172 set OutputArray electrons
173
174 source momentumResolutionVsP.tcl
175}
176
177
178###############################
179# Momentum resolution for muons
180###############################
181
182module MomentumSmearing MuonMomentumSmearing {
183 set InputArray MuonTrackingEfficiency/muons
184 set OutputArray muons
185
186 # TBC for just putting tracker resolution/ need to add improvement at high pT
187
188 source muonMomentumResolutionVsP.tcl
189}
190
191##############
192# Track merger
193##############
194
195module Merger TrackMerger {
196# add InputArray InputArray
197 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
198 add InputArray ElectronMomentumSmearing/electrons
199 add InputArray MuonMomentumSmearing/muons
200 set OutputArray tracks
201}
202
203#############
204# ECAL
205#############
206
207# TBC : calos seems ok, check eta max value though.
208
209module SimpleCalorimeter ECal {
210 set ParticleInputArray ParticlePropagator/stableParticles
211 set TrackInputArray TrackMerger/tracks
212
213 set TowerOutputArray ecalTowers
214 set EFlowTrackOutputArray eflowTracks
215 set EFlowTowerOutputArray eflowPhotons
216
217 set IsEcal true
218
219 set EnergyMin 0.5
220 set EnergySignificanceMin 2.0
221
222 set SmearTowerCenter true
223
224 set pi [expr {acos(-1)}]
225
226 # lists of the edges of each tower in eta and phi
227 # each list starts with the lower edge of the first tower
228 # the list ends with the higher edged of the last tower
229
230 # 0.5 degree towers (5x5 mm^2)
231 set PhiBins {}
232 for {set i -256} {$i <= 256} {incr i} {
233 add PhiBins [expr {$i * $pi/256.0}]
234 }
235
236 # TBC
237 foreach eta {-2.4875 -2.475 -2.4625 -2.45 -2.4375 -2.425 -2.4125 -2.4 -2.3875 -2.375 -2.3625 -2.35 -2.3375 -2.325 -2.3125 -2.3 -2.2875 -2.275 -2.2625 -2.25 -2.2375 -2.225 -2.2125 -2.2 -2.1875 -2.175 -2.1625 -2.15 -2.1375 -2.125 -2.1125 -2.1 -2.0875 -2.075 -2.0625 -2.05 -2.0375 -2.025 -2.0125 -2 -1.9875 -1.975 -1.9625 -1.95 -1.9375 -1.925 -1.9125 -1.9 -1.8875 -1.875 -1.8625 -1.85 -1.8375 -1.825 -1.8125 -1.8 -1.7875 -1.775 -1.7625 -1.75 -1.7375 -1.725 -1.7125 -1.7 -1.6875 -1.675 -1.6625 -1.65 -1.6375 -1.625 -1.6125 -1.6 -1.5875 -1.575 -1.5625 -1.55 -1.5375 -1.525 -1.5125 -1.5 -1.4875 -1.475 -1.4625 -1.45 -1.4375 -1.425 -1.4125 -1.4 -1.3875 -1.375 -1.3625 -1.35 -1.3375 -1.325 -1.3125 -1.3 -1.2875 -1.275 -1.2625 -1.25 -1.2375 -1.225 -1.2125 -1.2 -1.1875 -1.175 -1.1625 -1.15 -1.1375 -1.125 -1.1125 -1.1 -1.0875 -1.075 -1.0625 -1.05 -1.0375 -1.025 -1.0125 -1.0 -0.9875 -0.975 -0.9625 -0.95 -0.9375 -0.925 -0.9125 -0.9 -0.8875 -0.875 -0.8625 -0.85 -0.8375 -0.825 -0.8125 -0.8 -0.7875 -0.775 -0.7625 -0.75 -0.7375 -0.725 -0.7125 -0.7 -0.6875 -0.675 -0.6625 -0.65 -0.6375 -0.625 -0.6125 -0.6 -0.5875 -0.575 -0.5625 -0.55 -0.5375 -0.525 -0.5125 -0.5 -0.4875 -0.475 -0.4625 -0.45 -0.4375 -0.425 -0.4125 -0.4 -0.3875 -0.375 -0.3625 -0.35 -0.3375 -0.325 -0.3125 -0.3 -0.2875 -0.275 -0.2625 -0.25 -0.2375 -0.225 -0.2125 -0.2 -0.1875 -0.175 -0.1625 -0.15 -0.1375 -0.125 -0.1125 -0.1 -0.0875 -0.075 -0.0625 -0.05 -0.0375 -0.025 -0.0125 0 0.0125 0.025 0.0375 0.05 0.0625 0.075 0.0875 0.1 0.1125 0.125 0.1375 0.15 0.1625 0.175 0.1875 0.2 0.2125 0.225 0.2375 0.25 0.2625 0.275 0.2875 0.3 0.3125 0.325 0.3375 0.35 0.3625 0.375 0.3875 0.4 0.4125 0.425 0.4375 0.45 0.4625 0.475 0.4875 0.5 0.5125 0.525 0.5375 0.55 0.5625 0.575 0.5875 0.6 0.6125 0.625 0.6375 0.65 0.6625 0.675 0.6875 0.7 0.7125 0.725 0.7375 0.75 0.7625 0.775 0.7875 0.8 0.8125 0.825 0.8375 0.85 0.8625 0.875 0.8875 0.9 0.9125 0.925 0.9375 0.95 0.9625 0.975 0.9875 1.0 1.0125 1.025 1.0375 1.05 1.0625 1.075 1.0875 1.1 1.1125 1.125 1.1375 1.15 1.1625 1.175 1.1875 1.2 1.2125 1.225 1.2375 1.25 1.2625 1.275 1.2875 1.3 1.3125 1.325 1.3375 1.35 1.3625 1.375 1.3875 1.4 1.4125 1.425 1.4375 1.45 1.4625 1.475 1.4875 1.5 1.5125 1.525 1.5375 1.55 1.5625 1.575 1.5875 1.6 1.6125 1.625 1.6375 1.65 1.6625 1.675 1.6875 1.7 1.7125 1.725 1.7375 1.75 1.7625 1.775 1.7875 1.8 1.8125 1.825 1.8375 1.85 1.8625 1.875 1.8875 1.9 1.9125 1.925 1.9375 1.95 1.9625 1.975 1.9875 2 2.0125 2.025 2.0375 2.05 2.0625 2.075 2.0875 2.1 2.1125 2.125 2.1375 2.15 2.1625 2.175 2.1875 2.2 2.2125 2.225 2.2375 2.25 2.2625 2.275 2.2875 2.3 2.3125 2.325 2.3375 2.35 2.3625 2.375 2.3875 2.4 2.4125 2.425 2.4375 2.45 2.4625 2.475 2.4875 2.5} {
238 add EtaPhiBins $eta $PhiBins
239 }
240 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
241 set PhiBins {}
242 for {set i -128} {$i <= 128} {incr i} {
243 add PhiBins [expr {$i * $pi/128.0}]
244 }
245 foreach eta {-3.975 -3.95 -3.925 -3.9 -3.875 -3.85 -3.825 -3.8 -3.775 -3.75 -3.725 -3.7 -3.675 -3.65 -3.625 -3.6 -3.575 -3.55 -3.525 -3.5 -3.475 -3.45 -3.425 -3.4 -3.375 -3.35 -3.325 -3.3 -3.275 -3.25 -3.225 -3.2 -3.175 -3.15 -3.125 -3.1 -3.075 -3.05 -3.025 -3.0 -2.975 -2.95 -2.925 -2.9 -2.875 -2.85 -2.825 -2.8 -2.775 -2.75 -2.725 -2.7 -2.675 -2.65 -2.625 -2.6 -2.575 -2.55 -2.525 -2.5 2.525 2.55 2.575 2.6 2.625 2.65 2.675 2.7 2.725 2.75 2.775 2.8 2.825 2.85 2.875 2.9 2.925 2.95 2.975 3.0 3.025 3.05 3.075 3.1 3.125 3.15 3.175 3.2 3.225 3.25 3.275 3.3 3.325 3.35 3.375 3.4 3.425 3.45 3.475 3.5 3.525 3.55 3.575 3.6 3.625 3.65 3.675 3.7 3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 3.925 3.95 3.975 4.0} {
246 add EtaPhiBins $eta $PhiBins
247 }
248
249 # 0.05 x 0.05 for eta between 4.0 and 6.0
250 set PhiBins {}
251 for {set i -64} {$i <= 64} {incr i} {
252 add PhiBins [expr {$i * $pi/64.0}]
253 }
254 foreach eta {-6.0 -5.95 -5.9 -5.85 -5.8 -5.75 -5.7 -5.65 -5.6 -5.55 -5.5 -5.45 -5.4 -5.35 -5.3 -5.25 -5.2 -5.15 -5.1 -5.05 -5.0 -4.95 -4.9 -4.85 -4.8 -4.75 -4.7 -4.65 -4.6 -4.55 -4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 -4 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 4.95 5.0 5.05 5.1 5.15 5.2 5.25 5.3 5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6.0} {
255 add EtaPhiBins $eta $PhiBins
256 }
257 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
258
259 add EnergyFraction {0} {0.0}
260 # energy fractions for e, gamma and pi0
261 add EnergyFraction {11} {1.0}
262 add EnergyFraction {22} {1.0}
263 add EnergyFraction {111} {1.0}
264 # energy fractions for muon, neutrinos and neutralinos
265 add EnergyFraction {12} {0.0}
266 add EnergyFraction {13} {0.0}
267 add EnergyFraction {14} {0.0}
268 add EnergyFraction {16} {0.0}
269 add EnergyFraction {1000022} {0.0}
270 add EnergyFraction {1000023} {0.0}
271 add EnergyFraction {1000025} {0.0}
272 add EnergyFraction {1000035} {0.0}
273 add EnergyFraction {1000045} {0.0}
274 # energy fractions for K0short and Lambda
275 # add EnergyFraction {310} {0.3}
276 # add EnergyFraction {3122} {0.3}
277
278 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
279 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
280 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
281
282
283}
284
285#############
286# HCAL
287#############
288
289module SimpleCalorimeter HCal {
290 set ParticleInputArray ParticlePropagator/stableParticles
291 set TrackInputArray ECal/eflowTracks
292
293 set TowerOutputArray hcalTowers
294 set EFlowTrackOutputArray eflowTracks
295 set EFlowTowerOutputArray eflowNeutralHadrons
296
297 set IsEcal false
298
299 set EnergyMin 1.0
300 set EnergySignificanceMin 2.0
301
302 set SmearTowerCenter true
303
304 set pi [expr {acos(-1)}]
305
306 # lists of the edges of each tower in eta and phi
307 # each list starts with the lower edge of the first tower
308 # the list ends with the higher edged of the last tower
309
310 # 6 degree towers
311 set PhiBins {}
312 for {set i -64} {$i <= 64} {incr i} {
313 add PhiBins [expr {$i * $pi/64.0}]
314 }
315
316 # TBC
317 foreach eta {-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.65 -1.6 -1.55 -1.5 -1.45 -1.4 -1.35 -1.3 -1.25 -1.2 -1.15 -1.1 -1.05 -1.0 -0.95 -0.9 -0.85 -0.8 -0.75 -0.7 -0.65 -0.6 -0.55 -0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 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} {
318 add EtaPhiBins $eta $PhiBins
319 }
320
321 # 0.1 x 0.1 between eta 2.5 and 4.0
322 set PhiBins {}
323 for {set i -32} {$i <= 32} {incr i} {
324 add PhiBins [expr {$i * $pi/32.0}]
325 }
326 foreach eta {-3.9 -3.8 -3.7 -3.6 -3.5 -3.4 -3.3 -3.2 -3.1 -3.0 -2.9 -2.8 -2.7 -2.6 -2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0} {
327 add EtaPhiBins $eta $PhiBins
328 }
329
330 # 0.2 x 0.2 between eta 4.0 and 6.0
331 set PhiBins {}
332 for {set i -16} {$i <= 16} {incr i} {
333 add PhiBins [expr {$i * $pi/16.0}]
334 }
335 foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 -4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
336 add EtaPhiBins $eta $PhiBins
337 }
338
339 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
340 add EnergyFraction {0} {1.0}
341 # energy fractions for e, gamma and pi0
342 add EnergyFraction {11} {0.0}
343 add EnergyFraction {22} {0.0}
344 add EnergyFraction {111} {0.0}
345 # energy fractions for muon, neutrinos and neutralinos
346 add EnergyFraction {12} {0.0}
347 add EnergyFraction {13} {0.0}
348 add EnergyFraction {14} {0.0}
349 add EnergyFraction {16} {0.0}
350 add EnergyFraction {1000022} {0.0}
351 add EnergyFraction {1000023} {0.0}
352 add EnergyFraction {1000025} {0.0}
353 add EnergyFraction {1000035} {0.0}
354 add EnergyFraction {1000045} {0.0}
355 # energy fractions for K0short and Lambda
356 # add EnergyFraction {310} {0.7}
357 # add EnergyFraction {3122} {0.7}
358
359 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
360 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
361 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
362}
363
364#################
365# Electron filter
366#################
367
368module PdgCodeFilter ElectronFilter {
369 set InputArray HCal/eflowTracks
370 set OutputArray electrons
371 set Invert true
372 add PdgCode {11}
373 add PdgCode {-11}
374}
375
376
377######################
378# ChargedHadronFilter
379######################
380
381module PdgCodeFilter ChargedHadronFilter {
382 set InputArray HCal/eflowTracks
383 set OutputArray chargedHadrons
384
385 add PdgCode {11}
386 add PdgCode {-11}
387 add PdgCode {13}
388 add PdgCode {-13}
389}
390
391
392###################################################
393# Tower Merger (in case not using e-flow algorithm)
394###################################################
395
396module Merger Calorimeter {
397# add InputArray InputArray
398 add InputArray ECal/ecalTowers
399 add InputArray HCal/hcalTowers
400 set OutputArray towers
401}
402
403####################
404# Energy flow merger
405####################
406
407module Merger EFlowMerger {
408# add InputArray InputArray
409 add InputArray HCal/eflowTracks
410 add InputArray ECal/eflowPhotons
411 add InputArray HCal/eflowNeutralHadrons
412 set OutputArray eflow
413}
414
415######################
416# EFlowFilter
417######################
418
419module PdgCodeFilter EFlowFilter {
420 set InputArray EFlowMerger/eflow
421 set OutputArray eflow
422
423 add PdgCode {11}
424 add PdgCode {-11}
425 add PdgCode {13}
426 add PdgCode {-13}
427}
428
429
430###################
431# Missing ET merger
432###################
433
434module Merger MissingET {
435# add InputArray InputArray
436 add InputArray EFlowMerger/eflow
437 set MomentumOutputArray momentum
438}
439
440
441
442##################
443# Scalar HT merger
444##################
445
446module Merger ScalarHT {
447# add InputArray InputArray
448 add InputArray EFlowMerger/eflow
449 set EnergyOutputArray energy
450}
451
452#################
453# Neutrino Filter
454#################
455
456module PdgCodeFilter NeutrinoFilter {
457
458 set InputArray Delphes/stableParticles
459 set OutputArray filteredParticles
460
461 set PTMin 0.0
462
463 add PdgCode {12}
464 add PdgCode {14}
465 add PdgCode {16}
466 add PdgCode {-12}
467 add PdgCode {-14}
468 add PdgCode {-16}
469
470}
471
472
473#####################
474# MC truth jet finder
475#####################
476
477# TBC: is jet radius fine?
478
479module FastJetFinder GenJetFinder {
480# set InputArray NeutrinoFilter/filteredParticles
481 set InputArray Delphes/stableParticles
482
483 set OutputArray jets
484
485 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
486 set JetAlgorithm 6
487 set ParameterR 0.4
488
489 set JetPTMin 5.0
490}
491
492#########################
493# Gen Missing ET merger
494########################
495
496module Merger GenMissingET {
497
498# add InputArray InputArray
499 add InputArray NeutrinoFilter/filteredParticles
500 set MomentumOutputArray momentum
501}
502
503
504
505############
506# Jet finder
507############
508
509# TBC need to include jet substructure variables
510# TBC is jet radius fine?
511
512module FastJetFinder FastJetFinder {
513# set InputArray Calorimeter/towers
514 set InputArray EFlowMerger/eflow
515
516 set OutputArray jets
517
518 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
519 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
520
521 set JetAlgorithm 6
522 set ParameterR 0.4
523
524 set JetPTMin 30.0
525}
526
527##################
528# Fat Jet finder
529##################
530
531module FastJetFinder FatJetFinder {
532 set InputArray EFlowMerger/eflow
533
534 set OutputArray jets
535
536 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
537 set JetAlgorithm 6
538 set ParameterR 0.8
539
540 set ComputeNsubjettiness 1
541 set Beta 1.0
542 set AxisMode 4
543
544 set ComputeTrimming 1
545 set RTrim 0.2
546 set PtFracTrim 0.05
547
548 set ComputePruning 1
549 set ZcutPrun 0.1
550 set RcutPrun 0.5
551 set RPrun 0.8
552
553 set ComputeSoftDrop 1
554 set BetaSoftDrop 0.0
555 set SymmetryCutSoftDrop 0.1
556 set R0SoftDrop 0.8
557
558 set JetPTMin 200.0
559}
560
561
562
563##################
564# Jet Energy Scale
565##################
566
567module EnergyScale JetEnergyScale {
568 set InputArray FastJetFinder/jets
569 set OutputArray jets
570
571 # scale formula for jets
572 set ScaleFormula {1.00}
573}
574
575
576########################
577# Jet Flavor Association
578########################
579
580module JetFlavorAssociation JetFlavorAssociation {
581
582 set PartonInputArray Delphes/partons
583 set ParticleInputArray Delphes/allParticles
584 set ParticleLHEFInputArray Delphes/allParticlesLHEF
585 set JetInputArray JetEnergyScale/jets
586
587 set DeltaR 0.5
588 set PartonPTMin 5.0
589 set PartonEtaMax 6.0
590
591}
592
593###################
594# Photon efficiency
595###################
596
597module Efficiency PhotonEfficiency {
598 set InputArray ECal/eflowPhotons
599 set OutputArray photons
600
601 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
602
603 set EfficiencyFormula {
604 (pt <= 1.0) * (0.00) + \
605 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
606 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
607 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
608
609 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
610 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
611 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
612
613 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
614 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
615 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
616 (abs(eta) > 6.0) * (0.00)}
617
618}
619
620##################
621# Photon isolation
622##################
623
624# TBC: check values for iso cuts
625
626module Isolation PhotonIsolation {
627 set CandidateInputArray PhotonEfficiency/photons
628 set IsolationInputArray EFlowFilter/eflow
629
630 set OutputArray photons
631
632 set DeltaRMax 0.3
633
634 set PTMin 0.5
635
636 set PTRatioMax 0.1
637}
638
639
640####################
641# Electron isolation
642####################
643
644# TBC: check values for iso cuts
645
646module Isolation ElectronIsolation {
647 set CandidateInputArray ElectronFilter/electrons
648 set IsolationInputArray EFlowFilter/eflow
649
650 set OutputArray electrons
651
652 set DeltaRMax 0.3
653
654 set PTMin 0.5
655
656 set PTRatioMax 0.1
657}
658
659
660################
661# Muon isolation
662################
663
664# TBC: check values for iso cuts
665
666module Isolation MuonIsolation {
667 set CandidateInputArray MuonMomentumSmearing/muons
668 set IsolationInputArray EFlowFilter/eflow
669
670 set OutputArray muons
671
672 set DeltaRMax 0.3
673
674 set PTMin 0.5
675
676 set PTRatioMax 0.2
677}
678
679
680###########
681# b-tagging
682###########
683
684module BTagging BTagging {
685 set JetInputArray JetEnergyScale/jets
686
687 set BitNumber 0
688
689 add EfficiencyFormula {0} {
690
691 (pt <= 10.0) * (0.00) +
692 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
693 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
694 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
695 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
696 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
697 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
698 (abs(eta) > 4.0) * (0.00)}
699
700 add EfficiencyFormula {4} {
701
702 (pt <= 10.0) * (0.00) +
703 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
704 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
705 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
706 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
707 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
708 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
709 (abs(eta) > 4.0) * (0.00)}
710
711 add EfficiencyFormula {5} {
712
713 (pt <= 10.0) * (0.00) +
714 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
715 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
716 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
717 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
718 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
719 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
720 (abs(eta) >= 4.0) * (0.00)}
721
722}
723
724###########
725# c-tagging
726###########
727
728module BTagging CTagging {
729 set JetInputArray JetEnergyScale/jets
730
731 set BitNumber 1
732
733 add EfficiencyFormula {0} {
734
735 (pt <= 10.0) * (0.00) +
736 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
737 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
738
739 add EfficiencyFormula {4} {
740
741 (pt <= 10.0) * (0.00) +
742 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
743 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
744
745 add EfficiencyFormula {5} {
746
747 (pt <= 10.0) * (0.00) +
748 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
749 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
750
751}
752
753
754#############
755# tau-tagging
756#############
757
758
759module TauTagging TauTagging {
760 set ParticleInputArray Delphes/allParticles
761 set PartonInputArray Delphes/partons
762 set JetInputArray JetEnergyScale/jets
763
764 set DeltaR 0.5
765
766 set TauPTMin 1.0
767
768 set TauEtaMax 4.0
769
770 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
771 add EfficiencyFormula {0} {
772
773 (pt <= 10.0) * (0.00) +
774 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
775 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
776 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
777 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
778 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
779 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
780 (abs(eta) > 4.0) * (0.00)}
781
782 add EfficiencyFormula {11} {
783
784 (pt <= 10.0) * (0.00) +
785 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
786 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
787 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
788 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
789 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
790 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
791 (abs(eta) > 4.0) * (0.00)}
792
793 add EfficiencyFormula {15} {
794
795 (pt <= 10.0) * (0.00) +
796 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
797 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
798 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
799 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
800 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
801 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
802 (abs(eta) > 4.0) * (0.00)}
803
804}
805
806#####################################################
807# Find uniquely identified photons/electrons/tau/jets
808#####################################################
809
810module UniqueObjectFinder UniqueObjectFinder {
811# earlier arrays take precedence over later ones
812# add InputArray InputArray OutputArray
813 add InputArray PhotonIsolation/photons photons
814 add InputArray ElectronIsolation/electrons electrons
815 add InputArray MuonIsolation/muons muons
816 add InputArray JetEnergyScale/jets jets
817
818}
819
820
821##################
822# ROOT tree writer
823##################
824
825module TreeWriter TreeWriter {
826# add Branch InputArray BranchName BranchClass
827 add Branch Delphes/allParticles Particle GenParticle
828
829 add Branch GenJetFinder/jets GenJet Jet
830 add Branch GenMissingET/momentum GenMissingET MissingET
831
832 add Branch TrackMerger/tracks Track Track
833 add Branch Calorimeter/towers Tower Tower
834
835 add Branch HCal/eflowTracks EFlowTrack Track
836 add Branch ECal/eflowPhotons EFlowPhoton Tower
837 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
838
839 add Branch UniqueObjectFinder/photons Photon Photon
840 add Branch UniqueObjectFinder/electrons Electron Electron
841 add Branch UniqueObjectFinder/muons Muon Muon
842 add Branch UniqueObjectFinder/jets Jet Jet
843
844 add Branch FatJetFinder/jets FatJet Jet
845
846 add Branch MissingET/momentum MissingET MissingET
847 add Branch ScalarHT/energy ScalarHT ScalarHT
848}
849
Note: See TracBrowser for help on using the repository browser.