Changeset a2983ec in git for cards/delphes_card_ATLAS.tcl
- Timestamp:
- Dec 22, 2014, 11:38:33 AM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 8c72f09
- Parents:
- 4e09c3a
- File:
-
- 1 edited
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 }
Note:
See TracChangeset
for help on using the changeset viewer.