Changeset a2983ec in git
- Timestamp:
- Dec 22, 2014, 11:38:33 AM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 8c72f09
- Parents:
- 4e09c3a
- Location:
- cards
- Files:
-
- 7 edited
-
delphes_card_ATLAS.tcl (modified) (15 diffs)
-
delphes_card_ATLAS_PileUp.tcl (modified) (15 diffs)
-
delphes_card_CMS.tcl (modified) (14 diffs)
-
delphes_card_CMS_NoFastJet.tcl (modified) (9 diffs)
-
delphes_card_CMS_PileUp.tcl (modified) (14 diffs)
-
delphes_card_FCC_basic.tcl (modified) (22 diffs)
-
delphes_card_LHCb.tcl (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_ATLAS.tcl
r4e09c3a ra2983ec 77 77 78 78 # tracking efficiency formula for charged hadrons 79 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \80 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \81 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \83 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \79 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 80 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 81 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + 83 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + 84 84 (abs(eta) > 2.5) * (0.00)} 85 85 } … … 96 96 97 97 # tracking efficiency formula for electrons 98 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \100 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \101 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \98 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + 100 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + 101 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + 102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + 103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + 104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + 105 105 (abs(eta) > 2.5) * (0.00)} 106 106 } … … 117 117 118 118 # tracking efficiency formula for muons 119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + 124 124 (abs(eta) > 2.5) * (0.00)} 125 125 } … … 136 136 137 137 # resolution formula for charged hadrons 138 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \139 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \140 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \141 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \138 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + 139 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + 140 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + 141 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + 142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + 144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + 145 145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 146 146 } … … 156 156 # set ResolutionFormula {resolution formula as a function of eta and energy} 157 157 158 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \159 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \160 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \158 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + 159 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + 160 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + 161 161 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 162 162 … … 174 174 175 175 # resolution formula for muons 176 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \177 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \178 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \179 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \176 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 177 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + 178 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + 179 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + 180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + 181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + 182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + 183 183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)} 184 184 } … … 207 207 set PhotonOutputArray photons 208 208 209 set ECalEnergyMin 0.50210 set HCalEnergyMin 1.00211 212 set ECalEnergySignificanceMin 1.0213 set HCalEnergySignificanceMin 1.0214 215 209 set EFlowTrackOutputArray eflowTracks 216 210 set EFlowPhotonOutputArray eflowPhotons 217 211 set EFlowNeutralHadronOutputArray eflowNeutralHadrons 212 213 set ECalEnergyMin 0.5 214 set HCalEnergyMin 1.0 215 216 set ECalEnergySignificanceMin 1.0 217 set HCalEnergySignificanceMin 1.0 218 219 set SmearTowerCenter true 218 220 219 221 set pi [expr {acos(-1)}] … … 265 267 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf 266 268 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf 267 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \269 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + 268 270 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)} 269 271 … … 271 273 # http://arxiv.org/pdf/hep-ex/0004009v1 272 274 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf 273 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \274 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \275 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + 276 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + 275 277 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)} 276 278 } … … 299 301 300 302 # efficiency formula for photons 301 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \302 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \303 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \303 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 304 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 305 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 304 306 (abs(eta) > 2.5) * (0.00)} 305 307 } … … 333 335 334 336 # efficiency formula for electrons 335 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \336 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \337 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \337 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 338 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 339 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 338 340 (abs(eta) > 2.5) * (0.00)} 339 341 } … … 367 369 368 370 # efficiency formula for muons 369 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \370 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \371 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \371 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 372 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 373 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + 372 374 (abs(eta) > 2.7) * (0.00)} 373 375 } … … 419 421 420 422 module PdgCodeFilter NeutrinoFilter { 421 423 422 424 set InputArray Delphes/stableParticles 423 425 set OutputArray filteredParticles 424 426 425 427 set PTMin 0.0 426 428 427 429 add PdgCode {12} 428 430 add PdgCode {14} … … 475 477 set OutputArray jets 476 478 477 # scale formula for jets 479 # scale formula for jets 478 480 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) } 479 481 } … … 503 505 504 506 # efficiency formula for c-jets (misidentification rate) 505 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \506 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \507 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \507 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + 508 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + 509 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + 508 510 (abs(eta) > 2.5) * (0.000)} 509 511 510 512 # efficiency formula for b-jets 511 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \512 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \513 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \513 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + 514 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + 515 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + 514 516 (abs(eta) > 2.5) * (0.000)} 515 517 } -
cards/delphes_card_ATLAS_PileUp.tcl
r4e09c3a ra2983ec 79 79 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))} 80 80 81 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \82 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \83 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \81 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + 82 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + 83 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + 84 84 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)} 85 85 … … 119 119 120 120 # tracking efficiency formula for charged hadrons 121 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \122 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \123 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \124 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \121 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 122 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 123 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 124 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + 125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + 126 126 (abs(eta) > 2.5) * (0.00)} 127 127 } … … 138 138 139 139 # tracking efficiency formula for electrons 140 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \141 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \142 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \143 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \140 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 141 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + 142 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + 143 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + 144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + 145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + 146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + 147 147 (abs(eta) > 2.5) * (0.00)} 148 148 } … … 159 159 160 160 # tracking efficiency formula for muons 161 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \162 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \163 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \161 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 162 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 163 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + 166 166 (abs(eta) > 2.5) * (0.00)} 167 167 } … … 178 178 179 179 # resolution formula for charged hadrons 180 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \181 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \182 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \183 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \180 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + 181 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + 182 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + 183 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + 184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + 186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + 187 187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 188 188 } … … 198 198 # set ResolutionFormula {resolution formula as a function of eta and energy} 199 199 200 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \201 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \202 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \200 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + 201 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + 202 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + 203 203 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 204 204 … … 216 216 217 217 # resolution formula for muons 218 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \219 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \220 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \221 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \222 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \223 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \224 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \218 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 219 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + 220 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + 221 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + 222 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + 223 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + 224 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + 225 225 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)} 226 226 } … … 249 249 set PhotonOutputArray photons 250 250 251 set ECalEnergyMin 0.50252 set HCalEnergyMin 1.00253 254 set ECalEnergySignificanceMin 1.0255 set HCalEnergySignificanceMin 1.0256 257 251 set EFlowTrackOutputArray eflowTracks 258 252 set EFlowPhotonOutputArray eflowPhotons 259 253 set EFlowNeutralHadronOutputArray eflowNeutralHadrons 254 255 set ECalEnergyMin 0.5 256 set HCalEnergyMin 1.0 257 258 set ECalEnergySignificanceMin 1.0 259 set HCalEnergySignificanceMin 1.0 260 261 set SmearTowerCenter true 260 262 261 263 set pi [expr {acos(-1)}] … … 307 309 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf 308 310 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf 309 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \311 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + 310 312 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)} 311 313 … … 313 315 # http://arxiv.org/pdf/hep-ex/0004009v1 314 316 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf 315 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \316 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \317 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + 318 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + 317 319 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)} 318 320 } … … 482 484 483 485 # efficiency formula for photons 484 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \485 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \486 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \486 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 487 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 488 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 487 489 (abs(eta) > 2.5) * (0.00)} 488 490 } … … 517 519 518 520 # efficiency formula for electrons 519 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \520 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \521 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \521 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 522 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 523 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 522 524 (abs(eta) > 2.5) * (0.00)} 523 525 } … … 552 554 553 555 # efficiency formula for muons 554 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \555 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \556 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \556 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 557 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 558 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + 557 559 (abs(eta) > 2.7) * (0.00)} 558 560 } … … 624 626 625 627 # efficiency formula for c-jets (misidentification rate) 626 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \627 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \628 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \628 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + 629 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + 630 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + 629 631 (abs(eta) > 2.5) * (0.000)} 630 632 631 633 # efficiency formula for b-jets 632 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \633 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \634 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \634 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + 635 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + 636 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + 635 637 (abs(eta) > 2.5) * (0.000)} 636 638 } … … 698 700 699 701 } 700 701 702 -
cards/delphes_card_CMS.tcl
r4e09c3a ra2983ec 77 77 78 78 # tracking efficiency formula for charged hadrons 79 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \80 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \81 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \83 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \79 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 80 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 81 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + 83 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + 84 84 (abs(eta) > 2.5) * (0.00)} 85 85 } … … 96 96 97 97 # tracking efficiency formula for electrons 98 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \100 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \101 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \98 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + 100 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + 101 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + 102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + 103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + 104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + 105 105 (abs(eta) > 2.5) * (0.00)} 106 106 } … … 117 117 118 118 # tracking efficiency formula for muons 119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + 124 124 (abs(eta) > 2.5) * (0.00)} 125 125 } … … 136 136 137 137 # resolution formula for charged hadrons 138 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \139 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \140 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \141 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \138 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + 139 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + 140 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + 141 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + 142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + 144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + 145 145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 146 146 } … … 157 157 158 158 # resolution formula for electrons 159 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \160 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \161 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \159 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + 160 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 161 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 162 162 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 163 163 … … 175 175 176 176 # resolution formula for muons 177 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \178 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \179 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \180 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \181 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \182 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \183 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \184 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \177 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + 178 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + 179 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + 180 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 181 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + 182 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + 183 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + 184 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + 186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + 187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + 188 188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)} 189 189 } … … 212 212 set PhotonOutputArray photons 213 213 214 set ECalEnergyMin 0.5215 set HCalEnergyMin 1.0216 217 set ECalEnergySignificanceMin 1.0218 set HCalEnergySignificanceMin 1.0219 220 214 set EFlowTrackOutputArray eflowTracks 221 215 set EFlowPhotonOutputArray eflowPhotons 222 216 set EFlowNeutralHadronOutputArray eflowNeutralHadrons 217 218 set ECalEnergyMin 0.5 219 set HCalEnergyMin 1.0 220 221 set ECalEnergySignificanceMin 1.0 222 set HCalEnergySignificanceMin 1.0 223 224 set SmearTowerCenter true 223 225 224 226 set pi [expr {acos(-1)}] … … 276 278 277 279 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 278 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \280 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 279 281 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 280 282 281 283 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 282 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \284 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + 283 285 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)} 284 286 } … … 307 309 308 310 # efficiency formula for photons 309 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \310 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \311 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \311 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 312 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 313 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 312 314 (abs(eta) > 2.5) * (0.00)} 313 315 } … … 341 343 342 344 # efficiency formula for electrons 343 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \344 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \345 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \345 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 346 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 347 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 346 348 (abs(eta) > 2.5) * (0.00)} 347 349 } … … 375 377 376 378 # efficiency formula for muons 377 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \378 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \379 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \380 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \381 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \379 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 380 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + 381 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + 382 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + 383 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + 382 384 (abs(eta) > 2.4) * (0.00)} 383 385 } … … 429 431 430 432 module PdgCodeFilter NeutrinoFilter { 431 433 432 434 set InputArray Delphes/stableParticles 433 435 set OutputArray filteredParticles 434 436 435 437 set PTMin 0.0 436 438 437 439 add PdgCode {12} 438 440 add PdgCode {14} … … 518 520 519 521 # efficiency formula for c-jets (misidentification rate) 520 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \521 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \522 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \522 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + 523 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + 524 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + 523 525 (abs(eta) > 2.5) * (0.000)} 524 526 525 527 # efficiency formula for b-jets 526 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \527 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \528 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \528 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + 529 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + 530 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + 529 531 (abs(eta) > 2.5) * (0.000)} 530 532 } … … 589 591 add Branch ScalarHT/energy ScalarHT ScalarHT 590 592 } 591 -
cards/delphes_card_CMS_NoFastJet.tcl
r4e09c3a ra2983ec 51 51 52 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) + \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 58 (abs(eta) > 2.5) * (0.00)} 59 59 } … … 70 70 71 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) + \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 79 (abs(eta) > 2.5) * (0.00)} 80 80 } … … 91 91 92 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) + \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 98 (abs(eta) > 2.5) * (0.00)} 99 99 } … … 110 110 111 111 # resolution formula for charged hadrons 112 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \113 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \114 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \115 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \117 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \118 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \112 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + 113 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + 114 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + 115 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + 116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 117 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + 118 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + 119 119 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 120 120 } … … 131 131 132 132 # resolution formula for electrons 133 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \134 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \135 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \133 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + 134 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 135 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 136 136 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 137 137 … … 149 149 150 150 # resolution formula for muons 151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \152 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \153 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \154 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \155 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \156 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \157 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \158 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \159 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \160 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \161 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + 152 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + 153 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + 154 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 155 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + 156 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + 157 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + 158 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 159 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + 160 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + 161 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + 162 162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)} 163 163 } … … 186 186 set PhotonOutputArray photons 187 187 188 set ECalEnergyMin 0.5189 set HCalEnergyMin 1.0190 191 set ECalEnergySignificanceMin 1.0192 set HCalEnergySignificanceMin 1.0193 194 188 set EFlowTrackOutputArray eflowTracks 195 189 set EFlowPhotonOutputArray eflowPhotons 196 190 set EFlowNeutralHadronOutputArray eflowNeutralHadrons 191 192 set ECalEnergyMin 0.5 193 set HCalEnergyMin 1.0 194 195 set ECalEnergySignificanceMin 1.0 196 set HCalEnergySignificanceMin 1.0 197 198 set SmearTowerCenter true 197 199 198 200 set pi [expr {acos(-1)}] … … 250 252 251 253 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 252 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \254 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 253 255 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 254 256 255 257 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 256 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \258 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + 257 259 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)} 258 260 } … … 269 271 set OutputArray eflow 270 272 } 271 272 -
cards/delphes_card_CMS_PileUp.tcl
r4e09c3a ra2983ec 80 80 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))} 81 81 82 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \83 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \84 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \82 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + 83 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + 84 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + 85 85 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)} 86 86 … … 120 120 121 121 # tracking efficiency formula for charged hadrons 122 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \123 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \124 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \126 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \122 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 123 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 124 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + 126 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + 127 127 (abs(eta) > 2.5) * (0.00)} 128 128 } … … 139 139 140 140 # tracking efficiency formula for electrons 141 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \142 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \143 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \144 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \141 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 142 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + 143 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + 144 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + 145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + 146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + 147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + 148 148 (abs(eta) > 2.5) * (0.00)} 149 149 } … … 160 160 161 161 # tracking efficiency formula for muons 162 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \163 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \164 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \162 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 163 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 164 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + 166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + 167 167 (abs(eta) > 2.5) * (0.00)} 168 168 } … … 179 179 180 180 # resolution formula for charged hadrons 181 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \182 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \183 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \184 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \181 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + 182 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + 183 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + 184 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + 185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + 186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + 187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + 188 188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 189 189 } … … 200 200 201 201 # resolution formula for electrons 202 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \203 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \204 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \202 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + 203 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 204 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 205 205 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 206 206 … … 218 218 219 219 # resolution formula for muons 220 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \221 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \222 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \223 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \224 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \225 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \226 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \227 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \229 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \220 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + 221 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + 222 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + 223 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 224 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + 225 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + 226 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + 227 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + 228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + 229 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + 230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + 231 231 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)} 232 232 } … … 255 255 set PhotonOutputArray photons 256 256 257 set ECalEnergyMin 0.5258 set HCalEnergyMin 1.0259 260 set ECalEnergySignificanceMin 1.0261 set HCalEnergySignificanceMin 1.0262 263 257 set EFlowTrackOutputArray eflowTracks 264 258 set EFlowPhotonOutputArray eflowPhotons 265 259 set EFlowNeutralHadronOutputArray eflowNeutralHadrons 260 261 set ECalEnergyMin 0.5 262 set HCalEnergyMin 1.0 263 264 set ECalEnergySignificanceMin 1.0 265 set HCalEnergySignificanceMin 1.0 266 267 set SmearTowerCenter true 266 268 267 269 set pi [expr {acos(-1)}] … … 319 321 320 322 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 321 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \323 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + 322 324 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 323 325 324 326 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 325 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \327 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + 326 328 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)} 327 329 } … … 518 520 519 521 # efficiency formula for photons 520 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \521 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \522 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \522 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 523 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 524 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 523 525 (abs(eta) > 2.5) * (0.00)} 524 526 } … … 554 556 555 557 # efficiency formula for electrons 556 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \557 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \558 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \558 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 559 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + 560 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + 559 561 (abs(eta) > 2.5) * (0.00)} 560 562 } … … 589 591 590 592 # efficiency formula for muons 591 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \592 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \593 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \594 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \595 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \593 set EfficiencyFormula { (pt <= 10.0) * (0.00) + 594 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + 595 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + 596 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + 597 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + 596 598 (abs(eta) > 2.4) * (0.00)} 597 599 } … … 665 667 666 668 # efficiency formula for c-jets (misidentification rate) 667 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \668 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \669 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \669 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + 670 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + 671 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + 670 672 (abs(eta) > 2.5) * (0.000)} 671 673 672 674 # efficiency formula for b-jets 673 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \674 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \675 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \675 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + 676 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + 677 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + 676 678 (abs(eta) > 2.5) * (0.000)} 677 679 } … … 738 740 add Branch PileUpMerger/vertices Vertex Vertex 739 741 } 740 -
cards/delphes_card_FCC_basic.tcl
r4e09c3a ra2983ec 17 17 AngularSmearing 18 18 ImpactParameterSmearing 19 19 20 20 ECal 21 21 HCal 22 22 23 23 TowerMerger 24 24 EFlowMerger 25 25 26 26 MissingET 27 27 28 28 NeutrinoFilter 29 29 GenJetFinder … … 31 31 32 32 JetEnergyScale 33 33 34 34 TrackCountingBTagging 35 35 TauTagging … … 72 72 73 73 # tracking efficiency formula for charged hadrons 74 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \75 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \76 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \77 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) + \78 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.90) + \74 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 75 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 76 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 77 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) + 78 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.90) + 79 79 (abs(eta) > 4.0) * (0.00)} 80 80 81 81 } 82 82 … … 92 92 93 93 # tracking efficiency formula for electrons 94 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \95 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \96 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \97 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + \98 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + \94 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 95 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 96 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 97 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + 98 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + 99 99 (abs(eta) > 4.0) * (0.00)} 100 100 } … … 111 111 112 112 # tracking efficiency formula for muons 113 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \114 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \115 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \116 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + \117 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + \113 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 114 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 115 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 116 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + 117 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + 118 118 (abs(eta) > 4.0) * (0.00)} 119 119 } … … 130 130 131 131 # resolution formula for charged hadrons 132 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*2.e-5) + \132 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*2.e-5) + 133 133 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)} 134 134 … … 147 147 148 148 # resolution formula for electrons 149 set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) + \150 (abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) + \149 set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) + 150 (abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) + 151 151 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)} 152 152 … … 164 164 165 165 # resolution formula for muons 166 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*5.e-6) + \166 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*5.e-6) + 167 167 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)} 168 168 … … 210 210 211 211 # absolute impact parameter smearing formula (in mm) as a function of pt and eta 212 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \212 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + 213 213 (pt > 5.0) * (0.005)} 214 214 … … 225 225 set TowerOutputArray ecalTowers 226 226 set EFlowTowerOutputArray eflowPhotons 227 227 228 228 set EnergyMin 0.5 229 229 set EnergySignificanceMin 1.0 230 230 231 set SmearTowerCenter true 232 231 233 set pi [expr {acos(-1)}] 232 234 … … 234 236 # each list starts with the lower edge of the first tower 235 237 # the list ends with the higher edged of the last tower 236 238 237 239 # 0.5 degree towers 238 240 set PhiBins {} … … 246 248 add EtaPhiBins $eta $PhiBins 247 249 } 248 250 249 251 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL} 250 252 251 253 add EnergyFraction {0} {0.0} 252 254 # energy fractions for e, gamma and pi0 … … 269 271 270 272 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 271 273 272 274 # This is the CMS ECAL resolution, extended up eta = 6.0 273 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) + \275 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) + 274 276 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 275 277 276 278 277 279 } 278 280 … … 287 289 set TowerOutputArray hcalTowers 288 290 set EFlowTowerOutputArray eflowNeutralHadrons 289 291 290 292 set EnergyMin 1.0 291 293 set EnergySignificanceMin 1.0 292 294 295 set SmearTowerCenter true 296 293 297 set pi [expr {acos(-1)}] 294 298 … … 296 300 # each list starts with the lower edge of the first tower 297 301 # the list ends with the higher edged of the last tower 298 299 302 303 300 304 # 5 degree towers 301 305 set PhiBins {} … … 309 313 add EtaPhiBins $eta $PhiBins 310 314 } 311 312 315 316 313 317 # default energy fractions {abs(PDG code)} {Fecal Fhcal} 314 318 add EnergyFraction {0} {1.0} … … 332 336 333 337 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 334 338 335 339 # This is the ATLAS HCAL resolution, extended up eta = 6.0 336 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \337 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \340 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + 341 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + 338 342 (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)} 339 343 … … 373 377 set MomentumOutputArray momentum 374 378 } 375 379 376 380 377 381 ################## … … 390 394 391 395 module PdgCodeFilter NeutrinoFilter { 392 396 393 397 set InputArray Delphes/stableParticles 394 398 set OutputArray filteredParticles 395 399 396 400 set PTMin 0.0 397 401 398 402 add PdgCode {12} 399 403 add PdgCode {14} … … 463 467 # maximum distance between jet and track 464 468 set DeltaR 0.3 465 466 # minimum pt of tracks 469 470 # minimum pt of tracks 467 471 set TrackPTMin 1.0 468 472 469 473 # minimum transverse impact parameter (in mm) 470 474 set TrackIPMax 2.0 471 475 472 476 # minimum ip significance for the track to be counted 473 477 set SigMin 6.5 474 478 475 479 # minimum number of tracks (high efficiency n=2, high purity n=3) 476 480 set Ntracks 3 … … 510 514 add Branch Delphes/allParticles Particle GenParticle 511 515 add Branch GenJetFinder/jets GenJet Jet 512 516 513 517 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track 514 518 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower … … 521 525 add Branch ScalarHT/energy ScalarHT ScalarHT 522 526 } 523 -
cards/delphes_card_LHCb.tcl
r4e09c3a ra2983ec 8 8 9 9 set ExecutionPath { 10 10 11 11 ParticlePropagator 12 12 13 13 ChargedHadronMomentumSmearing 14 14 ElectronEnergySmearing … … 17 17 TrackMerger 18 18 ImpactParameterSmearing 19 19 20 20 IdentificationMap 21 21 22 22 ECal 23 23 HCal … … 41 41 42 42 # radius of the magnetic field coverage, in m 43 43 44 44 set Radius 3.31 45 46 45 46 47 47 # half-length of the magnetic field coverage, in m 48 48 set HalfLength 12.0 … … 53 53 # Need to veto anything with theta > 0.269 rad -> eta = 2 54 54 # theta < 0.0135 rad -> eta = 5 55 56 # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269) 57 58 } 59 60 55 56 # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269) 57 58 } 59 60 61 61 ######################################## 62 62 # Momentum resolution for charged tracks … … 84 84 85 85 # resolution formula for electrons 86 set ResolutionFormula { (eta > 2.0 && eta <= 5.0) * (energy > 0.1 && energy <= 8.0) * (energy*0.05) + \86 set ResolutionFormula { (eta > 2.0 && eta <= 5.0) * (energy > 0.1 && energy <= 8.0) * (energy*0.05) + 87 87 (eta > 2.0 && eta <= 5.0) * (energy > 8.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)} 88 88 } … … 139 139 set InputArray ImpactParameterSmearing/tracks 140 140 set OutputArray tracks 141 141 142 142 # {PID in} {PID out} {formula} 143 143 # make sure "PID in" and "PID out" have the same charge (e.g {-13} {211} or {-321} {211}) 144 144 # {211} {-13} is equivalent to {-211} {13} (and needs to be written once only...) 145 146 147 148 149 150 151 # --- pions --- 152 153 add EfficiencyFormula {211} {211} { (eta <= 2.0) * (0.00) + \154 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \155 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \145 146 147 148 149 150 151 # --- pions --- 152 153 add EfficiencyFormula {211} {211} { (eta <= 2.0) * (0.00) + 154 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 155 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + 156 156 (eta > 5.0) * (0.00)} 157 158 add EfficiencyFormula {211} {-13} { (eta <= 2.0) * (0.00) + \159 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \160 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) + \157 158 add EfficiencyFormula {211} {-13} { (eta <= 2.0) * (0.00) + 159 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 160 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) + 161 161 (eta > 5.0) * (0.00)} 162 162 163 163 164 164 # --- kaons --- 165 166 167 add EfficiencyFormula {321} {321} { (eta <= 2.0) * (0.00) + \168 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \169 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \165 166 167 add EfficiencyFormula {321} {321} { (eta <= 2.0) * (0.00) + 168 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 169 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + 170 170 (eta > 5.0) * (0.00)} 171 172 add EfficiencyFormula {321} {-13} { (eta <= 2.0) * (0.00) + \173 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \174 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) + \171 172 add EfficiencyFormula {321} {-13} { (eta <= 2.0) * (0.00) + 173 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 174 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) + 175 175 (eta > 5.0) * (0.00)} 176 176 177 177 178 178 # --- protons --- 179 180 181 add EfficiencyFormula {2212} {2212} { (eta <= 2.0) * (0.00) + \182 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \183 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \179 180 181 add EfficiencyFormula {2212} {2212} { (eta <= 2.0) * (0.00) + 182 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 183 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + 184 184 (eta > 5.0) * (0.00)} 185 186 add EfficiencyFormula {2212} {-13} { (eta <= 2.0) * (0.00) + \187 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \188 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.002) + \185 186 add EfficiencyFormula {2212} {-13} { (eta <= 2.0) * (0.00) + 187 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 188 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.002) + 189 189 (eta > 5.0) * (0.00)} 190 190 … … 195 195 196 196 197 add EfficiencyFormula {-13} {-13} { (eta <= 2.0) * (0.00) + \198 (eta > 2.0 && eta <= 5.0) * (pt < 0.8)* (0.00) + \199 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.97) + \197 add EfficiencyFormula {-13} {-13} { (eta <= 2.0) * (0.00) + 198 (eta > 2.0 && eta <= 5.0) * (pt < 0.8)* (0.00) + 199 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.97) + 200 200 (eta > 5.0) * (0.00)} 201 201 202 202 203 203 # efficiency for other charged particles (should be always {0} {0} {formula}) 204 205 add EfficiencyFormula {0} {0} { (eta <= 2.0) * (0.00) + \206 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \207 (eta > 2.0 && eta <= 5.0) * (pt > 0.8) * (0.95) + \204 205 add EfficiencyFormula {0} {0} { (eta <= 2.0) * (0.00) + 206 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + 207 (eta > 2.0 && eta <= 5.0) * (pt > 0.8) * (0.95) + 208 208 (eta > 5.0) * (0.00)} 209 209 … … 223 223 set TowerOutputArray ecalTowers 224 224 set EFlowTowerOutputArray eflowPhotons 225 225 226 set EnergyMin 0.0 227 set EnergySignificanceMin 0.0 228 229 set SmearTowerCenter true 230 226 231 set pi [expr {acos(-1)}] 227 232 … … 241 246 add EtaPhiBins $eta $PhiBins 242 247 } 243 248 244 249 # 1.25 degree towers 245 250 set PhiBins {} … … 285 290 add EnergyFraction {3122} {0.3} 286 291 287 292 set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)} 288 293 } 289 294 … … 298 303 set TowerOutputArray hcalTowers 299 304 set EFlowTowerOutputArray eflowNeutralHadrons 300 305 306 set EnergyMin 0.0 307 set EnergySignificanceMin 0.0 308 309 set SmearTowerCenter true 310 301 311 set pi [expr {acos(-1)}] 302 312 … … 330 340 331 341 332 342 333 343 # default energy fractions {abs(PDG code)} {Fecal Fhcal} 334 344 add EnergyFraction {0} {1.0} … … 352 362 353 363 set ResolutionFormula { (eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.05^2 + energy*0.80^2)} 354 364 } 355 365 356 366 … … 360 370 361 371 # tracks, towers and eflow objects are not stored by default in the output. 362 # if needed (for jet constituent or other studies), uncomment the relevant 363 # "add Branch ..." lines. 372 # if needed (for jet constituent or other studies), uncomment the relevant 373 # "add Branch ..." lines. 364 374 365 375 module TreeWriter TreeWriter { 366 376 # add Branch InputArray BranchName BranchClass 367 377 368 378 add Branch Delphes/allParticles Particle GenParticle 369 379 … … 371 381 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower 372 382 add Branch ECal/eflowPhotons Photon Photon 373 374 } 375 383 384 } 385
Note:
See TracChangeset
for help on using the changeset viewer.
