Fork me on GitHub

source: git/cards/delphes_card_CMS_NoFastJet.tcl@ 5198eed

ImprovedOutputFile Timing dual_readout llp 3.3.3pre16
Last change on this file since 5198eed was 934d037, checked in by Pavel Demin <pavel.demin@…>, 9 years ago

replace ElectronEnergySmearing with ElectronMomentumSmearing and add ElectronFilter

  • Property mode set to 100644
File size: 9.5 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 Calorimeter
18 EFlowMerger
19}
20
21#################################
22# Propagate particles in cylinder
23#################################
24
25module ParticlePropagator ParticlePropagator {
26 set InputArray Delphes/stableParticles
27
28 set OutputArray stableParticles
29 set ChargedHadronOutputArray chargedHadrons
30 set ElectronOutputArray electrons
31 set MuonOutputArray muons
32
33 # radius of the magnetic field coverage, in m
34 set Radius 1.29
35 # half-length of the magnetic field coverage, in m
36 set HalfLength 3.00
37
38 # magnetic field
39 set Bz 3.8
40}
41
42####################################
43# Charged hadron tracking efficiency
44####################################
45
46module Efficiency ChargedHadronTrackingEfficiency {
47 set InputArray ParticlePropagator/chargedHadrons
48 set OutputArray chargedHadrons
49
50 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
51
52 # tracking efficiency formula for charged hadrons
53 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
54 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
55 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
56 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
57 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
58 (abs(eta) > 2.5) * (0.00)}
59}
60
61##############################
62# Electron tracking efficiency
63##############################
64
65module Efficiency ElectronTrackingEfficiency {
66 set InputArray ParticlePropagator/electrons
67 set OutputArray electrons
68
69 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
70
71 # tracking efficiency formula for electrons
72 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
73 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
74 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
75 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
76 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
77 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
78 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
79 (abs(eta) > 2.5) * (0.00)}
80}
81
82##########################
83# Muon tracking efficiency
84##########################
85
86module Efficiency MuonTrackingEfficiency {
87 set InputArray ParticlePropagator/muons
88 set OutputArray muons
89
90 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
91
92 # tracking efficiency formula for muons
93 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
94 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
95 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
96 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
97 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
98 (abs(eta) > 2.5) * (0.00)}
99}
100
101########################################
102# Momentum resolution for charged tracks
103########################################
104
105module MomentumSmearing ChargedHadronMomentumSmearing {
106 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
107 set OutputArray chargedHadrons
108
109 # set ResolutionFormula {resolution formula as a function of eta and pt}
110
111 # resolution formula for charged hadrons
112 # based on arXiv:1405.6569
113 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
114 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
115 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
116}
117
118###################################
119# Momentum resolution for electrons
120###################################
121
122module MomentumSmearing ElectronMomentumSmearing {
123 set InputArray ElectronTrackingEfficiency/electrons
124 set OutputArray electrons
125
126 # set ResolutionFormula {resolution formula as a function of eta and energy}
127
128 # resolution formula for electrons
129 # based on arXiv:1405.6569
130 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
131 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
132 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
133}
134
135###############################
136# Momentum resolution for muons
137###############################
138
139module MomentumSmearing MuonMomentumSmearing {
140 set InputArray MuonTrackingEfficiency/muons
141 set OutputArray muons
142
143 # set ResolutionFormula {resolution formula as a function of eta and pt}
144
145 # resolution formula for muons
146 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
147 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
148 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
149}
150
151##############
152# Track merger
153##############
154
155module Merger TrackMerger {
156# add InputArray InputArray
157 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
158 add InputArray ElectronMomentumSmearing/electrons
159 add InputArray MuonMomentumSmearing/muons
160 set OutputArray tracks
161}
162
163#############
164# Calorimeter
165#############
166
167module Calorimeter Calorimeter {
168 set ParticleInputArray ParticlePropagator/stableParticles
169 set TrackInputArray TrackMerger/tracks
170
171 set TowerOutputArray towers
172 set PhotonOutputArray photons
173
174 set EFlowTrackOutputArray eflowTracks
175 set EFlowPhotonOutputArray eflowPhotons
176 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
177
178 set ECalEnergyMin 0.5
179 set HCalEnergyMin 1.0
180
181 set ECalEnergySignificanceMin 1.0
182 set HCalEnergySignificanceMin 1.0
183
184 set SmearTowerCenter true
185
186 set pi [expr {acos(-1)}]
187
188 # lists of the edges of each tower in eta and phi
189 # each list starts with the lower edge of the first tower
190 # the list ends with the higher edged of the last tower
191
192 # 5 degrees towers
193 set PhiBins {}
194 for {set i -36} {$i <= 36} {incr i} {
195 add PhiBins [expr {$i * $pi/36.0}]
196 }
197 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
198 add EtaPhiBins $eta $PhiBins
199 }
200
201 # 10 degrees towers
202 set PhiBins {}
203 for {set i -18} {$i <= 18} {incr i} {
204 add PhiBins [expr {$i * $pi/18.0}]
205 }
206 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
207 add EtaPhiBins $eta $PhiBins
208 }
209
210 # 20 degrees towers
211 set PhiBins {}
212 for {set i -9} {$i <= 9} {incr i} {
213 add PhiBins [expr {$i * $pi/9.0}]
214 }
215 foreach eta {-5 -4.7 -4.525 4.7 5} {
216 add EtaPhiBins $eta $PhiBins
217 }
218
219 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
220 add EnergyFraction {0} {0.0 1.0}
221 # energy fractions for e, gamma and pi0
222 add EnergyFraction {11} {1.0 0.0}
223 add EnergyFraction {22} {1.0 0.0}
224 add EnergyFraction {111} {1.0 0.0}
225 # energy fractions for muon, neutrinos and neutralinos
226 add EnergyFraction {12} {0.0 0.0}
227 add EnergyFraction {13} {0.0 0.0}
228 add EnergyFraction {14} {0.0 0.0}
229 add EnergyFraction {16} {0.0 0.0}
230 add EnergyFraction {1000022} {0.0 0.0}
231 add EnergyFraction {1000023} {0.0 0.0}
232 add EnergyFraction {1000025} {0.0 0.0}
233 add EnergyFraction {1000035} {0.0 0.0}
234 add EnergyFraction {1000045} {0.0 0.0}
235 # energy fractions for K0short and Lambda
236 add EnergyFraction {310} {0.3 0.7}
237 add EnergyFraction {3122} {0.3 0.7}
238
239 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
240 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
241 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
242
243 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
244 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
245 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
246}
247
248####################
249# Energy flow merger
250####################
251
252module Merger EFlowMerger {
253# add InputArray InputArray
254 add InputArray Calorimeter/eflowTracks
255 add InputArray Calorimeter/eflowPhotons
256 add InputArray Calorimeter/eflowNeutralHadrons
257 set OutputArray eflow
258}
Note: See TracBrowser for help on using the repository browser.