Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ 4c9835d

Last change on this file since 4c9835d was b8a6aa3, checked in by Michele Selvaggi <michele@…>, 4 years ago

removed duplication of IDEA cards

  • Property mode set to 100644
File size: 26.9 KB
Line 
1#################################################################### l
2# FCC-ee IDEA detector model
3#
4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
5# email: efontane@bo.infn.it,
6# lorenzo.pezzotti01@universitadipavia.it,
7# m.antonello@uninsubria.it,
8# michele.selvaggi@cern.ch
9#####################################################################
10
11set B 2.0
12
13#######################################
14# Order of execution of various modules
15#######################################
16
17set ExecutionPath {
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
24 TrackMergerPre
25 TrackSmearing
26
27 TrackMerger
28 Calorimeter
29 EFlowMerger
30
31 PhotonEfficiency
32 PhotonIsolation
33
34 MuonFilter
35
36 ElectronFilter
37 ElectronEfficiency
38 ElectronIsolation
39
40 MuonEfficiency
41 MuonIsolation
42
43 MissingET
44
45 NeutrinoFilter
46 GenJetFinder
47 GenMissingET
48
49 FastJetFinder
50
51 JetEnergyScale
52
53 JetFlavorAssociation
54
55 BTagging
56 TauTagging
57
58 UniqueObjectFinder
59
60 ScalarHT
61 TreeWriter
62}
63
64#################################
65# Propagate particles in cylinder
66#################################
67
68module ParticlePropagator ParticlePropagator {
69 set InputArray Delphes/stableParticles
70
71 set OutputArray stableParticles
72 set ChargedHadronOutputArray chargedHadrons
73 set ElectronOutputArray electrons
74 set MuonOutputArray muons
75
76 # inner radius of the solenoid, in m
77 set Radius 2.25
78
79 # half-length: z of the solenoid, in m
80 set HalfLength 2.5
81
82 # magnetic field, in T
83 set Bz $B
84}
85
86####################################
87# Charged hadron tracking efficiency
88####################################
89
90module Efficiency ChargedHadronTrackingEfficiency {
91 set InputArray ParticlePropagator/chargedHadrons
92 set OutputArray chargedHadrons
93 # We use only one efficiency, we set only 0 effincency out of eta bounds:
94
95 set EfficiencyFormula {
96 (abs(eta) > 3.0) * (0.000) +
97 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
98 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
99 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
100 }
101}
102
103# (pt <= 0.1) * (0.00) +
104# (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
105# (abs(eta) > 3) * (0.00)
106
107
108
109##############################
110# Electron tracking efficiency
111##############################
112
113module Efficiency ElectronTrackingEfficiency {
114 set InputArray ParticlePropagator/electrons
115 set OutputArray electrons
116
117
118 # Current full simulation with CLICdet provides for electrons:
119 set EfficiencyFormula {
120 (abs(eta) > 3.0) * (0.000) +
121 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
122 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
123 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
124 }
125}
126
127
128##########################
129# Muon tracking efficiency
130##########################
131
132module Efficiency MuonTrackingEfficiency {
133 set InputArray ParticlePropagator/muons
134 set OutputArray muons
135
136 # Current full simulation with CLICdet provides for muons:
137 set EfficiencyFormula {
138 (abs(eta) > 3.0) * (0.000) +
139 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
140 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
141 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
142 }
143}
144
145##############
146# Track merger
147##############
148
149module Merger TrackMergerPre {
150# add InputArray InputArray
151 add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
152 add InputArray ElectronTrackingEfficiency/electrons
153 add InputArray MuonTrackingEfficiency/muons
154 set OutputArray tracks
155}
156
157
158########################################
159# Smearing for charged tracks
160########################################
161
162module TrackCovariance TrackSmearing {
163 set InputArray TrackMergerPre/tracks
164 set OutputArray tracks
165
166
167 set InputArray TrackMergerPre/tracks
168 set OutputArray tracks
169
170 set Bz 2.0
171
172 ## minimum number of hits to accept a track
173 set NMinHits 6
174
175 ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
176 set DetectorGeometry {
177
178
179 # Layer type 1 = R (barrel) or 2 = z (forward/backward)
180 # Layer label
181 # Minimum dimension z for barrel or R for forward
182 # Maximum dimension z for barrel or R for forward
183 # R/z location of layer
184 # Thickness (meters)
185 # Radiation length (meters)
186 # Number of measurements in layers (1D or 2D)
187 # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Upper side
188 # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Lower side
189 # Resolution Upper side (meters) - 0 = no measurement
190 # Resolution Lower side (meters) - 0 = no measurement
191 # measurement flag = T, scattering only = F
192
193 # barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
194
195 # barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
196
197 1 PIPE -100 100 0.015 0.001655 0.2805 0 0 0 0 0 0
198 1 VTXLOW -0.12 0.12 0.017 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
199 1 VTXLOW -0.16 0.16 0.023 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
200 1 VTXLOW -0.16 0.16 0.031 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
201 1 VTXHIGH -1 1 0.32 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
202 1 VTXHIGH -1.05 1.05 0.34 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
203
204 # endcap name rmin rmax z w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
205
206 2 VTXDSK 0.141 0.3 -0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
207 2 VTXDSK 0.138 0.3 -0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
208 2 VTXDSK 0.065 0.3 -0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
209 2 VTXDSK 0.062 0.3 -0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
210 2 VTXDSK 0.062 0.3 0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
211 2 VTXDSK 0.065 0.3 0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
212 2 VTXDSK 0.138 0.3 0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
213 2 VTXDSK 0.141 0.3 0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
214
215 1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
216 1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
217 1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
218 1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
219 1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
220 1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
221 1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
222 1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
223 1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
224 1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
225 1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
226 1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
227 1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
228 1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
229 1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
230 1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
231 1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
232 1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
233 1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
234 1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
235 1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
236 1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
237 1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
238 1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
239 1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
240 1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
241 1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
242 1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
243 1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
244 1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
245 1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
246 1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
247 1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
248 1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
249 1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
250 1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
251 1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
252 1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
253 1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
254 1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
255 1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
256 1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
257 1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
258 1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
259 1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
260 1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
261 1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
262 1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
263 1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
264 1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
265 1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
266 1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
267 1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
268 1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
269 1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
270 1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
271 1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
272 1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
273 1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
274 1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
275 1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
276 1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
277 1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
278 1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
279 1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
280 1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
281 1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
282 1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
283 1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
284 1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
285 1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
286 1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
287 1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
288 1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
289 1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
290 1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
291 1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
292 1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
293 1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
294 1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
295 1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
296 1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
297 1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
298 1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
299 1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
300 1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
301 1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
302 1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
303 1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
304 1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
305 1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
306 1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
307 1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
308 1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
309 1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
310 1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
311 1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
312 1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
313 1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
314 1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
315 1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
316 1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
317 1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
318 1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
319 1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
320 1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
321 1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
322 1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
323 1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
324 1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
325 1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
326 1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
327 1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
328 1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
329 1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
330 1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
331 1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
332 1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
333 2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
334 2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
335 2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
336 2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
337 2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
338 2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
339 2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
340 2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
341 2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
342 2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
343 }
344
345 set Bz $B
346}
347
348##############
349# Track merger
350##############
351
352module Merger TrackMerger {
353# add InputArray InputArray
354 add InputArray TrackSmearing/tracks
355 set OutputArray tracks
356}
357
358
359#############
360# Calorimeter
361#############
362module DualReadoutCalorimeter Calorimeter {
363 set ParticleInputArray ParticlePropagator/stableParticles
364 set TrackInputArray TrackMerger/tracks
365
366 set TowerOutputArray towers
367 set PhotonOutputArray photons
368
369 set EFlowTrackOutputArray eflowTracks
370 set EFlowPhotonOutputArray eflowPhotons
371 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
372
373 set ECalEnergyMin 0.5
374 set HCalEnergyMin 0.5
375 set EnergyMin 0.5
376 set ECalEnergySignificanceMin 1.0
377 set HCalEnergySignificanceMin 1.0
378 set EnergySignificanceMin 1.0
379
380 set SmearTowerCenter true
381 set pi [expr {acos(-1)}]
382
383 # Lists of the edges of each tower in eta and phi;
384 # each list starts with the lower edge of the first tower;
385 # the list ends with the higher edged of the last tower.
386 # Barrel: deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
387 # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
388 # Cell size: about 6 cm x 6 cm
389
390 #barrel:
391 set PhiBins {}
392 for {set i -120} {$i <= 120} {incr i} {
393 add PhiBins [expr {$i * $pi/120}]
394 }
395 #deta=0.02 units for |eta| <= 0.88
396 for {set i -44} {$i < 45} {incr i} {
397 set eta [expr {$i * 0.02}]
398 add EtaPhiBins $eta $PhiBins
399 }
400
401 #endcaps:
402 set PhiBins {}
403 for {set i -120} {$i <= 120} {incr i} {
404 add PhiBins [expr {$i* $pi/120}]
405 }
406 #deta=0.02 units for 0.88 < |eta| <= 3.0
407 #first, from -3.0 to -0.88
408 for {set i 1} {$i <=106} {incr i} {
409 set eta [expr {-3.00 + $i * 0.02}]
410 add EtaPhiBins $eta $PhiBins
411 }
412 #same for 0.88 to 3.0
413 for {set i 1} {$i <=106} {incr i} {
414 set eta [expr {0.88 + $i * 0.02}]
415 add EtaPhiBins $eta $PhiBins
416 }
417
418 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
419 add EnergyFraction {0} {0.0 1.0}
420 # energy fractions for e, gamma and pi0
421 add EnergyFraction {11} {1.0 0.0}
422 add EnergyFraction {22} {1.0 0.0}
423 add EnergyFraction {111} {1.0 0.0}
424 # energy fractions for muon, neutrinos and neutralinos
425 add EnergyFraction {12} {0.0 0.0}
426 add EnergyFraction {13} {0.0 0.0}
427 add EnergyFraction {14} {0.0 0.0}
428 add EnergyFraction {16} {0.0 0.0}
429 add EnergyFraction {1000022} {0.0 0.0}
430 add EnergyFraction {1000023} {0.0 0.0}
431 add EnergyFraction {1000025} {0.0 0.0}
432 add EnergyFraction {1000035} {0.0 0.0}
433 add EnergyFraction {1000045} {0.0 0.0}
434 # energy fractions for K0short and Lambda
435 add EnergyFraction {310} {0.3 0.7}
436 add EnergyFraction {3122} {0.3 0.7}
437
438
439 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
440 set ECalResolutionFormula {
441 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.11^2)+
442 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.11^2)
443 }
444
445 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
446 set HCalResolutionFormula {
447 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.30^2)+
448 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.30^2)
449 }
450}
451
452####################
453# Energy flow merger
454####################
455
456module Merger EFlowMerger {
457# add InputArray InputArray
458 add InputArray Calorimeter/eflowTracks
459 add InputArray Calorimeter/eflowPhotons
460 add InputArray Calorimeter/eflowNeutralHadrons
461 set OutputArray eflow
462}
463
464###################
465# Photon efficiency
466###################
467
468module Efficiency PhotonEfficiency {
469 set InputArray Calorimeter/eflowPhotons
470 set OutputArray photons
471
472 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
473 # efficiency formula for photons
474 set EfficiencyFormula {
475 (energy < 2.0) * (0.000)+
476 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
477 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
478 (abs(eta) > 3.0) * (0.000)
479 }
480}
481
482##################
483# Photon isolation
484##################
485
486module Isolation PhotonIsolation {
487 set CandidateInputArray PhotonEfficiency/photons
488 set IsolationInputArray EFlowMerger/eflow
489
490 set OutputArray photons
491
492 set DeltaRMax 0.5
493
494 set PTMin 0.5
495
496 set PTRatioMax 999.
497}
498
499#################
500# Electron filter
501#################
502
503module PdgCodeFilter ElectronFilter {
504 set InputArray Calorimeter/eflowTracks
505 set OutputArray electrons
506 set Invert true
507 add PdgCode {11}
508 add PdgCode {-11}
509}
510
511#################
512# Muon filter
513#################
514
515module PdgCodeFilter MuonFilter {
516 set InputArray Calorimeter/eflowTracks
517 set OutputArray muons
518 set Invert true
519 add PdgCode {13}
520 add PdgCode {-13}
521}
522
523
524#####################
525# Electron efficiency
526#####################
527
528module Efficiency ElectronEfficiency {
529 set InputArray ElectronFilter/electrons
530 set OutputArray electrons
531
532 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
533
534 # efficiency formula for electrons
535 set EfficiencyFormula {
536 (energy < 2.0) * (0.000)+
537 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
538 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
539 (abs(eta) > 3.0) * (0.000)
540 }
541}
542
543####################
544# Electron isolation
545####################
546
547module Isolation ElectronIsolation {
548 set CandidateInputArray ElectronEfficiency/electrons
549 set IsolationInputArray EFlowMerger/eflow
550
551 set OutputArray electrons
552
553 set DeltaRMax 0.5
554
555 set PTMin 0.5
556
557 set PTRatioMax 0.12
558}
559
560#################
561# Muon efficiency
562#################
563
564module Efficiency MuonEfficiency {
565 set InputArray MuonFilter/muons
566 set OutputArray muons
567
568 # set EfficiencyFormula {efficiency as a function of eta and pt}
569
570 # efficiency formula for muons
571 set EfficiencyFormula {
572 (energy < 2.0) * (0.000)+
573 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
574 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
575 (abs(eta) > 3.0) * (0.000)
576 }
577}
578
579################
580# Muon isolation
581################
582
583module Isolation MuonIsolation {
584 set CandidateInputArray MuonEfficiency/muons
585 set IsolationInputArray EFlowMerger/eflow
586
587 set OutputArray muons
588
589 set DeltaRMax 0.5
590
591 set PTMin 0.5
592
593 set PTRatioMax 0.25
594}
595
596###################
597# Missing ET merger
598###################
599
600module Merger MissingET {
601# add InputArray InputArray
602 add InputArray EFlowMerger/eflow
603 set MomentumOutputArray momentum
604}
605
606##################
607# Scalar HT merger
608##################
609
610module Merger ScalarHT {
611# add InputArray InputArray
612 add InputArray UniqueObjectFinder/jets
613 add InputArray UniqueObjectFinder/electrons
614 add InputArray UniqueObjectFinder/photons
615 add InputArray UniqueObjectFinder/muons
616 set EnergyOutputArray energy
617}
618
619#####################
620# Neutrino Filter
621#####################
622
623module PdgCodeFilter NeutrinoFilter {
624
625 set InputArray Delphes/stableParticles
626 set OutputArray filteredParticles
627
628 set PTMin 0.0
629
630 add PdgCode {12}
631 add PdgCode {14}
632 add PdgCode {16}
633 add PdgCode {-12}
634 add PdgCode {-14}
635 add PdgCode {-16}
636}
637
638
639#####################
640# MC truth jet finder
641#####################
642
643module FastJetFinder GenJetFinder {
644 set InputArray NeutrinoFilter/filteredParticles
645
646 set OutputArray jets
647
648 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
649 set JetAlgorithm 6
650 set ParameterR 0.4
651 set JetPTMin 1.0
652}
653
654
655#########################
656# Gen Missing ET merger
657########################
658
659module Merger GenMissingET {
660# add InputArray InputArray
661 add InputArray NeutrinoFilter/filteredParticles
662 set MomentumOutputArray momentum
663}
664
665############
666# Jet finder
667############
668
669module FastJetFinder FastJetFinder {
670# set InputArray Calorimeter/towers
671 set InputArray EFlowMerger/eflow
672
673 set OutputArray jets
674
675 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
676 set JetAlgorithm 6
677 set ParameterR 0.4
678 set JetPTMin 1.0
679}
680
681##################
682# Jet Energy Scale
683##################
684
685module EnergyScale JetEnergyScale {
686 set InputArray FastJetFinder/jets
687 set OutputArray jets
688
689 # scale formula for jets
690 set ScaleFormula {1.08}
691}
692
693########################
694# Jet Flavor Association
695########################
696
697module JetFlavorAssociation JetFlavorAssociation {
698
699 set PartonInputArray Delphes/partons
700 set ParticleInputArray Delphes/allParticles
701 set ParticleLHEFInputArray Delphes/allParticlesLHEF
702 set JetInputArray JetEnergyScale/jets
703
704 set DeltaR 0.5
705 set PartonPTMin 1.0
706 set PartonEtaMax 3.0
707}
708
709###########
710# b-tagging
711###########
712
713module BTagging BTagging {
714 set JetInputArray JetEnergyScale/jets
715
716 set BitNumber 0
717
718 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
719
720 # default efficiency formula (misidentification rate)
721 add EfficiencyFormula {0} {0.01}
722
723 # efficiency formula for c-jets (misidentification rate)
724 add EfficiencyFormula {4} {0.10}
725
726 # efficiency formula for b-jets
727 add EfficiencyFormula {5} {0.80}
728}
729
730#############
731# tau-tagging
732#############
733
734module TauTagging TauTagging {
735 set ParticleInputArray Delphes/allParticles
736 set PartonInputArray Delphes/partons
737 set JetInputArray JetEnergyScale/jets
738
739 set DeltaR 0.5
740 set TauPTMin 1.0
741 set TauEtaMax 3.0
742
743 # default efficiency formula (misidentification rate)
744 add EfficiencyFormula {0} {0.001}
745 # efficiency formula for tau-jets
746 add EfficiencyFormula {15} {0.6}
747}
748
749
750#####################################################
751# Find uniquely identified photons/electrons/tau/jets
752#####################################################
753
754module UniqueObjectFinder UniqueObjectFinder {
755# earlier arrays take precedence over later ones
756# add InputArray InputArray OutputArray
757 add InputArray PhotonIsolation/photons photons
758 add InputArray ElectronIsolation/electrons electrons
759 add InputArray MuonIsolation/muons muons
760 add InputArray JetEnergyScale/jets jets
761}
762
763
764
765##################
766# ROOT tree writer
767##################
768
769# Tracks, towers and eflow objects are not stored by default in the output.
770# If needed (for jet constituent or other studies), uncomment the relevant
771# "add Branch ..." lines.
772
773module TreeWriter TreeWriter {
774 # add Branch InputArray BranchName BranchClass
775
776 add Branch Delphes/allParticles Particle GenParticle
777
778 add Branch TrackMerger/tracks Track Track
779 add Branch Calorimeter/towers Tower Tower
780
781 add Branch Calorimeter/eflowTracks EFlowTrack Track
782 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
783 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
784
785 add Branch Calorimeter/photons CaloPhoton Photon
786 add Branch PhotonEfficiency/photons PhotonEff Photon
787 add Branch PhotonIsolation/photons PhotonIso Photon
788
789 add Branch GenJetFinder/jets GenJet Jet
790 add Branch GenMissingET/momentum GenMissingET MissingET
791
792 add Branch UniqueObjectFinder/jets Jet Jet
793 add Branch UniqueObjectFinder/electrons Electron Electron
794 add Branch UniqueObjectFinder/photons Photon Photon
795 add Branch UniqueObjectFinder/muons Muon Muon
796
797 add Branch JetEnergyScale/jets AntiKtJet Jet
798
799 add Branch MissingET/momentum MissingET MissingET
800 add Branch ScalarHT/energy ScalarHT ScalarHT
801
802 # add Info InfoName InfoValue
803 add Info Bz $B
804}
Note: See TracBrowser for help on using the repository browser.