Changeset 1254 in svn for trunk/examples/delphes_card_CMS.tcl
- Timestamp:
- Aug 15, 2013, 1:32:50 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/delphes_card_CMS.tcl
r1099 r1254 32 32 FastJetFinder 33 33 34 ConstituentFilter34 JetEnergyScale 35 35 36 36 BTagging … … 135 135 136 136 # resolution formula for charged hadrons 137 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \ 138 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \ 139 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \ 140 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \ 141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \ 142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \ 143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \ 144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 137 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0e1) * (0.20) + \ 138 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \ 139 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.20) + \ 140 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.20) + \ 141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \ 142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.20)} 145 143 } 146 144 … … 155 153 # set ResolutionFormula {resolution formula as a function of eta and energy} 156 154 157 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \ 158 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \ 159 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \ 155 # resolution formula for electrons 156 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (0.0225)*energy + \ 157 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 158 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 160 159 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 161 160 … … 173 172 174 173 # resolution formula for muons 175 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \ 176 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \ 177 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \ 178 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \ 179 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \ 180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.03) + \ 181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \ 182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)} 174 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \ 175 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \ 176 (abs(eta) <= 0.5) * (pt > 1.0e2) * (0.03) + \ 177 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \ 178 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \ 179 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2) * (0.04) + \ 180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \ 181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \ 182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.05)} 183 183 } 184 184 … … 262 262 263 263 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 264 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.00 5^2 + energy*0.05^2 + 0.25^2)+ \264 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \ 265 265 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 266 266 267 267 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 268 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \ 269 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)} 268 set HCalResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.055^2 + energy*1.20^2 + 3.50^2) + \ 269 (abs(eta) > 1.5 && abs(eta) <= 3) * sqrt(energy^2*0.065^2 + energy*2.00^2 + 3.50^2) + \ 270 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*3.13^2)} 270 271 } 271 272 … … 403 404 add InputArray UniqueObjectFinder/electrons 404 405 add InputArray UniqueObjectFinder/photons 405 add InputArray MuonIsolation/muons406 add InputArray UniqueObjectFinder/muons 406 407 set EnergyOutputArray energy 407 408 } … … 418 419 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 419 420 set JetAlgorithm 6 420 set ParameterR 0.7 421 422 set ConeRadius 0.5 423 set SeedThreshold 1.0 424 set ConeAreaFraction 1.0 425 set AdjacencyCut 2.0 426 set OverlapThreshold 0.75 427 428 set MaxIterations 100 429 set MaxPairSize 2 430 set Iratch 1 421 set ParameterR 0.5 431 422 432 423 set JetPTMin 20.0 … … 445 436 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 446 437 set JetAlgorithm 6 447 set ParameterR 0.7 448 449 set ConeRadius 0.5 450 set SeedThreshold 1.0 451 set ConeAreaFraction 1.0 452 set AdjacencyCut 2.0 453 set OverlapThreshold 0.75 454 455 set MaxIterations 100 456 set MaxPairSize 2 457 set Iratch 1 438 set ParameterR 0.5 458 439 459 440 set JetPTMin 20.0 460 441 } 461 442 462 #################### 463 # Constituent filter 464 #################### 465 466 module ConstituentFilter ConstituentFilter { 467 468 # add JetInputArray InputArray 469 add JetInputArray GenJetFinder/jets 470 add JetInputArray FastJetFinder/jets 471 472 # add ConstituentInputArray InputArray OutputArray 473 add ConstituentInputArray Delphes/stableParticles stableParticles 474 add ConstituentInputArray Calorimeter/eflowTracks eflowTracks 475 add ConstituentInputArray Calorimeter/eflowTowers eflowTowers 476 add ConstituentInputArray MuonMomentumSmearing/muons muons 443 ################## 444 # Jet Energy Scale 445 ################## 446 447 module EnergyScale JetEnergyScale { 448 set InputArray FastJetFinder/jets 449 set OutputArray jets 450 451 # scale formula for jets 452 set ScaleFormula {1.08} 477 453 } 478 454 … … 483 459 module BTagging BTagging { 484 460 set PartonInputArray Delphes/partons 485 set JetInputArray FastJetFinder/jets461 set JetInputArray JetEnergyScale/jets 486 462 487 463 set BitNumber 0 … … 542 518 add InputArray PhotonIsolation/photons photons 543 519 add InputArray ElectronIsolation/electrons electrons 544 add InputArray FastJetFinder/jets jets 520 add InputArray MuonIsolation/muons muons 521 add InputArray JetEnergyScale/jets jets 545 522 } 546 523 … … 554 531 add Branch TrackMerger/tracks Track Track 555 532 add Branch Calorimeter/towers Tower Tower 556 add Branch C onstituentFilter/eflowTracks EFlowTrack Track557 add Branch C onstituentFilter/eflowTowers EFlowTower Tower558 add Branch ConstituentFilter/muons EFlowMuon Muon533 add Branch Calorimeter/eflowTracks EFlowTrack Track 534 add Branch Calorimeter/eflowTowers EFlowTower Tower 535 add Branch MuonMomentumSmearing/muons EFlowMuon Muon 559 536 add Branch GenJetFinder/jets GenJet Jet 560 537 add Branch UniqueObjectFinder/jets Jet Jet 561 538 add Branch UniqueObjectFinder/electrons Electron Electron 562 539 add Branch UniqueObjectFinder/photons Photon Photon 563 add Branch MuonIsolation/muons Muon Muon540 add Branch UniqueObjectFinder/muons Muon Muon 564 541 add Branch MissingET/momentum MissingET MissingET 565 542 add Branch ScalarHT/energy ScalarHT ScalarHT
Note:
See TracChangeset
for help on using the changeset viewer.