Fork me on GitHub

source: git/cards/delphes_card_FCCeh.tcl

Last change on this file was c856b3a, checked in by Michele Selvaggi <michele.selvaggi@…>, 4 years ago

added LHeC and FCCeh cards
  • Property mode set to 100644
File size: 29.9 KB
RevLine 
[c856b3a]1#######################################
2# Uta Klein
3# Order of execution of various modules
4# 19.12.16 add paramterisations according to Peter Kostka's design
5# muons not considered here - take CMS/FCC values
6# 22.12.16 new extensions for FCC-eh based on new LHeC cards
7# needs eflow and photons and JetFlavour Ass
8# 3.3.17 add more Filter 16.3.17 smaller constant CAL terms c.f. Max Klein
9# 17.3.17 add track resolutions as given by Peter Kostka
10#######################################
11
12set ExecutionPath {
13 ParticlePropagator
14
15 ChargedHadronTrackingEfficiency
16 ElectronTrackingEfficiency
17 MuonTrackingEfficiency
18
19 ChargedHadronMomentumSmearing
20 ElectronMomentumSmearing
21 MuonMomentumSmearing
22
23 TrackMerger
24 ImpactParameterSmearing
25 Calorimeter
26 EFlowMerger
27 EFlowFilter
28
29 PhotonEfficiency
30 PhotonIsolation
31
32 ElectronFilter
33 ElectronEfficiency
34 ElectronIsolation
35
36
37 MuonEfficiency
38 MuonIsolation
39
40 NeutrinoFilter
41
42 MissingET
43 GenMissingET
44 GenJetFinder
45 FastJetFinder
46 JetEnergyScale
47
48 GenJetFlavorAssociation
49 JetFlavorAssociation
50
51 GenBTagging
52 BTagging
53
54 GenCTagging
55 CTagging
56
57 GenTauTagging
58 TauTagging
59
60 UniqueObjectFinder
61
62 ScalarHT
63
64 TreeWriter
65}
66
67# exclude since it adds also an entry to Jet->BTag
68# TrackCountingBTagging
69
70#################################
71# Propagate particles in cylinder
72#################################
73
74module ParticlePropagator ParticlePropagator {
75 set InputArray Delphes/stableParticles
76
77 set OutputArray stableParticles
78 set ChargedHadronOutputArray chargedHadrons
79 set ElectronOutputArray electrons
80 set MuonOutputArray muons
81
82# FCC-eh
83# radius of the magnetic field coverage, in m
84 set Radius 1.36
85 # half-length of the magnetic field coverage, in m
86 set HalfLength 5.24
87
88 # magnetic field
89 set Bz 3.8
90
91}
92
93####################################
94# Charged hadron tracking efficiency
95####################################
96
97module Efficiency ChargedHadronTrackingEfficiency {
98 set InputArray ParticlePropagator/chargedHadrons
99 set OutputArray chargedHadrons
100
101# add EfficiencyFormula {efficiency formula as a function of eta and pt}
102
103# tracking efficiency formula for charged hadrons
104# set to 1 for full range (-0.2 eta already subtracted from edges)
105 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( eta <= 5.2 && eta >= -5.0) * (1.0) + \
106 (eta > 5.2 || eta < -5.0 ) * (0.00)}
107}
108##############################
109# Electron tracking efficiency
110##############################
111
112module Efficiency ElectronTrackingEfficiency {
113 set InputArray ParticlePropagator/electrons
114 set OutputArray electrons
115
116 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
117
118 # tracking efficiency formula for electrons
119# set to 1
120 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( (eta <=5.2 && eta >= -5.0) ) * (1.0) + \
121 ((eta > 5.2 || eta < -5.0)) * (0.00)}
122 }
123##########################
124# Muon tracking efficiency
125##########################
126
127module Efficiency MuonTrackingEfficiency {
128 set InputArray ParticlePropagator/muons
129 set OutputArray muons
130
131 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
132
133 # tracking efficiency formula for muons
134# set to 1
135 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *(abs(eta) <= 4) * (1.0) + \
136 (abs(eta) > 4) * (0.00)}
137}
138
139
140########################################
141# Momentum resolution for charged tracks
142########################################
143
144module MomentumSmearing ChargedHadronMomentumSmearing {
145 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
146 set OutputArray chargedHadrons
147
148 # set ResolutionFormula {resolution formula as a function of eta and pt}
149### 17.3.17 Peter's values
150 set ResolutionFormula { (pt <= 0.1 ) * (0.00) + ( pt > 0.1 && eta < 5.2 && eta > -5.0)*sqrt( pt^2*5.e-5^2) + \
151 (eta > 5.2 || eta < -5.0) * (0.00)} }
152
153 # resolution formula for charged hadrons
154# taken from FCC but adjust to ep
155#the abs(eta) <= 2.0 should be the "central best" and abs(eta) >
156#2.0 and <= 4.9 the fwd / <= 4.3 the bwd (LHeC)
157# and abs(eta) > 2.0 and <= 5.2 fwd / <= 5.0 the bwd (FCC)
158#
159# set ResolutionFormula { (abs(eta) <= 2.0) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
160# (abs(eta) <= 2.0) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
161# (abs(eta) <= 2.0) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
162# (abs(eta) <= 2.0) * (pt > 2.0e2) * (0.05) + \
163# (abs(eta) > 2.0 && ( eta <= 5.2 && eta >= -5.0 ) ) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
164# (abs(eta) > 2.0 && ( eta <= 5.2 && eta >= -5.0 ) ) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
165# (abs(eta) > 2.0 && ( eta <= 5.2 && eta >= -5.0 ) ) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
166# (abs(eta) > 2.0 && ( eta <= 5.2 && eta >= -5.0 ) ) * (pt > 2.0e2) * (0.05)}
167#}
168
169
170#################################
171# Energy resolution for electrons
172#################################
173
174module MomentumSmearing ElectronMomentumSmearing {
175 set InputArray ElectronTrackingEfficiency/electrons
176 set OutputArray electrons
177
178
179### 17.3.17 Peter's values
180 set ResolutionFormula { (pt <= 0.1 ) * (0.00) + ( pt > 0.1 && eta < 5.2 && eta > -5.0)*sqrt( pt^2*5.e-5^2) + \
181 (eta > 5.2 || eta < -5.0) * (0.00)} }
182
183# 22.12.16 Uta
184## take value from ECAL WITHOUT some energy ranges
185## usually there were energy ranges given
186
187# set ResolutionFormula { ( eta <= 5.5 && eta > 3.0 && energy > 0.1 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) + \
188# ( eta <= 3.0 && eta >= -2.7 && energy > 0.1 ) * sqrt(energy*0.09^2 + energy^2*0.02^2) +
189# \ ( eta < -2.7 && eta >= -5.2 && energy > 0.1 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) }
190#
191#}
192
193###############################
194# Momentum resolution for muons
195###############################
196
197module MomentumSmearing MuonMomentumSmearing {
198 set InputArray MuonTrackingEfficiency/muons
199 set OutputArray muons
200
201 # set ResolutionFormula {resolution formula as a function of eta and pt}
202
203### 17.3.17 Peter's values
204 set ResolutionFormula { (pt <= 0.1 ) * (0.00) + ( pt > 0.1 && eta < 5.2 && eta > -5.0)*sqrt( pt^2*5.e-5^2) + \
205 (eta > 5.2 || eta < -5.0) * (0.00)} }
206
207
208 # resolution formula for muons
209# taken from FCC
210# set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
211# (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
212# (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
213# (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
214# (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
215# (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
216# (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
217# (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
218# (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
219# (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
220# (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
221# (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
222#}
223
224
225##############
226# Track merger
227##############
228
229module Merger TrackMerger {
230# add InputArray InputArray
231 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
232# add InputArray ElectronEnergySmearing/electrons
233# change here to the momentum smearing
234 add InputArray ElectronMomentumSmearing/electrons
235 add InputArray MuonMomentumSmearing/muons
236 set OutputArray tracks
237}
238###################
239# Impact Parameters
240###################
241
242module ImpactParameterSmearing ImpactParameterSmearing {
243 set TrackInputArray TrackMerger/tracks
244 set OutputArray tracks
245
246# set ResolutionFormula {resolution formula as a function of eta and pt}
247# taken from FCC example - and about achieved at ATLAS Run-2
248# absolute impact parameter smearing formula (in mm) as a function of pt and eta
249#set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \
250# (pt > 5.0) * (0.005)}
251#
252#}
253# value from Peter 5mu m x 20 mu m /( p sin theta ^3/2)
254set ResolutionFormula { (pt > 0.5 ) * ( sqrt( (0.0001/pt)^2/(sin(2*atan(exp(-eta))))^3 ) ) }
255
256##
257
258##################
259# LHeC Calorimeter
260##################
261
262module Calorimeter Calorimeter {
263 set ParticleInputArray ParticlePropagator/stableParticles
264## 17.3.17 use also TrackMerger/tracks
265 set TrackInputArray TrackMerger/tracks
266# set TrackInputArray ImpactParameterSmearing/tracks
267
268 set TowerOutputArray towers
269 set PhotonOutputArray photons
270
271 set EFlowTrackOutputArray eflowTracks
272 set EFlowPhotonOutputArray eflowPhotons
273 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
274
275 set pi [expr {acos(-1)}]
276
277 # lists of the edges of each tower in eta and phi
278 # each list starts with the lower edge of the first tower
279 # the list ends with the higher edged of the last tower
280
281 # use a delta_phi=0.025 for 2*phi range == 252 cells
282 # use a delta_eta-0.025 for eta = -5 to 5 == 400 cells
283 # Uta 22.12.16 use eta range from Peter Kostka -4.7 - 5.2
284 set PhiBins {}
285 for {set i -126} {$i <= 126} {incr i} {
286 add PhiBins [expr {$i * $pi/126.0}]
287 }
288 foreach eta { -5.600 -5.575 -5.550 -5.525 -5.500 -5.475 -5.450 -5.425 -5.400 -5.375 -5.350 -5.325 -5.300 -5.275 -5.250 -5.225 -5.200 -5.175 -5.150 -5.125 -5.100 -5.075 -5.050 -5.025 -5.000 -4.975 -4.950 -4.925 -4.900 -4.875 -4.850 -4.825 -4.800 -4.775 -4.750 -4.725 -4.700 -4.675 -4.650 -4.625 -4.600 -4.575 -4.550 -4.525 -4.500 -4.475 -4.450 -4.425 -4.400 -4.375 -4.350 -4.325 -4.300 -4.275 -4.250 -4.225 -4.200 -4.175 -4.150 -4.125 -4.100 -4.075 -4.050 -4.025 -4.000 -3.975 -3.950 -3.925 -3.900 -3.875 -3.850 -3.825 -3.800 -3.775 -3.750 -3.725 -3.700 -3.675 -3.650 -3.625 -3.600 -3.575 -3.550 -3.525 -3.500 -3.475 -3.450 -3.425 -3.400 -3.375 -3.350 -3.325 -3.300 -3.275 -3.250 -3.225 -3.200 -3.175 -3.150 -3.125 -3.100 -3.075 -3.050 -3.025 -3.000 -2.975 -2.950 -2.925 -2.900 -2.875 -2.850 -2.825 -2.800 -2.775 -2.750 -2.725 -2.700 -2.675 -2.650 -2.625 -2.600 -2.575 -2.550 -2.525 -2.500 -2.475 -2.450 -2.425 -2.400 -2.375 -2.350 -2.325 -2.300 -2.275 -2.250 -2.225 -2.200 -2.175 -2.150 -2.125 -2.100 -2.075 -2.050 -2.025 -2.000 -1.975 -1.950 -1.925 -1.900 -1.875 -1.850 -1.825 -1.800 -1.775 -1.750 -1.725 -1.700 -1.675 -1.650 -1.625 -1.600 -1.575 -1.550 -1.525 -1.500 -1.475 -1.450 -1.425 -1.400 -1.375 -1.350 -1.325 -1.300 -1.275 -1.250 -1.225 -1.200 -1.175 -1.150 -1.125 -1.100 -1.075 -1.050 -1.025 -1.000 -0.975 -0.950 -0.925 -0.900 -0.875 -0.850 -0.825 -0.800 -0.775 -0.750 -0.725 -0.700 -0.675 -0.650 -0.625 -0.600 -0.575 -0.550 -0.525 -0.500 -0.475 -0.450 -0.425 -0.400 -0.375 -0.350 -0.325 -0.300 -0.275 -0.250 -0.225 -0.200 -0.175 -0.150 -0.125 -0.100 -0.075 -0.050 -0.025 0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 0.225 0.250 0.275 0.300 0.325 0.350 0.375 0.400 0.425 0.450 0.475 0.500 0.525 0.550 0.575 0.600 0.625 0.650 0.675 0.700 0.725 0.750 0.775 0.800 0.825 0.850 0.875 0.900 0.925 0.950 0.975 1.000 1.025 1.050 1.075 1.100 1.125 1.150 1.175 1.200 1.225 1.250 1.275 1.300 1.325 1.350 1.375 1.400 1.425 1.450 1.475 1.500 1.525 1.550 1.575 1.600 1.625 1.650 1.675 1.700 1.725 1.750 1.775 1.800 1.825 1.850 1.875 1.900 1.925 1.950 1.975 2.000 2.025 2.050 2.075 2.100 2.125 2.150 2.175 2.200 2.225 2.250 2.275 2.300 2.325 2.350 2.375 2.400 2.425 2.450 2.475 2.500 2.525 2.550 2.575 2.600 2.625 2.650 2.675 2.700 2.725 2.750 2.775 2.800 2.825 2.850 2.875 2.900 2.925 2.950 2.975 3.000 3.025 3.050 3.075 3.100 3.125 3.150 3.175 3.200 3.225 3.250 3.275 3.300 3.325 3.350 3.375 3.400 3.425 3.450 3.475 3.500 3.525 3.550 3.575 3.600 3.625 3.650 3.675 3.700 3.725 3.750 3.775 3.800 3.825 3.850 3.875 3.900 3.925 3.950 3.975 4.000 4.025 4.050 4.075 4.100 4.125 4.150 4.175 4.200 4.225 4.250 4.275 4.300 4.325 4.350 4.375 4.400 4.425 4.450 4.475 4.500 4.525 4.550 4.575 4.600 4.625 4.650 4.675 4.700 4.725 4.750 4.775 4.800 4.825 4.850 4.875 4.900 4.925 4.950 4.975 5.00 5.025 5.050 5.075 5.100 5.125 5.150 5.175 5.200 5.225 5.250 5.275 5.300 5.325 5.350 5.375 5.400 5.425 5.450 5.475 5.500 5.525 5.550 5.575 5.600 5.625 5.650 5.675 5.700 5.725 5.750 5.775 5.800 5.825 5.850 5.875 5.900 5.925 5.950 5.975 6.000
289 } {
290 add EtaPhiBins $eta $PhiBins
291 }
292
293
294 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
295 add EnergyFraction {0} {0.0 1.0}
296 # energy fractions for e, gamma and pi0
297 add EnergyFraction {11} {1.0 0.0}
298 add EnergyFraction {22} {1.0 0.0}
299 add EnergyFraction {111} {1.0 0.0}
300 # energy fractions for muon, neutrinos and neutralinos
301 add EnergyFraction {12} {0.0 0.0}
302 add EnergyFraction {13} {0.0 0.0}
303 add EnergyFraction {14} {0.0 0.0}
304 add EnergyFraction {16} {0.0 0.0}
305 add EnergyFraction {1000022} {0.0 0.0}
306 add EnergyFraction {1000023} {0.0 0.0}
307 add EnergyFraction {1000025} {0.0 0.0}
308 add EnergyFraction {1000035} {0.0 0.0}
309 add EnergyFraction {1000045} {0.0 0.0}
310 # energy fractions for K0short and Lambda
311 add EnergyFraction {310} {0.3 0.7}
312 add EnergyFraction {3122} {0.3 0.7}
313
314## no minimum p range given (Uta 22.12.16) - update with max 16.3.17 improved constant terms
315## sigma_E = E* (a/sqrt(E) + b) = sqrt ( E*a^2 E +E^2*b^2 )
316## a = sampling term, b = constant term
317# ECAL eta : 5.0 - 2.7 - -2.1 - -4.4 (subtract -0.1 in eta as safety distance to edges)
318 set ECalResolutionFormula { ( eta <= 5.5 && eta > 3.0 && energy > 0.1 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) + \
319 ( eta <= 3.0 && eta >= -2.7 && energy > 0.1 ) * sqrt(energy*0.09^2 + energy^2*0.01^2) +
320 \ ( eta < -2.7 && eta >= -5.2 && energy > 0.1 ) * sqrt(energy*0.1^2 + energy^2*0.01^2) }
321
322#### HCAL eta : 5.1 - 2.1 - -1.7 - -4.6 (subtract -0.1 in eta as safety distance to edges)
323#set HCalResolutionFormula {resolution formula as a function of eta and energy}
324
325 set HCalResolutionFormula { ( eta >= 5.9 && eta > 2.5 ) * sqrt(energy*0.3^2 + energy^2*0.01^2) + \
326 ( eta <=2.5 && eta >= -2.2 ) * sqrt(energy*0.4^2 + energy^2*0.01^2) + \
327 ( eta < -2.2 && eta >= -5.5 ) * sqrt(energy*0.4^2 + energy^2*0.02^2) }
328
329### HCal from Tanaka-san LHEC Hbb study
330#set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
331# (abs(eta) > 3.0 && abs(eta) <= 4.0) * sqrt(energy^2*0.02^2 + energy*0.2^2) +
332# \ (abs(eta) > 4.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.10^2 + energy*0.4^2)}
333
334# uta gues inspired by ILC particle flow resolution
335 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
336# set HCalResolutionFormula { (abs(eta) <= 2.0) * sqrt(energy^2*0.030^2 + energy*0.35^2) + \
337# (abs(eta) > 2.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.050^2 + energy*0.45^2)}
338
339
340
341}
342####################
343# Energy flow merger
344####################
345
346module Merger EFlowMerger {
347# add InputArray InputArray
348# exchange UK: 25.7.16 change back to all eflowtracks!
349# add InputArray ImpactParameterSmearing/tracks
350 add InputArray Calorimeter/eflowTracks
351 add InputArray Calorimeter/eflowPhotons
352 add InputArray Calorimeter/eflowNeutralHadrons
353 set OutputArray eflow
354}
355######################
356# EFlowFilter
357######################
358
359module PdgCodeFilter EFlowFilter {
360 set InputArray EFlowMerger/eflow
361 set OutputArray eflow
362
363 add PdgCode {11}
364 add PdgCode {-11}
365 add PdgCode {13}
366 add PdgCode {-13}
367}
368
369
370###################
371# Photon efficiency
372###################
373
374module Efficiency PhotonEfficiency {
375# set InputArray Calorimeter/eflowPhotons
376# 21.7.16 use photons
377 set InputArray Calorimeter/photons
378 set OutputArray photons
379
380 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
381
382 # efficiency formula for photons - same as for electrons
383# expand and set flat
384 # set to 1
385 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( (eta <=5.4 && eta >= -5.1) ) * (1.0) + \
386 ((eta > 5.4 || eta < -5.1)) * (0.00)}
387
388}
389
390##################
391# Photon isolation
392##################
393
394module Isolation PhotonIsolation {
395 set CandidateInputArray PhotonEfficiency/photons
396# set IsolationInputArray EFlowMerger/eflow
397 set IsolationInputArray EFlowFilter/eflow
398
399 set OutputArray photons
400
401 set DeltaRMax 0.3
402
403 set PTMin 0.5
404
405 set PTRatioMax 0.1
406}
407#################
408# Electron filter
409#################
410
411module PdgCodeFilter ElectronFilter {
412 set InputArray Calorimeter/eflowTracks
413 set OutputArray electrons
414 set Invert true
415 add PdgCode {11}
416 add PdgCode {-11}
417}
418
419
420
421#################
422# Muon filter
423#################
424
425module PdgCodeFilter MuonFilter {
426 set InputArray Calorimeter/eflowTracks
427 set OutputArray muons
428 set Invert true
429 add PdgCode {13}
430 add PdgCode {-13}
431}
432
433
434#####################
435# Electron efficiency
436#####################
437
438module Efficiency ElectronEfficiency {
439# set InputArray ElectronEnergySmearing/electrons
440 set InputArray ElectronFilter/electrons
441 set OutputArray electrons
442
443 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
444
445 # efficiency formula for electrons
446# set to 1 for eta<4.7 1 degree
447 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *( (eta <=5.4 && eta >= -5.1) ) * (1.0) + \
448 ((eta > 5.4 || eta < -5.1)) * (0.00)}
449# set EfficiencyFormula { (abs(eta) <= 4.7) * (1.0) + \
450# (abs(eta) > 4.7) * (0.00)}
451}
452
453####################
454# Electron isolation
455####################
456
457module Isolation ElectronIsolation {
458 set CandidateInputArray ElectronEfficiency/electrons
459# set IsolationInputArray EFlowMerger/eflow
460 set IsolationInputArray EFlowFilter/eflow
461
462 set OutputArray electrons
463
464 set DeltaRMax 0.3
465
466 set PTMin 0.5
467
468 set PTRatioMax 0.1
469}
470
471#################
472# Muon efficiency
473#################
474
475module Efficiency MuonEfficiency {
476 set InputArray MuonMomentumSmearing/muons
477 set OutputArray muons
478
479 # set EfficiencyFormula {efficiency as a function of eta and pt}
480
481 # efficiency formula for muons
482# set to 1 for eta<4
483 set EfficiencyFormula { (pt <= 0.5) * (0.00) + (pt > 0.5) *(abs(eta) <= 4) * (1.0) + \
484 (abs(eta) > 4) * (0.00)}
485}
486
487################
488# Muon isolation
489################
490
491module Isolation MuonIsolation {
492 set CandidateInputArray MuonEfficiency/muons
493# set IsolationInputArray EFlowMerger/eflow
494 set IsolationInputArray EFlowFilter/eflow
495
496 set OutputArray muons
497
498 set DeltaRMax 0.3
499
500 set PTMin 0.5
501
502 set PTRatioMax 0.1
503}
504
505###################
506# Missing ET merger
507###################
508
509module Merger MissingET {
510# add InputArray InputArray
511 add InputArray EFlowMerger/eflow
512 set MomentumOutputArray momentum
513}
514#########################
515# Gen Missing ET merger
516########################
517
518module Merger GenMissingET {
519
520# add InputArray InputArray
521 add InputArray NeutrinoFilter/filteredParticles
522 set MomentumOutputArray momentum
523}
524
525
526##################
527# Scalar HT merger
528##################
529
530module Merger ScalarHT {
531# add InputArray InputArray
532 add InputArray UniqueObjectFinder/jets
533 add InputArray UniqueObjectFinder/electrons
534 add InputArray UniqueObjectFinder/photons
535 add InputArray UniqueObjectFinder/muons
536 set EnergyOutputArray energy
537}
538#################
539# Neutrino Filter
540#################
541
542module PdgCodeFilter NeutrinoFilter {
543
544 set InputArray Delphes/stableParticles
545 set OutputArray filteredParticles
546
547 set PTMin 0.0
548
549 add PdgCode {12}
550 add PdgCode {14}
551 add PdgCode {16}
552 add PdgCode {-12}
553 add PdgCode {-14}
554 add PdgCode {-16}
555
556}
557
558
559#####################
560# MC truth jet finder
561#####################
562
563module FastJetFinder GenJetFinder {
564 set InputArray Delphes/stableParticles
565
566 set OutputArray jets
567
568 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
569 set JetAlgorithm 6
570 set ParameterR 0.4
571
572 set JetPTMin 3.0
573}
574
575############
576# Jet finder
577############
578
579module FastJetFinder FastJetFinder {
580# ATLAS uses towers
581# set InputArray Calorimeter/towers
582# CMS and we use eflow
583 set InputArray EFlowMerger/eflow
584
585 set OutputArray jets
586
587 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
588 set JetAlgorithm 6
589 set ParameterR 0.4
590
591 set JetPTMin 3.0
592}
593
594##################
595# Jet Energy Scale
596##################
597
598module EnergyScale JetEnergyScale {
599 set InputArray FastJetFinder/jets
600 set OutputArray jets
601
602 # scale formula for jets
603 set ScaleFormula {1.00}
604}
605########################
606# GenJet Flavor Association
607########################
608
609module JetFlavorAssociation GenJetFlavorAssociation {
610 set PartonInputArray Delphes/partons
611 set ParticleInputArray Delphes/allParticles
612 set ParticleLHEFInputArray Delphes/allParticlesLHEF
613 set JetInputArray GenJetFinder/jets
614# 18.4.17 enhance pt cut and eta range -> eta cut in tagging modules
615# change deltaR of 0.4 back to FCC-hh delta_R cut of 0.5
616 set DeltaR 0.5
617# set PartonPTMin 0.5
618 set PartonPTMin 1.0
619 # set PartonEtaMax 4.0
620 set PartonEtaMax 6.0
621
622}
623########################
624# Jet Flavor Association
625########################
626
627module JetFlavorAssociation JetFlavorAssociation {
628
629 set PartonInputArray Delphes/partons
630 set ParticleInputArray Delphes/allParticles
631 set ParticleLHEFInputArray Delphes/allParticlesLHEF
632 set JetInputArray JetEnergyScale/jets
633
634# 18.4.17 enhance pt cut and eta range -> eta cut in tagging modules
635 set DeltaR 0.5
636# set PartonPTMin 0.5
637 set PartonPTMin 1.0
638 # set PartonEtaMax 4.0
639 set PartonEtaMax 6.0
640
641}
642
643#################
644# Gen b-tagging
645#################
646module BTagging GenBTagging {
647# set PartonInputArray Delphes/partons
648 set JetInputArray GenJetFinder/jets
649
650# 18.4.17 taken from FCC hh but lower pt cut from 10 to 4
651 set BitNumber 0
652
653 add EfficiencyFormula {0} {
654
655 (pt <= 4.0) * (0.00) +
656 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.001) + \
657 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
658 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
659 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.00075) + \
660 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
661 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
662 (abs(eta) > 4.0) * (0.00)}
663
664 add EfficiencyFormula {4} {
665
666 (pt <= 4.0) * (0.00) +
667 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.04) + \
668 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
669 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
670 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.03) + \
671 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
672 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
673 (abs(eta) > 4.0) * (0.00)}
674
675 add EfficiencyFormula {5} {
676
677 (pt <= 4.0) * (0.00) +
678 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.85) +
679 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
680 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
681 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.64) +
682 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
683 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
684 (abs(eta) >= 4.0) * (0.00)}
685
686}
687
688#################
689# Rec b-tagging
690#################
691module BTagging BTagging {
692# set PartonInputArray Delphes/partons
693 set JetInputArray JetEnergyScale/jets
694
695# 18.4.17 taken from FCC hh but lower pt cut from 10 to 4
696 set BitNumber 0
697
698 add EfficiencyFormula {0} {
699
700 (pt <= 4.0) * (0.00) +
701 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.001) + \
702 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
703 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
704 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.00075) + \
705 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
706 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
707 (abs(eta) > 4.0) * (0.00)}
708
709 add EfficiencyFormula {4} {
710
711 (pt <= 4.0) * (0.00) +
712 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.04) + \
713 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
714 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
715 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.03) + \
716 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
717 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
718 (abs(eta) > 4.0) * (0.00)}
719
720 add EfficiencyFormula {5} {
721
722 (pt <= 4.0) * (0.00) +
723 (abs(eta) < 2.5) * (pt > 4.0 && pt < 500) * (0.85) +
724 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
725 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
726 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 4.0 && pt < 500) * (0.64) +
727 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
728 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
729 (abs(eta) >= 4.0) * (0.00)}
730
731}
732#################
733# Gen c-tagging
734#################
735
736module BTagging GenCTagging {
737# set PartonInputArray Delphes/partons
738 set JetInputArray GenJetFinder/jets
739
740
741 set BitNumber 1
742
743 add EfficiencyFormula {0} {
744
745 (pt <= 4.0) * (0.00) +
746 (abs(eta) < 4.0) * (pt > 4.0) * (0.01) + \
747 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
748
749 add EfficiencyFormula {4} {
750
751 (pt <= 4.0) * (0.00) +
752 (abs(eta) < 4.0) * (pt > 4.0) * (0.10) + \
753 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
754
755 add EfficiencyFormula {5} {
756
757 (pt <= 4.0) * (0.00) +
758 (abs(eta) < 4.0) * (pt > 4.0) * (0.25) + \
759 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
760
761}
762
763#################
764# Rec c-tagging
765#################
766
767module BTagging CTagging {
768# set PartonInputArray Delphes/partons
769 set JetInputArray JetEnergyScale/jets
770
771
772 set BitNumber 1
773
774 add EfficiencyFormula {0} {
775
776 (pt <= 4.0) * (0.00) +
777 (abs(eta) < 4.0) * (pt > 4.0) * (0.01) + \
778 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
779
780 add EfficiencyFormula {4} {
781
782 (pt <= 4.0) * (0.00) +
783 (abs(eta) < 4.0) * (pt > 4.0) * (0.10) + \
784 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
785
786 add EfficiencyFormula {5} {
787
788 (pt <= 4.0) * (0.00) +
789 (abs(eta) < 4.0) * (pt > 4.0) * (0.25) + \
790 (abs(eta) > 4.0) * (pt > 4.0) * (0.00)}
791
792}
793module TauTagging TauTagging {
794 set ParticleInputArray Delphes/allParticles
795 set PartonInputArray Delphes/partons
796 set JetInputArray JetEnergyScale/jets
797
798 set DeltaR 0.5
799
800 set TauPTMin 0.5
801
802 set TauEtaMax 6.0
803
804 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
805
806 # default efficiency formula (misidentification rate)
807 add EfficiencyFormula {0} {0.001}
808 # efficiency formula for tau-jets
809 add EfficiencyFormula {15} { (pt <= 0.5) * (0.00) + (pt > 0.5) * (0.4) }
810}
811
812module TauTagging GenTauTagging {
813 set ParticleInputArray Delphes/allParticles
814 set PartonInputArray Delphes/partons
815 set JetInputArray GenJetFinder/jets
816
817 set DeltaR 0.4
818
819 set TauPTMin 0.5
820
821 set TauEtaMax 6.0
822
823 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
824
825 # default efficiency formula (misidentification rate)
826 add EfficiencyFormula {0} {0.001}
827 # efficiency formula for tau-jets
828 add EfficiencyFormula {15} { (pt <= 0.5) * (0.00) + (pt > 0.5) * (0.4) }
829}
830
831
832
833#####################################################
834# Find uniquely identified photons/electrons/tau/jets
835#####################################################
836
837module UniqueObjectFinder UniqueObjectFinder {
838# earlier arrays take precedence over later ones
839# add InputArray InputArray OutputArray
840 add InputArray PhotonIsolation/photons photons
841 add InputArray ElectronIsolation/electrons electrons
842 add InputArray MuonIsolation/muons muons
843 add InputArray JetEnergyScale/jets jets
844}
845
846##################
847# ROOT tree writer
848##################
849
850# tracks, towers and eflow objects are not stored by default in the output.
851# if needed (for jet constituent or other studies), uncomment the relevant
852# "add Branch ..." lines.
853
854module TreeWriter TreeWriter {
855# add Branch InputArray BranchName BranchClass
856 add Branch Delphes/allParticles Particle GenParticle
857
858 add Branch TrackMerger/tracks Track Track
859 #add Branch Calorimeter/towers Tower Tower
860
861 #add Branch Calorimeter/eflowTracks EFlowTrack Track
862 #add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
863 #add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
864
865 add Branch GenJetFinder/jets GenJet Jet
866 add Branch UniqueObjectFinder/jets Jet Jet
867 add Branch UniqueObjectFinder/electrons Electron Electron
868 add Branch UniqueObjectFinder/photons Photon Photon
869 add Branch UniqueObjectFinder/muons Muon Muon
870 add Branch MissingET/momentum MissingET MissingET
871 add Branch ScalarHT/energy ScalarHT ScalarHT
872}
Note: See TracBrowser for help on using the repository browser.