Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ 85dc372

Last change on this file since 85dc372 was fa0b54e, checked in by michele <michele.selvaggi@…>, 3 years ago

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