Fork me on GitHub

Changeset 1258 in svn for trunk/examples


Ignore:
Timestamp:
Aug 15, 2013, 3:40:22 PM (11 years ago)
Author:
Pavel Demin
Message:

update ATLAS card

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/delphes_card_ATLAS.tcl

    r1099 r1258  
    3232  FastJetFinder
    3333
    34   ConstituentFilter
     34  JetEnergyScale
    3535
    3636  BTagging
    … …  
    135135
    136136  # 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)}
    145143}
    146144
    … …  
    174172  # resolution formula for muons
    175173  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) + \
     174                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) + \
     175                                           (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
     176                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) + \
    179177                         (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)}
     178                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) + \
     179                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
     180                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
    183181}
    184182
    … …  
    302300module Isolation PhotonIsolation {
    303301  set CandidateInputArray PhotonEfficiency/photons
    304   set IsolationInputArray EFlowMerger/eflow
     302  set IsolationInputArray TrackMerger/tracks
    305303
    306304  set OutputArray photons
    … …  
    336334module Isolation ElectronIsolation {
    337335  set CandidateInputArray ElectronEfficiency/electrons
    338   set IsolationInputArray EFlowMerger/eflow
     336  set IsolationInputArray TrackMerger/tracks
    339337
    340338  set OutputArray electrons
    … …  
    370368module Isolation MuonIsolation {
    371369  set CandidateInputArray MuonEfficiency/muons
    372   set IsolationInputArray EFlowMerger/eflow
     370  set IsolationInputArray TrackMerger/tracks
    373371
    374372  set OutputArray muons
    … …  
    400398  add InputArray UniqueObjectFinder/electrons
    401399  add InputArray UniqueObjectFinder/photons
    402   add InputArray MuonIsolation/muons
     400  add InputArray UniqueObjectFinder/muons
    403401  set EnergyOutputArray energy
    404402}
    … …  
    415413  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    416414  set JetAlgorithm 6
    417   set ParameterR 0.7
    418 
    419   set ConeRadius 0.5
    420   set SeedThreshold 1.0
    421   set ConeAreaFraction 1.0
    422   set AdjacencyCut 2.0
    423   set OverlapThreshold 0.75
    424 
    425   set MaxIterations 100
    426   set MaxPairSize 2
    427   set Iratch 1
     415  set ParameterR 0.6
    428416
    429417  set JetPTMin 20.0
    … …  
    435423
    436424module FastJetFinder FastJetFinder {
    437 #  set InputArray Calorimeter/towers
    438   set InputArray EFlowMerger/eflow
     425  set InputArray Calorimeter/towers
    439426
    440427  set OutputArray jets
    … …  
    442429  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    443430  set JetAlgorithm 6
    444   set ParameterR 0.7
    445 
    446   set ConeRadius 0.5
    447   set SeedThreshold 1.0
    448   set ConeAreaFraction 1.0
    449   set AdjacencyCut 2.0
    450   set OverlapThreshold 0.75
    451 
    452   set MaxIterations 100
    453   set MaxPairSize 2
    454   set Iratch 1
     431  set ParameterR 0.6
    455432
    456433  set JetPTMin 20.0
    457434}
    458435
    459 ####################
    460 # Constituent filter
    461 ####################
    462 
    463 module ConstituentFilter ConstituentFilter {
    464 
    465 # add JetInputArray InputArray
    466   add JetInputArray GenJetFinder/jets
    467   add JetInputArray FastJetFinder/jets
    468 
    469 # add ConstituentInputArray InputArray OutputArray
    470   add ConstituentInputArray Delphes/stableParticles stableParticles
    471   add ConstituentInputArray Calorimeter/eflowTracks eflowTracks
    472   add ConstituentInputArray Calorimeter/eflowTowers eflowTowers
    473   add ConstituentInputArray MuonMomentumSmearing/muons muons
     436##################
     437# Jet Energy Scale
     438##################
     439
     440module EnergyScale JetEnergyScale {
     441  set InputArray FastJetFinder/jets
     442  set OutputArray jets
     443
     444 # scale formula for jets
     445  set ScaleFormula {1.08}
    474446}
    475447
    … …  
    496468  # default efficiency formula (misidentification rate)
    497469  add EfficiencyFormula {0} {0.001}
     470
    498471  # efficiency formula for c-jets (misidentification rate)
    499   add EfficiencyFormula {4} {0.1}
     472  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
     473                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
     474                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     475                              (abs(eta) > 2.5)                                  * (0.000)}
     476
    500477  # efficiency formula for b-jets
    501   add EfficiencyFormula {5} {0.4}
     478  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
     479                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
     480                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     481                              (abs(eta) > 2.5)                                  * (0.000)}
    502482}
    503483
    … …  
    530510  add InputArray PhotonIsolation/photons photons
    531511  add InputArray ElectronIsolation/electrons electrons
     512  add InputArray MuonIsolation/muons muons
    532513  add InputArray FastJetFinder/jets jets
    533514}
    … …  
    542523  add Branch TrackMerger/tracks Track Track
    543524  add Branch Calorimeter/towers Tower Tower
    544   add Branch ConstituentFilter/eflowTracks EFlowTrack Track
    545   add Branch ConstituentFilter/eflowTowers EFlowTower Tower
    546   add Branch ConstituentFilter/muons EFlowMuon Muon
     525#  add Branch Calorimeter/eflowTracks EFlowTrack Track
     526#  add Branch Calorimeter/eflowTowers EFlowTower Tower
     527#  add Branch MuonMomentumSmearing/muons EFlowMuon Muon
    547528  add Branch GenJetFinder/jets GenJet Jet
    548529  add Branch UniqueObjectFinder/jets Jet Jet
    549530  add Branch UniqueObjectFinder/electrons Electron Electron
    550531  add Branch UniqueObjectFinder/photons Photon Photon
    551   add Branch MuonIsolation/muons Muon Muon
     532  add Branch UniqueObjectFinder/muons Muon Muon
    552533  add Branch MissingET/momentum MissingET MissingET
    553534  add Branch ScalarHT/energy ScalarHT ScalarHT
Note: See TracChangeset for help on using the changeset viewer.