Fork me on GitHub

Changeset 2aaf039 in git for cards/CMS_PhaseII_0PU.tcl


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

added muon id efficiencies

File:
1 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
Note: See TracChangeset for help on using the changeset viewer.