SFS: ATLAS_default.ma5

File ATLAS_default.ma5, 5.3 KB (added by jackaraz, 4 years ago)
Line 
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
14set main.fastsim.package = fastjet
15set main.fastsim.algorithm = antikt
16set main.fastsim.radius = 0.6
17set main.fastsim.bjet_id.matching_dr = 0.3
18set main.fastsim.bjet_id.exclusive = true
19set main.fastsim.exclusive_id = false
20set main.fastsim.ptmin = 20.0
21set main.fastsim.jetrecomode = jets
22
23# Jet Smearing
24
25## Tracker effects are approximated (optimized for dijet process)
26define reco_efficiency j 0.6 [abseta <= 1.5]
27define reco_efficiency j 0.975 [abseta > 1.5 and abseta <= 2.5]
28define reco_efficiency j 0.98 [abseta > 2.5 and abseta <= 3.5]
29define 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
41define smearer j with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
42define smearer j with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
43define 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
48define smearer e with PT sqrt(0.03^2 + pt^2*1.3e-3^2) [abseta <= 0.5 and pt > 0.1]
49define 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]
50define 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
58define smearer e with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
59define 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
63define reco_efficiency e 0.0 [pt <= 10.0 or abseta > 2.5]
64define reco_efficiency e 0.9 [abseta <=1.5 and pt > 10.0]
65define 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
73define smearer mu with PT sqrt(0.01^2 + pt^2*1.0e-4^2) [abseta <= 0.5 and pt > 0.1]
74define 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]
75define 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)
81define reco_efficiency mu 0.0 [abseta > 2.5 or pt <= 10.0]
82define reco_efficiency mu 0.95 [ABSETA <= 1.5 and pt > 10.0]
83define 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
95define smearer a with E sqrt(E^2*0.0017^2 + E*0.101^2) [abseta <= 3.2]
96define 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)
99define reco_efficiency a 0.0 [pt <= 10.0 or abseta > 2.5]
100define reco_efficiency a 0.95 [abseta <=1.5 and pt > 10.0]
101define 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
108define tagger j as b 0.002+7.3e-06*pt
109define tagger c as b 0.20*tanh(0.02*pt)*(1/(1+0.0034*pt)) [abseta < 2.5]
110define tagger c as b 0.0 [abseta >=2.5]
111define tagger b as b 0.80*tanh(0.003*pt)*(30/(1+0.086*pt)) [abseta < 2.5]
112define tagger b as b 0.0 [abseta >= 2.5]
113
114## Hadronic Taus
115
116## HCAL
117
118define smearer ta with E sqrt(E^2*0.0302^2 + E*0.5205^2 + 1.59^2) [abseta <= 1.7]
119define smearer ta with E sqrt(E^2*0.05^2 + E*0.706^2) [abseta > 1.7 and abseta <= 3.2]
120define smearer ta with E sqrt(E^2*0.0942^2 + E) [abseta > 3.2 and abseta <= 4.9]
121
122define tagger j as ta 0.01 [ntracks >= 2 and abseta <= 2.7]
123define tagger j as ta 0.02 [ntracks == 1 and abseta <= 2.7]
124define tagger j as ta 0.0 [abseta > 2.7]
125
126define tagger ta as ta 0.6 [abseta <= 2.5 and ntracks >= 2]
127define tagger ta as ta 0.7 [abseta <= 2.5 and ntracks == 1]
128define tagger ta as ta 0.0 [abseta > 2.5]
129
130# reconstruction efficiencies (optimized for ditau process)
131define reco_efficiency ta 0.0 [pt <= 20.]
132define reco_efficiency ta 0.65 [pt > 20. and abseta <= 2.5]
133define reco_efficiency ta 1. [pt > 20. and abseta > 2.5] # safety
134