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