1 | #<><><><><><><><><><><><><><><><><><><><><><><><><>
|
---|
2 | # <>
|
---|
3 | # ATLAS-Default <>
|
---|
4 | # Author: Jack Y. Araz, Benjamin Fuks <>
|
---|
5 | # <>
|
---|
6 | # Please cite arXiv:2006.09387 <>
|
---|
7 | # <>
|
---|
8 | # Based on MadAnalysis 5 v1.8 <>
|
---|
9 | # <>
|
---|
10 | #<><><><><><><><><><><><><><><><><><><><><><><><><>
|
---|
11 |
|
---|
12 | # Set FastJet variables
|
---|
13 |
|
---|
14 | set main.fastsim.package = fastjet
|
---|
15 | set main.fastsim.algorithm = antikt
|
---|
16 | set main.fastsim.radius = 0.6
|
---|
17 | set main.fastsim.bjet_id.matching_dr = 0.3
|
---|
18 | set main.fastsim.bjet_id.exclusive = true
|
---|
19 | set main.fastsim.exclusive_id = false
|
---|
20 | set main.fastsim.ptmin = 20.0
|
---|
21 | set main.fastsim.jetrecomode = jets
|
---|
22 |
|
---|
23 | # Jet Smearing
|
---|
24 |
|
---|
25 | ## Tracker effects are approximated (optimized for dijet process)
|
---|
26 | define reco_efficiency j 0.6 [abseta <= 1.5]
|
---|
27 | define reco_efficiency j 0.975 [abseta > 1.5 and abseta <= 2.5]
|
---|
28 | define reco_efficiency j 0.98 [abseta > 2.5 and abseta <= 3.5]
|
---|
29 | define reco_efficiency j 0.975 [abseta > 3.5 ]
|
---|
30 |
|
---|
31 |
|
---|
32 | #define smearer j with PT sqrt(0.06^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
|
---|
33 | #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]
|
---|
34 | #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]
|
---|
35 |
|
---|
36 | ## HCAL
|
---|
37 | # http://arxiv.org/pdf/hep-ex/0004009v1
|
---|
38 | # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
|
---|
39 | # set HCalResolutionFormula {resolution formula as a function of eta and energy}
|
---|
40 |
|
---|
41 | define smearer j with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
|
---|
42 | define smearer j with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
|
---|
43 | define smearer j with E sqrt(E^2*0.0942^2 + E) [abseta > 3.2 and abseta <= 4.9]
|
---|
44 |
|
---|
45 | ### Electron Smearing
|
---|
46 |
|
---|
47 | ## Momentum Smearing
|
---|
48 | define smearer e with PT sqrt(0.03^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
|
---|
49 | define smearer e with PT sqrt(0.05^2 + pt^2*1.7e-3^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
|
---|
50 | define smearer e with PT sqrt(0.15^2 + pt^2*3.1e-3^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]
|
---|
51 |
|
---|
52 | ## ECAL
|
---|
53 | # set ECalResolutionFormula {resolution formula as a function of eta and energy}
|
---|
54 | # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
|
---|
55 | # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
|
---|
56 | # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
|
---|
57 |
|
---|
58 | define smearer e with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
|
---|
59 | define smearer e with E sqrt(E^2*0.0350^2 + E*0.285^2) [abseta>3.2 and abseta<=4.9]
|
---|
60 |
|
---|
61 | ## Reconstruction efficiency (optimized for Drell Yan process)
|
---|
62 |
|
---|
63 | define reco_efficiency e 0.0 [pt <= 10.0 or abseta > 2.5]
|
---|
64 | define reco_efficiency e 0.9 [abseta <=1.5 and pt > 10.0]
|
---|
65 | define reco_efficiency e 0.7 [abseta > 1.5 and abseta <= 2.5 and pt > 10.0]
|
---|
66 |
|
---|
67 |
|
---|
68 | ## Muon Smearing
|
---|
69 |
|
---|
70 |
|
---|
71 | ## Momentum Smearing
|
---|
72 |
|
---|
73 | define smearer mu with PT sqrt(0.01^2 + pt^2*1.0e-4^2) [abseta <= 0.5 and pt > 0.1]
|
---|
74 | define smearer mu with PT sqrt(0.015^2 + pt^2*1.5e-4^2) [abseta > 0.5 and abseta <= 1.5 and pt > 0.1]
|
---|
75 | define smearer mu with PT sqrt(0.025^2 + pt^2*3.5e-4^2) [abseta > 1.5 and abseta <= 2.5 and pt > 0.1]
|
---|
76 |
|
---|
77 | ## Reconstruction efficiency
|
---|
78 |
|
---|
79 |
|
---|
80 | ## with Tracker (optimized for Drell Yan process)
|
---|
81 | define reco_efficiency mu 0.0 [abseta > 2.5 or pt <= 10.0]
|
---|
82 | define reco_efficiency mu 0.95 [ABSETA <= 1.5 and pt > 10.0]
|
---|
83 | define reco_efficiency mu 0.85 [ABSETA > 1.5 and abseta < 2.5 and pt > 10.0]
|
---|
84 |
|
---|
85 |
|
---|
86 |
|
---|
87 | ### Photon reconstruction
|
---|
88 |
|
---|
89 | ## ECAL
|
---|
90 | # set ECalResolutionFormula {resolution formula as a function of eta and energy}
|
---|
91 | # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
|
---|
92 | # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
|
---|
93 | # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
|
---|
94 |
|
---|
95 | define smearer a with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
|
---|
96 | define smearer a with E sqrt(E^2*0.0350^2 + E*0.285^2) [abseta>3.2 and abseta<=4.9]
|
---|
97 |
|
---|
98 | # Reconstruction efficiency (optimized for diphoton process)
|
---|
99 | define reco_efficiency a 0.0 [pt <= 10.0 or abseta > 2.5]
|
---|
100 | define reco_efficiency a 0.95 [abseta <=1.5 and pt > 10.0]
|
---|
101 | define reco_efficiency a 0.85 [abseta > 1.5 and abseta <= 2.5 and pt > 10.0]
|
---|
102 |
|
---|
103 |
|
---|
104 | ### Particle Identification
|
---|
105 |
|
---|
106 | ## B tagging
|
---|
107 | # based on ATL-PHYS-PUB-2015-022
|
---|
108 | define tagger j as b 0.002+7.3e-06*pt
|
---|
109 | define tagger c as b 0.20*tanh(0.02*pt)*(1/(1+0.0034*pt)) [abseta < 2.5]
|
---|
110 | define tagger c as b 0.0 [abseta >=2.5]
|
---|
111 | define tagger b as b 0.80*tanh(0.003*pt)*(30/(1+0.086*pt)) [abseta < 2.5]
|
---|
112 | define tagger b as b 0.0 [abseta >= 2.5]
|
---|
113 |
|
---|
114 | ## Hadronic Taus
|
---|
115 |
|
---|
116 | ## HCAL
|
---|
117 |
|
---|
118 | define smearer ta with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
|
---|
119 | define smearer ta with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
|
---|
120 | define smearer ta with E sqrt(E^2*0.0942^2 + E) [abseta > 3.2 and abseta <= 4.9]
|
---|
121 |
|
---|
122 | define tagger j as ta 0.01 [ntracks >= 2 and abseta <= 2.7]
|
---|
123 | define tagger j as ta 0.02 [ntracks == 1 and abseta <= 2.7]
|
---|
124 | define tagger j as ta 0.0 [abseta > 2.7]
|
---|
125 |
|
---|
126 | define tagger ta as ta 0.6 [abseta <= 2.5 and ntracks >= 2]
|
---|
127 | define tagger ta as ta 0.7 [abseta <= 2.5 and ntracks == 1]
|
---|
128 | define tagger ta as ta 0.0 [abseta > 2.5]
|
---|
129 |
|
---|
130 | # reconstruction efficiencies (optimized for ditau process)
|
---|
131 | define reco_efficiency ta 0.0 [pt <= 20.]
|
---|
132 | define reco_efficiency ta 0.65 [pt > 20. and abseta <= 2.5]
|
---|
133 | define reco_efficiency ta 1. [pt > 20. and abseta > 2.5] # safety
|
---|
134 |
|
---|