Fork me on GitHub

Changeset 129 in svn for trunk


Ignore:
Timestamp:
Jan 4, 2009, 7:40:49 PM (16 years ago)
Author:
Xavier Rouby
Message:

validation: premier jet terminé

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/paper/notes.tex

    r128 r129  
    3737\fi
    3838
    39 \title{\textsc{Delphes}, a framework for fast simulation \\of a general purpose LHC detector}
     39\title{\textsc{Delphes}, a framework for fast simulation \\of a general purpose \textsc{lhc} detector}
    4040\author{S. Ovyn and X. Rouby\thanks{Now in Physikalisches Institut, Albert-Ludwigs-Universit\"at Freiburg} \\
    4141        Center for Particle Physics and Phenomenology (CP3)\\ Universit\'e catholique de Louvain \\ B-1348 Louvain-la-Neuve, Belgium \\ \\
     
    5959
    6060\noindent
    61 \textit{Keywords:} \textsc{Delphes}, fast simulation, LHC, smearing, trigger, \textsc{FastJet}, \textsc{Hector}, \textsc{Frog}
     61\textit{Keywords:} \textsc{Delphes}, fast simulation, \textsc{lhc}, smearing, trigger, \textsc{FastJet}, \textsc{Hector}, \textsc{Frog}
    6262\vspace{1cm}
    6363\end{abstract}
     
    224224\end{figure}
    225225
    226 The calorimetric towers directly enter in the calculation of the missing transverse energy, and as input for the jet reconstruction algorithms. No longitudinal segmentation is available in the simulated calorimeters. No sharing between neighbouring towers is implemented when particles enter a tower very close to its geometrical edge.
     226The calorimetric towers directly enter in the calculation of the missing transverse energy (\textsc{met}), and as input for the jet reconstruction algorithms. No longitudinal segmentation is available in the simulated calorimeters. No sharing between neighbouring towers is implemented when particles enter a tower very close to its geometrical edge.
    227227
    228228\subsection{Very forward detectors simulation}
     
    331331The three following jet algorithms are safe for soft radiations (\textit{infrared}) and collinear splittings. They rely on recombination schemes where neighbouring calotower pairs are successively merged. The definitions of the jet algorithms are similar except for the definition of the \textit{distances} $d$ used during the merging procedure. Two such variables are defined: the distance $d_{ij}$ between each pair of towers $(i,j)$, and a variable $d_{iB}$ (\textit{beam distance}) depending on the transverse momentum of the tower $i$.
    332332
    333 The jet reconstruction algorithm browses the calotower list. It starts by finding the minimum value $d_\textrm{min}$ of all the distances $d_{ij}$ and $d_{iB}$. If $d_\textrm{min}$ is a $d_{ij}$, the towers $i$ and $j$ are merged into a single tower with a four-momentum $p^\mu = p^\mu (i) + p^\mu (j)$ (\textit{E-scheme recombination}). If $d_\textrm{min}$ is a $d_{iB}$, the tower is declared as a final jet and is removed from the input list. This procedure is repeated until no input towers are left. Further information on these jet algorithms is given here below, using $k_{ti}$, $y_{i}$ and $\phi_i$ as the transverse momentum, rapidity and azimuth of calotower $i$ and $\Delta R_{ij}= \sqrt{(y_i-y_j)^2+(\phi_i-\phi_j)^2}$ as the jet-radius parameter:
     333The jet reconstruction algorithm browses the calotower list. It starts by finding the minimum value $d_\textrm{min}$ of all the distances $d_{ij}$ and $d_{iB}$. If $d_\textrm{min}$ is a $d_{ij}$, the towers $i$ and $j$ are merged into a \textcolor{red}{single tower with a four-momentum $p^\mu = p^\mu (i) + p^\mu (j)$ (\textit{E-scheme recombination})}. If $d_\textrm{min}$ is a $d_{iB}$, the tower is declared as a final jet and is removed from the input list. This procedure is repeated until no towers are left in the input list. Further information on these jet algorithms is given here below, using $k_{ti}$, $y_{i}$ and $\phi_i$ as the transverse momentum, rapidity and azimuth of calotower $i$ and $\Delta R_{ij}= \sqrt{(y_i-y_j)^2+(\phi_i-\phi_j)^2}$ as the jet-radius parameter:
    334334 
    335335\begin{enumerate}[start=4]
     
    396396 $ \tau^- \rightarrow \mu^- \ \bar \nu_\mu  \ \nu_\tau$ & $17.36\%$ \\
    397397 \multicolumn{2}{l}{\textbf{Hadronic decays}}\\
    398  $ \tau^- \rightarrow h^-\ n\times h^\pm \ m\times h^0\  \nu_\tau$  & $64.79\%$ \\
     398 $ \tau^- \rightarrow h^-\ n\times h^\pm \ m\times h^0\  \nu_\tau$  & \textcolor{red}{$64.79\%$} \\
    399399 $ \tau^- \rightarrow h^-\ m\times h^0 \ \nu_\tau$  & $50.15\%$ \\
    400400 $ \tau^- \rightarrow h^-\ h^+ h^-  m\times h^0 \ \nu_\tau$  & $15.18\%$ \\
     
    475475Once both electromagnetic collimation and tracking isolation are applied, a threshold on the $p_T$ of the $\tau$-jet candidate is requested to purify the collection. This procedure selects $\tau$ leptons decaying hadronically with a typical efficiency of $60\%$.
    476476
    477 \subsection{Transverse missing energy}
    478 In an ideal detector, momentum conservation imposes the transverse momentum of the observed final state $p_T^\textrm{obs}$ to be equal to the $p_T$ sum of the invisible particles, written $p_T^\textrm{miss}$.
    479 \begin{equation}
    480  p_T^\textrm{miss} = - p_T^\textrm{obs}
    481 \end{equation}
    482 In a real experiment, energy is measured by calorimetry and any problem affecting the detector (dead channels, misalignment, noisy towers, cracks) worsens directly the measured missing transverse energy $E_T^\textrm{miss}$. In this document, missing transverse energy is based on the calorimetric towers and only muons and neutrinos are not taken into account for its evaluation:
    483 \begin{equation}
    484  E_T^\textrm{miss} = - \sum^\textrm{towers}_i E_T(i)
    485 \end{equation}
     477\subsection{Missing transverse energy}
     478In an ideal detector, momentum conservation imposes the transverse momentum of the observed final state $\overrightarrow{p_T}^\textrm{obs}$ to be equal to the $\overrightarrow{p_T}$ vector sum of the invisible particles, written $\overrightarrow{p_T}^\textrm{miss}$.
     479\begin{equation}
     480\overrightarrow{p_T} = \left(
     481\begin{array}{c}
     482p_x\\
     483p_y\\
     484\end{array}
     485\right)
     486~ \textrm{and} ~
     487\left\{
     488\begin{array}{l}
     489 p_x^\textrm{miss} = - p_x^\textrm{obs} \\
     490 p_y^\textrm{miss} = - p_y^\textrm{obs} \\
     491\end{array}
     492\right.
     493\end{equation}
     494The \textit{true} missing transverse energy, i.e. at generator-level, is calculated as the opposite of the vector sum of the transverse momenta of all visible particles -- or equivalently, to the vector sum of invisible particle transverse momenta.
     495In a real experiment, calorimeters measure energy and not momentum. Any problem affecting the detector (dead channels, misalignment, noisy towers, cracks) worsens directly the measured missing transverse energy $\overrightarrow {E_T}^\textrm{miss}$. In this document, \textsc{met} is based on the calorimetric towers and only muons and neutrinos are not taken into account for its evaluation:
     496\begin{equation}
     497\textcolor{red}{ \overrightarrow{E_T}^\textrm{miss} = - \sum^\textrm{towers}_i \overrightarrow{E_T}(i)}
     498\end{equation}
     499
    486500
    487501\section{Trigger emulation}
    488502
    489 New physics in collider experiment are often characterised by the phenomenology by low cross-section values. High statistics are required for their studies, which in turn imposes high luminosity collisions.
    490 
    491 On the other hand, due to the very high collision rate in recent collider ($40~\textrm{MHz}$ at the \textsc{lhc}) and the large total cross-section ($\mathcal{O}(110~\textrm{mb})$ at the \textsc{lhc}), the need for an online event selection is crucial in order to reject most of the event and keep
    492 
    493 \section{Validation}   
     503New physics in collider experiment are often characterised in phenomenology by low cross-section values, compared to the Standard Model (\textsc{sm}) processes. For instance at the \textsc{lhc} ($\sqrt{s}=14~\textrm{TeV}$), the cross-section of inclusive production of $b \bar b$ pairs is expected to be $10^7~\textrm{nb}$, or inclusive jets at $100~\textrm{nb}$ ($p_T > 200~\textrm{GeV}$), while \textcolor{red}{Higgs boson cross-section within the \textsc{sm} can be as small as $\ldots \times 10^{-6}~\textrm{nb}$}.
     504
     505High statistics are required for data analyses, consequently imposing high luminosity, i.e. a high collision rate.
     506As only a tiny fraction of the observed events can be stored for subsequent \textit{offline} analyses, a very large data rejection factor should be applied directly as the events are produced.
     507This data selection is supposed to reject only well-known \textsc{sm} events\footnote{However, some bandwidth is allocated to random triggers that stores a small fraction of the events without any selection criteria.}.
     508Dedicated algorithms of this \textit{online} selection, or \textit{trigger}, should be fast and very efficient for data rejection, in order to preserve the experiment output bandwidth. They must also be as inclusive as possible to avoid loosing interesting events.
     509
     510Most of the usual trigger algorithms select events containing objects (i.e. jets, particles, \textsc{met}) with an energy scale above some threshold. This is often expressed in terms of a cut on the transverse momentum of one or several objects of the measured event. Logical combinations of several conditions are also possible. For instance, a trigger path could select events containing at least one jet and one electron such as $p_T^\textrm{jet} > 100~\textrm{GeV}$ and $p_T^e > 50~\textrm{GeV}$.
     511
     512A trigger emulation is included in \textsc{Delphes}, using a fully parametrisable \textit{trigger table}\footnote{\texttt{[code] }The trigger card is the \texttt{data/trigger.dat} file.}. When enabled, this trigger is applied on analysis object data.
     513In a real experiment, the online selection is often divided into several steps (or \textit{levels}).
     514This splits the overall reduction factor into a product of smaller factors, corresponding to the different trigger levels.
     515This is related to the architecture of the experiment data acquisition chain, with limited electronic buffers requiring a quick decision for the first trigger level.
     516First level triggers are then fast and simple but based only on partial data as not all detector front-ends are readable within the decision latency.
     517Later levels are more complex, of finer-but-not-final quality and based on full detector data.
     518
     519Real triggers are thus intrinsically based on reconstructed data with a worse resolution than final analysis data.
     520On the contrary, same data are used in \textsc{Delphes} for trigger emulation and for final analyses.
     521
     522\section{Validation}
     523
     524\textsc{Delphes} performs a fast simulation of a collider experiment.
     525Its quality and validity are assessed by comparing to resolution of the reconstructed data to the \textsc{cms} detector expectations.
     526
     527Electrons and muons match by construction to the experiment designs, as the Gaussian smearing of their kinematical properties is defined according to the experiment resolution.
     528Similarly, the $b$-tagging efficiency (for real $b$-jets) and misidentification rates (for fake $b$-jets) are taken from the expected values of the experiment.
     529Unlike these simple objects, jets and missing transverse energy should be carefully cross-checked.
    494530
    495531\subsection{Jet resolution}
    496532 
    497 The majority of interesting processes at the \textsc{lhc} contain jets in the final state. The jet resolution obtained using \textsc{Delphes} is therefore a crucial point of the validation. While \textsc{Delphes} contains six jet reconstruction algorithms, only the jet clustering algorithm with $R=0.7$ is used to validate the jet collection. Cross-check has been made with the results obtained using the \textsc{cms} detector. This validation employs $pp \rightarrow gg$ events produced using \textsc{mg/me} and hadronized using \textsc{pythia}. The events were divided into 14 bins of $\hat{p_T}$ of the gluons. Each \textsc{Delphes} jet is matched to the closest {\it particle-level} jet using the spatial separation in $\eta - \phi$ between the two jet axis $\Delta R<0.25$, otherwise they are discarded. The particle-level jets are obtained by applying the same clustering algorithm to all particles considered as stable by \textsc{pythia}.
    498  
    499 For each $\hat{p}_T$ bin, the  \textsc{Delphes} jet transverse energy ($E_T^{rec}$) of all jets satisfying the matching criteria is compaired to the {\it particle level} transverse energy ($E_T^{MC}$). The obtained histograms of the $E_T^{rec}/E_T^{MC}$ response have been fitted with a Gaussian function in the interval $\pm 2.RMS$ centered around the mean value. The final jet resolution is obtained using the following formula:
    500  
    501 \begin{equation}
    502 \frac{\sigma(R_{jet})}{<R_{jet}>}=\frac{\sigma(\frac{E_T^{rec}}{E_T^{MC}})_{fit}}{<\frac{E_T^{rec}}{E_T^{MC}}>_{fit}}.
     533The majority of interesting processes at the \textsc{lhc} contain jets in the final state. The jet resolution obtained using \textsc{Delphes} is therefore a crucial point for its validation. Even if \textsc{Delphes} contains six algorithms for jet reconstruction, only the jet clustering algorithm (\textsc{jetclu}) with $R=0.7$ is used to validate the jet collection.
     534
     535This validation \textcolor{red}{employs} $pp \rightarrow gg$ events produced with \textsc{mg/me} and hadronised using \textsc{Pythia}~\cite{bib:mgme,bib:pythia}. The events were arranged in $14$ bins of gluon transverse momentum $\hat{p}_T$. In each $\hat{p}_T$ bin, every jet in \textsc{Delphes} is matched to the closest jet of generator-level particles, using the spatial separation between the two jet \textcolor{red}{axes}
     536\begin{equation}
     537\Delta R = \sqrt{ \big(\eta^\textrm{rec} - \eta^\textrm{MC} \big)^2 +  \big(\phi^\textrm{rec} - \phi^\textrm{MC} \big)^2}<0.25.
     538\end{equation}
     539The jets made of generator-level particles, or \textsc{mc} jets, are obtained by applying the same clustering algorithm to all particles considered as stable after hadronisation.
     540Jets produced by \textsc{Delphes} and satisfying the matching criterium are called hereafter \textit{reconstructed jets}.
     541
     542The ratio of the transverse energies of every reconstructed jet $E_T^\textrm{rec}$ and its corresponding \textsc{mc} jet $E_T^\textrm{MC}$ is calculated in each $\hat{p}_T$ bin.
     543The $E_T^\textrm{rec}/E_T^\textrm{MC}$ histogram is fitted with a Gaussian distribution in the interval \mbox{$\pm 2$~\textsc{rms}} centered around the mean value.
     544The resolution in each $\hat{p}_T$ bin is obtained by the fit mean $\langle x \rangle$ and variance $\sigma^2(x)$:
     545\begin{equation}
     546%\frac{\sigma(R_{jet})}{\langle R_{jet} \rangle }=
     547\frac{\sigma \Big (\frac{E_T^{rec}}{E_T^{MC}} \Big)_\textrm{fit}}{ \Big \langle \frac{E_T^{rec}}{E_T^{MC}} \Big \rangle_\textrm{fit}}~
     548\Big( \hat{p}_T(i) \Big)\textrm{, for all }i.
    503549\end{equation}
    504550 
     
    506552\begin{center}
    507553\includegraphics[width=\columnwidth]{resolutionJet}
    508 \caption{Distribution of the jet transverse energy resolution as a function of the {\it particle-level}  jet transverse energy. The maximum allowed separation between the \textsc{Delphes} and the {\it partile-level} jets is $\Delta R<0.25$.}
     554\caption{Resolution of the transverse energy of reconstructed jets $E_T^\textrm{rec}$ as a function of the transverse energy of the closest jet of generator-level particles $E_T^\textrm{MC}$. The maximum separation between the reconstructed and \textsc{mc} jets is $\Delta R= 0.25$. Pink line is the fit result for comparison to the \textsc{cms} resolution, in blue.}
    509555\label{fig:jetresol}
    510556\end{center}
    511557\end{figure}
    512558 
    513 The resulting jet resolution, plotted as a function of $E_T^{GEN}$ is shown in figure \ref{fig:jetresol}. The plots were then fitted with a function of the following form:
    514  
    515 \begin{equation}
    516 \frac{a}{E_T^{GEN}}\oplus \frac{b}{\sqrt{E_T^{GEN}}}\oplus c,
    517 \end{equation}
    518  
    519 where a, b, and c are the fit parameters. The obtained resolution is compared to the one obtained with a recent version of the simulation package of the CMS detector. Overall, the resolution curve of \textsc{Delphes} matches relatively well to those of \textsc{cms}.
    520  
    521 \subsection{$E_T^{mis}$ resolution}
    522  
    523 Because all major detectors at hadron colliders have been designed to be as mutch hermetic as possible in order to detect the presence of one or more neutrinos through apparent missing transverse energy, the resolution of the $E_T^{miss}$ obtained with \textsc{Delphes} is a crucial point. The samples used to study the transverse missing energy performance are identical to those used for the jet validation. The {\it particle-level} true transverse missing energy is calculated as the vector sum of the transverse momenta of all visible particles (or equivalently, to the vector sum of invisible particles). It should be noticed that the contribution to the transverse missing energy from muons is negligeable in the sample we are interested in.
    524  
    525 In order to obtain the x-component missing energy resolution ($E_x^{miss}$), the distribution of the difference between the \textsc{Delphes} and the {\it particle-level} $E_x^{miss}$ has been fitted with a Gaussian function. The resulting $E_x^{mis}$ is plotted in figure \ref{fig:resolETmis} as a function of the total visible transverse energy, defined as the scalar sum of transverse energy in all towers ($\Sigma E_T$).
    526  
    527 \begin{figure}[!h]
    528 \begin{center}
    529 \includegraphics[width=\columnwidth]{figures/resolutionETmis}
    530 \caption{$\sigma(E^{miss}_{x})$ as a function on the scalar sum of all towers ($\Sigma E_T$) for $pp \rightarrow gg$ events.}
     559The resulting jet resolution as a function of $E_T^\textrm{MC}$ is shown in Fig.~\ref{fig:jetresol}.
     560This distribution is fitted with a function of the following form:
     561\begin{equation}
     562\frac{a}{E_T^\textrm{MC}}\oplus \frac{b}{\sqrt{E_T^\textrm{MC}}}\oplus c,
     563\end{equation}
     564where $a$, $b$ and $c$ are the fit parameters.
     565It is then compared to the resolution obtained with a recent version of the simulation package of the \textsc{cms} detector~\cite{bib:cmsjetresolution}. The resolution curves from \textsc{Delphes} and \textsc{cms} are in good agreement.
     566 
     567\subsection{MET resolution}
     568 
     569All major detectors at hadron colliders have been designed to be as much hermetic as possible in order to detect the presence of one or more neutrinos through apparent missing transverse energy.
     570The resolution of the $\overrightarrow{E_T}^\textrm{miss}$ variable, as obtained with \textsc{Delphes}, is then crucial.
     571
     572The samples used to study the \textsc{met} performance are identical to those used for the jet validation.
     573It is worth noting that the contribution to $E_T^\textrm{miss}$ from muons is negligible in the studied sample.
     574\textcolor{red}{The\footnote{je n'ai pas tout compris. Ce que j'ai devin\'e est en rouge.} input samples are divided in five bins of scalar $E_T$ sums $(\Sigma E_T)$. This sum, called \textit{total visible transverse energy}, is defined as the scalar sum of transverse energy in all towers.}
     575The quality of the \textsc{met} reconstruction is checked via the resolution on its horizontal component $E_x^\textrm{miss}$.
     576
     577The $E_x^\textrm{miss}$ resolution is evaluated in the following way.
     578The distribution of the difference between $E_x^\textrm{miss}$ in \textsc{Delphes} and at generator-level is fitted with a Gaussian function \textcolor{red}{in each $(\Sigma E_T)$ bin. The fit mean gives the \textsc{met} bias in each bin.
     579The resulting value} is plotted in Fig.~\ref{fig:resolETmis} as a function of the total visible transverse
     580energy.\footnote{
     581\textcolor{red}{Entre nous, ca ressemble plus \`a un biais (= une diff\'erence entre le vrai et le simul\'e) plus qu'a une r\'esolution! Mais je suppose que c'est la definition que tu as trouv\'ee dans le CMS TDR.}}
     582 
     583\begin{figure}[!h]
     584\begin{center}
     585\includegraphics[width=\columnwidth]{resolutionETmis}
     586\caption{$\sigma(E^\textrm{miss}_{x})$ as a function on the scalar sum of all towers ($\Sigma E_T$) for $pp \rightarrow gg$ events.}
    531587\label{fig:resolETmis}
    532588\end{center}
    533589\end{figure}
    534590 
    535 The resolution is observed to follow the form
    536 \begin{equation}
    537 \sigma_X = \alpha ~\Sigma E_T ~\mathrm{GeV}^{1/2},
    538 \end{equation}
    539 whith $\alpha$ is depending on the resolution of the calorimeters. Knowing that the expected transverse missing energy resolution expected using the \textsc{cms} detector for similar events is $\sigma_X = (0.6-0.7) ~ \Sigma E_T ~ \mathrm{GeV}^{1/2}$ with no pile-up (no extra simultaneous $pp$ collision occuring at the same bunch crossing), we can conclude that the resolution obtained by \textsc{Delphes} ( $\sigma_X = 0.68~ \Sigma E_T ~\mathrm{GeV}^{1/2}$) is in excellent agreement with the expectations of a general purpose detector.
    540 
    541 \subsection{$tau$-jet efficiency}
    542 with an efficiciency of about $50\%$ for the $\tau$-jets in CMS~\cite{bib:cmstauresolution}.
     591The resolution $\sigma_x$ of the horizontal component of \textsc{met} is observed to behave like
     592\begin{equation}
     593\sigma_x = \alpha ~(\Sigma E_T) ~~~(\mathrm{GeV}^{1/2}),
     594\end{equation}
     595where the $\alpha$ parameter is depending on the resolution of the calorimeters.
     596
     597The \textsc{met} resolution expected for the \textsc{cms} detector for similar events is $\sigma_x = (0.6-0.7) ~ (\Sigma E_T) ~ \mathrm{GeV}^{1/2}$ with no pile-up\footnote{\textit{Pile-up} events are extra simultaneous $pp$ collision occuring at the same bunch crossing.}~\cite{bib:cmsjetresolution}.
     598The same quantity obtained by \textsc{Delphes} is in excellent agreement with the expectations of the general purpose detector, as $\alpha = 0.68$.
     599
     600\subsection{$\tau$-jet efficiency}
     601Due to the complexity of their reconstruction algorithm, $\tau$-jets have also to be checked.
     602Table~\ref{tab:taurecoefficiency} lists the reconstruction efficiencies for the hadronic $\tau$-jets in the \textsc{cms} experiment and in \textsc{Delphes}. Agreement is good enough to validate this reconstruction.
     603
     604~\cite{bib:cmstauresolution}.
     605
     606\begin{table}[!h]
     607\begin{center}
     608\caption{Reconstruction efficiencies of $\tau$-jets in decays from $Z$ or $H$ bosons.\vspace{0.5cm}}
     609\begin{tabular}{lll}
     610\hline
     611\multicolumn{2}{c}{\textsc{cms}} & \\
     612$Z \rightarrow \tau^+ \tau^-$                    & $38 \%$ &  \\
     613$H \rightarrow \tau^+ \tau^-$ & $36 \%$ & $m_H = 150~\textrm{GeV}$ \\
     614$H \rightarrow \tau^+ \tau^-$ & $47 \%$ & $m_H = 300~\textrm{GeV}$ \\
     615\multicolumn{2}{c}{\textsc{Delphes}} & \\
     616$H \rightarrow \tau^+ \tau^-$ &$42 \%$  & $m_H = 140~\textrm{GeV}$ \\
     617\hline
     618\end{tabular}
     619\label{tab:taurecoefficiency}
     620\end{center}
     621\end{table}
     622
    543623
    544624\section{Visualisation}
     
    587667\bibitem{bib:Hector} \textsc{Hector},
    588668\bibitem{bib:Frog} \textsc{Frog},
    589 \bibitem{bib:CMSresolution} CMS IN 2007/053
     669\bibitem{bib:cmsjetresolution} CMS IN 2007/053
    590670\bibitem{bib:Root} \textsc{Root} - An Object Oriented Data Analysis Framework, R. Brun and F. Rademakers, Nucl. Inst. \& Meth. in Phys. Res. A 389 (1997) 81-86, \url{http://root.cern.ch}
    591671\bibitem{bib:cmstaus} Tau reconstruction in CMS
    592672\bibitem{bib:whphotoproduction} WH photoproduction, S. Ovyn
    593 \bibitem{bib:mgme} Madgraph/Madevent
     673\bibitem{bib:mgme} Madgraph/Madevent version xx.yy
     674\bibitem{bib:pythia} \textsc{Pythia} version xx.yy
    594675\bibitem{bib:pdg} C. Amsler et al. (Particle Data Group), PL B667, 1 (2008) (URL: http://pdg.lbl.gov)
    595676\bibitem{bib:cmstauresolution} R. Kinnunen, \textit{Study of $\tau$-jet identification in CMS}, CMS NOTE 1997/002.
Note: See TracChangeset for help on using the changeset viewer.