SFS: sfs_atlas_susy_2016_07.ma5

#<><><><><><><><><><><><><><><><><><><><><><><><><>
#                                                <>
#   ATLAS-SUSY-2016-007                          <>
#   Author: Jack Y. Araz, Benjamin Fuks          <>
#                                                <>
#       Please cite arXiv:2006.09387,1712.02332  <>
#                                                <>
#   Based on MadAnalysis 5 v1.8                  <>
#                                                <>
#<><><><><><><><><><><><><><><><><><><><><><><><><>

set main.fastsim.package = fastjet
set main.fastsim.algorithm = antikt
set main.fastsim.radius = 0.4
set main.fastsim.bjet_id.matching_dr = 0.3
set main.fastsim.bjet_id.exclusive = true
set main.fastsim.exclusive_id = false
set main.fastsim.ptmin = 20.0
set main.fastsim.jetrecomode = jets

# Jets

define smearer j with PT sqrt(0.06^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
define smearer j with PT sqrt(0.10^2 + pt^2*1.7e-3^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
define smearer j with PT sqrt(0.25^2 + pt^2*3.1e-3^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]

define smearer j with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
define smearer j with E sqrt(E^2*0.05^2   + E*0.706^2)           [abseta > 1.7 and abseta <= 3.2]
define smearer j with E sqrt(E^2*0.0942^2 + E)                   [abseta > 3.2 and abseta <= 4.9]

# Electron

define smearer e with PT sqrt(0.06^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
define smearer e with PT sqrt(0.1^2  + pt^2*1.7e-3^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
define smearer e with PT sqrt(0.25^2 + pt^2*3.1e-3^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]

define smearer e with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
define smearer e with E sqrt(E^2*0.035^2  + E*0.285^2) [abseta>3.2 and abseta<=4.9]

define reco_efficiency e 0.0  [pt <= 10.0 or abseta > 2.5]
define reco_efficiency e 0.95 [abseta <= 1.5 and pt > 10.0]
define reco_efficiency e 0.85 [abseta > 1.5 and abseta <= 2.5 and pt > 10.0]

# Muon

define smearer mu with PT sqrt(0.02^2  + (pt^2)*2.0e-4^2) [abseta <= 0.5 and pt > 0.1]
define smearer mu with PT sqrt(0.03^2  + (pt^2)*3.0e-4^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
define smearer mu with PT sqrt(0.06^2  + (pt^2)*6.0e-4^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]

define reco_efficiency mu 0.0  [pt <= 10.0 or abseta > 2.7]
define reco_efficiency mu 0.95 [abseta <= 1.5 and pt > 10.0]
define reco_efficiency mu 0.85 [abseta > 1.5 and abseta <= 2.7 and pt > 10.0]

# Photon

define smearer a with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
define smearer a with E sqrt(E^2*0.0350^2 + E*0.285^2) [abseta>3.2 and abseta<=4.9]

define reco_efficiency a 0.0  [pt <= 10.0 or abseta > 2.5]
define reco_efficiency a 0.95 [abseta <= 1.5 and pt > 10.0]
define reco_efficiency a 0.85 [abseta > 1.5 and abseta <= 2.5 and pt > 10.0]

# B-tagging

define tagger j as b 0.002+7.3e-06*pt
define tagger c as b 0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))
define tagger b as b 0.80*tanh(0.003*pt)*(30/(1+0.086*pt)) [abseta < 2.5]
define tagger b as b 0.0                                   [abseta >= 2.5]


# Tau

define smearer ta with PT sqrt(0.06^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
define smearer ta with PT sqrt(0.10^2 + pt^2*1.7e-3^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
define smearer ta with PT sqrt(0.25^2 + pt^2*3.1e-3^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]

define smearer ta with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
define smearer ta with E sqrt(E^2*0.05^2   + E*0.706^2)           [abseta > 1.7 and abseta <= 3.2]
define smearer ta with E sqrt(E^2*0.0942^2 + E)                   [abseta > 3.2 and abseta <= 4.9]

define tagger j as ta  0.01
define tagger ta as ta 0.6 [abseta < 2.5]
define tagger ta as ta 0.0 [abseta >= 2.5]

define reco_efficiency ta 0.0 [pt <= 1.]
define reco_efficiency ta 1.0 [pt > 1.]