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
RevLine 
[7c9f3da]1#################################################################### l
[a0db751]2# FCC-ee IDEA detector model
3#
[b8a6aa3]4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
[7c9f3da]5# email: efontane@bo.infn.it,
[a0db751]6# lorenzo.pezzotti01@universitadipavia.it,
7# m.antonello@uninsubria.it,
[b8a6aa3]8# michele.selvaggi@cern.ch
[a0db751]9#####################################################################
[b8a6aa3]10
11set B 2.0
12
[ef8a06d]13#######################################
14# Order of execution of various modules
15#######################################
16
17set ExecutionPath {
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
[a0db751]24 TrackMergerPre
25 TrackSmearing
[ef8a06d]26
[a0db751]27 TrackMerger
[ef8a06d]28 Calorimeter
29 EFlowMerger
30
31 PhotonEfficiency
32 PhotonIsolation
33
[b8a6aa3]34 MuonFilter
35
[ef8a06d]36 ElectronFilter
37 ElectronEfficiency
38 ElectronIsolation
39
40 MuonEfficiency
41 MuonIsolation
42
43 MissingET
44
45 NeutrinoFilter
46 GenJetFinder
47 GenMissingET
[a0db751]48
[ef8a06d]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
[7c9f3da]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
[ef8a06d]81
[7c9f3da]82 # magnetic field, in T
[b8a6aa3]83 set Bz $B
[ef8a06d]84}
85
86####################################
87# Charged hadron tracking efficiency
88####################################
89
90module Efficiency ChargedHadronTrackingEfficiency {
[7c9f3da]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 }
[ef8a06d]101}
102
[a0db751]103# (pt <= 0.1) * (0.00) +
[7c9f3da]104# (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
105# (abs(eta) > 3) * (0.00)
106
107
108
[ef8a06d]109##############################
110# Electron tracking efficiency
111##############################
112
113module Efficiency ElectronTrackingEfficiency {
[7c9f3da]114 set InputArray ParticlePropagator/electrons
115 set OutputArray electrons
[ef8a06d]116
117
[7c9f3da]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)
[a0db751]124 }
[ef8a06d]125}
126
[7c9f3da]127
[ef8a06d]128##########################
129# Muon tracking efficiency
130##########################
131
132module Efficiency MuonTrackingEfficiency {
[7c9f3da]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 }
[ef8a06d]143}
144
[a0db751]145##############
146# Track merger
147##############
[7c9f3da]148
[a0db751]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}
[ef8a06d]156
157
[a0db751]158########################################
159# Smearing for charged tracks
160########################################
[ef8a06d]161
[a0db751]162module TrackCovariance TrackSmearing {
163 set InputArray TrackMergerPre/tracks
164 set OutputArray tracks
165
[b8a6aa3]166
167 set InputArray TrackMergerPre/tracks
168 set OutputArray tracks
169
[a0db751]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
[b8a6aa3]193 # barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
[a0db751]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
[b8a6aa3]197 1 PIPE -100 100 0.015 0.001655 0.2805 0 0 0 0 0 0
[a0db751]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
[7c9f3da]343 }
[ef8a06d]344
[b8a6aa3]345 set Bz $B
[ef8a06d]346}
347
348##############
349# Track merger
350##############
351
352module Merger TrackMerger {
353# add InputArray InputArray
[a0db751]354 add InputArray TrackSmearing/tracks
[ef8a06d]355 set OutputArray tracks
356}
357
358
[a0db751]359#############
360# Calorimeter
361#############
[ef8a06d]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
[7c9f3da]374 set HCalEnergyMin 0.5
[a1b19ea]375 set EnergyMin 0.5
[ef8a06d]376 set ECalEnergySignificanceMin 1.0
377 set HCalEnergySignificanceMin 1.0
378 set EnergySignificanceMin 1.0
379
380 set SmearTowerCenter true
[7c9f3da]381 set pi [expr {acos(-1)}]
382
[a0db751]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
[7c9f3da]389
[a0db751]390 #barrel:
[7c9f3da]391 set PhiBins {}
392 for {set i -120} {$i <= 120} {incr i} {
393 add PhiBins [expr {$i * $pi/120}]
394 }
[a0db751]395 #deta=0.02 units for |eta| <= 0.88
[7c9f3da]396 for {set i -44} {$i < 45} {incr i} {
397 set eta [expr {$i * 0.02}]
398 add EtaPhiBins $eta $PhiBins
399 }
400
[a0db751]401 #endcaps:
[7c9f3da]402 set PhiBins {}
403 for {set i -120} {$i <= 120} {incr i} {
404 add PhiBins [expr {$i* $pi/120}]
405 }
[a0db751]406 #deta=0.02 units for 0.88 < |eta| <= 3.0
407 #first, from -3.0 to -0.88
[7c9f3da]408 for {set i 1} {$i <=106} {incr i} {
409 set eta [expr {-3.00 + $i * 0.02}]
410 add EtaPhiBins $eta $PhiBins
411 }
[a0db751]412 #same for 0.88 to 3.0
[7c9f3da]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
[a0db751]418 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
[7c9f3da]419 add EnergyFraction {0} {0.0 1.0}
[a0db751]420 # energy fractions for e, gamma and pi0
[7c9f3da]421 add EnergyFraction {11} {1.0 0.0}
422 add EnergyFraction {22} {1.0 0.0}
423 add EnergyFraction {111} {1.0 0.0}
[a0db751]424 # energy fractions for muon, neutrinos and neutralinos
[7c9f3da]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}
[a0db751]434 # energy fractions for K0short and Lambda
[7c9f3da]435 add EnergyFraction {310} {0.3 0.7}
436 add EnergyFraction {3122} {0.3 0.7}
437
438
[a0db751]439 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]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
[a0db751]445 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]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 }
[ef8a06d]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
[7c9f3da]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 }
[ef8a06d]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
[7c9f3da]496 set PTRatioMax 999.
[ef8a06d]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
[a0db751]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
[ef8a06d]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
[a0db751]535 set EfficiencyFormula {
[7c9f3da]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 }
[ef8a06d]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 {
[a0db751]565 set InputArray MuonFilter/muons
[ef8a06d]566 set OutputArray muons
567
568 # set EfficiencyFormula {efficiency as a function of eta and pt}
569
570 # efficiency formula for muons
[a0db751]571 set EfficiencyFormula {
[7c9f3da]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 }
[ef8a06d]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
[7c9f3da]650 set ParameterR 0.4
651 set JetPTMin 1.0
[ef8a06d]652}
653
[7c9f3da]654
[ef8a06d]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
[7c9f3da]677 set ParameterR 0.4
678 set JetPTMin 1.0
[ef8a06d]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
[7c9f3da]690 set ScaleFormula {1.08}
[ef8a06d]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
[7c9f3da]706 set PartonEtaMax 3.0
[ef8a06d]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}
[a0db751]719
[ef8a06d]720 # default efficiency formula (misidentification rate)
[7c9f3da]721 add EfficiencyFormula {0} {0.01}
[ef8a06d]722
723 # efficiency formula for c-jets (misidentification rate)
[7c9f3da]724 add EfficiencyFormula {4} {0.10}
[ef8a06d]725
726 # efficiency formula for b-jets
[7c9f3da]727 add EfficiencyFormula {5} {0.80}
[ef8a06d]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
[7c9f3da]741 set TauEtaMax 3.0
[ef8a06d]742
743 # default efficiency formula (misidentification rate)
[7c9f3da]744 add EfficiencyFormula {0} {0.001}
[ef8a06d]745 # efficiency formula for tau-jets
746 add EfficiencyFormula {15} {0.6}
747}
748
[7c9f3da]749
[ef8a06d]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
[7c9f3da]763
764
[ef8a06d]765##################
766# ROOT tree writer
767##################
768
[7c9f3da]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
[ef8a06d]771# "add Branch ..." lines.
772
773module TreeWriter TreeWriter {
[7c9f3da]774 # add Branch InputArray BranchName BranchClass
[a0db751]775
[7c9f3da]776 add Branch Delphes/allParticles Particle GenParticle
777
778 add Branch TrackMerger/tracks Track Track
779 add Branch Calorimeter/towers Tower Tower
[a0db751]780
[7c9f3da]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
[a0db751]788
[7c9f3da]789 add Branch GenJetFinder/jets GenJet Jet
790 add Branch GenMissingET/momentum GenMissingET MissingET
[a0db751]791
[7c9f3da]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
[a0db751]796
797 add Branch JetEnergyScale/jets AntiKtJet Jet
798
[7c9f3da]799 add Branch MissingET/momentum MissingET MissingET
800 add Branch ScalarHT/energy ScalarHT ScalarHT
[b8a6aa3]801
802 # add Info InfoName InfoValue
803 add Info Bz $B
[ef8a06d]804}
Note: See TracBrowser for help on using the repository browser.