Context Navigation

Back to SFS

SFS: ATLAS_default.ma5

File ATLAS_default.ma5, 5.3 KB (added by jackaraz, 6 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
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^20.0302^2 + E0.5205^2 + 1.59^2) [abseta <= 1.7]
42	define smearer j with E sqrt(E^20.05^2 + E0.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^20.0017^2 + E0.101^2) [abseta <= 3.2]
59	define smearer e with E sqrt(E^20.0350^2 + E0.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^20.0017^2 + E0.101^2) [abseta <= 3.2]
96	define smearer a with E sqrt(E^20.0350^2 + E0.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.20tanh(0.02pt)(1/(1+0.0034pt)) [abseta < 2.5]
110	define tagger c as b 0.0 [abseta >=2.5]
111	define tagger b as b 0.80tanh(0.003pt)(30/(1+0.086pt)) [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^20.0302^2 + E0.5205^2 + 1.59^2) [abseta <= 1.7]
119	define smearer ta with E sqrt(E^20.05^2 + E0.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

Download in other formats:

Original Format