Fork me on GitHub

Changeset 2aaf039 in git for cards/CMS_PhaseII_200PU.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_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.