Fork me on GitHub

Ticket #833: myFCChh_DelphesCard_WithDipole_v01.tcl

File myFCChh_DelphesCard_WithDipole_v01.tcl, 29.7 KB (added by Michele Selvaggi, 9 years ago)
Line 
1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Heather Gray (CERN)
5# Filip Moortgat (CERN)
6#
7# Released on: Sep 15th, 2015
8#
9# Version: v00
10#
11# Configuration: FCC-hh generic detector with dipole region, no pile-up
12#
13
14#######################################
15# Order of execution of various modules
16#######################################
17
18set ExecutionPath {
19ParticlePropagator
20
21ChargedHadronTrackingEfficiency
22
23ElectronEfficiency
24MuonEfficiency
25
26ChargedHadronMomentumSmearing
27ElectronEnergySmearing
28MuonMomentumSmearing
29
30TrackMerger
31ImpactParameterSmearing
32
33Ecal
34Hcal
35
36TowerMerger
37EFlowMerger
38
39PhotonEfficiency
40
41MuonIsolation
42ElectronIsolation
43PhotonIsolation
44
45MissingET
46
47GenJetFinder
48FastJetFinder
49
50JetEnergyScale
51
52JetFlavorAssociation
53
54BTagging
55CTagging
56TauTagging
57
58UniqueObjectFinder
59
60ScalarHT
61
62TreeWriter
63}
64
65#################################
66# Propagate particles in cylinder
67#################################
68
69module ParticlePropagator ParticlePropagator {
70set InputArray Delphes/stableParticles
71
72set OutputArray stableParticles
73set ChargedHadronOutputArray chargedHadrons
74set ElectronOutputArray electrons
75set MuonOutputArray muons
76
77# radius of the magnetic field coverage, in m
78set Radius 6.00
79# half-length of the magnetic field coverage, in m
80set HalfLength 11.5
81
82# magnetic field
83set Bz 6.0
84}
85
86####################################
87# Charged hadron tracking efficiency
88####################################
89
90module Efficiency ChargedHadronTrackingEfficiency {
91set InputArray ParticlePropagator/chargedHadrons
92set OutputArray chargedHadrons
93
94# add EfficiencyFormula {efficiency formula as a function of eta and pt}
95
96# tracking efficiency formula for charged hadrons
97set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
98(abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
99(abs(eta) <= 1.5) * (pt > 1) * (0.95) + \
100(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
101(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
102(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
103(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
104(abs(eta) > 6.0) * (0.00)}
105}
106
107##############################
108# Electron tracking efficiency
109##############################
110
111module Efficiency ElectronEfficiency {
112set InputArray ParticlePropagator/electrons
113set OutputArray electrons
114
115# set EfficiencyFormula {efficiency formula as a function of eta and pt}
116
117# tracking efficiency formula for electrons
118# tracking efficiency formula for electrons
119set EfficiencyFormula { (pt <= 10) * (0.00) + \
120(abs(eta) <= 1.5) * (pt > 10 && pt <= 50) * (0.80) + \
121(abs(eta) <= 1.5) * (pt > 50) * (0.90) + \
122(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 10 && pt <= 50) * (0.80) + \
123(abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 50) * (0.90) + \
124(abs(eta) > 4 && abs(eta) <= 6) * (pt > 10 && pt <= 50) * (0.70) + \
125(abs(eta) > 4 && abs(eta) <= 6) * (pt > 50) * (0.80) + \
126(abs(eta) > 6) * (0.00)}
127}
128
129##########################
130# Muon tracking efficiency
131##########################
132
133module Efficiency MuonEfficiency {
134set InputArray ParticlePropagator/muons
135set OutputArray muons
136
137# set EfficiencyFormula {efficiency formula as a function of eta and pt}
138
139# tracking efficiency formula for muons
140set EfficiencyFormula { (pt <= 10) * (0.00) + \
141(abs(eta) <= 1.5) * (pt > 10 ) * (0.99) + \
142(abs(eta) <= 4 & abs(eta) > 1.5) * (pt > 10 ) * (0.99) + \
143(abs(eta) > 4) * (0.00)}
144}
145
146########################################
147# Momentum resolution for charged tracks
148########################################
149
150module MomentumSmearing ChargedHadronMomentumSmearing {
151set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
152set OutputArray chargedHadrons
153
154# set ResolutionFormula {resolution formula as a function of eta and pt}
155
156# resolution formula for charged hadrons
157#
158# Automatically generated tracker resolution formula for layout: FCChh_Option2 (by Z. Drasal, CERN)
159#
160# Central and Dipole
161#
162set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00194513) + \
163 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001945 + (pt-5.000000)* 0.000003) + \
164 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.001962 + (pt-10.000000)* 0.000006) + \
165 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.002547 + (pt-100.000000)* 0.000014) + \
166 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.014708 + (pt-1000.000000)* 0.000015) + \
167 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10000.0000) * (0.145553*pt/10000.000000) + \
168 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00247067) + \
169 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002471 + (pt-5.000000)* -0.000002) + \
170 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002462 + (pt-10.000000)* 0.000008) + \
171 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.003170 + (pt-100.000000)* 0.000015) + \
172 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.016485 + (pt-1000.000000)* 0.000016) + \
173 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10000.0000) * (0.162501*pt/10000.000000) + \
174 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00242524) + \
175 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002425 + (pt-5.000000)* 0.000013) + \
176 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002490 + (pt-10.000000)* 0.000042) + \
177 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.006273 + (pt-100.000000)* 0.000049) + \
178 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.050079 + (pt-1000.000000)* 0.000050) + \
179 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10000.0000) * (0.498773*pt/10000.000000) + \
180 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00168511) + \
181 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001685 + (pt-5.000000)* 0.000115) + \
182 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002262 + (pt-10.000000)* 0.000167) + \
183 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.017274 + (pt-100.000000)* 0.000172) + \
184 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.172095 + (pt-1000.000000)* 0.000172) + \
185 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10000.0000) * (1.720887*pt/10000.000000) + \
186 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00277124) + \
187 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002771 + (pt-5.000000)* 0.000430) + \
188 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.004921 + (pt-10.000000)* 0.000467) + \
189 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.046933 + (pt-100.000000)* 0.000469) + \
190 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.469096 + (pt-1000.000000)* 0.000469) + \
191 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10000.0000) * (4.690935*pt/10000.000000) + \
192 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00653887) + \
193 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.006539 + (pt-5.000000)* 0.001256) + \
194 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.012821 + (pt-10.000000)* 0.001272) + \
195 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.127345 + (pt-100.000000)* 0.001273) + \
196 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (1.273365 + (pt-1000.000000)* 0.001273) + \
197 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10000.0000) * (12.733637*pt/10000.000000)
198}
199}
200#################################
201# Energy resolution for electrons
202#################################
203
204module EnergySmearing ElectronEnergySmearing {
205set InputArray ElectronEfficiency/electrons
206set OutputArray electrons
207
208# set ResolutionFormula {resolution formula as a function of eta and energy}
209
210# check compared to compare to tracket resolution for low
211# resolution formula for electrons
212set ResolutionFormula {
213(abs(eta) <= 1.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
214(abs(eta) > 1.5 && abs(eta) <= 2.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
215(abs(eta) > 2.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
216(abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
217}
218
219
220#################################
221# Energy resolution for photons
222#################################
223#not used for the moment (has been removed from the ExecutionPath)
224
225module EnergySmearing PhotonEnergySmearing {
226set InputArray PhotonEfficiency/photons
227set OutputArray photons
228
229# set ResolutionFormula {resolution formula as a function of eta and energy}
230
231# check compared to compare to track resolution for low
232# resolution formula for photons
233set ResolutionFormula {
234(abs(eta) <= 1.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
235(abs(eta) > 1.5 && abs(eta) <= 2.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
236(abs(eta) > 2.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
237(abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
238}
239
240
241###############################
242# Momentum resolution for muons
243###############################
244
245module MomentumSmearing MuonMomentumSmearing {
246set InputArray MuonEfficiency/muons
247set OutputArray muons
248
249# set ResolutionFormula {resolution formula as a function of eta and pt}
250
251# resolution formula for muons
252# Automatically generated tracker resolution formula for layout: FCChh_Option2 (by Z. Drasal, CERN)
253# Central and Dipole
254#
255set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00194513) + \
256 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001945 + (pt-5.000000)* 0.000003) + \
257 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.001962 + (pt-10.000000)* 0.000006) + \
258 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.002547 + (pt-100.000000)* 0.000014) + \
259 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.014708 + (pt-1000.000000)* 0.000015) + \
260 (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10000.0000) * (0.145553*pt/10000.000000) + \
261 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00247067) + \
262 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002471 + (pt-5.000000)* -0.000002) + \
263 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002462 + (pt-10.000000)* 0.000008) + \
264 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.003170 + (pt-100.000000)* 0.000015) + \
265 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.016485 + (pt-1000.000000)* 0.000016) + \
266 (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10000.0000) * (0.162501*pt/10000.000000) + \
267 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00242524) + \
268 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002425 + (pt-5.000000)* 0.000013) + \
269 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002490 + (pt-10.000000)* 0.000042) + \
270 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.006273 + (pt-100.000000)* 0.000049) + \
271 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.050079 + (pt-1000.000000)* 0.000050) + \
272 (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10000.0000) * (0.498773*pt/10000.000000) + \
273 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00168511) + \
274 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001685 + (pt-5.000000)* 0.000115) + \
275 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002262 + (pt-10.000000)* 0.000167) + \
276 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.017274 + (pt-100.000000)* 0.000172) + \
277 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.172095 + (pt-1000.000000)* 0.000172) + \
278 (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10000.0000) * (1.720887*pt/10000.000000) + \
279 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00277124) + \
280 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002771 + (pt-5.000000)* 0.000430) + \
281 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.004921 + (pt-10.000000)* 0.000467) + \
282 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.046933 + (pt-100.000000)* 0.000469) + \
283 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.469096 + (pt-1000.000000)* 0.000469) + \
284 (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10000.0000) * (4.690935*pt/10000.000000) + \
285 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00653887) + \
286 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.006539 + (pt-5.000000)* 0.001256) + \
287 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.012821 + (pt-10.000000)* 0.001272) + \
288 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.127345 + (pt-100.000000)* 0.001273) + \
289 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (1.273365 + (pt-1000.000000)* 0.001273) + \
290 (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10000.0000) * (12.733637*pt/10000.000000)
291}
292}
293####################
294# Electron isolation
295####################
296
297module Isolation ElectronIsolation {
298set CandidateInputArray ElectronEnergySmearing/electrons
299set IsolationInputArray EFlowMerger/eflow
300
301set OutputArray electrons
302
303set DeltaRMax 0.3
304
305set PTMin 0.5
306
307set PTRatioMax 0.1
308}
309
310################
311# Muon isolation
312################
313
314module Isolation MuonIsolation {
315set CandidateInputArray MuonMomentumSmearing/muons
316set IsolationInputArray EFlowMerger/eflow
317
318set OutputArray muons
319
320set DeltaRMax 0.3
321
322set PTMin 0.5
323
324set PTRatioMax 0.1
325}
326
327###################
328# Photon efficiency
329###################
330
331module Efficiency PhotonEfficiency {
332set InputArray Ecal/eflowPhotons
333set OutputArray photons
334
335# set EfficiencyFormula {efficiency formula as a function of eta and pt}
336# efficiency formula for photons
337set EfficiencyFormula {
338(pt <= 10.0) * (0.00) + \
339(abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
340(abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) + \
341(abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
342(abs(eta) > 6.0) * (0.00)}
343}
344
345##################
346# Photon isolation
347##################
348
349module Isolation PhotonIsolation {
350set CandidateInputArray PhotonEfficiency/photons
351set IsolationInputArray EFlowMerger/eflow
352
353set OutputArray photons
354
355set DeltaRMax 0.3
356
357set PTMin 0.5
358
359set PTRatioMax 0.1
360}
361
362##############
363# Track merger
364##############
365
366module Merger TrackMerger {
367# add InputArray InputArray
368add InputArray ChargedHadronMomentumSmearing/chargedHadrons
369add InputArray ElectronEnergySmearing/electrons
370add InputArray MuonMomentumSmearing/muons
371set OutputArray tracks
372}
373
374################################
375# Track impact parameter smearing
376################################
377
378#not used at the moment
379
380module ImpactParameterSmearing ImpactParameterSmearing {
381set InputArray TrackMerger/tracks
382set OutputArray tracks
383
384
385# absolute impact parameter smearing formula (in mm) as a function of pt and eta
386set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \
387(pt > 5.0) * (0.005)}
388
389}
390
391#############
392# ECAL
393#############
394
395module SimpleCalorimeter Ecal {
396 set ParticleInputArray ParticlePropagator/stableParticles
397 set TrackInputArray ImpactParameterSmearing/tracks
398
399 set TowerOutputArray ecalTowers
400 set EFlowTowerOutputArray eflowPhotons
401
402 set pi [expr {acos(-1)}]
403
404 # lists of the edges of each tower in eta and phi
405 # each list starts with the lower edge of the first tower
406 # the list ends with the higher edged of the last tower
407 # 0.0125 x 0.0125 for eta between 0 and 2.5
408 set PhiBins {}
409 for {set i -256} {$i <= 256} {incr i} {
410 add PhiBins [expr {$i * $pi/256.0}]
411 }
412 foreach eta {-2.5 -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} {
413 add EtaPhiBins $eta $PhiBins
414 }
415 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
416 set PhiBins {}
417 for {set i -128} {$i <= 128} {incr i} {
418 add PhiBins [expr {$i * $pi/128.0}]
419 }
420 foreach eta {-4.0 -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.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} {
421 add EtaPhiBins $eta $PhiBins
422 }
423
424 # 0.05 x 0.05 for eta between 4.0 and 6.0
425 set PhiBins {}
426 for {set i -64} {$i <= 64} {incr i} {
427 add PhiBins [expr {$i * $pi/64.0}]
428 }
429 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.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} {
430 add EtaPhiBins $eta $PhiBins
431 }
432 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
433
434 add EnergyFraction {0} {0.4}
435 # energy fractions for e, gamma and pi0
436 add EnergyFraction {11} {1.0}
437 add EnergyFraction {22} {1.0}
438 add EnergyFraction {111} {1.0}
439 # energy fractions for muon, neutrinos and neutralinos
440 add EnergyFraction {12} {0.0}
441 add EnergyFraction {13} {0.0}
442 add EnergyFraction {14} {0.0}
443 add EnergyFraction {16} {0.0}
444 add EnergyFraction {1000022} {0.0}
445 add EnergyFraction {1000023} {0.0}
446 add EnergyFraction {1000025} {0.0}
447 add EnergyFraction {1000035} {0.0}
448 add EnergyFraction {1000045} {0.0}
449 # energy fractions for K0short and Lambda
450 # add EnergyFraction {310} {0.3}
451 # add EnergyFraction {3122} {0.3}
452
453 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
454 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
455 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
456
457
458}
459
460
461#############
462# HCAL
463#############
464
465module SimpleCalorimeter Hcal {
466 set ParticleInputArray ParticlePropagator/stableParticles
467 set TrackInputArray ImpactParameterSmearing/tracks
468
469 set TowerOutputArray hcalTowers
470 set EFlowTowerOutputArray eflowNeutralHadrons
471
472 set pi [expr {acos(-1)}]
473
474 # lists of the edges of each tower in eta and phi
475 # each list starts with the lower edge of the first tower
476 # the list ends with the higher edged of the last tower
477
478 # 0.05 x 0.05 between eta 0 and 2.5
479 set PhiBins {}
480 for {set i -128} {$i <= 128} {incr i} {
481 add PhiBins [expr {$i * $pi/128.0}]
482 }
483 foreach eta {-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.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} {
484 add EtaPhiBins $eta $PhiBins
485 }
486
487 # 0.1 x 0.1 between eta 2.5 and 4.0
488 set PhiBins {}
489 for {set i -64} {$i <= 64} {incr i} {
490 add PhiBins [expr {$i * $pi/64.0}]
491 }
492 foreach eta {-4.0 -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.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} {
493 add EtaPhiBins $eta $PhiBins
494 }
495
496 # 0.2 x 0.2 between eta 4.0 and 6.0
497 set PhiBins {}
498 for {set i -32} {$i <= 32} {incr i} {
499 add PhiBins [expr {$i * $pi/32.0}]
500 }
501 foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
502 add EtaPhiBins $eta $PhiBins
503 }
504
505 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
506 add EnergyFraction {0} {0.6}
507 # energy fractions for e, gamma and pi0
508 add EnergyFraction {11} {0.0}
509 add EnergyFraction {22} {0.0}
510 add EnergyFraction {111} {0.0}
511 # energy fractions for muon, neutrinos and neutralinos
512 add EnergyFraction {12} {0.0}
513 add EnergyFraction {13} {0.0}
514 add EnergyFraction {14} {0.0}
515 add EnergyFraction {16} {0.0}
516 add EnergyFraction {1000022} {0.0}
517 add EnergyFraction {1000023} {0.0}
518 add EnergyFraction {1000025} {0.0}
519 add EnergyFraction {1000035} {0.0}
520 add EnergyFraction {1000045} {0.0}
521 # energy fractions for K0short and Lambda
522 # add EnergyFraction {310} {0.7}
523 # add EnergyFraction {3122} {0.7}
524
525 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
526 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
527 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
528}
529
530
531
532
533####################
534# Tower Merger (in case not using e-flow algorithm)
535####################
536
537module Merger TowerMerger {
538# add InputArray InputArray
539add InputArray Ecal/ecalTowers
540add InputArray Hcal/hcalTowers
541set OutputArray towers
542}
543
544####################
545# Energy flow merger
546####################
547
548module Merger EFlowMerger {
549# add InputArray InputArray
550add InputArray ImpactParameterSmearing/tracks
551add InputArray Ecal/eflowPhotons
552add InputArray Hcal/eflowNeutralHadrons
553set OutputArray eflow
554}
555
556
557###################
558# Missing ET merger
559###################
560
561module Merger MissingET {
562# add InputArray InputArray
563add InputArray EFlowMerger/eflow
564set MomentumOutputArray momentum
565}
566
567
568##################
569# Scalar HT merger
570##################
571
572module Merger ScalarHT {
573# add InputArray InputArray
574add InputArray EFlowMerger/eflow
575set EnergyOutputArray energy
576}
577########################
578# Jet Flavor Association
579########################
580
581module JetFlavorAssociation JetFlavorAssociation {
582
583 set PartonInputArray Delphes/partons
584 set ParticleInputArray Delphes/allParticles
585 set ParticleLHEFInputArray Delphes/allParticlesLHEF
586 set JetInputArray JetEnergyScale/jets
587
588 set DeltaR 0.5
589 set PartonPTMin 1.0
590 set PartonEtaMax 6.0
591
592}
593#####################
594# MC truth jet finder
595#####################
596
597module FastJetFinder GenJetFinder {
598set InputArray Delphes/stableParticles
599
600set OutputArray jets
601
602# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 7 antikt wta algorithm
603set JetAlgorithm 7
604set ParameterR 0.4
605
606set JetPTMin 5.0
607}
608
609############
610# Jet finder
611############
612
613module FastJetFinder FastJetFinder {
614# set InputArray Calorimeter/towers
615set InputArray EFlowMerger/eflow
616
617set OutputArray jets
618
619# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 antikt wta algorithm
620set JetAlgorithm 7
621set ParameterR 0.4
622
623set JetPTMin 30.0
624}
625
626##################
627# Jet Energy Scale
628##################
629# not used at the moment
630
631module EnergyScale JetEnergyScale {
632set InputArray FastJetFinder/jets
633set OutputArray jets
634
635# scale formula for jets
636set ScaleFormula {1.00}
637}
638
639##########################
640# tau-tagging
641##########################
642
643
644module TauTagging TauTagging {
645set ParticleInputArray Delphes/allParticles
646set PartonInputArray Delphes/partons
647set JetInputArray JetEnergyScale/jets
648set DeltaR 0.4
649set TauPTMin 15.0
650set TauEtaMax 6
651# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
652add EfficiencyFormula {0} {0.01}
653add EfficiencyFormula {11} {0.005}
654add EfficiencyFormula {15} {0.6}
655}
656
657
658##########################
659# b-tagging
660##########################
661
662module BTagging BTagging { # need to check that this module exists or write it
663set JetInputArray JetEnergyScale/jets
664
665# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
666
667# default efficiency formula (misidentification rate)
668add EfficiencyFormula {0} {0.001}
669add EfficiencyFormula {4} {0.04}
670add EfficiencyFormula {5} {0.75}
671}
672
673##########################
674# c-tagging
675##########################
676
677
678module BTagging CTagging {
679set ParticleInputArray Delphes/allParticles
680set PartonInputArray Delphes/partons
681set JetInputArray JetEnergyScale/jets
682
683set DeltaR 0.4
684set BPTMin 30.0
685set BEtaMax 6.0
686
687# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
688add EfficiencyFormula {0} {0.01}
689add EfficiencyFormula {5} {0.10}
690add EfficiencyFormula {4} {0.25}
691}
692
693#####################################################
694# Find uniquely identified photons/electrons/tau/jets
695#####################################################
696
697module UniqueObjectFinder UniqueObjectFinder {
698# earlier arrays take precedence over later ones
699# add InputArray InputArray OutputArray
700add InputArray MuonIsolation/muons muons
701add InputArray ElectronIsolation/electrons electrons
702add InputArray PhotonIsolation/photons photons
703add InputArray JetEnergyScale/jets jets
704}
705
706##################
707# ROOT tree writer
708##################
709
710module TreeWriter TreeWriter {
711# add Branch InputArray BranchName BranchClass
712add Branch Delphes/allParticles Particle GenParticle
713add Branch GenJetFinder/jets GenJet Jet
714
715add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
716add Branch Hcal/eflowNeutralHadrons NeutralHadron Tower
717add Branch UniqueObjectFinder/photons Photon Photon
718
719add Branch UniqueObjectFinder/electrons Electron Electron
720add Branch UniqueObjectFinder/muons Muon Muon
721add Branch UniqueObjectFinder/jets Jet Jet
722#add Branch JetEnergyScale/jets Jet Jet
723add Branch MissingET/momentum MissingET MissingET
724add Branch ScalarHT/energy ScalarHT ScalarHT
725}
726
727
728