| 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]
|
|---|