Changeset a2983ec in git for cards/delphes_card_ATLAS_PileUp.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
-
cards/delphes_card_ATLAS_PileUp.tcl (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
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
Note:
See TracChangeset
for help on using the changeset viewer.
