From 2242acdc491de40c93e94316702aba4abf13642f Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Mon, 27 Apr 2020 20:17:57 +0900 Subject: [PATCH 01/11] By CT: modify validation note following Jack's comments Dear All, I have modified our validation note based on Jack's comments in the attached files. (I also attached the one from Jack.) Please read them carefully ! All modifications are marked in BLUE, but still one question (marked in RED) is left. On the other hand, I don't know how to modify the text on the page 3 to avoid them like a shopping list. If you have a better way to deal with them, please directly refine them. Thanks ! I will wait for your feedback, and then reply to Jack. Thanks a lot ! Best regards, Chih-Ting --- .../validation_note_ATLAS_SUSY_2018_04.tex | 88 ++++++++++--------- 1 file changed, 47 insertions(+), 41 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 7f13495..f6f8f4e 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -11,7 +11,8 @@ %% \documentclass{ws-mpla} -\usepackage[super]{cite} +\usepackage{cite} +%\usepackage[super]{cite} \usepackage{graphicx} \usepackage{color,url} \begin{document} @@ -45,8 +46,8 @@ \begin{abstract} We present the MADANALYSIS 5 implementation and validation of the ATLAS-SUSY-2018-04 search. -This ATLAS analysis targets the search for direct stau production in events with two hadronic tau leptons and probes 139 fb$^{-1}$ of LHC proton-proton collisions at a center-of-mass energy of 13 TeV. -The validation of our reimplementation relies on a comparison of our cutflow predictions with the auxiliary table of official ATLAS results in the context of two supersymmetry-inspired simplified benchmark models in which the Standard Model is extended by a neutralino and a stau decaying into a tau lepton and a neutralino. +This ATLAS analysis targets the search for direct stau production in events with two hadronic tau leptons and {\color{blue}uses a dataset of LHC proton-proton collisions with an integrated luninosity of 139 fb$^{-1}$} at a center-of-mass energy of 13 TeV. +The validation of our reimplementation relies on a comparison of our cutflow predictions with the official ATLAS results in the context of two supersymmetry-inspired simplified {\color{blue}model benchmarks} in which the Standard Model is extended by a neutralino and a stau decaying into a tau lepton and a neutralino. \keywords{supersymmetry; stau; hadronic tau lepton.} \end{abstract} @@ -77,28 +78,29 @@ \section{Introduction} \section{Description of the analysis} This analysis targets a final state containing two hadronic tau leptons with a certain amount of missing transverse energy. -The kinematics of di-$\tau +E^{miss}_T$ system is used to reduce the contributions from Standard Model backgrounds. +The kinematics of {\color{blue}hadronic} di-$\tau +E^{miss}_T$ system is used to reduce the contributions from Standard Model backgrounds. First, all the objects are reconstructed and defined. Then a sequence of event selections for the signal final state is applied. \subsection{Object definitions} Jets are reconstructed by means of the anti-$k_t$ algorithm~\cite{Cacciari:2008gp} with a radius parameter set to $R=0.4$. This analysis focuses on jets whose transverse momentum $p^j_T$ and pseudorapidity $\eta^j$ fullfill \begin{equation} -p^j_T > 20 \textrm{GeV}\quad \textrm{and}\quad |\eta^j| < 2.8. +p^j_T > 20\textrm{ GeV}\quad \textrm{and}\quad |\eta^j| < 2.8. \end{equation} -Moreover, the selected jets are tagged as originating from the fragmentation of a $b$-quark with +Moreover, {\color{blue}the selected jets which are required to satisfy \begin{equation} -p^b_T > 20 \textrm{GeV}\quad \textrm{and}\quad |\eta^b| < 2.5. +p^b_T > 20 \textrm{ GeV}\quad \textrm{and}\quad |\eta^b| < 2.5, \end{equation} +are tagged as originating from the fragmentation of $b$-quarks.} A working point with the average $b$-tagging efficiency of $77\%$ is used. This working point corresponds to a $c$-jet and light-jet rejection of $4.9$ and $110$, respectively. Electron candidates are required to have a transverse momentum $p^e_T$ and pseudorapidity $\eta^e$ obeying \begin{equation} -p^e_T > 17 \textrm{GeV}\quad \textrm{and}\quad |\eta^e| < 2.47. +p^e_T > 17 \textrm{ GeV}\quad \textrm{and}\quad |\eta^e| < 2.47. \end{equation} Furthermore, all electron candidates are required to have both track and calorimeter isolations. The condition of the track isolation is \begin{equation} -\sum p_{T,\textrm{tracks}}/p^e_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{GeV}/p^e_T,0.2), +\sum p_{T,\textrm{tracks}}/p^e_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{ GeV}/p^e_T,0.2), \end{equation} the condition of the calorimeter isolation is \begin{equation} @@ -106,16 +108,16 @@ \subsection{Object definitions} \end{equation} and for high transverse momentum electron, \begin{equation} -\sum E_{T,\textrm{tracks}} < max(0.015\times p^e_T,3.5\textrm{GeV})\quad \textrm{with}\quad \Delta R=0.2\quad \textrm{if}\quad p^e_T > 200\textrm{GeV}. +\sum E_{T,\textrm{tracks}} < max(0.015\times p^e_T,3.5\textrm{ GeV})\quad \textrm{with}\quad \Delta R=0.2\quad \textrm{if}\quad p^e_T > 200\textrm{ GeV}. \end{equation} Muon candidate definition is similar, although with slightly looser thresholds, \begin{equation} -p^{\mu}_T > 14 \textrm{GeV}\quad \textrm{and}\quad |\eta^{\mu}| < 2.7, +p^{\mu}_T > 14 \textrm{ GeV}\quad \textrm{and}\quad |\eta^{\mu}| < 2.7. \end{equation} The condition of the track isolation is \begin{equation} -\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{GeV}/p^{\mu}_T,0.3), +\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{ GeV}/p^{\mu}_T,0.3), \end{equation} and the condition of the calorimeter isolation is \begin{equation} @@ -124,14 +126,16 @@ \subsection{Object definitions} Tau lepton candidates are reconstructed with one or three associated charged pion tracks (prongs). %$\sum_i e_i (\textrm{tracks}) = \pm 1$. %For 1-prong (3-prong) $\tau$ lepton candidates, the signal efficiencies are $75\%$($60\%$) and $60\%$($45\%$) for the \textit{medium} and \textit{tight} working points, respectively. -For 1-prong (3-prong) $\tau$ lepton candidates, the signal efficiencies are $75\%$ and $60\%$ for the \textit{medium} working points. Due to the technical restriction, the \textit{tight} working point efficiencies are directly taken from the official ATLAS cutflow table. %Should add more details +For 1-prong (3-prong) $\tau$ lepton candidates, the signal efficiencies are $75\%$ and $60\%$ for the \textit{medium} working points. +{\color{blue}The \textit{tight} working point efficiencies are directly taken from the official ATLAS cutflow table in this recasting.} +%Due to the technical restriction, the \textit{tight} working point efficiencies are directly taken from the official ATLAS cutflow table. %Should add more details The baseline tau lepton candidates are required to have \begin{equation} -p^{\tau}_T > 50(40) \textrm{GeV}\quad \textrm{and}\quad |\eta^{\tau}| < 2.5 +p^{\tau}_T > 50(40) \textrm{ GeV}\quad \textrm{and}\quad |\eta^{\tau}| < 2.5 \end{equation} for the leading (subleading) ones and the transition region between the barrel and endcap calorimeters ($ 1.37 < |\eta^{\tau}| < 1.52 $) is excluded. -Finally, some overlap removal conditions are in order {\color{blue}which are consistent with the analysis code provided in HEPData\cite{hepdata}. Electron is removed if $\Delta R(e,e) < 0.05$}. Tau lepton is removed if $\Delta R(\tau,e/\mu) < 0.2$. {\color{blue}Electron is removed if $\Delta R(e,\mu) < 0.01$. Jet is removed if $\Delta R(j,e/\mu) < 0.2$, and then electron or muon is removed if $\Delta R(e/\mu,j) < 0.4$, and jet is removed if $\Delta R(j,\tau) < 0.4$.} +Finally, some overlap removal conditions are in order which are consistent with the analysis code provided in HEPData\cite{hepdata}. Electron is removed if $\Delta R(e,e) < 0.05$. Tau lepton is removed if $\Delta R(\tau,e/\mu) < 0.2$. Electron is removed if $\Delta R(e,\mu) < 0.01$. Jet is removed if $\Delta R(j,e/\mu) < 0.2$, and then electron or muon is removed if $\Delta R(e/\mu,j) < 0.4$, and jet is removed if $\Delta R(j,\tau) < 0.4$. \subsection{Event selection} @@ -145,9 +149,9 @@ \subsection{Event selection} The events with exactly two \textit{medium} tau lepton candidates with opposite-sign electric charge (OS) are selected. Then, $b$-jet veto is applied to reject events from top quark associated processes. Also, the events with additional light lepton (muon or electron) are rejected. -The reconstructed invariant mass of the two leading tau lepton candidates, $m(\tau_1,\tau_2)$, larger than $120$ GeV is required for removing tau lepton pair from low-mass resonances, $Z$ boson, and Higgs boson events ($Z/H$ veto). +The reconstructed invariant mass of the two leading {\color{blue}hadronic taus}, $m(\tau_1,\tau_2)$, larger than $120$ GeV is required for removing tau lepton pair from low-mass {\color{blue}$Z/H$ resonances}. -All events are required to pass either an \textit{asymmetric di-$\tau$} trigger for the low stau mass region (SR-lowMass) or a combined \textit{di-$\tau +E^{miss}_T$} ($E^{miss}_T > 150$ GeV) trigger for the high stau mass region (SR-highMass) (\textbf{Trigger and offline cuts}). +All events are required to pass either an \textit{asymmetric di-$\tau$} trigger for the low stau mass region or a combined \textit{di-$\tau +E^{miss}_T$} ($E^{miss}_T > 150$ GeV) trigger for the high stau mass region (Trigger and offline cuts). The trigger efficiencies about $80\%$ are applied in our recasting with the following offline $p_T$ thresholds for the leading (subleading) tau lepton candidates in Table~\ref{tab:trig-eff}. \begin{table}[h!] @@ -161,14 +165,15 @@ \subsection{Event selection} \end{tabular}\label{tab:trig-eff} } \end{table} -In SR-lowMass region, values of $75 < E^{miss}_T < 150$ GeV are required for SR-lowMass to increase signal sensitivity. +In {\color{blue}low mass} region, values of $75 < E^{miss}_T < 150$ GeV are required to increase signal sensitivity. Also, two selected tau leptons are required to be tight tagged. -The selection efficiency of two taus passing the \textit{tight} working point on top of two medium tagged tau leptons is taken from official ATLAS cutflow table as a ratio of raw number of event before and after applying cut, and the tight tagging efficiency is then applied as a probability per tau lepton by a square root of the efficiency. -On the other hand, in SR-highMass region, the tight tagging efficiency is applied in the same manner as SR-lowMass region, but allowing at least one of two tau leptons passing the tight selection. +The selection efficiency of two taus passing the \textit{tight} working point on top of two {\color{blue}tagged tau leptons with the \textit{medium} working point are} taken from official ATLAS cutflow table as a ratio of raw number of event before and after applying cut. {\color{blue}The tau tagging efficiency of \textit{tight} working point} is then applied as a probability per tau lepton by a square root of the efficiency. +On the other hand, in {\color{blue}high mass} region, the tight tagging efficiency is applied in the same manner as {\color{blue}low mass} region, but allowing at least one of two tau leptons passing the tight selection. -The \textit{stransverse mass} $m_{T2}$ variable is defined as\footnote{ -Notice $m_{T2}$ calculation is done with MADANALYSIS 5 function ($PHYSICS\rightarrow Transverse\rightarrow MT2(vec1,vec2,E^{miss}_T,m_{invisible})$) -} +The \textit{stransverse mass} $m_{T2}$ variable is defined as +%\footnote{ +%Notice $m_{T2}$ calculation is done with MADANALYSIS 5 function ($PHYSICS\rightarrow Transverse\rightarrow MT2(vec1,vec2,E^{miss}_T,m_{invisible})$) +%} \begin{equation} m_{T2} =min_{\mathbf{q}_T} \left[ @@ -203,24 +208,24 @@ \section{Validation} \subsection{Event generation} In order to validate our analysis, we rely on the MSSM UFO model file~\cite{Duhr:2011se} from Feynrules model database~\cite{Alloul:2013bka}. -Two benchmark points with masses $ m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $GeV and $ (280,1) $GeV are used in this note to illustrate the validation of our reimplementation. -We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the NNPDF23LO~\cite{Martin:2009iq} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. +Two benchmark points with masses $ m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ GeV and $ (280,1) $ GeV are used in this note to illustrate the validation of our reimplementation. +We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the {\color{blue}LO set of} NNPDF23LO~\cite{Martin:2009iq} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. The PYTHIA8 version 8.244~\cite{Sjostrand:2007gs} with $A14$ tune has been used for the simulation of the parton showering and hadronization. The simulation of the detector response has been performed by using DELPHES-3.4.2~\cite{deFavereau:2013fsa}, that relies on FASTJET~\cite{Cacciari:2011ma} for object reconstruction. The modified delphes card has been used with an appropriate tuned detector card. -For example, the loosened isolation criteria are applied to cover all offline object definitions. -And the radius parameter of jet and minimum transverse momentum are lowered to 0.4 and 15 GeV with updating $b$ and tau tagging efficiencies. -Also, UniqueObjectFinder is disabled for overlap removal which is done in MADANALYSIS5. -Finally, we have used the MADANALYSIS5 reimplementation to calculate the signal selection efficiencies. +{\color{red}For example, the loosened isolation criteria are applied to cover all offline object definitions.(I don't understand Jack's comment here)} +%And the radius parameter of jet and minimum transverse momentum are lowered to 0.4 and 15 GeV with updating $b$ and tau tagging efficiencies. +%Also, UniqueObjectFinder is disabled for overlap removal which is done in MADANALYSIS5. +Finally, we have used the {\color{blue}MADANALYSIS 5 framework for recasting} the signal selection efficiencies. \subsection{Comparison with the official results} -In Table~\ref{tab:120GeV} and~\ref{tab:280GeV}, we compare the results obtained with our implememtation to the official raw event numbers in the auxiliary tables provided by the ATLAS collaboration for the benchmark points with masses $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ and $(280,1)$GeV, respectively. +In Table~\ref{tab:120GeV} and~\ref{tab:280GeV}, we compare the results obtained with our implememtation to the official raw event numbers in the auxiliary tables provided by the ATLAS collaboration for the benchmark points with masses $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ and $(280,1)$ GeV, respectively. For each cut, we have calculated the related efficiency defined as \begin{equation} \epsilon_i =\frac{n_i}{n_{i-1}} \end{equation} -where $ n_i $ and $ n_{i-1} $ mean the event number after and before the considered cut, respectively. +where $ n_i $ and $ n_{i-1} $ {\color{blue}refer to} the event number after and before the considered cut, respectively. % On the other hand, we have also calculated the differences between $ \epsilon_i (MA5)$ and $ \epsilon_i (ATLAS)$ with the definition as \begin{equation} @@ -242,7 +247,7 @@ \subsection{Comparison with the official results} Light lepton veto & 22109 & 99.82 & 25414 & 99.94 & 0.12 \\ \hline $Z/H$-veto & 18188 & 82.27 & 21082 & 82.95 & 0.83 \\ \hline % -\multicolumn{5}{c}{ \textbf{SR-lowMass} }\\\hline +\multicolumn{5}{c}{ \textbf{low mass region} }\\\hline % Trigger and offline cuts & 6512 & 35.80 & 8144 & 38.63 & 7.91 \\ \hline $ 75 < E^{miss}_T < 150 $ GeV & 2228 & 34.21 & 2840 & 34.87 & 1.93 \\ \hline @@ -253,7 +258,7 @@ \subsection{Comparison with the official results} All & & 1.24 & & 2.06 & 66.13 \\ \hline % \hline -\multicolumn{5}{c}{ \textbf{SR-highMass} }\\\hline +\multicolumn{5}{c}{ \textbf{high mass region} }\\\hline % Trigger and offline cuts & 1272 & 6.99 & 1649 & 7.82 & 11.87 \\ \hline $ \geq 1 $ tight $\tau$ & 1249 & 98.19 & 1590 & 96.42 & $-1.80$ \\ \hline @@ -279,7 +284,7 @@ \subsection{Comparison with the official results} Light lepton veto & 3888 & 99.62 & 34308 & 99.91 & 0.29 \\ \hline $Z/H$-veto & 3382 & 86.99 & 30356 & 88.48 & 1.71 \\ \hline % -\multicolumn{5}{c}{ \textbf{SR-lowMass} }\\\hline +\multicolumn{5}{c}{ \textbf{low mass region} }\\\hline % Trigger and offline cuts & 1920 & 56.77 & 17468 & 57.54 & 1.36 \\ \hline $ 75 < E^{miss}_T < 150 $ GeV & 738 & 38.44 & 6329 & 36.23 & $-5.75$ \\ \hline @@ -290,7 +295,7 @@ \subsection{Comparison with the official results} All & & 6.98 & & 6.53 & $-6.45$ \\ \hline % \hline -\multicolumn{5}{c}{ \textbf{SR-highMass} }\\\hline +\multicolumn{5}{c}{ \textbf{high mass region} }\\\hline % Trigger and offline cuts & 1096 & 32.41 & 11305 & 37.24 & 14.90 \\ \hline $ \geq 1 $ tight $\tau$ & 1076 & 98.18 & 11110 & 98.28 & $-0.10$ \\ \hline @@ -303,20 +308,21 @@ \subsection{Comparison with the official results} \end{table} -We observe that the disagreement on a cut-by-cut basis, is {\color{blue}$8.66\%$ and $7.55\%$} before the $m_{T2}$ cut for SR-lowMass and SR-highMass in Table~\ref{tab:120GeV} for $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $GeV. The major part of the disagreement comes from the \textbf{Trigger and offline cuts} step. By lack of more public experimental information, we have not been able to validate these two steps more precisely. +We observe that the disagreement on a cut-by-cut basis, is $8.66\%$ and $7.55\%$ before the $m_{T2}$ cut for {\color{blue}low mass} and {\color{blue}high mass regions} in Table~\ref{tab:120GeV} for $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ GeV. The major part of the disagreement comes from the Trigger and offline cuts step. By lack of more public experimental information, we have not been able to validate these two steps more precisely. \begin{figure}[h!] \centerline{\includegraphics[width=2.0in]{m120_norm_1}\includegraphics[width=2.0in]{m120_norm_2}} \vspace*{8pt} - \caption{The $m_{T2}$ distributions for $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $GeV.\protect\label{fig:m120_norm}} + \caption{The $m_{T2}$ distributions for $m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1)$ GeV.\protect\label{fig:m120_norm}} \end{figure} \begin{figure}[t] \centerline{\includegraphics[width=2.0in]{m280_norm_1}\includegraphics[width=2.0in]{m280_norm_2}} \vspace*{8pt} - \caption{The $m_{T2}$ distributions for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1) $GeV.\protect\label{fig:m280_norm}} + \caption{The $m_{T2}$ distributions for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1)$ GeV.\protect\label{fig:m280_norm}} \end{figure} -In Fig.~\ref{fig:m120_norm}, we realize the $m_{T2}$ distributions from ATLAS analysis are sofer than our results. This causes the $m_{T2} > 70$GeV cut looser in our reimplementation than the original ATLAS results which is taken from HEPData\cite{hepdata}. -Similarly, the disagreement on a cut-by-cut basis is {\color{blue}$-6.45\%$ and $26.15\%$} with all cuts for SR-lowMass and SR-highMass in Table~\ref{tab:280GeV} for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1) $GeV. The major parts of the disagreement come from \textbf{Trigger and offline cuts}, tight $\tau$ selection, and $m_{T2}$ cut steps. +In Fig.~\ref{fig:m120_norm}, we {\color{blue}observed} the $m_{T2}$ distributions from ATLAS analysis are sofer than our results. This causes the $m_{T2} > 70$ GeV cut looser in our reimplementation than the original ATLAS results which is taken from HEPData\cite{hepdata}. +Similarly, the disagreement on a cut-by-cut basis is $-6.45\%$ and $26.15\%$ with all cuts +for low mass and high mass regions in Table~\ref{tab:280GeV} for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1)$ GeV. The major parts of the disagreement come from Trigger and offline cuts, tight $\tau$ selection, and $m_{T2}$ cut steps. The $m_{T2}$ distributions are shown in Fig.~\ref{fig:m280_norm} for the comparison of our results with ATLAS analysis. %\begin{figure}[t] @@ -328,7 +334,7 @@ \subsection{Comparison with the official results} \section{Conclusions} We have implemented the ATLAS-SUSY-2018-04 search in the MADANALYSIS 5 framework. Our analysis has been validated in the context of a supersymmetry-inspired simplified benchmark model in which the Standard Model is extended by a neutralino and a stau decaying into a tau lepton and a neutralino, employing two different benchmark points in the parameter space. -By comparing our predictions for the cutflow with the official one provided by ATLAS in Ref.\cite{Aad:2019byo}, we have found an agreement for each step in Table~\ref{tab:120GeV} and~\ref{tab:280GeV} except for the ones from \textbf{Trigger and offline cuts}, tight $\tau$ selection, and $m_{T2}$ cut. Due to the lack of more information, we have not been able to validate these steps more precisely. +By comparing our predictions for the cutflow with the official one provided by ATLAS in Ref.\cite{Aad:2019byo}, we have found an agreement for each step in Table~\ref{tab:120GeV} and~\ref{tab:280GeV} except for the ones from Trigger and offline cuts, tight $\tau$ selection, and $m_{T2}$ cut. Due to the lack of more information, we have not been able to validate these steps more precisely. \section*{Acknowledgments} From 3f2645a1c472d19f8379668572e4e2f4e0585695 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 10:42:44 +0900 Subject: [PATCH 02/11] Correct typos --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index f6f8f4e..58376e2 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -180,7 +180,7 @@ \subsection{Event selection} max(m_{T,\tau_1}(\mathbf{p}_{T,\tau_1},\mathbf{q}_T),m_{T,\tau_2}(\mathbf{p}_{T,\tau_2},\mathbf{p}^{miss}_T -\mathbf{q}_T)) \right], \end{equation} -where $\mathbf{p}_{T,\tau_1}$ and $\mathbf{p}_{T,\tau_2}$ are the transverse momenta of the two tau lepton candidates, and the transverse momentum vector of one of the invisible particle, $\mathbf{q}_T$, is choosen to minimize the larger of the two transverse mass $m_{T,\tau_1}$ and $m_{T,\tau_2}$. The transverse mass $m_T$ is defined by +where $\mathbf{p}_{T,\tau_1}$ and $\mathbf{p}_{T,\tau_2}$ are the transverse momenta of the two tau lepton candidates, and the transverse momentum vector of one of the invisible particle, $\mathbf{q}_T$, is chosen to minimize the larger of the two transverse mass $m_{T,\tau_1}$ and $m_{T,\tau_2}$. The transverse mass $m_T$ is defined by \begin{equation} m_{T}(\mathbf{p}_T,\mathbf{q}_T) = \sqrt{2(p_T q_T -\mathbf{p}_T\cdot\mathbf{q}_T)}. \end{equation} @@ -320,7 +320,7 @@ \subsection{Comparison with the official results} \vspace*{8pt} \caption{The $m_{T2}$ distributions for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1)$ GeV.\protect\label{fig:m280_norm}} \end{figure} -In Fig.~\ref{fig:m120_norm}, we {\color{blue}observed} the $m_{T2}$ distributions from ATLAS analysis are sofer than our results. This causes the $m_{T2} > 70$ GeV cut looser in our reimplementation than the original ATLAS results which is taken from HEPData\cite{hepdata}. +In Fig.~\ref{fig:m120_norm}, we {\color{blue}observed} the $m_{T2}$ distributions from ATLAS analysis are softer than our results. This causes the $m_{T2} > 70$ GeV cut looser in our reimplementation than the original ATLAS results which is taken from HEPData\cite{hepdata}. Similarly, the disagreement on a cut-by-cut basis is $-6.45\%$ and $26.15\%$ with all cuts for low mass and high mass regions in Table~\ref{tab:280GeV} for $m(\tilde{\tau},\tilde{\chi}^0_1)=(280,1)$ GeV. The major parts of the disagreement come from Trigger and offline cuts, tight $\tau$ selection, and $m_{T2}$ cut steps. The $m_{T2}$ distributions are shown in Fig.~\ref{fig:m280_norm} for the comparison of our results with ATLAS analysis. From 69f0832c24a3da461976294f3130277a4ef54e96 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 10:50:14 +0900 Subject: [PATCH 03/11] Use LaTeX commands for log-like functions --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 58376e2..002ba61 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -100,7 +100,7 @@ \subsection{Object definitions} \end{equation} Furthermore, all electron candidates are required to have both track and calorimeter isolations. The condition of the track isolation is \begin{equation} -\sum p_{T,\textrm{tracks}}/p^e_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{ GeV}/p^e_T,0.2), +\sum p_{T,\textrm{tracks}}/p^e_T < 0.15\quad \textrm{with}\quad \Delta R=\min(10\textrm{ GeV}/p^e_T,0.2), \end{equation} the condition of the calorimeter isolation is \begin{equation} @@ -108,7 +108,7 @@ \subsection{Object definitions} \end{equation} and for high transverse momentum electron, \begin{equation} -\sum E_{T,\textrm{tracks}} < max(0.015\times p^e_T,3.5\textrm{ GeV})\quad \textrm{with}\quad \Delta R=0.2\quad \textrm{if}\quad p^e_T > 200\textrm{ GeV}. +\sum E_{T,\textrm{tracks}} < \max(0.015\times p^e_T,3.5\textrm{ GeV})\quad \textrm{with}\quad \Delta R=0.2\quad \textrm{if}\quad p^e_T > 200\textrm{ GeV}. \end{equation} Muon candidate definition is similar, although with slightly looser thresholds, @@ -117,7 +117,7 @@ \subsection{Object definitions} \end{equation} The condition of the track isolation is \begin{equation} -\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15\quad \textrm{with}\quad \Delta R=min(10\textrm{ GeV}/p^{\mu}_T,0.3), +\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15\quad \textrm{with}\quad \Delta R=\min(10\textrm{ GeV}/p^{\mu}_T,0.3), \end{equation} and the condition of the calorimeter isolation is \begin{equation} From 70f3473b4fa1e4e1102d42a13800d9b7bcb54921 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 11:07:25 +0900 Subject: [PATCH 04/11] Minor changes --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 002ba61..1393c48 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -63,7 +63,7 @@ \section{Introduction} In this note, we describe the validation of the implementation, in MADANALYSIS 5 framework~\cite{Conte:2018vmg,Dumont:2014tja,Conte:2014zja,Conte:2012fm}, of the ATLAS-SUSY-2018-04 search~\cite{Aad:2019byo} for direct stau production in two hadronic $\tau +E^{miss}_T$ events. This process is illustrated by the representative Feynman diagram of Fig.~\ref{fig:fig_01}. -This analysis focuses on LHC proton-proton collisions at a center-of-mass energy of 13 TeV and an integrated luminosity of $139 fb^{-1}$. +This analysis is based on a data sample of LHC proton-proton collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of $139 fb^{-1}$. For the validation of our reimplementation, we have focused on the sector of sparticles with only electroweak interactions. The lightest neutralino ($\tilde{\chi}^0_1$) is taken as the lightest supersymmetric particle (LSP). @@ -146,7 +146,7 @@ \subsection{Event selection} %Furthermore, $b$-jet veto is applied to reject events from top quark associated processes. The reconstructed invariant mass of the two leading tau lepton candidates, $m(\tau_1,\tau_2)$, larger than $120$ GeV is required for removing tau lepton pair from low-mass resonances, $Z$ boson, and Higgs boson events ($Z/H$ veto). %Here mostly re-arranged sentences from above paragraphs -The events with exactly two \textit{medium} tau lepton candidates with opposite-sign electric charge (OS) are selected. +The events with exactly two \textit{medium} tau lepton candidates with opposite-sign (OS) electric charge are selected. Then, $b$-jet veto is applied to reject events from top quark associated processes. Also, the events with additional light lepton (muon or electron) are rejected. The reconstructed invariant mass of the two leading {\color{blue}hadronic taus}, $m(\tau_1,\tau_2)$, larger than $120$ GeV is required for removing tau lepton pair from low-mass {\color{blue}$Z/H$ resonances}. @@ -333,7 +333,7 @@ \subsection{Comparison with the official results} \section{Conclusions} -We have implemented the ATLAS-SUSY-2018-04 search in the MADANALYSIS 5 framework. Our analysis has been validated in the context of a supersymmetry-inspired simplified benchmark model in which the Standard Model is extended by a neutralino and a stau decaying into a tau lepton and a neutralino, employing two different benchmark points in the parameter space. +We have implemented the ATLAS-SUSY-2018-04 search in the MADANALYSIS 5 framework. Our analysis has been validated in the context of a supersymmetry-inspired simplified benchmark model in which the Standard Model is extended by a neutralino and a stau decaying into a tau lepton and the neutralino, employing two different benchmark points in the parameter space. By comparing our predictions for the cutflow with the official one provided by ATLAS in Ref.\cite{Aad:2019byo}, we have found an agreement for each step in Table~\ref{tab:120GeV} and~\ref{tab:280GeV} except for the ones from Trigger and offline cuts, tight $\tau$ selection, and $m_{T2}$ cut. Due to the lack of more information, we have not been able to validate these steps more precisely. From e8ff7cc7e2f2f40a1dfd5e2500453ca35a3f58a0 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 13:20:43 +0900 Subject: [PATCH 05/11] Implement Jack's correction to NNPDF23LO --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 1393c48..10a1c0d 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -209,7 +209,7 @@ \subsection{Event generation} In order to validate our analysis, we rely on the MSSM UFO model file~\cite{Duhr:2011se} from Feynrules model database~\cite{Alloul:2013bka}. Two benchmark points with masses $ m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ GeV and $ (280,1) $ GeV are used in this note to illustrate the validation of our reimplementation. -We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the {\color{blue}LO set of} NNPDF23LO~\cite{Martin:2009iq} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. +We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the {\color{blue}LO set of} NNPDF2.3~\cite{Martin:2009iq} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. The PYTHIA8 version 8.244~\cite{Sjostrand:2007gs} with $A14$ tune has been used for the simulation of the parton showering and hadronization. The simulation of the detector response has been performed by using DELPHES-3.4.2~\cite{deFavereau:2013fsa}, that relies on FASTJET~\cite{Cacciari:2011ma} for object reconstruction. The modified delphes card has been used with an appropriate tuned detector card. {\color{red}For example, the loosened isolation criteria are applied to cover all offline object definitions.(I don't understand Jack's comment here)} From fb5269ca603558f0253fcd6b84293a24f9fec301 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 13:42:08 +0900 Subject: [PATCH 06/11] Fix reference for NNPDF2.3 See http://nnpdf.mi.infn.it/for-users/unpolarized-pdf-sets . --- .../validation_note_ATLAS_SUSY_2018_04.tex | 12 +++++++++++- 1 file changed, 11 insertions(+), 1 deletion(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 10a1c0d..07e8201 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -209,7 +209,7 @@ \subsection{Event generation} In order to validate our analysis, we rely on the MSSM UFO model file~\cite{Duhr:2011se} from Feynrules model database~\cite{Alloul:2013bka}. Two benchmark points with masses $ m(\tilde{\tau},\tilde{\chi}^0_1)=(120,1) $ GeV and $ (280,1) $ GeV are used in this note to illustrate the validation of our reimplementation. -We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the {\color{blue}LO set of} NNPDF2.3~\cite{Martin:2009iq} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. +We make use of MADGRAPH5 aMC@NLO version 2.6.7~\cite{Alwall:2014hca} for hard-scattering event generation in which leading-order matrix elements are convoluted with the {\color{blue}LO set of} NNPDF2.3~\cite{Ball:2012cx} parton distribution function (PDF) set. The signal includes the emission of up to two additional partons. We apply the MLM scheme~\cite{Mangano:2006rw,Alwall:2008qv} of the ME-PS matching with $xqcut = m_{\tilde{\tau}}/4$. The PYTHIA8 version 8.244~\cite{Sjostrand:2007gs} with $A14$ tune has been used for the simulation of the parton showering and hadronization. The simulation of the detector response has been performed by using DELPHES-3.4.2~\cite{deFavereau:2013fsa}, that relies on FASTJET~\cite{Cacciari:2011ma} for object reconstruction. The modified delphes card has been used with an appropriate tuned detector card. {\color{red}For example, the loosened isolation criteria are applied to cover all offline object definitions.(I don't understand Jack's comment here)} @@ -418,6 +418,16 @@ \section*{Acknowledgments} [arXiv:1405.0301 [hep-ph]]. %%CITATION = doi:10.1007/JHEP07(2014)079;%% %4502 citations counted in INSPIRE as of 27 Mar 2020 + +% see http://nnpdf.mi.infn.it/for-users/unpolarized-pdf-sets +%\cite{Ball:2012cx} +\bibitem{Ball:2012cx} +R.~D.~Ball, V.~Bertone, S.~Carrazza, C.~S.~Deans, L.~Del Debbio, S.~Forte, A.~Guffanti, N.~P.~Hartland, J.~I.~Latorre, J.~Rojo and M.~Ubiali, +%``Parton distributions with LHC data,'' +Nucl. Phys. B \textbf{867}, 244-289 (2013) +doi:10.1016/j.nuclphysb.2012.10.003 +[arXiv:1207.1303 [hep-ph]]. +%1584 citations counted in INSPIRE as of 28 Apr 2020 %\cite{Martin:2009iq} \bibitem{Martin:2009iq} From 7c6a04f05ad7fcc4feb9eb6a76dc985962878268 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 13:52:27 +0900 Subject: [PATCH 07/11] Fix number agreement The selection efficiency ... are ... -> The selection efficiency ... is ... --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 07e8201..9b30150 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -167,7 +167,7 @@ \subsection{Event selection} In {\color{blue}low mass} region, values of $75 < E^{miss}_T < 150$ GeV are required to increase signal sensitivity. Also, two selected tau leptons are required to be tight tagged. -The selection efficiency of two taus passing the \textit{tight} working point on top of two {\color{blue}tagged tau leptons with the \textit{medium} working point are} taken from official ATLAS cutflow table as a ratio of raw number of event before and after applying cut. {\color{blue}The tau tagging efficiency of \textit{tight} working point} is then applied as a probability per tau lepton by a square root of the efficiency. +The selection efficiency of two taus passing the \textit{tight} working point on top of two {\color{blue}tagged tau leptons with the \textit{medium} working point is} taken from official ATLAS cutflow table as a ratio of raw number of event before and after applying cut. {\color{blue}The tau tagging efficiency of \textit{tight} working point} is then applied as a probability per tau lepton by a square root of the efficiency. On the other hand, in {\color{blue}high mass} region, the tight tagging efficiency is applied in the same manner as {\color{blue}low mass} region, but allowing at least one of two tau leptons passing the tight selection. The \textit{stransverse mass} $m_{T2}$ variable is defined as From 9fd5c8ce83355958125b43d6a36714c217ae88f1 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Wed, 29 Apr 2020 13:57:36 +0900 Subject: [PATCH 08/11] Add my interpretation of Jack's comment on detector card tuning --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 1 + 1 file changed, 1 insertion(+) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 9b30150..e715e16 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -213,6 +213,7 @@ \subsection{Event generation} The PYTHIA8 version 8.244~\cite{Sjostrand:2007gs} with $A14$ tune has been used for the simulation of the parton showering and hadronization. The simulation of the detector response has been performed by using DELPHES-3.4.2~\cite{deFavereau:2013fsa}, that relies on FASTJET~\cite{Cacciari:2011ma} for object reconstruction. The modified delphes card has been used with an appropriate tuned detector card. {\color{red}For example, the loosened isolation criteria are applied to cover all offline object definitions.(I don't understand Jack's comment here)} +\textcolor{magenta}{It might be unclear what you mean by ``the loosened isolation criteria''. I do not find its definition in the text either.} %And the radius parameter of jet and minimum transverse momentum are lowered to 0.4 and 15 GeV with updating $b$ and tau tagging efficiencies. %Also, UniqueObjectFinder is disabled for overlap removal which is done in MADANALYSIS5. Finally, we have used the {\color{blue}MADANALYSIS 5 framework for recasting} the signal selection efficiencies. From 7dd7d2fda8800b8c1282eec2165199730deae1e8 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Mon, 11 May 2020 21:36:02 +0900 Subject: [PATCH 09/11] Restyle paragraph on electron candidate definition as an attempt to address Jack's comment on it: "This sounds like a shopping list". --- .../validation_note_ATLAS_SUSY_2018_04.tex | 31 +++++++++++++++---- 1 file changed, 25 insertions(+), 6 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index e715e16..b2a5d3f 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -98,18 +98,37 @@ \subsection{Object definitions} \begin{equation} p^e_T > 17 \textrm{ GeV}\quad \textrm{and}\quad |\eta^e| < 2.47. \end{equation} -Furthermore, all electron candidates are required to have both track and calorimeter isolations. The condition of the track isolation is +Furthermore, all electron candidates are required to have both track and calorimeter isolations +to reduce the number of jets misidentified as charged leptons. +For the track isolation, the scalar sum of the $p_T$ of tracks inside a variable-size cone around the electron (excluding its own track) is required to satisfy \begin{equation} -\sum p_{T,\textrm{tracks}}/p^e_T < 0.15\quad \textrm{with}\quad \Delta R=\min(10\textrm{ GeV}/p^e_T,0.2), +\sum p_{T,\textrm{tracks}}/p^e_T < 0.15 , \end{equation} -the condition of the calorimeter isolation is +with the cone size \begin{equation} -\sum E_{T,\textrm{calorimeter}}/p^e_T < 0.2\quad \textrm{with}\quad \Delta R=0.2, +\Delta R=\min(10\textrm{ GeV}/p^e_T,0.2) . \end{equation} -and for high transverse momentum electron, +For the calorimeter isolation, +the sum of the transverse energy of the calorimeter energy clusters is required to satisfy \begin{equation} -\sum E_{T,\textrm{tracks}} < \max(0.015\times p^e_T,3.5\textrm{ GeV})\quad \textrm{with}\quad \Delta R=0.2\quad \textrm{if}\quad p^e_T > 200\textrm{ GeV}. +\sum E_{T,\textrm{calorimeter}}/p^e_T < 0.2 , \end{equation} +in a cone of size +\begin{equation} +\Delta R=0.2 +\end{equation} +around the electron, excluding the energy from the electron itself. +For electrons with a high transverse momentum +$ p^e_T > 200\textrm{ GeV} $, +only an upper limit instead is imposed on the transverse energy of the calorimeter energy clusters so that +\begin{equation} +\sum E_{T,\textrm{tracks}} < \max(0.015\times p^e_T,3.5\textrm{ GeV}) , +\end{equation} +in a cone of +\begin{equation} +\Delta R=0.2 +\end{equation} +around the electron. Muon candidate definition is similar, although with slightly looser thresholds, \begin{equation} From 97c0f412c2ded151d9beefe1bd4255c9a31aefc4 Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Mon, 11 May 2020 21:41:27 +0900 Subject: [PATCH 10/11] Restyle paragraph on muon candidate definition as an attempt to address Jack's comment on it: "This sounds like a shopping list". --- .../validation_note_ATLAS_SUSY_2018_04.tex | 21 ++++++++----------- 1 file changed, 9 insertions(+), 12 deletions(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index b2a5d3f..4e5b432 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -130,18 +130,15 @@ \subsection{Object definitions} \end{equation} around the electron. -Muon candidate definition is similar, although with slightly looser thresholds, -\begin{equation} -p^{\mu}_T > 14 \textrm{ GeV}\quad \textrm{and}\quad |\eta^{\mu}| < 2.7. -\end{equation} -The condition of the track isolation is -\begin{equation} -\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15\quad \textrm{with}\quad \Delta R=\min(10\textrm{ GeV}/p^{\mu}_T,0.3), -\end{equation} -and the condition of the calorimeter isolation is -\begin{equation} -\sum E_{T,\textrm{tracks}}/p^{\mu}_T < 0.3\quad \textrm{with}\quad \Delta R=0.2. -\end{equation} +Muon candidate definition is similar to that for electrons, although with slightly looser thresholds, +$p^{\mu}_T > 14 \textrm{ GeV}$ and $|\eta^{\mu}| < 2.7$. +The condition of the track isolation is the same, i.e.\ +$\sum p_{T,\textrm{tracks}}/p^{\mu}_T < 0.15$, +with an altered cone size +$\Delta R=\min(10\textrm{ GeV}/p^{\mu}_T,0.3)$. +The condition of the calorimeter isolation is changed to +$\sum E_{T,\textrm{tracks}}/p^{\mu}_T < 0.3$ +with the same cone size $\Delta R=0.2$. Tau lepton candidates are reconstructed with one or three associated charged pion tracks (prongs). %$\sum_i e_i (\textrm{tracks}) = \pm 1$. %For 1-prong (3-prong) $\tau$ lepton candidates, the signal efficiencies are $75\%$($60\%$) and $60\%$($45\%$) for the \textit{medium} and \textit{tight} working points, respectively. From 7888d84212d10894998850550e51c3dd071b84ad Mon Sep 17 00:00:00 2001 From: Jae-hyeon Park Date: Mon, 11 May 2020 21:46:46 +0900 Subject: [PATCH 11/11] Restyle paragraph on overlap removal procedure as an attempt to address Jack's comment on it: "This sounds like a shopping list". --- validation_note/validation_note_ATLAS_SUSY_2018_04.tex | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex index 4e5b432..e325848 100644 --- a/validation_note/validation_note_ATLAS_SUSY_2018_04.tex +++ b/validation_note/validation_note_ATLAS_SUSY_2018_04.tex @@ -151,7 +151,7 @@ \subsection{Object definitions} \end{equation} for the leading (subleading) ones and the transition region between the barrel and endcap calorimeters ($ 1.37 < |\eta^{\tau}| < 1.52 $) is excluded. -Finally, some overlap removal conditions are in order which are consistent with the analysis code provided in HEPData\cite{hepdata}. Electron is removed if $\Delta R(e,e) < 0.05$. Tau lepton is removed if $\Delta R(\tau,e/\mu) < 0.2$. Electron is removed if $\Delta R(e,\mu) < 0.01$. Jet is removed if $\Delta R(j,e/\mu) < 0.2$, and then electron or muon is removed if $\Delta R(e/\mu,j) < 0.4$, and jet is removed if $\Delta R(j,\tau) < 0.4$. +Finally, some overlap removal conditions are in order which are consistent with the analysis code provided in HEPData\cite{hepdata}. Electrons are removed if they overlap so closely with another electron that $\Delta R(e,e) < 0.05$. Tau candidates are discarded if they overlap with a light lepton, i.e.\ $\Delta R(\tau,e/\mu) < 0.2$. If an electron overlaps with a muon such that $\Delta R(e,\mu) < 0.01$, the muon is retained. If a jet is found to overlap with a light lepton leading to $\Delta R(j,e/\mu) < 0.2$, the jet is discarded. For overlapping light leptons and jets with $\Delta R(e/\mu,j) < 0.4$, the jet is kept. Lastly, jets are removed if they overlap with a tau such that $\Delta R(j,\tau) < 0.4$. \subsection{Event selection}