- Timestamp:
- Jan 4, 2009, 7:40:49 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/paper/notes.tex
r128 r129 37 37 \fi 38 38 39 \title{\textsc{Delphes}, a framework for fast simulation \\of a general purpose LHCdetector}39 \title{\textsc{Delphes}, a framework for fast simulation \\of a general purpose \textsc{lhc} detector} 40 40 \author{S. Ovyn and X. Rouby\thanks{Now in Physikalisches Institut, Albert-Ludwigs-Universit\"at Freiburg} \\ 41 41 Center for Particle Physics and Phenomenology (CP3)\\ Universit\'e catholique de Louvain \\ B-1348 Louvain-la-Neuve, Belgium \\ \\ … … 59 59 60 60 \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} 62 62 \vspace{1cm} 63 63 \end{abstract} … … 224 224 \end{figure} 225 225 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.226 The 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. 227 227 228 228 \subsection{Very forward detectors simulation} … … 331 331 The 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$. 332 332 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: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 \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: 334 334 335 335 \begin{enumerate}[start=4] … … 396 396 $ \tau^- \rightarrow \mu^- \ \bar \nu_\mu \ \nu_\tau$ & $17.36\%$ \\ 397 397 \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\%$} \\ 399 399 $ \tau^- \rightarrow h^-\ m\times h^0 \ \nu_\tau$ & $50.15\%$ \\ 400 400 $ \tau^- \rightarrow h^-\ h^+ h^- m\times h^0 \ \nu_\tau$ & $15.18\%$ \\ … … 475 475 Once 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\%$. 476 476 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} 478 In 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} 482 p_x\\ 483 p_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} 494 The \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. 495 In 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 486 500 487 501 \section{Trigger emulation} 488 502 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} 503 New 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 505 High statistics are required for data analyses, consequently imposing high luminosity, i.e. a high collision rate. 506 As 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. 507 This 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.}. 508 Dedicated 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 510 Most 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 512 A 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. 513 In a real experiment, the online selection is often divided into several steps (or \textit{levels}). 514 This splits the overall reduction factor into a product of smaller factors, corresponding to the different trigger levels. 515 This is related to the architecture of the experiment data acquisition chain, with limited electronic buffers requiring a quick decision for the first trigger level. 516 First 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. 517 Later levels are more complex, of finer-but-not-final quality and based on full detector data. 518 519 Real triggers are thus intrinsically based on reconstructed data with a worse resolution than final analysis data. 520 On 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. 525 Its quality and validity are assessed by comparing to resolution of the reconstructed data to the \textsc{cms} detector expectations. 526 527 Electrons and muons match by construction to the experiment designs, as the Gaussian smearing of their kinematical properties is defined according to the experiment resolution. 528 Similarly, the $b$-tagging efficiency (for real $b$-jets) and misidentification rates (for fake $b$-jets) are taken from the expected values of the experiment. 529 Unlike these simple objects, jets and missing transverse energy should be carefully cross-checked. 494 530 495 531 \subsection{Jet resolution} 496 532 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}}. 533 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 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 535 This 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} 539 The 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. 540 Jets produced by \textsc{Delphes} and satisfying the matching criterium are called hereafter \textit{reconstructed jets}. 541 542 The 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. 543 The $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. 544 The 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. 503 549 \end{equation} 504 550 … … 506 552 \begin{center} 507 553 \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.} 509 555 \label{fig:jetresol} 510 556 \end{center} 511 557 \end{figure} 512 558 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.} 559 The resulting jet resolution as a function of $E_T^\textrm{MC}$ is shown in Fig.~\ref{fig:jetresol}. 560 This 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} 564 where $a$, $b$ and $c$ are the fit parameters. 565 It 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 569 All 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. 570 The resolution of the $\overrightarrow{E_T}^\textrm{miss}$ variable, as obtained with \textsc{Delphes}, is then crucial. 571 572 The samples used to study the \textsc{met} performance are identical to those used for the jet validation. 573 It 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.} 575 The quality of the \textsc{met} reconstruction is checked via the resolution on its horizontal component $E_x^\textrm{miss}$. 576 577 The $E_x^\textrm{miss}$ resolution is evaluated in the following way. 578 The 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. 579 The resulting value} is plotted in Fig.~\ref{fig:resolETmis} as a function of the total visible transverse 580 energy.\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.} 531 587 \label{fig:resolETmis} 532 588 \end{center} 533 589 \end{figure} 534 590 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}. 591 The 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} 595 where the $\alpha$ parameter is depending on the resolution of the calorimeters. 596 597 The \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}. 598 The 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} 601 Due to the complexity of their reconstruction algorithm, $\tau$-jets have also to be checked. 602 Table~\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 543 623 544 624 \section{Visualisation} … … 587 667 \bibitem{bib:Hector} \textsc{Hector}, 588 668 \bibitem{bib:Frog} \textsc{Frog}, 589 \bibitem{bib: CMSresolution} CMS IN 2007/053669 \bibitem{bib:cmsjetresolution} CMS IN 2007/053 590 670 \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} 591 671 \bibitem{bib:cmstaus} Tau reconstruction in CMS 592 672 \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 594 675 \bibitem{bib:pdg} C. Amsler et al. (Particle Data Group), PL B667, 1 (2008) (URL: http://pdg.lbl.gov) 595 676 \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.