- Timestamp:
- Aug 15, 2013, 2:11:25 PM (11 years ago)
- File:
-
- 1 edited
-
trunk/examples/delphes_card_CMS_PileUp.tcl (modified) (14 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/delphes_card_CMS_PileUp.tcl
r1139 r1256 4 4 5 5 set ExecutionPath { 6 7 Weighter8 6 9 7 PileUpMerger … … 29 27 JetPileUpSubtractor 30 28 31 ConstituentFilter29 JetEnergyScale 32 30 33 31 PhotonEfficiency … … 50 48 51 49 TreeWriter 52 }53 54 ##########55 # Weighter56 ##########57 58 module Weighter Weighter {59 set InputArray Delphes/allParticles60 61 set OutputArray weight62 63 # add Weight {PID1 PID2} {weight}64 # default weight65 add Weight {0} {1.0}66 67 add Weight {23} {0.5}68 add Weight {24} {0.5}69 add Weight {-24} {0.5}70 add Weight {23 24} {0.25}71 add Weight {23 -24} {0.25}72 50 } 73 51 … … 181 159 182 160 # resolution formula for charged hadrons 183 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \ 184 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \ 185 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \ 186 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \ 187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \ 188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \ 189 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \ 190 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 161 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0e1) * (0.20) + \ 162 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \ 163 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.20) + \ 164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.20) + \ 165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \ 166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.20)} 191 167 } 192 168 … … 201 177 # set ResolutionFormula {resolution formula as a function of eta and energy} 202 178 203 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \ 204 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \ 205 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \ 179 # resolution formula for electrons 180 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (0.0225)*energy + \ 181 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 182 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 206 183 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 207 184 … … 219 196 220 197 # resolution formula for muons 221 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \ 222 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \ 223 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \ 224 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \ 225 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \ 226 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.03) + \ 227 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \ 228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 198 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \ 199 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \ 200 (abs(eta) <= 0.5) * (pt > 1.0e2) * (0.03) + \ 201 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \ 202 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \ 203 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2) * (0.04) + \ 204 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \ 205 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \ 206 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.05)} 229 207 } 230 208 … … 308 286 309 287 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 310 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.00 5^2 + energy*0.05^2 + 0.25^2)+ \288 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 311 289 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 312 290 313 291 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 314 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \ 315 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)} 292 set HCalResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.055^2 + energy*1.20^2 + 3.50^2) + \ 293 (abs(eta) > 1.5 && abs(eta) <= 3) * sqrt(energy^2*0.065^2 + energy*2.00^2 + 3.50^2) + \ 294 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*3.13^2)} 316 295 } 317 296 … … 402 381 } 403 382 404 ####################405 # Constituent filter406 ####################407 408 module ConstituentFilter ConstituentFilter {409 410 # add JetInputArray InputArray411 add JetInputArray GenJetFinder/jets412 add JetInputArray FastJetFinder/jets413 414 # add ConstituentInputArray InputArray OutputArray415 add ConstituentInputArray Delphes/stableParticles stableParticles416 add ConstituentInputArray TrackPileUpSubtractor/eflowTracks eflowTracks417 add ConstituentInputArray Calorimeter/eflowTowers eflowTowers418 add ConstituentInputArray MuonMomentumSmearing/muons muons419 }420 421 383 ########################### 422 384 # Jet Pile-Up Subtraction … … 430 392 431 393 set JetPTMin 20.0 394 } 395 396 ################## 397 # Jet Energy Scale 398 ################## 399 400 module EnergyScale JetEnergyScale { 401 set InputArray JetPileUpSubtractor/jets 402 set OutputArray jets 403 404 # scale formula for jets 405 set ScaleFormula {1.0} 432 406 } 433 407 … … 556 530 add InputArray UniqueObjectFinder/electrons 557 531 add InputArray UniqueObjectFinder/photons 558 add InputArray MuonIsolation/muons532 add InputArray UniqueObjectFinder/muons 559 533 set EnergyOutputArray energy 560 534 } … … 567 541 module BTagging BTagging { 568 542 set PartonInputArray Delphes/partons 569 set JetInputArray Jet PileUpSubtractor/jets543 set JetInputArray JetEnergyScale/jets 570 544 571 545 set BitNumber 0 … … 601 575 set ParticleInputArray Delphes/allParticles 602 576 set PartonInputArray Delphes/partons 603 set JetInputArray Jet PileUpSubtractor/jets577 set JetInputArray JetEnergyScale/jets 604 578 605 579 set DeltaR 0.5 … … 626 600 add InputArray PhotonIsolation/photons photons 627 601 add InputArray ElectronIsolation/electrons electrons 628 add InputArray JetPileUpSubtractor/jets jets 602 add InputArray MuonIsolation/muons muons 603 add InputArray JetEnergyScale/jets jets 629 604 } 630 605 … … 645 620 add Branch UniqueObjectFinder/electrons Electron Electron 646 621 add Branch UniqueObjectFinder/photons Photon Photon 647 add Branch MuonIsolation/muons Muon Muon622 add Branch UniqueObjectFinder/muons Muon Muon 648 623 add Branch MissingET/momentum MissingET MissingET 649 624 add Branch ScalarHT/energy ScalarHT ScalarHT 650 625 add Branch Rho/rho Rho Rho 651 add Branch Weighter/weight Weight Weight 652 } 653 626 } 627
Note:
See TracChangeset
for help on using the changeset viewer.
