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
|
---|
12 | set main.fastsim.package = fastjet
|
---|
13 | set main.fastsim.algorithm = antikt
|
---|
14 | set main.fastsim.radius = 0.4
|
---|
15 | set main.fastsim.bjet_id.matching_dr = 0.3
|
---|
16 | set main.fastsim.bjet_id.exclusive = true
|
---|
17 | set main.fastsim.exclusive_id = false
|
---|
18 | set main.fastsim.ptmin = 20.0
|
---|
19 |
|
---|
20 | # SFS Configuration
|
---|
21 | set 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
|
---|
26 | define smearer j with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
|
---|
27 | define smearer j with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
|
---|
28 | define 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
|
---|
31 | define jes 1.08
|
---|
32 |
|
---|
33 | ### Electron Smearing
|
---|
34 | ## Momentum Smearing
|
---|
35 | define smearer e with PT pt*0.015 [pt > 0.1 and pt <= 25 and abseta <= 2.5]
|
---|
36 | define smearer e with PT sqrt((pt^2)*0.005^2 + pt*0.05^2 + 0.25^2) [pt > 25 and abseta <= 2.5]
|
---|
37 | define smearer e with PT sqrt((pt^2)*0.005^2 + pt*0.05^2 + 0.25^2) [abseta > 2.5 and abseta <= 3]
|
---|
38 | define 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
|
---|
44 | define smearer e with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
|
---|
45 | define 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
|
---|
48 | define reco_efficiency e 0.0 [pt <= 10.0 or abseta > 2.5 or (abseta > 1.37 and abseta < 1.52)]
|
---|
49 | define reco_efficiency e 0.45 [abseta <=1.37 and pt > 10.0]
|
---|
50 | define reco_efficiency e 0.35 [abseta >= 1.52 and abseta <= 2.5 and pt > 10.0]
|
---|
51 |
|
---|
52 |
|
---|
53 | ## Muon Smearing
|
---|
54 | ## Momentum Smearing
|
---|
55 | define smearer mu with pt 0.03 [abseta <= 1.5 and pt > 0.1 and pt <= 1]
|
---|
56 | define smearer mu with pt 0.03 [abseta <= 1.5 and pt > 1 and pt <= 50]
|
---|
57 | define smearer mu with pt 0.04 [abseta <= 1.5 and pt > 50 and pt <= 100]
|
---|
58 | define smearer mu with pt 0.07 [abseta <= 1.5 and pt > 100]
|
---|
59 | define smearer mu with pt 0.04 [abseta > 1.5 and abseta <= 2.5 and pt > 0.1 and pt <= 1]
|
---|
60 | define smearer mu with pt 0.04 [abseta > 1.5 and abseta <= 2.5 and pt > 1 and pt <= 50]
|
---|
61 | define smearer mu with pt 0.05 [abseta > 1.5 and abseta <= 2.5 and pt > 50 and pt <= 100]
|
---|
62 | define smearer mu with pt 0.1 [abseta > 1.5 and abseta <= 2.5 and pt > 100]
|
---|
63 |
|
---|
64 |
|
---|
65 | ## Reconstruction efficiency
|
---|
66 | define reco_efficiency mu 0.0 [abseta > 2.5 or pt <= 10.0]
|
---|
67 | define reco_efficiency mu 0.45 [ABSETA <= 1.5 and pt > 10.0]
|
---|
68 | define reco_efficiency mu 0.55 [ABSETA > 1.5 and abseta < 2.5 and pt > 10.0]
|
---|
69 |
|
---|
70 | ### Photon reconstruction
|
---|
71 | ## ECAL
|
---|
72 | define smearer a with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
|
---|
73 | define 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)
|
---|
76 | define reco_efficiency a 0.0 [pt <= 10.0 or abseta > 2.5]
|
---|
77 | define reco_efficiency a 0.95 [abseta <=1.5 and pt > 10.0]
|
---|
78 | define 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
|
---|
84 | define tagger j as b 0.12 [pt>=20 and pt<40]
|
---|
85 | define tagger j as b 0.04 [pt>=40 and pt<60]
|
---|
86 | define tagger j as b 0.03 [pt>=60 and pt<90]
|
---|
87 | define tagger j as b 0.02 [pt>=90 and pt<150]
|
---|
88 | define tagger j as b 0.02 [pt>=150 and pt<200]
|
---|
89 | define tagger j as b 0.04 [pt>=200 and pt<300]
|
---|
90 | define tagger j as b 0.08 [pt>=300 and pt<500]
|
---|
91 | define 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
|
---|
94 | define tagger c as b 0.31 [pt>=20 and pt<40 and abseta<2.5]
|
---|
95 | define tagger c as b 0.21 [pt>=40 and pt<60 and abseta<2.5]
|
---|
96 | define tagger c as b 0.17 [pt>=60 and pt<90 and abseta<2.5]
|
---|
97 | define tagger c as b 0.14 [pt>=90 and pt<150 and abseta<2.5]
|
---|
98 | define tagger c as b 0.16 [pt>=150 and pt<200 and abseta<2.5]
|
---|
99 | define tagger c as b 0.18 [pt>=200 and pt<300 and abseta<2.5]
|
---|
100 | define tagger c as b 0.24 [pt>=300 and pt<500 and abseta<2.5]
|
---|
101 | define tagger c as b 0.31 [pt>=500 and abseta<2.5]
|
---|
102 | define 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
|
---|
105 | define tagger b as b 0.6 [pt>=20 and pt<40 and abseta<2.5]
|
---|
106 | define tagger b as b 0.7 [pt>=40 and pt<60 and abseta<2.5]
|
---|
107 | define tagger b as b 0.8 [pt>=60 and pt<200 and abseta<2.5]
|
---|
108 | define tagger b as b 0.7 [pt>=200 and pt<500 and abseta<2.5]
|
---|
109 | define tagger b as b 0.6 [pt>=500 and abseta<2.5]
|
---|
110 | define tagger b as b 0.0 [abseta >= 2.5 or pt < 20]
|
---|
111 |
|
---|
112 | ## Hadronic Taus
|
---|
113 |
|
---|
114 | ## HCAL
|
---|
115 | define smearer ta with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
|
---|
116 | define smearer ta with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
|
---|
117 | define smearer ta with E sqrt(E^2*0.0942^2 + E) [abseta > 3.2 and abseta <= 4.9]
|
---|
118 |
|
---|
119 | define tagger j as ta 0.01 [abseta <= 2.7]
|
---|
120 | define tagger j as ta 0.0 [abseta > 2.5]
|
---|
121 |
|
---|
122 | define tagger ta as ta 0.6 [abseta <= 2.5]
|
---|
123 | define tagger ta as ta 0.0 [abseta > 2.5]
|
---|