Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ 8562e4e

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 8562e4e was 8562e4e, checked in by Michele Selvaggi <michele.selvaggi@…>, 7 years ago

first commit HLLHC generic card
  • Property mode set to 100644
File size: 58.3 KB
Line 
1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Michele Selvaggi (CERN)
5#
6# Released on: Apr. 6th, 2017
7#
8# Configuration: FCC-hh baseline detector
9#
10#######################################
11# Order of execution of various modules
12#######################################
13
14set ExecutionPath {
15
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
23 ElectronMomentumSmearing
24 MuonMomentumSmearing
25
26 DenseProp
27 DenseMergeTracks
28 TrackMergerProp
29 TrackMerger
30
31 ECal
32 HCal
33
34 Calorimeter
35 EFlowMerger
36 EFlowFilter
37
38 PhotonIsolation
39 PhotonEfficiency
40
41 ElectronFilter
42 ElectronIsolation
43 ElectronEfficiency
44
45 ChargedHadronFilter
46
47 MuonIsolation
48 MuonEfficiency
49
50 MissingET
51
52 NeutrinoFilter
53 GenJetFinder
54 GenMissingET
55
56 FastJetFinder
57 FatJetFinder
58
59 JetEnergyScale
60
61 JetFlavorAssociation
62
63 BTagging
64 TauTagging
65
66 UniqueObjectFinder
67
68 ScalarHT
69
70 TreeWriter
71}
72
73
74#################################
75# Propagate particles in cylinder
76#################################
77
78module ParticlePropagator ParticlePropagator {
79 set InputArray Delphes/stableParticles
80
81 set OutputArray stableParticles
82 set ChargedHadronOutputArray chargedHadrons
83 set ElectronOutputArray electrons
84 set MuonOutputArray muons
85
86 # radius of the magnetic field coverage, in m
87 set Radius 1.2
88 # half-length of the magnetic field coverage, in m
89 set HalfLength 3.25
90
91 # magnetic field
92 set Bz 3.0
93}
94
95####################################
96# Charged hadron tracking efficiency
97####################################
98
99module Efficiency ChargedHadronTrackingEfficiency {
100 set InputArray ParticlePropagator/chargedHadrons
101 set OutputArray chargedHadrons
102
103 # TBC (which eta_max ? which pT min?)
104
105 # tracking efficiency formula for charged hadrons
106
107 set EfficiencyFormula {
108 (pt <= 0.5) * (0.00) + \
109 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
110 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
111 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
112 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
113 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
114 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
115 (abs(eta) > 4.0) * (0.00)
116 }
117
118}
119
120##############################
121# Electron tracking efficiency
122##############################
123
124module Efficiency ElectronTrackingEfficiency {
125 set InputArray ParticlePropagator/electrons
126 set OutputArray electrons
127
128# TBC (which eta_max ?)
129# putting same as charged hadrons for now...
130
131 set EfficiencyFormula {
132 (pt <= 0.5) * (0.00) + \
133 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
134 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
135 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
136 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
137 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
138 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
139 (abs(eta) > 4.0) * (0.00)
140 }
141
142}
143##########################
144# Muon tracking efficiency
145##########################
146
147module Efficiency MuonTrackingEfficiency {
148 set InputArray ParticlePropagator/muons
149 set OutputArray muons
150
151 set EfficiencyFormula {
152 (pt <= 0.5) * (0.00) + \
153 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 1.00) + \
154 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
155 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.5 && pt <= 1.0) * (pt*1.00) + \
156 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
157 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.95) + \
158 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
159 (abs(eta) > 4.0) * (0.00)
160
161 }
162
163}
164
165########################################
166# Momentum resolution for charged tracks
167########################################
168
169module MomentumSmearing ChargedHadronMomentumSmearing {
170 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
171 set OutputArray chargedHadrons
172
173 # --- CMS resolution Phase II ---
174
175 set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
176 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
177 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
178 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
179 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
180 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
181 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
182 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
183 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
184 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
185 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
186 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
187 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
188 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
189 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
190 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
191 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
192 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
193 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
194 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
195 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
196 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
197 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
198 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
199 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
200 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
201 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
202 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
203 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
204 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
205 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
206 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
207 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
208 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
209 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
210 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
211 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
212 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
213 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
214 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
215 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
216 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
217 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
218 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
219 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
220 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
221 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
222 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
223 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
224 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
225 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
226 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
227 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
228 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
229 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
230 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
231 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
232 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
233 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
234 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
235 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
236 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
237 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
238 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
239 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
240 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
241 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
242 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
243 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
244 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
245 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
246 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
247 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
248 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
249 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
250 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
251 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
252 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
253 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
254 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
255 }
256
257}
258
259
260###################################
261# Momentum resolution for electrons
262###################################
263
264module MomentumSmearing ElectronMomentumSmearing {
265 set InputArray ElectronTrackingEfficiency/electrons
266 set OutputArray electrons
267
268 set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
269 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
270 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
271 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
272 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
273 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
274 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
275 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
276 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
277 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
278 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
279 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
280 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
281 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
282 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
283 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
284 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
285 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
286 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
287 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
288 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
289 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
290 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
291 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
292 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
293 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
294 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
295 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
296 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
297 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
298 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
299 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
300 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
301 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
302 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
303 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
304 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
305 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
306 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
307 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
308 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
309 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
310 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
311 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
312 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
313 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
314 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
315 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
316 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
317 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
318 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
319 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
320 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
321 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
322 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
323 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
324 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
325 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
326 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
327 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
328 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
329 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
330 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
331 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
332 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
333 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
334 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
335 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
336 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
337 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
338 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
339 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
340 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
341 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
342 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
343 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
344 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
345 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
346 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
347 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
348 }
349}
350
351
352###############################
353# Momentum resolution for muons
354###############################
355
356module MomentumSmearing MuonMomentumSmearing {
357 set InputArray MuonTrackingEfficiency/muons
358 set OutputArray muons
359
360 # Equal to Tk resolution up to pT = 100 GeV, then tk slope/5.
361
362 set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
363 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
364 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
365 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244/5.*pt/100.000000) + \
366 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
367 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
368 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
369 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172/5.*pt/100.000000) + \
370 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
371 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
372 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
373 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077/5.*pt/100.000000) + \
374 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
375 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
376 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
377 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177/5.*pt/100.000000) + \
378 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
379 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
380 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
381 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429/5.*pt/100.000000) + \
382 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
383 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
384 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
385 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867/5.*pt/100.000000) + \
386 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
387 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
388 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
389 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800/5.*pt/100.000000) + \
390 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
391 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
392 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
393 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161/5.*pt/100.000000) + \
394 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
395 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
396 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
397 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039/5.*pt/100.000000) + \
398 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
399 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
400 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
401 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098/5.*pt/100.000000) + \
402 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
403 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
404 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
405 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289/5.*pt/100.000000) + \
406 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
407 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
408 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
409 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524/5.*pt/100.000000) + \
410 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
411 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
412 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
413 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477/5.*pt/100.000000) + \
414 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
415 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
416 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
417 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392/5.*pt/100.000000) + \
418 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
419 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
420 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
421 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445/5.*pt/100.000000) + \
422 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
423 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
424 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
425 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053/5.*pt/100.000000) + \
426 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
427 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
428 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
429 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710/5.*pt/100.000000) + \
430 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
431 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
432 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
433 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629/5.*pt/100.000000) + \
434 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
435 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
436 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
437 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345/5.*pt/100.000000) + \
438 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
439 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
440 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
441 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617/5.*pt/100.000000))
442 }
443
444}
445
446##############
447# Track merger
448##############
449
450module Merger TrackMergerProp {
451# add InputArray InputArray
452 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
453 add InputArray ElectronMomentumSmearing/electrons
454 add InputArray MuonMomentumSmearing/muons
455 set OutputArray tracks
456}
457
458
459#################################
460# Dense Track propagation
461#################################
462
463module ParticlePropagator DenseProp {
464 set InputArray Delphes/stableParticles
465
466 set OutputArray stableParticles
467 set NeutralOutputArray neutralParticles
468 set ChargedHadronOutputArray chargedHadrons
469 set ElectronOutputArray electrons
470 set MuonOutputArray muons
471
472 # radius of the magnetic field coverage, in m
473 set Radius 0.3
474 set RadiusMax 1.2
475 # half-length of the magnetic field coverage, in m
476 set HalfLength 0.7
477 set HalfLengthMax 3.25
478
479 # magnetic field
480 set Bz 3.0
481}
482
483
484####################
485# Dense Track merger
486###################
487
488module Merger DenseMergeTracks {
489# add InputArray InputArray
490 add InputArray DenseProp/chargedHadrons
491 add InputArray DenseProp/electrons
492 add InputArray DenseProp/muons
493 set OutputArray tracks
494}
495
496
497######################
498# Dense Track Filter
499######################
500
501module DenseTrackFilter TrackMerger {
502 set TrackInputArray TrackMergerProp/tracks
503 set DenseChargedInputArray DenseMergeTracks/tracks
504
505 set TrackOutputArray tracks
506
507 set EtaPhiRes 0.003
508 set EtaMax 4.0
509
510 set pi [expr {acos(-1)}]
511
512 set nbins_phi [expr {$pi/$EtaPhiRes} ]
513 set nbins_phi [expr {int($nbins_phi)} ]
514
515 set PhiBins {}
516 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
517 add PhiBins [expr {$i * $pi/$nbins_phi}]
518 }
519
520 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
521 set nbins_eta [expr {int($nbins_eta)} ]
522
523 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
524 set eta [expr {$i * $EtaPhiRes}]
525 add EtaPhiBins $eta $PhiBins
526 }
527}
528
529
530#############
531# ECAL
532#############
533
534module SimpleCalorimeter ECal {
535 set ParticleInputArray ParticlePropagator/stableParticles
536 set TrackInputArray TrackMerger/tracks
537
538 set TowerOutputArray ecalTowers
539 set EFlowTrackOutputArray eflowTracks
540 set EFlowTowerOutputArray eflowPhotons
541
542 set IsEcal true
543
544 set EnergyMin 0.5
545 set EnergySignificanceMin 1.0
546
547 set SmearTowerCenter true
548
549 set pi [expr {acos(-1)}]
550
551 # lists of the edges of each tower in eta and phi
552 # each list starts with the lower edge of the first tower
553 # the list ends with the higher edged of the last tower
554
555 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
556
557 set PhiBins {}
558 for {set i -180} {$i <= 180} {incr i} {
559 add PhiBins [expr {$i * $pi/180.0}]
560 }
561
562 # 0.02 unit in eta up to eta = 1.5 (barrel)
563 for {set i -85} {$i <= 86} {incr i} {
564 set eta [expr {$i * 0.0174}]
565 add EtaPhiBins $eta $PhiBins
566 }
567
568 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
569
570 set PhiBins {}
571 for {set i -180} {$i <= 180} {incr i} {
572 add PhiBins [expr {$i * $pi/180.0}]
573 }
574
575 # 0.02 unit in eta up to eta = 3
576 for {set i 1} {$i <= 84} {incr i} {
577 set eta [expr { -2.958 + $i * 0.0174}]
578 add EtaPhiBins $eta $PhiBins
579 }
580
581 for {set i 1} {$i <= 84} {incr i} {
582 set eta [expr { 1.4964 + $i * 0.0174}]
583 add EtaPhiBins $eta $PhiBins
584 }
585
586 # take present CMS granularity for HF
587
588 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
589 set PhiBins {}
590 for {set i -18} {$i <= 18} {incr i} {
591 add PhiBins [expr {$i * $pi/18.0}]
592 }
593
594 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
595 add EtaPhiBins $eta $PhiBins
596 }
597
598
599 add EnergyFraction {0} {0.0}
600 # energy fractions for e, gamma and pi0
601 add EnergyFraction {11} {1.0}
602 add EnergyFraction {22} {1.0}
603 add EnergyFraction {111} {1.0}
604 # energy fractions for muon, neutrinos and neutralinos
605 add EnergyFraction {12} {0.0}
606 add EnergyFraction {13} {0.0}
607 add EnergyFraction {14} {0.0}
608 add EnergyFraction {16} {0.0}
609 add EnergyFraction {1000022} {0.0}
610 add EnergyFraction {1000023} {0.0}
611 add EnergyFraction {1000025} {0.0}
612 add EnergyFraction {1000035} {0.0}
613 add EnergyFraction {1000045} {0.0}
614 # energy fractions for K0short and Lambda
615 add EnergyFraction {310} {0.3}
616 add EnergyFraction {3122} {0.3}
617
618
619 # ECAL is CMS
620 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
621 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
622
623 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
624 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
625 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
626 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
627 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
628
629
630}
631
632#############
633# HCAL
634#############
635module SimpleCalorimeter HCal {
636 set ParticleInputArray ParticlePropagator/stableParticles
637 set TrackInputArray ECal/eflowTracks
638
639 set TowerOutputArray hcalTowers
640 set EFlowTrackOutputArray eflowTracks
641 set EFlowTowerOutputArray eflowNeutralHadrons
642
643 set IsEcal false
644
645 set EnergyMin 1.0
646 set EnergySignificanceMin 2.0
647
648 set SmearTowerCenter true
649
650 set pi [expr {acos(-1)}]
651
652 # lists of the edges of each tower in eta and phi
653 # each list starts with the lower edge of the first tower
654 # the list ends with the higher edged of the last tower
655
656 # 10 degrees towers
657 set PhiBins {}
658 for {set i -18} {$i <= 18} {incr i} {
659 add PhiBins [expr {$i * $pi/18.0}]
660 }
661 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
662 add EtaPhiBins $eta $PhiBins
663 }
664
665 # 20 degrees towers
666 set PhiBins {}
667 for {set i -9} {$i <= 9} {incr i} {
668 add PhiBins [expr {$i * $pi/9.0}]
669 }
670 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
671 add EtaPhiBins $eta $PhiBins
672 }
673
674 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
675 add EnergyFraction {0} {1.0}
676 # energy fractions for e, gamma and pi0
677 add EnergyFraction {11} {0.0}
678 add EnergyFraction {22} {0.0}
679 add EnergyFraction {111} {0.0}
680 # energy fractions for muon, neutrinos and neutralinos
681 add EnergyFraction {12} {0.0}
682 add EnergyFraction {13} {0.0}
683 add EnergyFraction {14} {0.0}
684 add EnergyFraction {16} {0.0}
685 add EnergyFraction {1000022} {0.0}
686 add EnergyFraction {1000023} {0.0}
687 add EnergyFraction {1000025} {0.0}
688 add EnergyFraction {1000035} {0.0}
689 add EnergyFraction {1000045} {0.0}
690 # energy fractions for K0short and Lambda
691 add EnergyFraction {310} {0.7}
692 add EnergyFraction {3122} {0.7}
693
694 # HCAL is ALTAS
695 # http://arxiv.org/pdf/hep-ex/0004009v1
696 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
697 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
698 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
699 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
700 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
701}
702
703
704#################
705# Electron filter
706#################
707
708module PdgCodeFilter ElectronFilter {
709 set InputArray HCal/eflowTracks
710 set OutputArray electrons
711 set Invert true
712 add PdgCode {11}
713 add PdgCode {-11}
714}
715
716
717######################
718# ChargedHadronFilter
719######################
720
721module PdgCodeFilter ChargedHadronFilter {
722 set InputArray HCal/eflowTracks
723 set OutputArray chargedHadrons
724
725 add PdgCode {11}
726 add PdgCode {-11}
727 add PdgCode {13}
728 add PdgCode {-13}
729}
730
731
732###################################################
733# Tower Merger (in case not using e-flow algorithm)
734###################################################
735
736module Merger Calorimeter {
737# add InputArray InputArray
738 add InputArray ECal/ecalTowers
739 add InputArray HCal/hcalTowers
740 set OutputArray towers
741}
742
743####################
744# Energy flow merger
745####################
746
747module Merger EFlowMerger {
748# add InputArray InputArray
749 add InputArray HCal/eflowTracks
750 add InputArray ECal/eflowPhotons
751 add InputArray HCal/eflowNeutralHadrons
752 set OutputArray eflow
753}
754
755######################
756# EFlowFilter
757######################
758
759module PdgCodeFilter EFlowFilter {
760 set InputArray EFlowMerger/eflow
761 set OutputArray eflow
762
763 add PdgCode {11}
764 add PdgCode {-11}
765 add PdgCode {13}
766 add PdgCode {-13}
767}
768
769
770###################
771# Missing ET merger
772###################
773
774module Merger MissingET {
775# add InputArray InputArray
776 add InputArray EFlowMerger/eflow
777 set MomentumOutputArray momentum
778}
779
780
781
782##################
783# Scalar HT merger
784##################
785
786module Merger ScalarHT {
787# add InputArray InputArray
788 add InputArray EFlowMerger/eflow
789 set EnergyOutputArray energy
790}
791
792#################
793# Neutrino Filter
794#################
795
796module PdgCodeFilter NeutrinoFilter {
797
798 set InputArray Delphes/stableParticles
799 set OutputArray filteredParticles
800
801 set PTMin 0.0
802
803 add PdgCode {12}
804 add PdgCode {14}
805 add PdgCode {16}
806 add PdgCode {-12}
807 add PdgCode {-14}
808 add PdgCode {-16}
809
810}
811
812
813#####################
814# MC truth jet finder
815#####################
816
817module FastJetFinder GenJetFinder {
818 set InputArray NeutrinoFilter/filteredParticles
819
820 set OutputArray jets
821
822 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
823 set JetAlgorithm 6
824 set ParameterR 0.4
825
826 set JetPTMin 15.0
827}
828
829#########################
830# Gen Missing ET merger
831########################
832
833module Merger GenMissingET {
834
835# add InputArray InputArray
836 add InputArray NeutrinoFilter/filteredParticles
837 set MomentumOutputArray momentum
838}
839
840
841############
842# Jet finder
843############
844
845module FastJetFinder FastJetFinder {
846# set InputArray Calorimeter/towers
847 set InputArray EFlowMerger/eflow
848
849 set OutputArray jets
850
851 set JetAlgorithm 6
852 set ParameterR 0.4
853
854 set ComputeNsubjettiness 1
855 set Beta 1.0
856 set AxisMode 4
857
858 set ComputeSoftDrop 1
859 set BetaSoftDrop 0.0
860 set SymmetryCutSoftDrop 0.1
861 set R0SoftDrop 0.4
862
863 set JetPTMin 30.0
864}
865
866##################
867# Fat Jet finder
868##################
869
870module FastJetFinder FatJetFinder {
871 set InputArray EFlowMerger/eflow
872
873 set OutputArray jets
874
875 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
876 set JetAlgorithm 6
877 set ParameterR 0.8
878
879 set ComputeNsubjettiness 1
880 set Beta 1.0
881 set AxisMode 4
882
883 set ComputeSoftDrop 1
884 set BetaSoftDrop 0.0
885 set SymmetryCutSoftDrop 0.1
886 set R0SoftDrop 0.8
887
888 set JetPTMin 30.0
889}
890
891
892
893##################
894# Jet Energy Scale
895##################
896
897module EnergyScale JetEnergyScale {
898 set InputArray FastJetFinder/jets
899 set OutputArray jets
900
901 # scale formula for jets
902 set ScaleFormula {1.00}
903}
904
905
906########################
907# Jet Flavor Association
908########################
909
910module JetFlavorAssociation JetFlavorAssociation {
911
912 set PartonInputArray Delphes/partons
913 set ParticleInputArray Delphes/allParticles
914 set ParticleLHEFInputArray Delphes/allParticlesLHEF
915 set JetInputArray JetEnergyScale/jets
916
917 set DeltaR 0.5
918 set PartonPTMin 5.0
919 set PartonEtaMax 4.0
920
921}
922
923##################
924# Photon isolation
925##################
926
927# TBC: check values for iso cuts
928
929module Isolation PhotonIsolation {
930 set CandidateInputArray ECal/eflowPhotons
931 set IsolationInputArray EFlowFilter/eflow
932
933 set OutputArray photons
934
935 set DeltaRMax 0.3
936 set PTMin 0.0
937 set PTRatioMax 0.1
938}
939
940
941
942###################
943# Photon efficiency
944###################
945
946module Efficiency PhotonEfficiency {
947 set InputArray PhotonIsolation/photons
948 set OutputArray photons
949
950 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
951
952 set EfficiencyFormula {
953 (pt <= 1.0) * (0.00) + \
954 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
955 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
956 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
957
958 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
959 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
960 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
961
962 (abs(eta) > 4.0) * (0.00)}
963
964}
965
966
967
968####################
969# Electron isolation
970####################
971
972# TBC: check values for iso cuts
973
974module Isolation ElectronIsolation {
975 set CandidateInputArray ElectronFilter/electrons
976 set IsolationInputArray EFlowFilter/eflow
977
978 set OutputArray electrons
979
980 set DeltaRMax 0.3
981
982 set PTMin 0.0
983
984 set PTRatioMax 0.1
985}
986
987
988#######################
989# Electron efficiency #
990#######################
991
992module Efficiency ElectronEfficiency {
993
994 set InputArray ElectronIsolation/electrons
995 set OutputArray electrons
996
997 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
998 # efficiency formula for electrons (taken from CMS PhaseII card)
999 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
1000 set EfficiencyFormula {
1001 (pt <= 2.0) * (0.00) +
1002 (abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1003 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1004 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1005 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1006 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1007 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1008 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1009 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1010 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1011 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1012 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1013 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1014 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1015 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1016 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1017 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1018 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1019 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1020 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1021 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1022 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1023 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1024 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1025 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1026 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1027 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1028 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1029 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1030 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1031 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1032 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1033 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1034 (abs(eta) > 4.0) * (0.00)
1035
1036 }
1037
1038}
1039
1040
1041
1042################
1043# Muon isolation
1044################
1045
1046# TBC: check values for iso cuts
1047
1048module Isolation MuonIsolation {
1049 set CandidateInputArray MuonMomentumSmearing/muons
1050 set IsolationInputArray EFlowFilter/eflow
1051
1052 set OutputArray muons
1053
1054 set DeltaRMax 0.3
1055
1056 set PTMin 0.0
1057
1058 set PTRatioMax 0.2
1059}
1060
1061####################
1062# Muon Efficiency #
1063####################
1064
1065module Efficiency MuonEfficiency {
1066 set InputArray MuonIsolation/muons
1067 set OutputArray muons
1068 # tracking + TightID efficiency formula for muons
1069 set EfficiencyFormula {
1070 (pt <= 2.0) * (0.00) +
1071 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1072 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1073 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1074 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1075 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1076 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1077 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1078 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1079 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1080 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1081 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1082 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
1083 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1084 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
1085 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
1086 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1087 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
1088 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
1089 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1090 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
1091 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1092 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1093 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1094 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1095 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
1096 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1097 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
1098 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
1099 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
1100 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1101 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1102 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1103 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1104 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
1105 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
1106 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1107 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1108 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1109 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1110 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
1111 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1112 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1113 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1114 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1115 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1116 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1117 (abs(eta) > 4.0) * (0.00)
1118 }
1119}
1120
1121
1122
1123###########
1124# b-tagging
1125###########
1126
1127module BTagging BTagging {
1128 set JetInputArray JetEnergyScale/jets
1129
1130 set BitNumber 0
1131
1132 add EfficiencyFormula {0} {
1133
1134 (pt <= 10.0) * (0.00) +
1135 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
1136 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
1137 (abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
1138 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
1139 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
1140 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1141 (abs(eta) > 4.0) * (0.00)}
1142
1143 add EfficiencyFormula {4} {
1144
1145 (pt <= 10.0) * (0.00) +
1146 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
1147 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
1148 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1149 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
1150 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
1151 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1152 (abs(eta) > 4.0) * (0.00)}
1153
1154 add EfficiencyFormula {5} {
1155
1156 (pt <= 10.0) * (0.00) +
1157 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
1158 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
1159 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
1160 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
1161 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
1162 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
1163 (abs(eta) >= 4.0) * (0.00)}
1164
1165}
1166
1167
1168#############
1169# tau-tagging
1170#############
1171
1172
1173module TauTagging TauTagging {
1174 set ParticleInputArray Delphes/allParticles
1175 set PartonInputArray Delphes/partons
1176 set JetInputArray JetEnergyScale/jets
1177
1178 set DeltaR 0.5
1179
1180 set TauPTMin 1.0
1181
1182 set TauEtaMax 4.0
1183
1184 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1185 add EfficiencyFormula {0} {
1186
1187 (pt <= 10.0) * (0.00) +
1188 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
1189 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
1190 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1191 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
1192 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
1193 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1194 (abs(eta) > 4.0) * (0.00)}
1195
1196 add EfficiencyFormula {11} {
1197
1198 (pt <= 10.0) * (0.00) +
1199 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
1200 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
1201 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1202 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
1203 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
1204 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1205 (abs(eta) > 4.0) * (0.00)}
1206
1207 add EfficiencyFormula {15} {
1208
1209 (pt <= 10.0) * (0.00) +
1210 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
1211 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
1212 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1213 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
1214 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
1215 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1216 (abs(eta) > 4.0) * (0.00)}
1217
1218}
1219
1220#####################################################
1221# Find uniquely identified photons/electrons/tau/jets
1222#####################################################
1223
1224module UniqueObjectFinder UniqueObjectFinder {
1225# earlier arrays take precedence over later ones
1226# add InputArray InputArray OutputArray
1227 add InputArray PhotonEfficiency/photons photons
1228 add InputArray ElectronEfficiency/electrons electrons
1229 add InputArray MuonEfficiency/muons muons
1230 add InputArray JetEnergyScale/jets jets
1231
1232}
1233
1234
1235##################
1236# ROOT tree writer
1237##################
1238
1239module TreeWriter TreeWriter {
1240# add Branch InputArray BranchName BranchClass
1241 add Branch Delphes/allParticles Particle GenParticle
1242
1243 add Branch GenJetFinder/jets GenJet Jet
1244 add Branch GenMissingET/momentum GenMissingET MissingET
1245
1246 add Branch TrackMerger/tracks Track Track
1247 add Branch Calorimeter/towers Tower Tower
1248
1249 add Branch HCal/eflowTracks EFlowTrack Track
1250 add Branch ECal/eflowPhotons EFlowPhoton Tower
1251 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1252
1253 add Branch UniqueObjectFinder/photons Photon Photon
1254 add Branch UniqueObjectFinder/electrons Electron Electron
1255 add Branch UniqueObjectFinder/muons Muon Muon
1256 add Branch UniqueObjectFinder/jets Jet Jet
1257
1258 add Branch FatJetFinder/jets FatJet Jet
1259
1260 add Branch MissingET/momentum MissingET MissingET
1261 add Branch ScalarHT/energy ScalarHT ScalarHT
1262}
Note: See TracBrowser for help on using the repository browser.