Fork me on GitHub

Changeset 2aaf039 in git for cards


Ignore:
Timestamp:
Apr 25, 2016, 10:52:28 AM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
bd225dd
Parents:
576428a
Message:

added muon id efficiencies

Location:
cards
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII_0PU.tcl

    r576428a r2aaf039  
    4040 
    4141  MuonIsolation
    42   MuonEfficiency
     42  MuonLooseIdEfficiency
     43  MuonTightIdEfficiency
    4344 
    4445  NeutrinoFilter
    … …  
    142143  set EfficiencyFormula {
    143144      (pt <= 0.2) * (0.00) + \
    144           (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
    145           (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
    146           (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
    147           (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
    148           (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
    149           (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
     145          (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
     146          (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
     147          (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
     148          (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
     149          (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
     150          (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
    150151          (abs(eta) > 4.0) * (0.00)
    151152         
    … …  
    800801}
    801802
    802 
    803 ###################
    804 # Muon efficiency #
    805 ###################
    806 
    807 module Efficiency MuonEfficiency {
    808  
    809   set InputArray MuonIsolation/muons
    810  
    811   set OutputArray muons
    812   # set EfficiencyFormula {efficiency as a function of eta and pt}
    813   # efficiency formula for muons
    814   set EfficiencyFormula { (pt <= 2.0)  * (0.00) + \
    815                           (abs(eta) <= 4.00) * (pt >  2.0 && pt <= 3.0)  * (0.51) + \
    816                           (abs(eta) <= 4.00) * (pt >  3.0 && pt <= 4.0)  * (0.85) + \
    817                           (abs(eta) <= 4.00) * (pt >  4.0 && pt <= 11.0) * (0.93) + \
    818                           (abs(eta) <= 4.00) * (pt >  11. && pt <= 50.)  * (0.96) + \
    819                           (abs(eta) <= 4.00) * (pt >  50. && pt <= 70.)  * (0.98) + \
    820                           (abs(eta) <= 4.00) * (pt > 70.0 )  * (1.00) + \
    821                           (abs(eta) > 4.00)  * (0.00)
    822  }
    823 
    824 }
    825 
    826 
    827 
    828 
     803##################
     804# Muon Loose Id  #
     805##################
     806
     807module Efficiency MuonLooseIdEfficiency {
     808    set InputArray MuonIsolation/muons
     809    set OutputArray muons
     810    # tracking + LooseID efficiency formula for muons
     811    set EfficiencyFormula {
     812       
     813        (pt <= 5.0) * (0.00) + \
     814                (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.967 + \
     815                (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.968 + \
     816                (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
     817                (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.982 + \
     818                (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.972 + \
     819                (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.974 + \
     820                (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.969 + \
     821                (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.988 + \
     822                (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.989 + \
     823                (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.992 + \
     824                (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.979 + \
     825                (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.976 + \
     826                (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.882 + \
     827                (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.883 + \
     828                (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.851 + \
     829                (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.828 + \
     830
     831                (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.989 + \
     832                (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
     833                (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
     834                (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.995 + \
     835                (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.987 + \
     836                (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
     837                (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
     838                (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.99 + \
     839                (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.997 + \
     840                (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
     841                (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.989 + \
     842                (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.995 + \
     843                (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.888 + \
     844                (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.894 + \
     845                (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.886 + \
     846                (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.883 + \
     847
     848        (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.981 + \
     849                (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.994 + \
     850                (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
     851                (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.997 + \
     852                (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.986 + \
     853                (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.991 + \
     854                (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.987 + \
     855                (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.992 + \
     856                (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.993 + \
     857                (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.999 + \
     858                (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.994 + \
     859                (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.997 + \
     860                (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.894 + \
     861                (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.879 + \
     862                (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.871 + \
     863                (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.864 + \
     864
     865        (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.986 + \
     866                (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
     867                (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.996 + \
     868                (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.994 + \
     869                (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.991 + \
     870                (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
     871                (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.985 + \
     872                (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.993 + \
     873                (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.996 + \
     874                (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
     875                (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.99 + \
     876                (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.993 + \
     877                (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.905 + \
     878                (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.86 + \
     879                (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.845 + \
     880                (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.855 + \
     881
     882        (abs(eta) > 2.8) * (0.00)
     883
     884    }
     885}
     886
     887
     888##################
     889# Muon Tight Id  #
     890##################
     891
     892module Efficiency MuonTightIdEfficiency {
     893    set InputArray MuonIsolation/muons
     894    set OutputArray muons
     895    # tracking + TightID efficiency formula for muons
     896    set EfficiencyFormula {
     897
     898        (pt <= 5.0) * (0.00) + \
     899                (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.949 + \
     900                (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.894 + \
     901                (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.966 + \
     902                (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.964 + \
     903                (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.923 + \
     904                (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.94 + \
     905                (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.959 + \
     906                (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.972 + \
     907                (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.962 + \
     908                (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.972 + \
     909                (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.947 + \
     910                (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.911 + \
     911                (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.825 + \
     912                (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.821 + \
     913                (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.783 + \
     914                (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.746 + \
     915
     916                (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
     917                (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
     918                (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
     919                (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.976 + \
     920                (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.954 + \
     921                (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.978 + \
     922                (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.975 + \
     923                (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.98 + \
     924                (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.982 + \
     925                (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
     926                (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.962 + \
     927                (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.946 + \
     928                (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.839 + \
     929                (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.844 + \
     930                (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.803 + \
     931                (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.791 + \
     932
     933        (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.965 + \
     934                (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.938 + \
     935                (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.99 + \
     936                (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.977 + \
     937                (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.959 + \
     938                (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.976 + \
     939                (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.984 + \
     940                (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.981 + \
     941                (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.978 + \
     942                (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
     943                (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.969 + \
     944                (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.964 + \
     945                (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.863 + \
     946                (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.828 + \
     947                (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.807 + \
     948                (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.764 + \
     949
     950        (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
     951                (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
     952                (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.991 + \
     953                (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.983 + \
     954                (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.957 + \
     955                (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.979 + \
     956                (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
     957                (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.983 + \
     958                (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.981 + \
     959                (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.984 + \
     960                (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.964 + \
     961                (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.956 + \
     962                (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.855 + \
     963                (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.817 + \
     964                (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.764 + \
     965                (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.752 + \
     966
     967        (abs(eta) > 2.8) * (0.00)
     968
     969    }
     970}
    829971
    830972########################
    … …  
    14831625  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    14841626 
    1485   add Branch PhotonIsolation/photons Photon Photon
    1486   add Branch ElectronIsolation/electrons Electron Electron
    1487   add Branch MuonIsolation/muons Muon Muon
     1627  add Branch PhotonEfficiency/photons Photon Photon
     1628  add Branch ElectronEfficiency/electrons Electron Electron
     1629  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
     1630  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    14881631 
    14891632  add Branch JetEnergyScale/jets Jet Jet

  • cards/CMS_PhaseII_200PU.tcl

    r576428a r2aaf039  
    191191  set EfficiencyFormula {
    192192      (pt <= 0.2) * (0.00) + \
    193           (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
    194           (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
    195           (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
    196           (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
    197           (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
    198           (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
     193          (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
     194          (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
     195          (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
     196          (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
     197          (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
     198          (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
    199199          (abs(eta) > 4.0) * (0.00)
    200200         
    … …  
    935935
    936936
    937 ###################
    938 # Muon efficiency #
    939 ###################
    940 
    941 module Efficiency MuonEfficiency {
    942  
    943   set InputArray MuonIsolation/muons
    944  
    945   set OutputArray muons
    946   # set EfficiencyFormula {efficiency as a function of eta and pt}
    947   # efficiency formula for muons
    948   set EfficiencyFormula { (pt <= 2.0)  * (0.00) + \
    949                           (abs(eta) <= 4.00) * (pt >  2.0 && pt <= 3.0)  * (0.51) + \
    950                           (abs(eta) <= 4.00) * (pt >  3.0 && pt <= 4.0)  * (0.85) + \
    951                           (abs(eta) <= 4.00) * (pt >  4.0 && pt <= 11.0) * (0.93) + \
    952                           (abs(eta) <= 4.00) * (pt >  11. && pt <= 50.)  * (0.96) + \
    953                           (abs(eta) <= 4.00) * (pt >  50. && pt <= 70.)  * (0.98) + \
    954                           (abs(eta) <= 4.00) * (pt > 70.0 )  * (1.00) + \
    955                           (abs(eta) > 4.00)  * (0.00)
    956  }
    957 
    958 }
    959 
    960 
     937##################
     938# Muon Loose Id  #
     939##################
     940
     941module Efficiency MuonLooseIdEfficiency {
     942    set InputArray MuonIsolation/muons
     943    set OutputArray muons
     944    # tracking + LooseID efficiency formula for muons
     945    set EfficiencyFormula {
     946       
     947        (pt <= 5.0) * (0.00) + \
     948                (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.967 + \
     949                (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.968 + \
     950                (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
     951                (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.982 + \
     952                (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.972 + \
     953                (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.974 + \
     954                (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.969 + \
     955                (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.988 + \
     956                (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.989 + \
     957                (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.992 + \
     958                (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.979 + \
     959                (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.976 + \
     960                (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.882 + \
     961                (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.883 + \
     962                (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.851 + \
     963                (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.828 + \
     964
     965                (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.989 + \
     966                (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
     967                (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
     968                (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.995 + \
     969                (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.987 + \
     970                (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
     971                (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
     972                (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.99 + \
     973                (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.997 + \
     974                (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
     975                (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.989 + \
     976                (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.995 + \
     977                (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.888 + \
     978                (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.894 + \
     979                (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.886 + \
     980                (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.883 + \
     981
     982        (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.981 + \
     983                (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.994 + \
     984                (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
     985                (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.997 + \
     986                (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.986 + \
     987                (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.991 + \
     988                (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.987 + \
     989                (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.992 + \
     990                (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.993 + \
     991                (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.999 + \
     992                (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.994 + \
     993                (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.997 + \
     994                (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.894 + \
     995                (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.879 + \
     996                (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.871 + \
     997                (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.864 + \
     998
     999        (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.986 + \
     1000                (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
     1001                (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.996 + \
     1002                (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.994 + \
     1003                (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.991 + \
     1004                (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
     1005                (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.985 + \
     1006                (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.993 + \
     1007                (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.996 + \
     1008                (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
     1009                (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.99 + \
     1010                (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.993 + \
     1011                (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.905 + \
     1012                (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.86 + \
     1013                (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.845 + \
     1014                (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.855 + \
     1015
     1016        (abs(eta) > 2.8) * (0.00)
     1017
     1018    }
     1019}
     1020
     1021
     1022##################
     1023# Muon Tight Id  #
     1024##################
     1025
     1026module Efficiency MuonTightIdEfficiency {
     1027    set InputArray MuonIsolation/muons
     1028    set OutputArray muons
     1029    # tracking + TightID efficiency formula for muons
     1030    set EfficiencyFormula {
     1031
     1032        (pt <= 5.0) * (0.00) + \
     1033                (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.949 + \
     1034                (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.894 + \
     1035                (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.966 + \
     1036                (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.964 + \
     1037                (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.923 + \
     1038                (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.94 + \
     1039                (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.959 + \
     1040                (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.972 + \
     1041                (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.962 + \
     1042                (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.972 + \
     1043                (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.947 + \
     1044                (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.911 + \
     1045                (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.825 + \
     1046                (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.821 + \
     1047                (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.783 + \
     1048                (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.746 + \
     1049
     1050                (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
     1051                (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
     1052                (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
     1053                (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.976 + \
     1054                (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.954 + \
     1055                (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.978 + \
     1056                (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.975 + \
     1057                (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.98 + \
     1058                (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.982 + \
     1059                (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
     1060                (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.962 + \
     1061                (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.946 + \
     1062                (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.839 + \
     1063                (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.844 + \
     1064                (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.803 + \
     1065                (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.791 + \
     1066
     1067        (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.965 + \
     1068                (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.938 + \
     1069                (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.99 + \
     1070                (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.977 + \
     1071                (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.959 + \
     1072                (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.976 + \
     1073                (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.984 + \
     1074                (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.981 + \
     1075                (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.978 + \
     1076                (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
     1077                (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.969 + \
     1078                (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.964 + \
     1079                (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.863 + \
     1080                (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.828 + \
     1081                (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.807 + \
     1082                (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.764 + \
     1083
     1084        (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
     1085                (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
     1086                (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.991 + \
     1087                (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.983 + \
     1088                (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.957 + \
     1089                (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.979 + \
     1090                (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
     1091                (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.983 + \
     1092                (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.981 + \
     1093                (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.984 + \
     1094                (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.964 + \
     1095                (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.956 + \
     1096                (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.855 + \
     1097                (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.817 + \
     1098                (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.764 + \
     1099                (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.752 + \
     1100
     1101        (abs(eta) > 2.8) * (0.00)
     1102
     1103    }
     1104}
    9611105
    9621106
    … …  
    16301774#  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    16311775 
    1632   add Branch PhotonIsolation/photons Photon Photon
    1633   add Branch ElectronIsolation/electrons Electron Electron
    1634   add Branch MuonIsolation/muons Muon Muon
    1635  
     1776  add Branch PhotonEfficiency/photons Photon Photon
     1777  add Branch ElectronEfficiency/electrons Electron Electron
     1778  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
     1779  add Branch MuonTightIdEfficiency/muons MuonTight Muon
     1780   
    16361781  add Branch JetEnergyScale/jets Jet Jet
    16371782 
Note: See TracChangeset for help on using the changeset viewer.