SFS: sfs_atlas_susy_2018_031.ma5

File sfs_atlas_susy_2018_031.ma5, 5.7 KB (added by Benjamin Fuks, 4 years ago)
Line 
1#<><><><><><><><><><><><><><><><><><><><><><><><><>
2# <>
3# ATLAS-SUSY-2018-031 <>
4# Author: Jack Y. Araz <>
5# Please cite arXiv:2006.09387 <>
6# arXiv:1908.03122 <>
7# Based on MadAnalysis 5 v1.9.4+ <>
8# <>
9#<><><><><><><><><><><><><><><><><><><><><><><><><>
10
11# FastJet Configuration
12set main.fastsim.package = fastjet
13set main.fastsim.algorithm = antikt
14set main.fastsim.radius = 0.4
15set main.fastsim.bjet_id.matching_dr = 0.3
16set main.fastsim.bjet_id.exclusive = true
17set main.fastsim.exclusive_id = false
18set main.fastsim.ptmin = 20.0
19
20# SFS Configuration
21set main.fastsim.jetrecomode = jets
22
23## HCAL
24# http://arxiv.org/pdf/hep-ex/0004009v1
25# http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
26define smearer j with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
27define smearer j with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
28define smearer j with E sqrt(E^2*0.0942^2 + E*0.075^2) [abseta > 3.2 and abseta <= 4.9]
29
30# Jet Energy Scale this feature only available for v1.9.4 and beyond
31define jes 1.08
32
33### Electron Smearing
34## Momentum Smearing
35define smearer e with PT pt*0.015 [pt > 0.1 and pt <= 25 and abseta <= 2.5]
36define smearer e with PT sqrt((pt^2)*0.005^2 + pt*0.05^2 + 0.25^2) [pt > 25 and abseta <= 2.5]
37define smearer e with PT sqrt((pt^2)*0.005^2 + pt*0.05^2 + 0.25^2) [abseta > 2.5 and abseta <= 3]
38define smearer e with PT sqrt((pt^2)*0.107^2 + pt*2.08^2) [abseta > 3 and abseta <= 5]
39
40## ECAL
41# http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
42# http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
43# http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
44define smearer e with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
45define smearer e with E sqrt(E^2*0.0350^2 + E*0.285^2) [abseta>3.2 and abseta<=4.9]
46
47# the electrons are disallowed in the transition region between the barrel and the endcap
48define reco_efficiency e 0.0 [pt <= 10.0 or abseta > 2.5 or (abseta > 1.37 and abseta < 1.52)]
49define reco_efficiency e 0.45 [abseta <=1.37 and pt > 10.0]
50define reco_efficiency e 0.35 [abseta >= 1.52 and abseta <= 2.5 and pt > 10.0]
51
52
53## Muon Smearing
54## Momentum Smearing
55define smearer mu with pt 0.03 [abseta <= 1.5 and pt > 0.1 and pt <= 1]
56define smearer mu with pt 0.03 [abseta <= 1.5 and pt > 1 and pt <= 50]
57define smearer mu with pt 0.04 [abseta <= 1.5 and pt > 50 and pt <= 100]
58define smearer mu with pt 0.07 [abseta <= 1.5 and pt > 100]
59define smearer mu with pt 0.04 [abseta > 1.5 and abseta <= 2.5 and pt > 0.1 and pt <= 1]
60define smearer mu with pt 0.04 [abseta > 1.5 and abseta <= 2.5 and pt > 1 and pt <= 50]
61define smearer mu with pt 0.05 [abseta > 1.5 and abseta <= 2.5 and pt > 50 and pt <= 100]
62define smearer mu with pt 0.1 [abseta > 1.5 and abseta <= 2.5 and pt > 100]
63
64
65## Reconstruction efficiency
66define reco_efficiency mu 0.0 [abseta > 2.5 or pt <= 10.0]
67define reco_efficiency mu 0.45 [ABSETA <= 1.5 and pt > 10.0]
68define reco_efficiency mu 0.55 [ABSETA > 1.5 and abseta < 2.5 and pt > 10.0]
69
70### Photon reconstruction
71## ECAL
72define smearer a with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
73define smearer a with E sqrt(E^2*0.0350^2 + E*0.285^2) [abseta>3.2 and abseta<=4.9]
74
75# Reconstruction efficiency (optimized for diphoton process)
76define reco_efficiency a 0.0 [pt <= 10.0 or abseta > 2.5]
77define reco_efficiency a 0.95 [abseta <=1.5 and pt > 10.0]
78define reco_efficiency a 0.85 [abseta > 1.5 and abseta <= 2.5 and pt > 10.0]
79
80### Particle Identification
81
82## B tagging
83# extracted from https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/FTAG-2018-01/figaux_03a.png
84define tagger j as b 0.12 [pt>=20 and pt<40]
85define tagger j as b 0.04 [pt>=40 and pt<60]
86define tagger j as b 0.03 [pt>=60 and pt<90]
87define tagger j as b 0.02 [pt>=90 and pt<150]
88define tagger j as b 0.02 [pt>=150 and pt<200]
89define tagger j as b 0.04 [pt>=200 and pt<300]
90define tagger j as b 0.08 [pt>=300 and pt<500]
91define tagger j as b 0.21 [pt>=500]
92
93# extracted from https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/FTAG-2018-01/figaux_03b.png
94define tagger c as b 0.31 [pt>=20 and pt<40 and abseta<2.5]
95define tagger c as b 0.21 [pt>=40 and pt<60 and abseta<2.5]
96define tagger c as b 0.17 [pt>=60 and pt<90 and abseta<2.5]
97define tagger c as b 0.14 [pt>=90 and pt<150 and abseta<2.5]
98define tagger c as b 0.16 [pt>=150 and pt<200 and abseta<2.5]
99define tagger c as b 0.18 [pt>=200 and pt<300 and abseta<2.5]
100define tagger c as b 0.24 [pt>=300 and pt<500 and abseta<2.5]
101define tagger c as b 0.31 [pt>=500 and abseta<2.5]
102define tagger c as b 0.0 [abseta >=2.5 or pt<20]
103
104# extracted from https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/FTAG-2018-01/figaux_03c.png
105define tagger b as b 0.6 [pt>=20 and pt<40 and abseta<2.5]
106define tagger b as b 0.7 [pt>=40 and pt<60 and abseta<2.5]
107define tagger b as b 0.8 [pt>=60 and pt<200 and abseta<2.5]
108define tagger b as b 0.7 [pt>=200 and pt<500 and abseta<2.5]
109define tagger b as b 0.6 [pt>=500 and abseta<2.5]
110define tagger b as b 0.0 [abseta >= 2.5 or pt < 20]
111
112## Hadronic Taus
113
114## HCAL
115define smearer ta with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
116define smearer ta with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
117define smearer ta with E sqrt(E^2*0.0942^2 + E) [abseta > 3.2 and abseta <= 4.9]
118
119define tagger j as ta 0.01 [abseta <= 2.7]
120define tagger j as ta 0.0 [abseta > 2.5]
121
122define tagger ta as ta 0.6 [abseta <= 2.5]
123define tagger ta as ta 0.0 [abseta > 2.5]