Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ b4ec6ac

ImprovedOutputFile Timing dual_readout llp
Last change on this file since b4ec6ac was 5348ed9, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

first commit of new generation of FCC detector cards
  • Property mode set to 100644
File size: 21.3 KB
RevLine 
[5348ed9]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: Oct 02th, 2016
9#
10# Configuration: FCC-hh generic detector without dipole in the forward region, no pile-up
11#
12#######################################
13# Order of execution of various modules
14#######################################
15
16set ExecutionPath {
17 ParticlePropagator
18
19 ChargedHadronTrackingEfficiency
20 ElectronTrackingEfficiency
21 MuonTrackingEfficiency
22
23 ChargedHadronMomentumSmearing
24 ElectronMomentumSmearing
25 MuonMomentumSmearing
26
27 TrackMerger
28
29 ECal
30 HCal
31
32 Calorimeter
33 EFlowMerger
34
35 PhotonEfficiency
36 PhotonIsolation
37
38 ElectronFilter
39 ElectronIsolation
40
41 MuonIsolation
42
43 NeutrinoFilter
44
45 MissingET
46 GenMissingET
47
48 GenJetFinder
49 FastJetFinder
50
51 JetEnergyScale
52
53 JetFlavorAssociation
54
55 BTagging
56 CTagging
57 TauTagging
58
59 ScalarHT
60
61 UniqueObjectFinder
62
63 TreeWriter
64}
65
66#################################
67# Propagate particles in cylinder
68#################################
69
70module ParticlePropagator ParticlePropagator {
71 set InputArray Delphes/stableParticles
72
73 set OutputArray stableParticles
74 set ChargedHadronOutputArray chargedHadrons
75 set ElectronOutputArray electrons
76 set MuonOutputArray muons
77
78 # radius of the magnetic field coverage, in m
79 set Radius 1.5
80 # half-length of the magnetic field coverage, in m
81 set HalfLength 5
82
83 # magnetic field
84 set Bz 4.0
85}
86
87####################################
88# Charged hadron tracking efficiency
89####################################
90
91module Efficiency ChargedHadronTrackingEfficiency {
92 set InputArray ParticlePropagator/chargedHadrons
93 set OutputArray chargedHadrons
94
95 # TBC (which eta_max ? which pT min?)
96
97 # tracking efficiency formula for charged hadrons
98
99 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
100(abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
101(abs(eta) <= 1.5) * (pt > 1) * (0.95) + \
102(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
103(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
104(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
105(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
106(abs(eta) > 6.0) * (0.00)}
107
108}
109
110##############################
111# Electron tracking efficiency
112##############################
113
114module Efficiency ElectronTrackingEfficiency {
115 set InputArray ParticlePropagator/electrons
116 set OutputArray electrons
117
118# TBC (which eta_max ?)
119# putting same as charged hadrons for now...
120
121 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
122 (abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
123 (abs(eta) <= 1.5) * (pt > 1) * (0.95) + \
124 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
125 (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
126 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
127 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
128 (abs(eta) > 6.0) * (0.00)}
129
130}
131##########################
132# Muon tracking efficiency
133##########################
134
135module Efficiency MuonTrackingEfficiency {
136 set InputArray ParticlePropagator/muons
137 set OutputArray muons
138
139# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
140# what about high pT ?
141 # tracking efficiency formula for muons
142 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
143 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
144 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
145 (abs(eta) > 6.0) * (0.00)}
146
147}
148
149########################################
150# Momentum resolution for charged tracks
151########################################
152
153module MomentumSmearing ChargedHadronMomentumSmearing {
154 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
155 set OutputArray chargedHadrons
156
157 source momentumResolutionVsP.tcl
158}
159
160
161###################################
162# Momentum resolution for electrons
163###################################
164
165module MomentumSmearing ElectronMomentumSmearing {
166 set InputArray ElectronTrackingEfficiency/electrons
167 set OutputArray electrons
168
169 source momentumResolutionVsP.tcl
170}
171
172
173###############################
174# Momentum resolution for muons
175###############################
176
177module MomentumSmearing MuonMomentumSmearing {
178 set InputArray MuonTrackingEfficiency/muons
179 set OutputArray muons
180
181 # TBC for just putting tracker resolution/ need to add improvement at high pT
182
183 source muonMomentumResolutionVsP.tcl
184}
185
186##############
187# Track merger
188##############
189
190module Merger TrackMerger {
191# add InputArray InputArray
192 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
193 add InputArray ElectronMomentumSmearing/electrons
194 add InputArray MuonMomentumSmearing/muons
195 set OutputArray tracks
196}
197
198#############
199# ECAL
200#############
201
202# TBC : calos seems ok, check eta max value though.
203
204module SimpleCalorimeter ECal {
205 set ParticleInputArray ParticlePropagator/stableParticles
206 set TrackInputArray TrackMerger/tracks
207
208 set TowerOutputArray ecalTowers
209 set EFlowTrackOutputArray eflowTracks
210 set EFlowTowerOutputArray eflowPhotons
211
212 set IsEcal true
213
214 set EnergyMin 0.5
215 set EnergySignificanceMin 2.0
216
217 set SmearTowerCenter true
218
219 set pi [expr {acos(-1)}]
220
221 # lists of the edges of each tower in eta and phi
222 # each list starts with the lower edge of the first tower
223 # the list ends with the higher edged of the last tower
224
225 # 0.5 degree towers (5x5 mm^2)
226 set PhiBins {}
227 for {set i -256} {$i <= 256} {incr i} {
228 add PhiBins [expr {$i * $pi/256.0}]
229 }
230
231 # TBC
232 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} {
233 add EtaPhiBins $eta $PhiBins
234 }
235 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
236 set PhiBins {}
237 for {set i -128} {$i <= 128} {incr i} {
238 add PhiBins [expr {$i * $pi/128.0}]
239 }
240 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} {
241 add EtaPhiBins $eta $PhiBins
242 }
243
244 # 0.05 x 0.05 for eta between 4.0 and 6.0
245 set PhiBins {}
246 for {set i -64} {$i <= 64} {incr i} {
247 add PhiBins [expr {$i * $pi/64.0}]
248 }
249 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} {
250 add EtaPhiBins $eta $PhiBins
251 }
252 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
253
254 add EnergyFraction {0} {0.0}
255 # energy fractions for e, gamma and pi0
256 add EnergyFraction {11} {1.0}
257 add EnergyFraction {22} {1.0}
258 add EnergyFraction {111} {1.0}
259 # energy fractions for muon, neutrinos and neutralinos
260 add EnergyFraction {12} {0.0}
261 add EnergyFraction {13} {0.0}
262 add EnergyFraction {14} {0.0}
263 add EnergyFraction {16} {0.0}
264 add EnergyFraction {1000022} {0.0}
265 add EnergyFraction {1000023} {0.0}
266 add EnergyFraction {1000025} {0.0}
267 add EnergyFraction {1000035} {0.0}
268 add EnergyFraction {1000045} {0.0}
269 # energy fractions for K0short and Lambda
270 # add EnergyFraction {310} {0.3}
271 # add EnergyFraction {3122} {0.3}
272
273 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
274 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
275 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
276
277
278}
279
280#############
281# HCAL
282#############
283
284module SimpleCalorimeter HCal {
285 set ParticleInputArray ParticlePropagator/stableParticles
286 set TrackInputArray ECal/eflowTracks
287
288 set TowerOutputArray hcalTowers
289 set EFlowTrackOutputArray eflowTracks
290 set EFlowTowerOutputArray eflowNeutralHadrons
291
292 set IsEcal false
293
294 set EnergyMin 1.0
295 set EnergySignificanceMin 2.0
296
297 set SmearTowerCenter true
298
299 set pi [expr {acos(-1)}]
300
301 # lists of the edges of each tower in eta and phi
302 # each list starts with the lower edge of the first tower
303 # the list ends with the higher edged of the last tower
304
305 # 6 degree towers
306 set PhiBins {}
307 for {set i -64} {$i <= 64} {incr i} {
308 add PhiBins [expr {$i * $pi/64.0}]
309 }
310
311 # TBC
312 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} {
313 add EtaPhiBins $eta $PhiBins
314 }
315
316 # 0.1 x 0.1 between eta 2.5 and 4.0
317 set PhiBins {}
318 for {set i -32} {$i <= 32} {incr i} {
319 add PhiBins [expr {$i * $pi/32.0}]
320 }
321 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} {
322 add EtaPhiBins $eta $PhiBins
323 }
324
325 # 0.2 x 0.2 between eta 4.0 and 6.0
326 set PhiBins {}
327 for {set i -16} {$i <= 16} {incr i} {
328 add PhiBins [expr {$i * $pi/16.0}]
329 }
330 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} {
331 add EtaPhiBins $eta $PhiBins
332 }
333
334 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
335 add EnergyFraction {0} {1.0}
336 # energy fractions for e, gamma and pi0
337 add EnergyFraction {11} {0.0}
338 add EnergyFraction {22} {0.0}
339 add EnergyFraction {111} {0.0}
340 # energy fractions for muon, neutrinos and neutralinos
341 add EnergyFraction {12} {0.0}
342 add EnergyFraction {13} {0.0}
343 add EnergyFraction {14} {0.0}
344 add EnergyFraction {16} {0.0}
345 add EnergyFraction {1000022} {0.0}
346 add EnergyFraction {1000023} {0.0}
347 add EnergyFraction {1000025} {0.0}
348 add EnergyFraction {1000035} {0.0}
349 add EnergyFraction {1000045} {0.0}
350 # energy fractions for K0short and Lambda
351 # add EnergyFraction {310} {0.7}
352 # add EnergyFraction {3122} {0.7}
353
354 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
355 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
356 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
357}
358
359#################
360# Electron filter
361#################
362
363module PdgCodeFilter ElectronFilter {
364 set InputArray HCal/eflowTracks
365 set OutputArray electrons
366 set Invert true
367 add PdgCode {11}
368 add PdgCode {-11}
369}
370
371###################################################
372# Tower Merger (in case not using e-flow algorithm)
373###################################################
374
375module Merger Calorimeter {
376# add InputArray InputArray
377 add InputArray ECal/ecalTowers
378 add InputArray HCal/hcalTowers
379 set OutputArray towers
380}
381
382####################
383# Energy flow merger
384####################
385
386module Merger EFlowMerger {
387# add InputArray InputArray
388 add InputArray HCal/eflowTracks
389 add InputArray ECal/eflowPhotons
390 add InputArray HCal/eflowNeutralHadrons
391 set OutputArray eflow
392}
393
394
395###################
396# Missing ET merger
397###################
398
399module Merger MissingET {
400# add InputArray InputArray
401 add InputArray EFlowMerger/eflow
402 set MomentumOutputArray momentum
403}
404
405
406
407##################
408# Scalar HT merger
409##################
410
411module Merger ScalarHT {
412# add InputArray InputArray
413 add InputArray EFlowMerger/eflow
414 set EnergyOutputArray energy
415}
416
417#################
418# Neutrino Filter
419#################
420
421module PdgCodeFilter NeutrinoFilter {
422
423 set InputArray Delphes/stableParticles
424 set OutputArray filteredParticles
425
426 set PTMin 0.0
427
428 add PdgCode {12}
429 add PdgCode {14}
430 add PdgCode {16}
431 add PdgCode {-12}
432 add PdgCode {-14}
433 add PdgCode {-16}
434
435}
436
437
438#####################
439# MC truth jet finder
440#####################
441
442# TBC: is jet radius fine?
443
444module FastJetFinder GenJetFinder {
445# set InputArray NeutrinoFilter/filteredParticles
446 set InputArray Delphes/stableParticles
447
448 set OutputArray jets
449
450 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
451 set JetAlgorithm 6
452 set ParameterR 0.4
453
454 set JetPTMin 5.0
455}
456
457#########################
458# Gen Missing ET merger
459########################
460
461module Merger GenMissingET {
462
463# add InputArray InputArray
464 add InputArray NeutrinoFilter/filteredParticles
465 set MomentumOutputArray momentum
466}
467
468
469
470############
471# Jet finder
472############
473
474# TBC need to include jet substructure variables
475# TBC is jet radius fine?
476
477module FastJetFinder FastJetFinder {
478# set InputArray Calorimeter/towers
479 set InputArray EFlowMerger/eflow
480
481 set OutputArray jets
482
483 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
484 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
485
486 set JetAlgorithm 6
487 set ParameterR 0.4
488
489 set JetPTMin 30.0
490}
491
492##################
493# Jet Energy Scale
494##################
495
496module EnergyScale JetEnergyScale {
497 set InputArray FastJetFinder/jets
498 set OutputArray jets
499
500 # scale formula for jets
501 set ScaleFormula {1.00}
502}
503
504
505########################
506# Jet Flavor Association
507########################
508
509module JetFlavorAssociation JetFlavorAssociation {
510
511 set PartonInputArray Delphes/partons
512 set ParticleInputArray Delphes/allParticles
513 set ParticleLHEFInputArray Delphes/allParticlesLHEF
514 set JetInputArray JetEnergyScale/jets
515
516 set DeltaR 0.5
517 set PartonPTMin 5.0
518 set PartonEtaMax 6.0
519
520}
521
522###################
523# Photon efficiency
524###################
525
526module Efficiency PhotonEfficiency {
527 set InputArray ECal/eflowPhotons
528 set OutputArray photons
529
530 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
531
532 set EfficiencyFormula {
533 (pt <= 10.0) * (0.00) + \
534 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
535 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) + \
536 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
537 (abs(eta) > 6.0) * (0.00)}
538
539}
540
541##################
542# Photon isolation
543##################
544
545# TBC: check values for iso cuts
546
547module Isolation PhotonIsolation {
548 set CandidateInputArray PhotonEfficiency/photons
549 set IsolationInputArray EFlowMerger/eflow
550
551 set OutputArray photons
552
553 set DeltaRMax 0.3
554
555 set PTMin 0.5
556
557 set PTRatioMax 0.1
558}
559
560
561####################
562# Electron isolation
563####################
564
565# TBC: check values for iso cuts
566
567module Isolation ElectronIsolation {
568 set CandidateInputArray ElectronFilter/electrons
569 set IsolationInputArray EFlowMerger/eflow
570
571 set OutputArray electrons
572
573 set DeltaRMax 0.3
574
575 set PTMin 0.5
576
577 set PTRatioMax 0.1
578}
579
580
581################
582# Muon isolation
583################
584
585# TBC: check values for iso cuts
586
587module Isolation MuonIsolation {
588 set CandidateInputArray MuonMomentumSmearing/muons
589 set IsolationInputArray EFlowMerger/eflow
590
591 set OutputArray muons
592
593 set DeltaRMax 0.3
594
595 set PTMin 0.5
596
597 set PTRatioMax 0.2
598}
599
600
601###########
602# b-tagging
603###########
604
605module BTagging BTagging {
606 set JetInputArray JetEnergyScale/jets
607
608 set BitNumber 0
609
610 add EfficiencyFormula {0} {
611
612 (pt <= 10.0) * (0.00) +
613 (abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.001) + \
614 (abs(eta) < 4.0) * (pt > 500.0) * (0.001*exp(0.5 - pt*0.001)) + \
615 (abs(eta) > 4.0) * (0.00)}
616
617 add EfficiencyFormula {4} {
618
619 (pt <= 10.0) * (0.00) +
620 (abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.04) + \
621 (abs(eta) < 4.0) * (pt > 500.0) * (0.04*exp(0.5 - pt*0.001)) + \
622 (abs(eta) > 4.0) * (0.00)}
623
624 add EfficiencyFormula {5} {
625
626 (pt <= 10.0) * (0.00) +
627 (abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.80) + \
628 (abs(eta) < 4.0) * (pt > 500.0) * (0.85*exp(0.5 - pt*0.001)) + \
629 (abs(eta) > 4.0) * (0.00)}
630
631}
632
633###########
634# c-tagging
635###########
636
637module BTagging CTagging {
638 set JetInputArray JetEnergyScale/jets
639
640 set BitNumber 1
641
642 add EfficiencyFormula {0} {
643
644 (pt <= 10.0) * (0.00) +
645 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
646 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
647
648 add EfficiencyFormula {4} {
649
650 (pt <= 10.0) * (0.00) +
651 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
652 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
653
654 add EfficiencyFormula {5} {
655
656 (pt <= 10.0) * (0.00) +
657 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
658 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
659
660}
661
662
663#############
664# tau-tagging
665#############
666
667
668module TauTagging TauTagging {
669 set ParticleInputArray Delphes/allParticles
670 set PartonInputArray Delphes/partons
671 set JetInputArray JetEnergyScale/jets
672
673 set DeltaR 0.5
674
675 set TauPTMin 1.0
676
677 set TauEtaMax 4.0
678
679 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
680 add EfficiencyFormula {0} {
681
682 (pt <= 10.0) * (0.00) +
683 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
684 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
685
686 add EfficiencyFormula {11} {
687
688 (pt <= 10.0) * (0.00) +
689 (abs(eta) < 4.0) * (pt > 10.0) * (0.005) + \
690 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
691
692 add EfficiencyFormula {15} {
693
694 (pt <= 10.0) * (0.00) +
695 (abs(eta) < 4.0) * (pt > 10.0) * (0.6) + \
696 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
697
698}
699
700#####################################################
701# Find uniquely identified photons/electrons/tau/jets
702#####################################################
703
704module UniqueObjectFinder UniqueObjectFinder {
705# earlier arrays take precedence over later ones
706# add InputArray InputArray OutputArray
707 add InputArray PhotonIsolation/photons photons
708 add InputArray ElectronIsolation/electrons electrons
709 add InputArray MuonIsolation/muons muons
710 add InputArray JetEnergyScale/jets jets
711
712}
713
714
715##################
716# ROOT tree writer
717##################
718
719module TreeWriter TreeWriter {
720# add Branch InputArray BranchName BranchClass
721 add Branch Delphes/allParticles Particle GenParticle
722
723 add Branch GenJetFinder/jets GenJet Jet
724 add Branch GenMissingET/momentum GenMissingET MissingET
725
726 add Branch TrackMerger/tracks Track Track
727 add Branch Calorimeter/towers Tower Tower
728
729 add Branch HCal/eflowTracks EFlowTrack Track
730 add Branch ECal/eflowPhotons EFlowPhoton Tower
731 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
732
733 add Branch UniqueObjectFinder/photons Photon Photon
734 add Branch UniqueObjectFinder/electrons Electron Electron
735 add Branch UniqueObjectFinder/muons Muon Muon
736 add Branch UniqueObjectFinder/jets Jet Jet
737
738 add Branch MissingET/momentum MissingET MissingET
739 add Branch ScalarHT/energy ScalarHT ScalarHT
740}
741
Note: See TracBrowser for help on using the repository browser.