Changeset a2983ec in git for cards/delphes_card_FCC_basic.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_FCC_basic.tcl (modified) (22 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_FCC_basic.tcl
r4e09c3a ra2983ec 17 17 AngularSmearing 18 18 ImpactParameterSmearing 19 19 20 20 ECal 21 21 HCal 22 22 23 23 TowerMerger 24 24 EFlowMerger 25 25 26 26 MissingET 27 27 28 28 NeutrinoFilter 29 29 GenJetFinder … … 31 31 32 32 JetEnergyScale 33 33 34 34 TrackCountingBTagging 35 35 TauTagging … … 72 72 73 73 # tracking efficiency formula for charged hadrons 74 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \75 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \76 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \77 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) + \78 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.90) + \74 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 75 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 76 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + 77 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) + 78 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.90) + 79 79 (abs(eta) > 4.0) * (0.00)} 80 80 81 81 } 82 82 … … 92 92 93 93 # tracking efficiency formula for electrons 94 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \95 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \96 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \97 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + \98 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + \94 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 95 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 96 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 97 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + 98 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + 99 99 (abs(eta) > 4.0) * (0.00)} 100 100 } … … 111 111 112 112 # tracking efficiency formula for muons 113 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \114 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \115 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \116 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + \117 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + \113 set EfficiencyFormula { (pt <= 0.1) * (0.00) + 114 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + 115 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + 116 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + 117 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) + 118 118 (abs(eta) > 4.0) * (0.00)} 119 119 } … … 130 130 131 131 # resolution formula for charged hadrons 132 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*2.e-5) + \132 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*2.e-5) + 133 133 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)} 134 134 … … 147 147 148 148 # resolution formula for electrons 149 set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) + \150 (abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) + \149 set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) + 150 (abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) + 151 151 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)} 152 152 … … 164 164 165 165 # resolution formula for muons 166 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*5.e-6) + \166 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*5.e-6) + 167 167 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)} 168 168 … … 210 210 211 211 # absolute impact parameter smearing formula (in mm) as a function of pt and eta 212 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \212 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + 213 213 (pt > 5.0) * (0.005)} 214 214 … … 225 225 set TowerOutputArray ecalTowers 226 226 set EFlowTowerOutputArray eflowPhotons 227 227 228 228 set EnergyMin 0.5 229 229 set EnergySignificanceMin 1.0 230 230 231 set SmearTowerCenter true 232 231 233 set pi [expr {acos(-1)}] 232 234 … … 234 236 # each list starts with the lower edge of the first tower 235 237 # the list ends with the higher edged of the last tower 236 238 237 239 # 0.5 degree towers 238 240 set PhiBins {} … … 246 248 add EtaPhiBins $eta $PhiBins 247 249 } 248 250 249 251 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL} 250 252 251 253 add EnergyFraction {0} {0.0} 252 254 # energy fractions for e, gamma and pi0 … … 269 271 270 272 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 271 273 272 274 # This is the CMS ECAL resolution, extended up eta = 6.0 273 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) + \275 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) + 274 276 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)} 275 277 276 278 277 279 } 278 280 … … 287 289 set TowerOutputArray hcalTowers 288 290 set EFlowTowerOutputArray eflowNeutralHadrons 289 291 290 292 set EnergyMin 1.0 291 293 set EnergySignificanceMin 1.0 292 294 295 set SmearTowerCenter true 296 293 297 set pi [expr {acos(-1)}] 294 298 … … 296 300 # each list starts with the lower edge of the first tower 297 301 # the list ends with the higher edged of the last tower 298 299 302 303 300 304 # 5 degree towers 301 305 set PhiBins {} … … 309 313 add EtaPhiBins $eta $PhiBins 310 314 } 311 312 315 316 313 317 # default energy fractions {abs(PDG code)} {Fecal Fhcal} 314 318 add EnergyFraction {0} {1.0} … … 332 336 333 337 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 334 338 335 339 # This is the ATLAS HCAL resolution, extended up eta = 6.0 336 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \337 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \340 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + 341 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + 338 342 (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)} 339 343 … … 373 377 set MomentumOutputArray momentum 374 378 } 375 379 376 380 377 381 ################## … … 390 394 391 395 module PdgCodeFilter NeutrinoFilter { 392 396 393 397 set InputArray Delphes/stableParticles 394 398 set OutputArray filteredParticles 395 399 396 400 set PTMin 0.0 397 401 398 402 add PdgCode {12} 399 403 add PdgCode {14} … … 463 467 # maximum distance between jet and track 464 468 set DeltaR 0.3 465 466 # minimum pt of tracks 469 470 # minimum pt of tracks 467 471 set TrackPTMin 1.0 468 472 469 473 # minimum transverse impact parameter (in mm) 470 474 set TrackIPMax 2.0 471 475 472 476 # minimum ip significance for the track to be counted 473 477 set SigMin 6.5 474 478 475 479 # minimum number of tracks (high efficiency n=2, high purity n=3) 476 480 set Ntracks 3 … … 510 514 add Branch Delphes/allParticles Particle GenParticle 511 515 add Branch GenJetFinder/jets GenJet Jet 512 516 513 517 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track 514 518 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower … … 521 525 add Branch ScalarHT/energy ScalarHT ScalarHT 522 526 } 523
Note:
See TracChangeset
for help on using the changeset viewer.
