Fork me on GitHub

source: svn/trunk/paper/notes.tex@ 255

Last change on this file since 255 was 208, checked in by Xavier Rouby, 16 years ago

updating bibliography

File size: 82.5 KB
Line 
1\documentclass[a4paper,11pt,oneside,twocolumn]{article}
2%\usepackage[english]{babel}
3\usepackage[ansinew]{inputenc}
4\usepackage{abstract}
5
6\usepackage{amsmath}
7\usepackage{epic}
8\usepackage{wrapfig}
9\usepackage{eepic}
10\usepackage{color}
11\usepackage{latexsym}
12\usepackage{array}
13\usepackage{multicol}
14
15\usepackage{fancyhdr}
16\usepackage{verbatim}
17\addtolength{\textwidth}{1cm} \addtolength{\hoffset}{-0.5cm}
18\usepackage[colorlinks=true, pdfstartview=FitV, linkcolor=blue, citecolor=blue, urlcolor=blue, unicode]{hyperref}
19\usepackage{ifpdf}
20\usepackage{cite}
21
22\usepackage{enumitem}
23
24\newcommand{\dollar}{\$}
25
26\ifpdf
27 \usepackage[pdftex]{graphicx}
28 \graphicspath{{all_png/}}
29 \pdfinfo{
30 /Author (S. Ovyn, X. Rouby)
31 /Title (Delphes, a framework for fast simulation of a general purpose LHC detector)
32 /Subject ()
33 /Keywords (Delphes, Fast simulation, LHC, FROG, Hector, Smearing, FastJet)}
34\else
35 \usepackage[dvips]{graphicx}
36 \graphicspath{{figures/}}
37\fi
38
39%\title{\textsc{Delphes}, a framework for fast simulation \\of a general purpose \textsc{lhc} detector}
40\title{\textsc{Delphes}, a framework for fast simulation \\of a generic collider experiment}
41\author{S. Ovyn and X. Rouby$^\textrm{a}$\\
42 \small{Center for Particle Physics and Phenomenology (CP3)}\\
43 \small{Universit\'e catholique de Louvain}\\
44 \small{B-1348 Louvain-la-Neuve, Belgium}\\ \\
45 \texttt{severine.ovyn@uclouvain.be, xavier.rouby@cern.ch} \\
46}
47\date{}
48
49
50\begin{document}
51
52\twocolumn[
53\maketitle
54
55\begin{center}
56\includegraphics{DelphesLogoSml}
57\end{center}
58
59
60\begin{abstract}
61It is always delicate to know whether theoretical predictions are visible and measurable in a high energy experiment due to the complexity of the related detectors, data acquisition chain and software.
62%Knowing whether theoretical predictions are visible and measurable in a high energy experiment is always delicate due to the complexity of the related detectors, data acquisition chain and software.
63We introduce here a new framework, \textsc{Delphes}, for fast simulation of
64a general purpose experiment. The simulation includes a tracking system, embedded into a magnetic field, calorimetry and a muon
65system, and possible very forward detectors arranged along the beamline.
66The framework is interfaced to standard file formats (e.g. Les Houches Event File) and outputs observable objects for analysis, like missing transverse energy and collections of electrons or jets.
67The simulation of detector response takes into account the detector resolution, and usual reconstruction algorithms, like \textsc{FastJet}. A simplified preselection can also be applied on processed data for trigger emulation. Detection of very forward scattered particles relies on the transport in beamlines with the \textsc{Hector} software. Finally, the \textsc{Frog} 2D/3D event display is used for visualisation of the collision final states.
68An overview of \textsc{Delphes} is given as well as a few use-cases for illustration.
69\vspace{0.5cm}
70
71\noindent
72\textit{Keywords:} \textsc{Delphes}, fast simulation, \textsc{lhc}, smearing, trigger, \textsc{FastJet}, \textsc{Hector}, \textsc{Frog}\\
73\href{http://www.fynu.ucl.ac.be/delphes.html}{http://www.fynu.ucl.ac.be/delphes.html}\\
74\textit{Preprint:} \texttt{CP3-09-01}
75\vspace{1.5cm}
76
77\end{abstract}
78\small{$^\textrm{a}$ Now in Physikalisches Institut, Albert-Ludwigs-Universit\"at Freiburg}
79]
80%\saythanks
81
82\section{Introduction}
83% Motiver l'utilisation d'un simulateur rapide
84% - 1) rapide VS lent
85% - 2) relativement bonne prédiction en premiÚre approximation
86% - 3) permet de comparer
87
88Experiments at high energy colliders are very complex systems for several reasons. Firstly, in terms of the various detector subsystems, including tracking, central calorimetry, forward calorimetry, and muon chambers. These detectors differ with their principles, technologies, geometries and sensitivities. Secondly, due to the requirement of a highly effective online selection (i.e. a \textit{trigger}), subdivided into several levels for an optimal reduction factor, but based only on partially processed data. Finally, in terms of the experiment software, with different data formats (like \textit{raw} or \textit{reconstructed} data), many reconstruction algorithms and particle identification schemes.
89
90This complexity is handled by large collaborations of thousands of people, but the data and the expertise are only available to their members. Real data analyses require a full detector simulation, including the various detector inefficiencies, the dead material, the imperfections and the geometrical details. Moreover, detector calibration and alignment are crucial. Such simulation is very complicated, technical and requires a large \texttt{CPU} power. On the other hand, phenomenological studies, looking for the observability of given signals, may require only fast but realistic estimates of the observables.
91
92A new framework, called \textsc{Delphes}~\cite{bib:Delphes}, is introduced here, for the fast simulation of a general purpose collider experiment.
93Using the framework, observables can be estimated for specific signal and background channels, as well as their production and measurement rates.
94Starting from the output of event generators, the simulation of the detector response takes into account the subdetector resolutions, by smearing the kinematic properties of the final state particles\footnote{throughout the paper, final state particles refer as particles considered as stable by the event generator.}. Tracks of charged particles and deposits of energy in calorimetric cells (or \textit{calotowers}) are then created.
95
96\textsc{Delphes} includes the most crucial experimental features, like
97\begin{enumerate}
98\item the geometry of both central or forward detectors,
99\item reconstruction of photons, leptons, jets, $b$-jets, $\tau$-jets and missing transverse energy,
100\item lepton isolation,
101\item trigger emulation,
102\item an event display (Fig.~\ref{fig:FlowChart}, at the end).
103\end{enumerate}
104
105\begin{figure*}[t]
106\begin{center}
107%\includegraphics[width=0.9\textwidth]{FlowDelphes}
108\includegraphics[scale=0.78]{FlowDelphes}
109\caption{Flow chart describing the principles behind \textsc{Delphes}. Event files coming from external Monte Carlo generators are read by a converter stage.
110The kinematics variables of the final state particles are then smeared according to the subdetector resolutions.
111Tracks are reconstructed in a simulated dipolar magnetic field and calorimetric towers sample the energy deposits. Based on these, dedicated algorithms are applied for particle identification, isolation and reconstruction.
112The transport of very forward particle to the near-beam detectors is also simulated.
113Finally, an output file is written, including generator level and analysis object data. If requested, a fully parametrisable trigger can be emulated. Optionally, the geometry and visualisation files for the 3D event display can also be produced.
114All user parameters are set in the \textit{Smearing Card} and the \textit{Trigger Card}. }
115\label{fig:FlowChart}
116\end{center}
117\end{figure*}
118
119Although this kind of approach yields much realistic results than a simple ``parton-level" analysis, a fast simulation comes with some limitations. Detector geometry is idealised, being uniform, symmetric around the beam axis, and having no cracks nor dead material. Secondary interactions, multiple scatterings, photon conversion and bremsstrahlung are also neglected.
120
121%The simulation package proceeds in two stages. The first part is executed on the generated events. ``Particle-level" informations are read from input files and stored in a {\it \textsc{gen}} \textsc{root} tree.
122
123Three formats of input files can be used as input in \textsc{Delphes}\footnote{\texttt{[code] }See the \texttt{HEPEVTConverter}, \texttt{LHEFConverter} and \texttt{STDHEPConverter} classes.}. In order to process events from many different generators, the standard Monte Carlo event structure \texttt{StdHEP}~\cite{bib:stdhep} can be used as an input. Besides, \textsc{Delphes} can also provide detector response for events read in ``Les Houches Event Format'' (\textsc{lhef}~\cite{bib:lhe}) and \textsc{root} files obtained using the \texttt{h2root} utility from the \textsc{root} framework~\cite{bib:Root}.
124%Afterwards, \textsc{Delphes} performs a simple trigger simulation and reconstruct "high-level objects". These informations are organised in classes and each objects are ordered with respect to the transverse momentum.
125
126\textsc{Delphes} uses the \texttt{ExRootAnalysis} utility~\cite{bib:ExRootAnalysis} to create output data in a \texttt{*.root} ntuple.
127This output contains a copy of the generator level data (\textsc{gen} tree), the analysis data objects after reconstruction (\mbox{\textsc{A}nalysis} tree), and possibly the results of the trigger emulation (\mbox{\textsc{T}rigger} tree). The program is driven by input cards. The detector card (\texttt{data/DataCardDet.dat}) allows a large spectrum of running conditions by modifying basic detector parameters, including calorimeter and tracking coverage and resolution, thresholds or jet algorithm parameters. The trigger card (\texttt{data/trigger.dat}) lists the user algorithms for the simplified online preselection.\\
128
129
130\section{Detector simulation}
131
132The overall layout of the general purpose detector simulated by \textsc{Delphes} is shown in Fig.~\ref{fig:GenDet3}.
133A central tracking system (\textsc{tracker}) is surrounded by an electromagnetic and a hadron calorimeters (\textsc{ecal} and \textsc{hcal}, resp.). Two forward calorimeters (\textsc{fcal}) ensure a larger geometric coverage for the measurement of the missing transverse energy. Finally, a muon system (\textsc{muon}) encloses the central detector volume
134The fast simulation of the detector response takes into account geometrical acceptance of sub-detectors and their finite resolution, as defined in the smearing data card\footnote{\texttt{[code] }See the \texttt{RESOLution} class.}.
135If no such file is provided, predefined values are used. The coverage of the various subsystems used in the default configuration are summarised in Tab.~\ref{tab:defEta}.
136
137\begin{table*}[t]
138\begin{center}
139\caption{Default extension in pseudorapidity $\eta$ of the different subdetectors.
140The corresponding parameter name, in the smearing card, is given. \vspace{0.5cm}}
141\begin{tabular}{lll}
142\hline
143\textsc{tracker} & {\verb CEN_max_tracker } & $0.0 \leq |\eta| \leq 2.5$\\
144\textsc{ecal}, \textsc{hcal} & {\verb CEN_max_calo_cen } & $0.0 \leq |\eta| \leq 3.0$\\
145\textsc{fcal} & {\verb CEN_max_calo_fwd } & $3.0 \leq |\eta| \leq5.0$\\
146\textsc{muon} & {\verb CEN_max_mu } & $0.0 \leq |\eta| \leq 2.4$\\\hline
147\end{tabular}
148\label{tab:defEta}
149\end{center}
150\end{table*}
151
152\begin{figure}[!h]
153\begin{center}
154\includegraphics[width=\columnwidth]{Detector_Delphes_3}
155\caption{
156Profile of layout of the generic detector geometry assumed in \textsc{Delphes}. The innermost layer, close to the interaction point, is a central tracking system (pink).
157It is surrounded by a central calorimeter volume (green) with both electromagnetic and hadronic sections.
158The outer layer of the central system (red) consist of a muon system. In addition, two end-cap calorimeters (blue) extend the pseudorapidity coverage of the central detector.
159The detector parameters are defined in the user-configuration card. The extension of the various subdetectors, as defined in Tab.~\ref{tab:defEta}, are clearly visible. The detector is assumed to be strictly symmetric around the beam axis (black line). Additional forward detectors are not depicted.
160}
161\label{fig:GenDet3}
162\end{center}
163\end{figure}
164
165
166\subsubsection*{Magnetic field}
167In addition to the subdetectors, the effects of a dipolar magnetic field is simulated for the charged particles\footnote{\texttt{[code] }See the \texttt{TrackPropagation} class.}. This affects the position at which charged particles enter the calorimeters.
168
169
170
171\subsection{Tracks reconstruction}
172Every stable charged particle with a transverse momentum above some threshold and lying inside the detector volume covered by the tracker provides a track.
173By default, a track is assumed to be reconstructed with $90\%$ probability\footnote{\texttt{[code]} The reconstruction efficiency is defined in the smearing datacard by the \texttt{TRACKING\_EFF} term.} if its transverse momentum $p_T$ is higher than $0.9~\textrm{GeV}$ and if its pseudorapidity $|\eta| \leq 2.5$.
174
175
176\subsection{Simulation of calorimeters}
177
178The energy of each particle considered as stable in the generator particle list is smeared, with a Gaussian distribution depending on the calorimeter resolution. This resolution varies with the sub-calorimeter (\textsc{ecal}, \textsc{hcal}, \textsc{fcal}) measuring the particle.
179The response of each sub-calorimeter is parametrised as a function of the energy:
180\begin{equation}
181\frac{\sigma}{E} = \frac{S}{\sqrt{E}} \oplus \frac{N}{E} \oplus C,
182\label{eq:caloresolution}
183\end{equation}
184where $S$, $N$ and $C$ are the \textit{stochastic}, \textit{noise} and \textit{constant} terms, respectively.\\
185
186
187The particle four-momentum $p^\mu$ are smeared with a parametrisation directly derived from typical detector technical designs\footnote{\texttt{[code] }~\cite{bib:cmsjetresolution,bib:ATLASresolution}. The response of the detector is applied to the electromagnetic and the hadronic particles through the \texttt{SmearElectron} and \texttt{SmearHadron} functions.}.
188In the default parametrisation, the calorimeter is assumed to cover the pseudorapidity range $|\eta|<3$ and consists in an electromagnetic and an hadronic part. Coverage between pseudorapidities of $3.0$ and $5.0$ is provided by forward calorimeters, with different response to electromagnetic objects ($e^\pm, \gamma$) or hadrons.
189Muons and neutrinos are assumed not to interact with the calorimeters\footnote{In the current \textsc{Delphes} version, particles other than electrons ($e^\pm$), photons ($\gamma$), muons ($\mu^\pm$) and neutrinos ($\nu_e$, $\nu_\mu$ and $\nu_\tau$) are simulated as hadrons for their interactions with the calorimeters. The simulation of stable particles beyond the Standard Model should subsequently be handled with care.}.
190The default values of the stochastic, noisy and constant terms are given in Tab.~\ref{tab:defResol}.\\
191
192\begin{table}[!h]
193\begin{center}
194\caption{Default values for the resolution of the central and forward calorimeters. Resolution is parametrised by the \textit{stochastic} ($S$), \textit{noise} ($N$) and \textit{constant} ($C$) terms (Eq.~\ref{eq:caloresolution}).
195The corresponding parameter name, in the smearing card, is given. \vspace{0.5cm}}
196\begin{tabular}[!h]{lclc}
197\hline
198\multicolumn{2}{c}{Resolution Term} & Card flag & Value\\\hline
199 \multicolumn{4}{l}{\textsc{ecal}} \\
200 & $S$ & {\verb ELG_Scen } & $0.05$ \\
201 & $N$ & {\verb ELG_Ncen } & $0.25$ \\
202 & $C$ & {\verb ELG_Ccen } & $0.0055$ \\
203 \multicolumn{4}{l}{\textsc{fcal}, electromagnetic part} \\
204 & $S$ & {\verb ELG_Sfwd } & $2.084$ \\
205 & $N$ & {\verb ELG_Nfwd } & $0$ \\
206 & $C$ & {\verb ELG_Cfwd } & $0.107$ \\
207 \multicolumn{4}{l}{\textsc{hcal}} \\
208 & $S$ & {\verb HAD_Shcal } & $1.5$ \\
209 & $N$ & {\verb HAD_Nhcal } & $0$\\
210 & $C$ & {\verb HAD_Chcal } & $0.05$\\
211 \multicolumn{4}{l}{\textsc{fcal}, hadronic part} \\
212 & $S$ & {\verb HAD_Shf } & $2.7$\\
213 & $N$ & {\verb HAD_Nhf } & $0$. \\
214 & $C$ & {\verb HAD_Chf } & $0.13$\\
215\hline
216\end{tabular}
217\label{tab:defResol}
218\end{center}
219\end{table}
220
221The energy of electrons and photons found in the particle list are smeared using the \textsc{ecal} resolution terms. Charged and neutral final state hadrons interact with the \textsc{ecal}, \textsc{hcal} and \textsc{fcal}.
222Some long-living particles, such as the $K^0_s$, possessing lifetime $c\tau$ smaller than $10~\textrm{mm}$ are considered as stable particles although they decay before the calorimeters. The energy smearing of such particles is performed using the expected fraction of the energy, determined according to their decay products, that would be deposited into the \textsc{ecal} ($E_{\textsc{ecal}}$) and into the \textsc{hcal} ($E_{\textsc{hcal}}$). Defining $F$ as the fraction of the energy leading to a \textsc{hcal} deposit, the two energy values are given by
223\begin{equation}
224\left\{
225\begin{array}{l}
226E_{\textsc{hcal}} = E \times F \\
227E_{\textsc{ecal}} = E \times (1-F) \\
228\end{array}
229\right.
230\end{equation}
231where $0 \leq F \leq 1$. The electromagnetic part is handled as the same way as the electrons. The resulting energy measurement given after the application of the smearing is then $E = E_{\textsc{hcal}} + E_{\textsc{ecal}}$. For $K_S^0$ and $\Lambda$ hadrons, the energy fraction is $F$ is assumed to be worth $0.7$.\\
232
233\subsection{Calorimetric towers}
234
235The smallest unit for geometrical sampling of the calorimeters is a \textit{tower}; it segments the $(\eta,\phi)$ plane for the energy measurement. No longitudinal segmentation is available in the simulated calorimeters. All undecayed particles, except muons and neutrinos produce a calorimetric tower, either in \textsc{ecal}, in \textsc{hcal} or \textsc{fcal}.
236As the detector is assumed to be cylindical (e.g. symmetric in $\phi$ and with respect to the $\eta=0$ plane), the smearing card stores the number of calorimetric towers with $\phi=0$ and $\eta>0$ (default: $40$ towers). For a given $\eta$, the size of the $\phi$ segmentation is also specified. Fig.~\ref{fig:calosegmentation} illustrates the default segmentation of the $(\eta,\phi)$ plane.
237
238\begin{figure}[!h]
239\begin{center}
240\includegraphics[width=\columnwidth]{calosegmentation}
241\caption{Default segmentation of the calorimeters in the $(\eta,\phi)$ plane. Only the central detectors (\textsc{ecal}, \textsc{hcal} and \textsc{fcal}) are considered.}
242\label{fig:calosegmentation}
243\end{center}
244\end{figure}
245
246The calorimetric towers directly enter in the calculation of the missing transverse energy (\textsc{met}), and as input for the jet reconstruction algorithms. No sharing between neighbouring towers is implemented when particles enter a tower very close to its geometrical edge.
247
248\subsection{Very forward detectors simulation}
249
250Most of the recent experiments in beam colliders have additional instrumentation along the beamline. These extend the $\eta$ coverage to higher values, for the detection of very forward final-state particles.
251Zero Degree Calorimeters (\textsc{zdc}) are located at zero angle, i.e. are aligned with the beamline axis at the interaction point, and placed at the distance where the paths of incoming and outgoing beams separate (Fig.~\ref{fig:fdets}). These allow the measurement of stable neutral particles ($\gamma$ and $n$) coming from the interaction point, with large pseudorapidities (e.g. $|\eta_{\textrm{n,}\gamma}| > 8.3$ in \textsc{cms}).
252Forward taggers (called here \textsc{rp220} and \textsc{fp420} as at the \textsc{lhc}) are meant for the measurement of particles following very closely the beam path. To be able to reach these detectors, such particles must have a charge identical to the beam particles, and a momentum very close to the nominal value for the beam. These taggers are near-beam detectors located a few millimetres from the true beam trajectory and this distance defines their acceptance (Tab.~\ref{tab:fdetacceptance}).
253
254\begin{figure}[!h]
255\begin{center}
256\includegraphics[width=\columnwidth]{fdets}
257\caption{Default location of the very forward detectors, including \textsc{zdc}, \textsc{rp220} and \textsc{fp420} in the \textsc{lhc} beamline.
258Incoming (red) and outgoing (black) beams on one side of the interaction point ($s=0~\textrm{m}$).
259The Zero Degree Calorimeter is located in perfect alignment with the beamline axis at the interaction point, at $140~\textrm{m}$, where the beam paths are separated. The forward taggers are near-beam detectors located at $220~\textrm{m}$ and $420~\textrm{m}$. Beamline simulation with \textsc{Hector}~\cite{bib:Hector}.}
260\label{fig:fdets}
261\end{center}
262\end{figure}
263
264\begin{table*}[t]
265\begin{center}
266\caption{Default parameters for the forward detectors: distance from the interaction point and detector acceptance. The \textsc{lhc} beamline is assumed around the fifth interaction point. For the \textsc{zdc}, the acceptance depends only on the pseudorapidity $\eta$ of the particle, which should be neutral and stable.
267The tagger acceptance is fully determined by the distance in the transverse plane of the detector to the real beam position~\cite{bib:Hector}. It is expressed in terms of the particle energy.
268\vspace{0.5cm}}
269\begin{tabular}{llcl}
270\hline
271Detector & Distance & Acceptance & \\ \hline
272\textsc{zdc} & $140$ m & $|\eta|> 8.3$ & for $n$ and $\gamma$\\
273\textsc{rp220} & $220$ m & $E \in [6100 ; 6880]$ (GeV) & at $2~\textrm{mm}$\\
274\textsc{fp420} & $420$ m & $E \in [6880 ; 6980]$ (GeV) & at $4~\textrm{mm}$\\
275\hline
276\end{tabular}
277\label{tab:fdetacceptance}
278\end{center}
279\end{table*}
280
281
282While neutral particles propagate along a straight line to the \textsc{zdc}, a dedicated simulation of the transport of charged particles is needed for \textsc{rp220} and \textsc{fp420}. This fast simulation uses the \textsc{Hector} software~\cite{bib:Hector}, which includes the chromaticity effects and the geometrical aperture of the beamline elements.
283
284Some subdetectors have the ability to measure the time of flight of the particle.
285This corresponds to the delay after which the particle is observed in the detector, after the bunch crossing. The time of flight measurement of \textsc{zdc} and \textsc{fp420} detector is implemented here. For the \textsc{zdc}, the formula is simply
286\begin{equation}
287 t = t_0 + \frac{1}{v} \times \Big( \frac{s-z}{\cos \theta}\Big),
288\end{equation}
289where $t$ is the time of flight, $t_0$ is the true time coordinate of the vertex from which the particle originates, $v$ the particle velocity, $s$ is the \textsc{zdc} distance to the interaction point, $z$ is the longitudinal coordinate of the vertex from which the particle comes from, $\theta$ is the particle emission angle. This assumes that the neutral particle observed in the \textsc{zdc} is highly relativistic, i.e. travelling at the speed of light $c$. We also assume that $\cos \theta = 1$, i.e. $\theta \approx 0$ or equivalently $\eta$ is large. As an example, $\eta = 5$ leads to $\theta = 0.013$ and $1 - \cos \theta < 10^{-4}$.
290The formula then reduces to
291\begin{equation}
292 t = \frac{1}{c} \times (s-z)
293\end{equation}
294For example, a photon takes $0.47~\mu\textrm{s}$ to reach a \textsc{zdc} located at $s=140~\textrm{m}$, neglecting $z$ and $\theta$, and assuming that $v=c$. Only neutrons and photons are currently assumed to be able to reach the \textsc{zdc}. All other particles are neglected in the \textsc{zdc}.
295
296\section{High-level object reconstruction}
297
298Analysis object data contain the final collections of particles ($e^\pm$, $\mu^\pm$, $\gamma$) or objects (light jets, $b$-jets, $\tau$-jets, $E_T^\textrm{miss}$) and are stored\footnote{\texttt{[code] }All these processed data are located under the \texttt{Analysis} tree.} in the output file created by \textsc{Delphes}.
299In addition, some detector data are added: tracks, calorimetric towers and hits in \textsc{zdc}, \textsc{rp220} and \textsc{fp420}.
300While electrons, muons and photons are easily identified, some other objects are more difficult to measure, like jets or missing energy due to invisible particles.
301
302For most of these objects, their four-momentum and related quantities are directly accessible in \textsc{Delphes} output ($E$, $\vec{p}$, $p_T$, $\eta$ and $\phi$). Additional properties are available for specific objects (like the charge and the isolation status for $e^\pm$ and $\mu^\pm$, the result of application of $b$-tag for jets and time-of-flight for some detector hits).
303
304
305
306\subsection{Photon and charged lepton reconstruction}
307From here onwards, \textit{electrons} refer to both positrons ($e^+$) and electrons ($e^-$), and $\textit{charged leptons}$ refer to electrons and muons ($\mu^\pm$), leaving out the $\tau^\pm$ leptons as they decay before being detected.
308\subsubsection*{Electrons and photons}
309Photon and electron ($e^\pm$) candidates are reconstructed if they fall into the acceptance of the tracking system and have a transverse momentum above a threshold (default $p_T > 10~\textrm{GeV}$). A calorimetric tower will be seen in the detector, an electrons leave in addition a track. Consequently, electrons and photons creates as usual a candidate in the jet collection.
310
311\subsubsection*{Muons}
312
313Generator level muons entering the detector acceptance are considered as candidates for the analysis level.
314The acceptance is defined in terms of a transverse momentum threshold to be overpassed that should be computed using the chosen geometry of the detector and the magnetic field considered. (default : $p_T > 10~\textrm{GeV}$) and of the pseudorapidity coverage of the muon system of the detector (default: $-2.4 \leq \eta \leq 2.4$).
315The application of the detector resolution on the muon momentum depends on a Gaussian smearing of the $p_T$ variable\footnote{\texttt{[code]} See the \texttt{SmearMuon} method.}. Neither $\eta$ nor $\phi$ variables are modified beyond the calorimeters: no additional magnetic field is applied. In addition, multiple scattering is also neglected. This implies that low energy muons have in \textsc{Delphes} a better resolution than in a real detector. Moreover, muons leave no deposit in calorimeters.
316
317\subsubsection*{Charged lepton isolation}
318
319To improve the quality of the contents of the charged lepton collections, additional criteria can be applied to impose some isolation. This requires that electron or muon candidates are isolated in the detector from any other particle, within a small cone. In \textsc{Delphes}, charged lepton isolation demands that there is no other charged particle with $p_T>2~\textrm{GeV}$ within a cone of $\Delta R = \sqrt{\Delta \eta^2 + \Delta \phi^2} <0.5$ around the lepton. The result (i.e. \textit{isolated} or \textit{not}) is added to the charged lepton measured properties\footnote{\texttt{[code] }See the \texttt{IsolFlag} output of the \texttt{Electron} or \texttt{Muon} collections in the \texttt{Analysis} tree.}. No calorimetric isolation is applied. \\
320
321
322
323
324
325
326
327\subsection{Jet reconstruction}
328
329A realistic analysis requires a correct treatment of particles which have hadronised. Therefore, the most widely currently used jet algorithms have been integrated into the \textsc{Delphes} framework using the \textsc{FastJet} tools~\cite{bib:FastJet}.
330Six different jet reconstruction schemes are available\footnote{\texttt{[code] }The choice is done by allocating the \texttt{JET\_jetalgo } input parameter in the smearing card.}. The first three belong to the cone algorithm class while the last three are using a sequential recombination scheme. For all of them, the towers are used as input for the jet clustering. Jet algorithms also differ in their sensitivity to soft particles or collinear splittings, and with their computing speed performance.
331
332\subsubsection*{Cone algorithms}
333
334\begin{enumerate}
335
336\item {\it CDF Jet Clusters}~\cite{bib:jetclu}: Algorithm forming jets by associating together towers lying within a circle (default radius $\Delta R=0.7$) in the $(\eta$, $\phi)$ space.
337This so-called \textsc{Jetclu} cone jet algorithm is used by the \textsc{cdf} experiment in Run II.
338All towers with a transverse energy $E_T$ higher than a given threshold (default: $E_T > 1~\textrm{GeV}$) are used to seed the jet candidates.
339The existing \textsc{FastJet} code has been modified to allow easy modification of the tower pattern in $\eta$, $\phi$ space.
340In the following versions of \textsc{Delphes}, a new dedicated plug-in will be created on this purpose\footnote{\texttt{[code] }\texttt{JET\_coneradius} and \texttt{JET\_seed} variables in the smearing card.}.
341
342\item {\it CDF MidPoint}~\cite{bib:midpoint}: Algorithm developed for the \textsc{cdf} Run II to reduce infrared and collinear sensitivity compared to purely seed-based cone by adding `midpoints' (energy barycentres) in the list of cone seeds.
343
344\item {\it Seedless Infrared Safe Cone}~\cite{bib:SIScone}: The \textsc{SISCone} algorithm is simultaneously insensitive to additional soft particles and collinear splittings, and fast enough to be used in experimental analysis.
345
346\end{enumerate}
347
348\subsubsection*{Recombination algorithms}
349
350The three following jet algorithms are safe for soft radiations (\textit{infrared}) and collinear splittings. They rely on recombination schemes where calorimeter tower 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$.
351
352The 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 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:
353
354\begin{enumerate}[start=4]
355
356\item {\it Longitudinally invariant $k_t$ jet}~\cite{bib:ktjet}:
357\begin{equation}
358\begin{array}{l}
359 d_{ij} = \min(k_{ti}^2,k_{tj}^2)\Delta R_{ij}^2/R^2 \\
360 d_{iB}=k_{ti}^2 \\
361\end{array}
362\end{equation}
363
364\item {\it Cambridge/Aachen jet}~\cite{bib:aachen}:
365\begin{equation}
366\begin{array}{l}
367d_{ij} = \Delta R_{ij}^2/R^2\\
368d_{iB}=1 \\
369\end{array}
370\end{equation}
371
372\item {\it Anti $k_t$ jet}~\cite{bib:antikt}: where hard jets are exactly circular in the $(y,\phi)$ plane
373\begin{equation}
374\begin{array}{l}
375d_{ij} = \min(1/k_{ti}^2,1/k_{tj}^2)\Delta R_{ij}^2/R^2 \\
376d_{iB}=1/k_{ti}^2 \\
377\end{array}
378\end{equation}
379\end{enumerate}
380
381By default, reconstruction uses a cone algorithm with $\Delta R=0.7$. Jets are stored if their transverse energy is higher\footnote{\texttt{[code] PTCUT\_jet }variable in the smearing card.} than $20~\textrm{GeV}$.
382
383
384\subsection{$b$-tagging}
385
386A jet is tagged as $b$-jets if its direction lies in the acceptance of the tracker and if it is associated to a parent $b$-quark. A $b$-tagging efficiency of $40\%$ is assumed if the jet has a parent $b$ quark. For $c$-jets and light jets (i.e. originating in $u$,$d$,$s$ quarks or in gluons), a fake $b$-tagging efficiency of $10 \%$ and $1 \%$ respectively is assumed\footnote{\texttt{[code] }Corresponding to the \texttt{TAGGING\_B}, \texttt{MISTAGGING\_C} and \texttt{MISTAGGING\_L} constants, for (respectively) the efficiency of tagging of a $b$-jet, the efficiency of mistagging a $c$-jet as a $b$-jet, and the efficiency of mistagging a light jet ($u$,$d$,$s$,$g$) as a $b$-jet.}
387%(Fig.~\ref{fig:btag})
388.
389The (mis)tagging relies on the true particle identity (\textsc{pid}) of the most energetic particle within a cone around the observed $(\eta,\phi)$ region, with a radius equal to the one used to reconstruct the jet (default: $\Delta R$ of $0.7$).
390
391%\begin{figure}[!h]
392%\begin{center}
393%\includegraphics[width=0.6\columnwidth]{btag}
394%\caption{Default efficiency of $b$-tag for jets coming from $b$ quarks, $c$ quarks and from other particles (jets from gluons or $u$, $d$ and $s$ quarks).}
395%\label{fig:btag}
396%\end{center}
397%\end{figure}
398
399
400\subsection{\texorpdfstring{$\tau$}{\texttau} identification}
401
402Jets originating from $\tau$-decays are identified using an identification procedure consistent with the one applied in a full detector simulation~\cite{bib:cmsjetresolution}.
403The tagging relies on two properties of the $\tau$ lepton. First, $77\%$ of the $\tau$ hadronic decays contain only one charged hadron associated to a few neutrals (Tab.~\ref{tab:taudecay}). Tracks are useful for this criterion. Secondly, the particles arisen from the $\tau$ lepton produce narrow jets in the calorimeter (this is defined as the jet \textit{collimation}).
404
405\begin{table}[!h]
406\begin{center}
407\caption{ Branching rations for $\tau^-$ lepton~\cite{bib:pdg}. $h^\pm$ and $h^0$ refer to charged and neutral hadrons, respectively. $n \geq 0$ and $m \geq 0$ are integers.
408\vspace{0.5cm} }
409\begin{tabular}[!h]{ll}
410\hline
411 \multicolumn{2}{l}{\textbf{Leptonic decays}}\\
412 $ \tau^- \rightarrow e^- \ \bar \nu_e \ \nu_\tau$ & $17.9\% $ \\
413 $ \tau^- \rightarrow \mu^- \ \bar \nu_\mu \ \nu_\tau$ & $17.4\%$ \\
414 \multicolumn{2}{l}{\textbf{Hadronic decays}}\\
415 $ \tau^- \rightarrow h^-\ n\times h^\pm \ m\times h^0\ \nu_\tau$ & $64.7\%$ \\
416 $ \tau^- \rightarrow h^-\ m\times h^0 \ \nu_\tau$ & $50.1\%$ \\
417 $ \tau^- \rightarrow h^-\ h^+ h^- m\times h^0 \ \nu_\tau$ & $14.6\%$ \\
418\hline
419\end{tabular}
420\label{tab:taudecay}
421\end{center}
422\end{table}
423
424
425%\begin{wrapfigure}{l}{0.3\columnwidth}
426\begin{figure}[!h]
427\begin{center}
428\includegraphics[width=0.6\columnwidth]{Tau}
429\caption{Illustration of the identification of $\tau$-jets. The jet cone is narrow and contains only one track. The small cone shown as the red one is used for the \textit{electromagnetic collimation}, while the green cone is the cone radius used to reconstruct the jet originating from the $\tau$-decay.}
430\label{h_WW_ss_cut1}
431\end{center}
432\end{figure}
433%\end{wrapfigure}
434
435
436\subsubsection*{Electromagnetic collimation}
437
438To use the narrowness of the $\tau$-jet, the \textit{electromagnetic collimation} $C_{\tau}^{em}$ is defined as the sum of the energy of towers in a small cone of radius $R^\textrm{em}$ around the jet axis, divided by the energy of the reconstructed jet.
439To be taken into account, a calorimeter tower should have a transverse energy $E_T^\textrm{tower}$ above a given threshold.
440A large fraction of the jet energy is expected in this small cone. This fraction, or collimation factor, is represented in Fig.~\ref{fig:tau2} for the default values (see Tab.~\ref{tab:tauRef}).
441
442\begin{figure}[!h]
443\begin{center}
444\includegraphics[width=\columnwidth]{Tau2}
445\caption{Distribution of the electromagnetic collimation $C_\tau$ variable for true $\tau$-jets, normalised to unity. This distribution is shown for associated $WH$ photoproduction~\cite{bib:whphotoproduction}, where the Higgs boson decays into a $W^+ W^-$ pair. Each $W$ boson decays into a $\ell \nu_\ell$ pair, where $\ell = e, \mu, \tau$.
446Events generated with \textsc{MadGraph/MadEvent}~\cite{bib:mgme}.
447Final state hadronisation is performed by \textsc{Pythia}~\cite{bib:pythia}.
448Histogram entries correspond to true $\tau$-jets, matched with generator level data. }
449\label{fig:tau2}
450\end{center}
451\end{figure}
452
453\subsubsection*{Tracking isolation}
454
455The tracking isolation for the $\tau$ identification requires that the number of tracks associated to a particle with a significant transverse momentum is one and only one in a cone of radius $R^\textrm{tracks}$ (3-prong $\tau$s are dropped).
456This cone should be entirely incorporated into the tracker to be taken into account. Default values of these parameters are given in Tab.~\ref{tab:tauRef}.
457
458
459
460\begin{figure}[!h]
461\begin{center}
462\includegraphics[width=\columnwidth]{Tau1}
463\caption{Distribution of the number of tracks $N^\textrm{tracks}$ within a small jet cone for true $\tau$-jets, normalised to unity. Photoproduced $WH$ events, where $W$ bosons decay leptonically ($e,\mu,\tau$), as in Fig.~\ref{fig:tau2}.
464Histogram entries correspond to true $\tau$-jets, matched with generator level data.}
465\label{fig:tau1}
466\end{center}
467\end{figure}
468
469
470\begin{table}[!h]
471\begin{center}
472\caption{Default values for parameters used in $\tau$-jet reconstruction algorithm. Electromagnetic collimation requirements involve the inner \textit{small} cone radius $R^\textrm{em}$, the minimum transverse energy for calotowers $E_T^\textrm{tower}$ and the collimation factor $C_\tau$. Tracking isolation constrains the number of tracks with a significant transverse momentum $p_T^\textrm{tracks}$ in a cone of radius $R^\textrm{tracks}$. Finally, the $\tau$-jet collection is purified by the application of a cut on the $p_T$ of $\tau$-jet candidates.
473\vspace{0.5cm} }
474\begin{tabular}[!h]{lll}
475\hline
476Parameter & Card flag & Value\\\hline
477\multicolumn{3}{l}{\textbf{Electromagnetic collimation}} \\
478$R^\textrm{em}$ & \texttt{TAU\_energy\_scone } & $0.15$\\
479min $E_{T}^\textrm{tower}$ & {\verb JET_M_seed } & $1.0$~GeV\\
480$C_{\tau}$ & \texttt{TAU\_energy\_frac} & $0.95$\\
481\multicolumn{3}{l}{\textbf{Tracking isolation}} \\
482$R^\textrm{tracks}$ & \texttt{TAU\_track\_scone} & $0.4$\\
483min $p_T^{tracks}$ & \texttt{PTAU\_track\_pt } & $2$ GeV\\
484\multicolumn{3}{l}{\textbf{$\tau$-jet candidate}} \\
485$\min p_T$ & \texttt{TAUJET\_pt} & $10$ GeV\\
486\hline
487\end{tabular}
488\label{tab:tauRef}
489\end{center}
490\end{table}
491
492\subsubsection*{Purity}
493Once 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\%$.
494
495\subsection{Missing transverse energy}
496In 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}$.
497\begin{equation}
498\overrightarrow{p_T} = \left(
499\begin{array}{c}
500p_x\\
501p_y\\
502\end{array}
503\right)
504~ \textrm{and} ~
505\left\{
506\begin{array}{l}
507 p_x^\textrm{miss} = - p_x^\textrm{obs} \\
508 p_y^\textrm{miss} = - p_y^\textrm{obs} \\
509\end{array}
510\right.
511\end{equation}
512The \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.
513In 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:
514\begin{equation}
515\overrightarrow{E_T}^\textrm{miss} = - \sum^\textrm{towers}_i \overrightarrow{E_T}(i)
516\end{equation}
517
518
519\section{Trigger emulation}
520
521New 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 Higgs boson cross-section within the \textsc{sm} can be as small as $2 \times 10^{-3}~\textrm{nb}$ ($pp \rightarrow WH$, $m_H=115~\textrm{GeV}$).
522
523%High statistics are required for data analyses, consequently imposing high luminosity, i.e. a high collision rate.
524As 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.
525This 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.}.
526Dedicated 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.
527
528Most 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}$.
529
530A 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.
531In a real experiment, the online selection is often divided into several steps (or \textit{levels}).
532This splits the overall reduction factor into a product of smaller factors, corresponding to the different trigger levels.
533This is related to the architecture of the experiment data acquisition chain, with limited electronic buffers requiring a quick decision for the first trigger level.
534First 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.
535Later levels are more complex, of finer-but-not-final quality and based on full detector data.
536
537Real triggers are thus intrinsically based on reconstructed data with a worse resolution than final analysis data.
538On the contrary, same data are used in \textsc{Delphes} for trigger emulation and for final analyses.
539
540\section{Validation}
541
542\textsc{Delphes} performs a fast simulation of a collider experiment.
543Its quality and validity are assessed by comparing to resolution of the reconstructed data to the \textsc{cms} detector expectations.
544
545Electrons and muons are by construction equal to the experiment designs, as the Gaussian smearing of their kinematics properties is defined according to their resolutions.
546Similarly, the $b$-tagging efficiency (for real $b$-jets) and misidentification rates (for fake $b$-jets) are taken from the expected values of the experiment.
547Unlike these simple objects, jets and missing transverse energy should be carefully cross-checked.
548
549\subsection{Jet resolution}
550
551The 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.
552
553This validation is based on $pp \rightarrow gg$ events produced with \textsc{MadGraph/MadEvent} 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 axes
554\begin{equation}
555\Delta R = \sqrt{ \big(\eta^\textrm{rec} - \eta^\textrm{MC} \big)^2 + \big(\phi^\textrm{rec} - \phi^\textrm{MC} \big)^2}<0.25.
556\end{equation}
557The jets made of generator-level particles, here refered as \textit{MC jets}, are obtained by applying the same clustering algorithm to all particles considered as stable after hadronisation.
558Jets produced by \textsc{Delphes} and satisfying the matching criterion are called hereafter \textit{reconstructed jets}.
559
560The 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.
561The $E_T^\textrm{rec}/E_T^\textrm{MC}$ histogram is fitted with a Gaussian distribution in the interval \mbox{$\pm 2$~\textsc{rms}} centred around the mean value.
562The resolution in each $\hat{p}_T$ bin is obtained by the fit mean $\langle x \rangle$ and variance $\sigma^2(x)$:
563\begin{equation}
564%\frac{\sigma(R_{jet})}{\langle R_{jet} \rangle }=
565\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}}~
566\Big( \hat{p}_T(i) \Big)\textrm{, for all }i.
567\end{equation}
568
569\begin{figure}[!h]
570\begin{center}
571\includegraphics[width=\columnwidth]{resolutionJet}
572\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~\cite{bib:cmsjetresolution}, in blue.}
573\label{fig:jetresol}
574\end{center}
575\end{figure}
576
577The resulting jet resolution as a function of $E_T^\textrm{MC}$ is shown in Fig.~\ref{fig:jetresol}.
578This distribution is fitted with a function of the following form:
579\begin{equation}
580\frac{a}{E_T^\textrm{MC}}\oplus \frac{b}{\sqrt{E_T^\textrm{MC}}}\oplus c,
581\end{equation}
582where $a$, $b$ and $c$ are the fit parameters.
583It is then compared to the resolution published by the \textsc{cms} collaboration~\cite{bib:cmsjetresolution}. The resolution curves from \textsc{Delphes} and \textsc{cms} are in good agreement.
584
585\subsection{MET resolution}
586
587All 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.
588The resolution of the $\overrightarrow{E_T}^\textrm{miss}$ variable, as obtained with \textsc{Delphes}, is then crucial.
589
590The samples used to study the \textsc{met} performance are identical to those used for the jet validation.
591It is worth noting that the contribution to $E_T^\textrm{miss}$ from muons is negligible in the studied sample.
592The 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.
593The quality of the \textsc{met} reconstruction is checked via the resolution on its horizontal component $E_x^\textrm{miss}$.
594
595The $E_x^\textrm{miss}$ resolution is evaluated in the following way.
596The distribution of the difference between $E_x^\textrm{miss}$ in \textsc{Delphes} and at generator-level is fitted with a Gaussian function in each $(\Sigma E_T)$ bin. The fit \textsc{rms} gives the \textsc{met} resolution in each bin.
597The resulting value is plotted in Fig.~\ref{fig:resolETmis} as a function of the total visible transverse
598energy.
599
600\begin{figure}[!h]
601\begin{center}
602\includegraphics[width=\columnwidth]{resolutionETmis}
603\caption{$\sigma(E^\textrm{miss}_{x})$ as a function on the scalar sum of all towers ($\Sigma E_T$) for $pp \rightarrow gg$ events.}
604\label{fig:resolETmis}
605\end{center}
606\end{figure}
607
608The resolution $\sigma_x$ of the horizontal component of \textsc{met} is observed to behave like
609\begin{equation}
610\sigma_x = \alpha ~\sqrt(E_T) ~~~(\mathrm{GeV}^{1/2}),
611\end{equation}
612where the $\alpha$ parameter is depending on the resolution of the calorimeters.
613
614The \textsc{met} resolution expected for the \textsc{cms} detector for similar events is $\sigma_x = (0.6-0.7) ~ \sqrt(E_T) ~ \mathrm{GeV}^{1/2}$ with no pile-up\footnote{\textit{Pile-up} events are extra simultaneous $pp$ collision occurring at the same bunch crossing.}~\cite{bib:cmsjetresolution}, which compares very well with the $\alpha = 0.68$ obtained with \textsc{Delphes}.
615
616\subsection{\texorpdfstring{$\tau$}{\texttau}-jet efficiency}
617Due to the complexity of their reconstruction algorithm, $\tau$-jets have also to be checked.
618Tab.~\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.
619
620\begin{table}[!h]
621\begin{center}
622\caption{Reconstruction efficiencies of $\tau$-jets in decays from $Z$ or $H$ bosons, in \textsc{Delphes} and in the \textsc{cms} experiment~\cite{bib:cmstauresolution}.\vspace{0.5cm}}
623\begin{tabular}{lll}
624\hline
625\multicolumn{2}{c}{\textsc{cms}} & \\
626$Z \rightarrow \tau^+ \tau^-$ & $38 \%$ & \\
627$H \rightarrow \tau^+ \tau^-$ & $36 \%$ & $m_H = 150~\textrm{GeV}$ \\
628$H \rightarrow \tau^+ \tau^-$ & $47 \%$ & $m_H = 300~\textrm{GeV}$ \\
629\multicolumn{2}{c}{\textsc{Delphes}} & \\
630$H \rightarrow \tau^+ \tau^-$ &$42 \%$ & $m_H = 140~\textrm{GeV}$ \\
631\hline
632\end{tabular}
633\label{tab:taurecoefficiency}
634\end{center}
635\end{table}
636
637
638\section{Visualisation}
639
640When performing an event analysis, a visualisation tool is useful to convey information about the detector layout and the event topology in a simple way. The \textit{Fast and Realistic OpenGL Displayer} \textsc{frog}~\cite{bib:Frog} has been interfaced in \textsc{Delphes}, allowing an easy display of the defined detector configuration\footnote{\texttt{[code] } To prepare the visualisation, the \texttt{FLAG\_frog} parameter should be equal to $1$.}.
641
642% \begin{figure}[!h]
643% \begin{center}
644% \includegraphics[width=\columnwidth]{Detector_Delphes_1}
645% \caption{Layout of the generic detector geometry assumed in \textsc{Delphes}. The innermost layer, close to the interaction point, is a central tracking system (pink), embedded into a solenoidal magnetic field.
646% It is surrounded by a central calorimeter volume (green) with both electromagnetic and hadronic sections.
647% The outer layer of the central system (red) consist of a muon system.
648% In addition, two end-cap calorimeters (blue) extend the pseudorapidity coverage of the central detector.
649% The actual detector granularity and extension is defined in the smearing card.
650% The detector is assumed to be strictly symmetric around the beam axis (black line).
651% Additional forward detectors are not depicted.}
652% \label{fig:GenDet}
653% \end{center}
654% \end{figure}
655
656Two and three-dimensional representations of the detector configuration can be used for communication purpose, as it clearly shows the geometric coverage of the different detector subsystems. As an illustration, the generic detector geometry assumed in this paper is shown in Fig.~\ref{fig:GenDet3}
657%, \ref{fig:GenDet}
658 and~\ref{fig:GenDet2}.
659As pointed before, the detector is assumed to be strictly symmetric around the beam axis.
660The extensions of the central tracking system, the central calorimeters and both forward calorimeters are visible.
661Nevertheless, it should be noticed that only the geometrical coverage is depicted and that the calorimeter segmentation is not taken into account in the drawing of the detector. Moreover, both the radius and the length of each sub-detectors are just display parameters and are insignificant for the physics simulation.
662
663\begin{figure}[!h]
664\begin{center}
665\includegraphics[width=\columnwidth]{Detector_Delphes_2b}
666\caption{Layout of the generic detector geometry assumed in \textsc{Delphes}. Open 3D-view of the detector with solid volumes. Same colour codes as for Fig.~\ref{fig:GenDet3} are applied. Additional forward detectors are not depicted.}
667\label{fig:GenDet2}
668\end{center}
669\end{figure}
670
671Deeper understanding of interesting physics processes is possible by displaying the events themselves.
672The visibility of each set of objects ($e^\pm$, $\mu^\pm$, $\tau^\pm$, jets, transverse missing energy) is enhanced by a colour coding.
673Moreover, kinematics information of each object is visible by a simple mouse action.
674As an illustration, an associated photoproduction of a $W$ boson and a $t$ quark is shown in Fig.~\ref{fig:wt}.
675This corresponds to a $pp(\gamma p \rightarrow Wt)pX$ process, where the $Wt$ couple is induced by an incoming photon emitted by one interacting proton~\cite{bib:wtphotoproduction}.
676This leading proton survives from the photon emission and subsequently from the $pp$ interaction, and is present in the final state.
677As the energy and virtuality of the emitted photon are low, the surviving proton does not leave the beam and escapes from the central detector without being detected.
678The experimental signature is a lack of hadronic activity in one forward hemisphere, where the surviving proton escapes.
679The $t$ quark decays into a $W$ boson and a $b$ quark.
680Both $W$ bosons decay into leptons ($W \rightarrow \mu \nu_\mu$ and $W \rightarrow e \nu_e$).
681The balance between the missing transverse energy and the charged lepton pair is clear, as well as the presence of an empty forward region. It is interesting to notice that the reconstruction algorithms build a fake $\tau$-jet around the electron.
682
683\begin{figure}[!h]
684\begin{center}
685%\includegraphics[width=\columnwidth]{Events_Delphes_1}
686\includegraphics[width=\columnwidth]{DisplayWt}
687\caption{Example of $pp(\gamma p \rightarrow Wt)pY$ event, with $t \rightarrow Wb$.
688One $W$ boson decays into a $\mu \nu_\mu$ pair and the second one into a $e \nu_e$ pair.
689The surviving proton leaves a forward hemisphere with no hadronic activity.
690The isolated muon is shown as the dark blue vector.
691Around the electron, in red, is reconstructed a fake $\tau$-jet (green vector surrounded by a blue cone), while the reconstructed missing energy (in grey) is very small. One jet is visible in one forward region, along the beamline axis, opposite to the direction of the escaping proton.}
692\label{fig:wt}
693\end{center}
694\end{figure}
695
696For the comparison, Fig.~\ref{fig:gg} depicts an inclusive gluon pair production $pp \rightarrow ggX$.
697The event final state contains more jets, in particular along the beam axis, which is expected as the interacting protons are destroyed by the collision. Two muon candidates and large missing transverse energy are also visible.
698
699\begin{figure}[!h]
700\begin{center}
701%\includegraphics[width=\columnwidth]{Events_Delphes_1}
702\includegraphics[width=\columnwidth]{Displayppgg}
703\caption{Example of inclusive gluon pair production $pp \rightarrow ggX$. Many jets are visible in the event, in particular along the beam axis. Two muons (in blue) are produced and the missing transverse energy is significant in this event (grey vector).}
704\label{fig:gg}
705\end{center}
706\end{figure}
707
708
709\section{Conclusion and perspectives}
710
711% \subsection{version 1}
712% We have described here the major features of the \textsc{Delphes} framework, introduced for the fast simulation of a collider experiment.
713% It has already been used for several phenomenological studies, in particular in photon interactions at the \textsc{lhc}.
714%
715% \textsc{Delphes} takes the output of event generators, in various formats, and yields analysis object data.
716% The simulation applies the resolutions of central and forward detectors by smearing the kinematical properties of final state particles.
717% It yields tracks in a solenoidal magnetic field and calorimetric towers.
718% Realistic reconstruction algorithms are run, including the \textsc{FastJet} package, to produce collections of $e^\pm$, $\mu^\pm$, jets and $\tau$-jets. $b$-tag and missing transverse energy are also evaluated.
719% The output is validated by comparing to the \textsc{cms} expected performances.
720% A trigger stage can be emulated on the output data.
721% At last, event visualisation is possible through the \textsc{Frog} 3D event display.
722%
723%
724% \textsc{Delphes} has been developped using the parameters of the \textsc{cms} experiment but can be easily extended to \textsc{atlas} and other non-\textsc{lhc} experiments, as at Tevatron or at the \textsc{ilc}. Further developments include a more flexible design for the subdetector assembly and possibly the implementation of an event mixing module for pile-up event simulation.
725%
726%
727% \subsection{version 2}
728We have described here the major features of the \textsc{Delphes} framework, introduced for the fast simulation of a collider experiment. This framework is a tool meant for feasibility studies in phenomenology, probing the observability of models in collider experiments. It has already been used for several analyses, in particular in photon interactions at the \textsc{lhc}~\cite{bib:wtphotoproduction, bib:papierquisortirajamais, bib:papiersimon}.
729
730\textsc{Delphes} takes the output of event generators and yields analysis object data.
731The simulation includes central and forward detectors to produce realistic observables using standard reconstruction algorithms.
732Moreover, the framework allows trigger emulation and 3D event visualisation.
733
734\textsc{Delphes} has been developed using the parameters of the \textsc{cms} experiment but can be easily extended to \textsc{atlas} and other non-\textsc{lhc} experiments, as at Tevatron or at the \textsc{ilc}. Further developments include a more flexible design for the subdetector assembly and possibly the implementation of an event mixing module for pile-up event simulation.
735
736
737
738\section*{Acknowledgements}
739\addcontentsline{toc}{section}{Acknowledgements}
740The authors would like to thank Jer\^ome de Favereau, Christophe Delaere, Vincent Lema\^itre, Muriel Vander Donckt and David d'Enterria for useful discussions and comments, and Loic Quertenmont for support in interfacing \textsc{Frog}. We are also really grateful to Alice Dechambre and Simon de Visscher for being beta testers of the complete package.
741Part of this work was supported by the Belgian Federal Office for Scientific, Technical and Cultural Affairs through the Interuniversity Attraction Pole P6/11.
742
743
744\begin{thebibliography}{99}
745\addcontentsline{toc}{section}{References}
746
747\bibitem{bib:Delphes} \textsc{Delphes}, \href{http://www.fynu.ucl.ac.be/delphes.html}{www.fynu.ucl.ac.be/delphes.html}
748%hepforge:
749\bibitem{bib:stdhep} L.A. Garren, M. Fischler, \\ \href{http://cepa.fnal.gov/psm/stdhep/c++}{cepa.fnal.gov/psm/stdhep/c++}
750\bibitem{bib:lhe} J. Alwall, et al., \textbf{Comput. Phys. Commun.} \href{http://dx.doi.org/10.1016/j.cpc.2006.11.010}{176:300-304,2007}.
751\bibitem{bib:Root} %\textsc{Root}, \textit{An Object Oriented Data Analysis Framework},
752R. Brun, F. Rademakers, Nucl. Inst. \& Meth. in \textbf{Phys. Res. A} \href{http://dx.doi.org/10.1016/S0168-9002(97)00048-X}{389 (1997) 81-86}.
753\bibitem{bib:ExRootAnalysis} %\textit{The} \textsc{ExRootAnalysis} \textit{analysis steering utility},
754P. Demin, (2006), unpublished. Now part of \textsc{MadGraph/MadEvent}.
755\bibitem{bib:cmsjetresolution} The CMS Collaboration, \textbf{CERN/LHCC} \\ \href{http://documents.cern.ch/cgi-bin/setlink?base=lhcc&categ=public&id=lhcc-2006-001}{2006-001}.
756\bibitem{bib:ATLASresolution} The ATLAS Collaboration, \\ arXiv:\href{http://arxiv.org/abs/arxiv:0901.0512}{0901.0512v1}[hep-ex].
757\bibitem{bib:Hector} %\textsc{Hector}, \textit{a fast simulator for the transport of particles in beamlines},
758X. Rouby, J. de Favereau, K. Piotrzkowski, \textbf{JINST} \href{http://www.iop.org/EJ/abstract/1748-0221/2/09/P09005}{2 P09005 (2007)}.
759\bibitem{bib:FastJet} %\textit{The} \textsc{FastJet} \textit{package},
760M. Cacciari, G.P. Salam, \textbf{Phys. Lett. B} \href{http://dx.doi.org/10.1016/j.physletb.2006.08.037}{641 (2006) 57}.
761\bibitem{bib:jetclu} %\textsc{cdf} Run I legacy algorithm,
762F. Abe et al. (CDF Coll.), \textbf{Phys. Rev. D} \href{http://link.aps.org/doi/10.1103/PhysRevD.45.1448}{45 (1992) 1448}.
763\bibitem{bib:midpoint} %Run II Jet Physics: Proceedings of the Run II QCD and Weak Boson Physics Workshop,
764G.C. Blazey, et al., arXiv:\href{http://arxiv.org/abs/hep-ex/0005012}{0005012}[hep-ex].
765\bibitem{bib:SIScone} %\textsc{SIScone}, \textit{A practical Seedless Infrared-Safe Cone jet algorithm},
766G.P. Salam, G. Soyez, \textbf{JHEP} \href{http://dx.doi.org/10.1088/1126-6708/2007/05/086}{05 (2007) 086}.
767\bibitem{bib:ktjet} S. Catani, Y.L. Dokshitzer, M.H. Seymour, B.R. Webber, \textbf{Nucl. Phys. B} \href{http://dx.doi.org/10.1016/0550-3213(93)90166-M}{406 (1993) 187}; S.D. Ellis, D.E. Soper, \textbf{Phys. Rev. D} \href{http://link.aps.org/doi/10.1103/PhysRevD.48.3160}{48 (1993) 3160}.
768\bibitem{bib:aachen} Y.L. Dokshitzer, G.D. Leder, S. Moretti, B.R. Webber, \textbf{JHEP} \href{http://dx.doi.org/10.1088/1126-6708/1998/01/011}{08} \href{http://dx.doi.org/10.1088/1126-6708/1998/01/011}{(1997) 001}; M. Wobisch, T. Wengler, arXiv:\href{http://arxiv.org/abs/hep-ph/9907280}{9907280}[hep-ph].
769\bibitem{bib:antikt} %\textit{The anti-kt jet clustering algorithm},
770M. Cacciari, G.P. Salam, G. Soyez, \textbf{JHEP} \href{http://dx.doi.org/10.1088/1126-6708/2008/04/063}{04 (2008) 063}.
771\bibitem{bib:pdg} C. Amsler et al. (Particle Data Group), \textbf{Phys. Lett. B} \href{http://dx.doi.org/10.1016/j.physletb.2008.07.018}{667 (2008) 1}.
772\bibitem{bib:whphotoproduction} S. Ovyn, \textbf{Nucl. Phys. Proc. Suppl.} \href{http://dx.doi.org/10.1016/j.nuclphysbps.2008.07.034}{179-180 (2008) 269-276}.
773\bibitem{bib:mgme} %\textsc{MadGraph/MadEvent v4}, \textit{The New Web Generation},
774J. Alwall, et al., \textbf{JHEP} \href{http://dx.doi.org/10.1088/1126-6708/2007/09/028}{09 (2007) 028}.
775\bibitem{bib:pythia} %\textsc{Pythia 6.4}, \textit{Physics and Manual},
776T. Sjostrand, S. Mrenna, P. Skands, \textbf{JHEP} \href{http://dx.doi.org/10.1088/1126-6708/2006/05/026}{05 (2006) 026}.
777\bibitem{bib:cmstauresolution} %\textit{Study of $\tau$-jet identification in CMS},
778R. Kinnunen, A.N. Nikitenko, \textbf{CMS NOTE} \href{http://cdsweb.cern.ch/record/687274}{1997/002}.
779\bibitem{bib:Frog} L. Quertenmont, V. Roberfroid, \textbf{CMS CR} \href{http://cms.cern.ch/iCMS/jsp/openfile.jsp?type=CR&year=2009&files=CR2009_028.pdf}{2009/028}, arXiv:\href{http://arxiv.org/abs/0901.2718}{0901.2718v1}[hep-ex].
780\bibitem{bib:wtphotoproduction} J. de Favereau de Jeneret, S. Ovyn, \textbf{Nucl. Phys. Proc. Suppl.} \href{http://dx.doi.org/10.1016/j.nuclphysbps.2008.07.040}{179-180 (2008)} \href{http://dx.doi.org/10.1016/j.nuclphysbps.2008.07.040}{277-284}; S. Ovyn, J. de Favereau de Jeneret, \href{http://dx.doi.org/10.1393/ncb/i2008-10684-5}{Nuovo Cimento B}, arXiv:0806.4841[hep-ph].
781
782\bibitem{bib:papierquisortirajamais}J. de Favereau~et~al, \textbf{CP3-08-04} (2008), to be published in EPJ.
783
784\bibitem{bib:papiersimon} S. de Visscher, M. Herquet, to be published.
785
786\bibitem{bib:mcfio} P. Lebrun, L. Garren, Copyright (c) 1994-1995 Universities Research Association, Inc.
787
788
789\end{thebibliography}
790
791\onecolumn
792\appendix
793
794\section{User manual}
795
796The available code is a zipped tar file which comes with everything needed to run the \textsc{Delphes} package, assuming a running \textsc{root} installation.
797The package includes \texttt{ExRootAnalysis}~\cite{bib:ExRootAnalysis}, \textsc{Hector}~\cite{bib:Hector},
798\textsc{FastJet}~\cite{bib:FastJet}, and \textsc{Frog}~\cite{bib:Frog}, as well as the conversion codes to read standard \mbox{\textsc{s}td\textsc{hep}} input files (\texttt{mcfio} and \texttt{stdhep})~\cite{bib:mcfio}.
799Nevertheless in order to visualise the events with the \textsc{Frog} software, some external libraries may be required, as explained in \href{http://projects.hepforge.org/frog/}{http://projects.hepforge.org/frog/}.
800
801\subsection{Getting started}
802
803In order to run \textsc{Delphes} on your system, first download its sources and compile it:\\
804\texttt{wget http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes\_V\_*.tar.gz}\\
805Replace the \texttt{*} symbol by the proper version number\footnote{Refer to the download page on the \textsc{Delphes} website \href{http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/download.html}{http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/download.html}.}.
806
807\begin{quote}
808\begin{verbatim}
809me@mylaptop:~$ tar -xvf Delphes_V_*.tar.gz
810me@mylaptop:~$ cd Delphes_V_*.*
811me@mylaptop:~$ ./genMakefile.tcl > Makefile
812me@mylaptop:~$ make
813\end{verbatim}
814\end{quote}
815Due to the large number of external utilities, the number of printed lines during the compilation can be high.
816The user should not pay attention to possible warning messages.
817When compilation is completed, the following message is printed:
818\begin{quote}
819\begin{verbatim}
820me@mylaptop:~$ Delphes has been compiled
821me@mylaptop:~$ Ready to run
822\end{verbatim}
823\end{quote}
824
825\subsection{Running \textsc{Delphes} on your events}
826
827In this chapter, we will explain how to use \textsc{Delphes} to perform a fast simulation of a general purpose detector on your event files. The first step to use \textsc{Delphes} is to create the list of input event files (e.g. {\verb inputlist.list }). As an important comment, don't forget that all the files comprised in the list file should have the same type (\texttt{*.hep}, \texttt{*.lhe} or \texttt{*.root}). In the simplest way of running \textsc{Delphes}, you need this input file and you need to specify the name of the output file that will contain the generator-level data (\texttt{GEN} tree), the analysis data objects after reconstruction (\texttt{Analysis} tree), and the results of the trigger emulation (\texttt{Trigger} tree).
828
829\begin{quote}
830\begin{verbatim}
831me@mylaptop:~$ ./Delphes inputlist.list OutputRootFileName.root
832\end{verbatim}
833\end{quote}
834
835\subsubsection{Setting up the configuration}
836
837The program is driven by two datacards (default cards are {\verb data/DataCardDet.dat } and {\verb data/trigger.dat }) which allow a large spectrum of running conditions.
838Please note that either the user provides these two datacards, either the running will be done using the default parameters defined in the constructor of the class \texttt{RESOLution}. If you choose a different detector or running configuration, you will need to edit the datacards accordingly.
839
840\begin{enumerate}
841
842\item{\bf The smearing card }
843
844The \textit{smearing} or \textit{run} card is by default \texttt{data/DataCard.dat}.
845It contains all pieces of information needed to run \textsc{Delphes}:
846\begin{itemize}
847 \item detector parameters, including calorimeter and tracking coverage and resolution, transverse energy thresholds for object reconstruction and jet algorithm parameters.
848 \item four flags ({\verb FLAG_bfield }, {\verb FLAG_vfd }, {\verb FLAG_trigger } and {\verb FLAG_frog }), which should be assigned if the magnetic field propagation, the very forward detectors simulation, the trigger selection and the preparation for \textsc{Frog} display (respectively) have to be run by \textsc{Delphes}.
849 \end{itemize}
850
851If no datacard is provided by the user, the default smearing and running parameters are used:
852\begin{quote}
853\begin{verbatim}
854# Detector extension, in pseudorapidity units
855CEN_max_tracker 2.5 // Maximum tracker coverage
856CEN_max_calo_cen 3.0 // central calorimeter coverage
857CEN_max_calo_fwd 5.0 // forward calorimeter pseudorapidity coverage
858CEN_max_mu 2.4 // muon chambers pseudorapidity coverage
859
860# Energy resolution for electron/photon
861# \sigma/E = C + N/E + S/\sqrt{E}, E in GeV
862ELG_Scen 0.05 // S term for central ECAL
863ELG_Ncen 0.25 // N term for central ECAL
864ELG_Ccen 0.005 // C term for central ECAL
865ELG_Cfwd 0.107 // S term for FCAL
866ELG_Sfwd 2.084 // C term for FCAL
867ELG_Nfwd 0.0 // N term for FCAL
868
869# Energy resolution for hadrons in ecal/hcal/hf
870# \sigma/E = C + N/E + S/\sqrt{E}, E in GeV
871HAD_Shcal 1.5 // S term for central HCAL
872HAD_Nhcal 0. // N term for central HCAL
873HAD_Chcal 0.05 // C term for central HCAL
874HAD_Shf 2.7 // S term for FCAL
875HAD_Nhf 0. // N term for FCAL
876HAD_Chf 0.13 // C term for FCAL
877
878# Muon smearing
879MU_SmearPt 0.01 // transverse momentum Pt in GeV
880
881# Tracking efficiencies
882TRACK_ptmin 0.9 // minimal pT
883TRACK_eff 100 // efficiency associated to the tracking (%)
884
885
886\end{verbatim}
887\end{quote}
888
889\begin{quote}
890\begin{verbatim}
891# Calorimetric towers
892TOWER_number 40
893### list of the edges of each tower in eta for eta>0 assuming
894###a symmetric detector in eta<0
895### the list starts with the lower edge of the most central tower
896### the list ends with the higher edged of the most forward tower
897### there should be NTOWER+1 values
898TOWER_eta_edges 0. 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783
899 0.870 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653
900 1.740 1.830 1.930 2.043 2.172 2.322 2.500 2.650 2.868 2.950
901 3.125 3.300 3.475 3.650 3.825 4.000 4.175 4.350 4.525 4.700
902 5.000
903
904### list of the tower size in phi (in degrees), assuming that all
905### towers are similar in phi for a given eta value
906### the list starts with the phi-size of the most central tower (eta=0)
907### the list ends with the phi-size of the most forward tower
908### there should be NTOWER values
909#TOWER_dphi 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 10
910 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 20 20
911
912# Thresholds for reconstructed objects, in GeV
913PTCUT_elec 10.0
914PTCUT_muon 10.0
915PTCUT_jet 20.0
916PTCUT_gamma 10.0
917PTCUT_taujet 10.0
918
919# General jet variable
920JET_coneradius 0.7 // generic jet radius
921JET_jetalgo 1 // Jet algorithm selection
922JET_seed 1.0 // minimum seed to start jet reconstruction, in GeV
923
924# Tagging definition
925BTAG_b 40 // b-tag efficiency (%)
926BTAG_mistag_c 10 // mistagging (%)
927BTAG_mistag_l 1 // mistagging (%)
928
929# FLAGS
930FLAG_bfield 0 // 1 to run the bfield propagation else 0
931FLAG_vfd 1 // 1 to run the very forward detectors else 0
932FLAG_trigger 1 // 1 to run the trigger selection else 0
933FLAG_frog 1 // 1 to run the FROG event display
934
935# In case BField propagation allowed
936TRACK_radius 129 // radius of the BField coverage, in cm
937TRACK_length 300 // length of the BField coverage, in cm
938TRACK_bfield_x 0 // X composant of the BField, in T
939TRACK_bfield_y 0 // Y composant of the BField, in T
940TRACK_bfield_z 3.8 // Z composant of the BFieldn in T
941
942# Very forward detector extension, in pseudorapidity
943# if allowed
944VFD_min_calo_vfd 5.2 // very forward calorimeter (if any) like CASTOR
945VFD_max_calo_vfd 6.6
946VFD_min_zdc 8.3
947VFD_s_zdc 140 // distance of the ZDC, from the IP, in [m]
948
949RP_220_s 220 // distance of the RP to the IP, in meters
950RP_220_x 0.002 // distance of the RP to the beam, in meters
951RP_420_s 420 // distance of the RP to the IP, in meters
952RP_420_x 0.004 // distance of the RP to the beam, in meters
953
954# In case FROG event display allowed
955NEvents_Frog 100
956
957\end{verbatim}
958\end{quote}
959In general, energies and momenta are expressed in GeV, and magnetic fields in T.
960Geometrical extension are often referred in terms of pseudorapidity $\eta$, as the detectors are supposed to be symmetric in $\phi$.
961
962\item{\bf The trigger card }
963
964This card contains the definitions of all trigger bits. Cuts can be applied on the transverse momentum $p_T$ of electrons, muons, jets, $\tau$-jets, photons and the missing transverse energy. The following codes should be used so that \textsc{Delphes} can correctly translate the input list of trigger bits into selection algorithms:
965
966\begin{quote}
967\begin{tabular}{ll}
968{\it Trigger code} & {\it Corresponding object}\\
969{\verb ELEC_PT } & electron \\
970{\verb MUON_PT } & muon \\
971{\verb JET_PT } & jet \\
972{\verb TAUJET_PT } & $\tau$-jet \\
973{\verb ETMIS_PT } & missing transverse energy \\
974{\verb GAMMA_PT } & photon \\
975\end{tabular}
976\end{quote}
977
978Each line in the trigger datacard is allocated to exactly one trigger bit and starts with the name of the corresponding trigger.
979Logical combination of several conditions is also possible.
980If the trigger bit requires the presence of multiple identical objects, the order of their $p_T$ thresholds is very important: they must be defined in \textit{decreasing} order. Finally, the different requirements on the objects must be separated by a {\verb && } flag.
981The default trigger card can be found in the data repository of \textsc{Delphes} (\texttt{data/trigger.dat}).
982An example of trigger table consistent with the previous rules is given here:
983\begin{quote}
984\begin{verbatim}
985SingleJet >> JET_PT: '200'
986DoubleElec >> ELEC_PT: '20' && ELEC_PT: '10'
987SingleElec and Single Muon >> ELEC_PT: '20' && MUON_PT: '15'
988\end{verbatim}
989\end{quote}
990\end{enumerate}
991
992\subsubsection{Running the code}
993
994First, create the smearing and trigger cards (\texttt{data/mydetector.dat} and \texttt{data/mytrigger.dat}). \\
995Then, create a text file containing the list of input files that will be used by \textsc{Delphes} (with extension \texttt{*.lhe}, \texttt{*.root} or \texttt{*.hep}).
996To run the code, type the following command (in one line)
997\begin{quote}
998\begin{verbatim}
999me@mylaptop:~$ ./Delphes inputlist.list OutputRootFileName.root
1000 data/mydetector.dat data/mytrigger.dat
1001\end{verbatim}
1002\end{quote}
1003As a reminder, typing the \texttt{./Delphes} command simply displays the correct usage:
1004
1005\begin{quote}
1006\begin{verbatim}
1007me@mylaptop:~$ ./Delphes
1008 Usage: ./Delphes input_file output_file [detector_card] [trigger_card]
1009 input_list - list of files in Ntpl, StdHep of LHEF format,
1010 output_file - output file.
1011 detector_card - Card containing resolution variables for detector simulation (optional)
1012 trigger_card - Card containing the trigger algorithms (optional)
1013\end{verbatim}
1014\end{quote}
1015
1016
1017\subsection{Getting the \textsc{Delphes} information}
1018
1019\subsubsection{Contents of the \textsc{Delphes} ROOT trees}
1020
1021The \textsc{Delphes} output file (\texttt{*.root}) is subdivided into three \textit{trees}, corresponding to generator-level data, analysis object data and trigger output. These \textit{trees} are structures that organise the output data into \textit{branches} containing data (or \textit{leaves}) related with each others, like the kinematics properties ($E$, $p_x$, $\eta$, $\ldots$) of a given particle.
1022
1023Here is the exhaustive list of \textit{branches} availables in these \textit{trees}, together with their corresponding physical objet and \texttt{ExRootAnalysis} class:
1024\begin{quote}
1025\begin{tabular}{lll}
1026{\bf GEN \textsc{tree}} & &\\
1027~~~Particle & generator particles from \textsc{hepevt} & {\verb TRootGenParticle }\\
1028{\bf Trigger } & &\\
1029~~~TrigResult & Acceptance of different trigger bits & {\verb TRootTrigger }\\
1030\end{tabular}
1031\end{quote}
1032\begin{quote}
1033\begin{tabular}{lll}
1034{\bf Analysis \textsc{tree}} & & \\
1035~~~Tracks & Collection of tracks & {\verb TRootTracks }\\
1036~~~CaloTower & Calorimetric towers & {\verb TRootCalo }\\
1037~~~Electron & Collection of electrons & {\verb TRootElectron }\\
1038~~~Photon & Collection of photons & {\verb TRootPhoton }\\
1039~~~Muon & Collection of muons & {\verb TRootMuon }\\
1040~~~Jet & Collection of jets & {\verb TRootJet }\\
1041~~~TauJet & Collection of jets tagged as $\tau$-jets & {\verb TRootTauJet }\\
1042~~~ETmis & Transverse missing energy information & {\verb TRootETmis }\\
1043~~~ZDChits & Hits in the Zero Degree Calorimeters & {\verb TRootZdcHits }\\
1044~~~RP220hits & Hits in the first proton taggers & {\verb TRootRomanPotHits }\\
1045~~~FP420hits & Hits in the next proton taggers & {\verb TRootRomanPotHits }\\
1046\end{tabular}
1047\end{quote}
1048The third column shows the names of the corresponding classes to be written in a \textsc{root} tree.
1049All classes except \texttt{TRootTrigger}, \texttt{TRootETmis} and \texttt{TRootRomanPotHits} inherit from the class \texttt{TRootParticle} which includes the following data members (stored as \textit{leaves} in \textit{branches} of the \textit{trees}):
1050\begin{quote}
1051\begin{tabular}{ll}
1052\multicolumn{2}{l}{\textbf{Most common leaves}}\\
1053 \texttt{~~~float E; }&\texttt{ // particle energy in GeV }\\
1054 \texttt{~~~float Px; }&\texttt{ // particle momentum vector (x component) in GeV }\\
1055 \texttt{~~~float Py; }&\texttt{ // particle momentum vector (y component) in GeV }\\
1056 \texttt{~~~float Pz; }&\texttt{ // particle momentum vector (z component) in GeV }\\
1057 \texttt{~~~float PT; }&\texttt{ // particle transverse momentum in GeV }\\
1058 \texttt{~~~float Eta; }&\texttt{ // particle pseudorapidity }\\
1059 \texttt{~~~float Phi; }&\texttt{ // particle azimuthal angle in rad }\\
1060\end{tabular}
1061\end{quote}
1062
1063In addition to their kinematics, some additional properties are available for specific objects:
1064\begin{quote}
1065\begin{tabular}{ll}
1066\multicolumn{2}{l}{{\bf Leaves in the \texttt{Particle} branch}} \\
1067 \texttt{~~~int PID; }&\texttt{ // particle HEP ID number }\\
1068 \texttt{~~~int Status; }&\texttt{ // particle status }\\
1069 \texttt{~~~int M1; }&\texttt{ // particle 1st mother }\\
1070 \texttt{~~~int M2; }&\texttt{ // particle 2nd mother }\\
1071 \texttt{~~~int D1; }&\texttt{ // particle 1st daughter }\\
1072 \texttt{~~~int D2; }&\texttt{ // particle 2nd daughter }\\
1073 \texttt{~~~float Charge; }&\texttt{ // electrical charge in units of e}\\
1074 \texttt{~~~float T; }&\texttt{ // particle vertex position (t component, in mm/c) }\\
1075 \texttt{~~~float X; }&\texttt{ // particle vertex position (x component, in mm) }\\
1076 \texttt{~~~float Y; }&\texttt{ // particle vertex position (y component, in mm) }\\
1077 \texttt{~~~float Z; }&\texttt{ // particle vertex position (z component, in mm) }\\
1078 \texttt{~~~float M; }&\texttt{ // particle mass in GeV}\\
1079\end{tabular}
1080\end{quote}
1081\begin{quote}
1082\begin{tabular}{ll}
1083\multicolumn{2}{l}{\textbf{Additional leaves in \texttt{Electron} and \texttt{Muon} branches}} \\
1084 \texttt{~~~int Charge } &\\
1085 \texttt{~~~bool IsolFlag } &\\
1086& \\
1087\multicolumn{2}{l}{\textbf{Additional leaf in the \texttt{Jet} branch}} \\
1088 \texttt{~~~bool Btag } &\\
1089& \\
1090\multicolumn{2}{l}{\textbf{Additional leaves in the \texttt{ZDChits} branch}}\\
1091 \texttt{~~~float T; }&\texttt{ // time of flight in s }\\
1092 \texttt{~~~int side; }&\texttt{ // -1 or +1 }\\
1093\end{tabular}
1094\end{quote}
1095
1096\subsection{Running an analysis on your \textsc{Delphes} events}
1097
1098To analyse the \textsc{root} ntuple produced by \textsc{Delphes}, the simplest way is to use the {\verb Analysis_Ex.cpp } code which is coming in the {\verb Examples } repository of \textsc{Delphes}. Note that all of this is optional and done to facilitate the analyses, as the output from \textsc{Delphes} is viewable with the standard \textsc{root} \texttt{TBrowser} and can be analysed using the \texttt{MakeClass} facility.
1099As an example, here is a simple overview of a \texttt{myoutput.root} file created by \textsc{Delphes}:
1100\begin{quote}
1101\begin{verbatim}
1102me@mylaptop:~$ root -l myoutput.root
1103root [0]
1104Attaching file myoutput.root as _file0...
1105root [1] .ls
1106TFile** myoutput.root
1107 TFile* myoutput.root
1108 KEY: TTree GEN;1 Analysis tree
1109 KEY: TTree Analysis;1 Analysis tree
1110 KEY: TTree Trigger;1 Analysis tree
1111root [2] TBrowser t;
1112root [3] Analysis->GetEntries()
1113(const Long64_t)200
1114root [4] GEN->GetListOfBranches()->ls()
1115OBJ: TBranchElement Event Event_ : 0 at: 0x9108f30
1116OBJ: TBranch Event_size Event_size/I : 0 at: 0x910cfd0
1117OBJ: TBranchElement Particle Particle_ : 0 at: 0x910c6b0
1118OBJ: TBranch Particle_size Particle_size/I : 0 at: 0x9111c58
1119root [5] Trigger->GetListOfLeaves()->ls()
1120OBJ: TLeafElement TrigResult_ TrigResult_ : 0 at: 0x90f90a0
1121OBJ: TLeafElement TrigResult.Accepted Accepted[TrigResult_] : 0 at: 0x90f9000
1122OBJ: TLeafI TrigResult_size TrigResult_size : 0 at: 0x90fb860
1123\end{verbatim}
1124\end{quote}
1125The \texttt{.ls} command lists the current keys available and in particular the three \textit{tree} names.
1126\mbox{\texttt{TBrowser t}} launches a browser and the \texttt{GetEntries()} method outputs the number of data in the corresponding \textit{tree}.
1127The list of \textit{branches} or \textit{leaves} can be displayed with the \texttt{GetListOfBranches()} and \texttt{GetListOfLeaves()} methods, pointing to the \texttt{ls()} one. In particular, it is possible to shown only parts of the output, using wildcard characters (\texttt{*}):
1128\begin{quote}
1129\begin{verbatim}
1130root [6] Analysis->GetListOfLeaves()->ls("*.E")
1131OBJ: TLeafElement Jet.E E[Jet_] : 0 at: 0xa08bc68
1132OBJ: TLeafElement TauJet.E E[TauJet_] : 0 at: 0xa148910
1133OBJ: TLeafElement Electron.E E[Electron_] : 0 at: 0xa1d8a50
1134OBJ: TLeafElement Muon.E E[Muon_] : 0 at: 0xa28ac80
1135OBJ: TLeafElement Photon.E E[Photon_] : 0 at: 0xa33cd88
1136OBJ: TLeafElement Tracks.E E[Tracks_] : 0 at: 0xa3cced0
1137OBJ: TLeafElement CaloTower.E E[CaloTower_] : 0 at: 0xa4ba188
1138OBJ: TLeafElement ZDChits.E E[ZDChits_] : 0 at: 0xa54a3c8
1139OBJ: TLeafElement RP220hits.E E[RP220hits_] : 0 at: 0xa61e648
1140OBJ: TLeafElement FP420hits.E E[FP420hits_] : 0 at: 0xa6d0920
1141\end{verbatim}
1142\end{quote}
1143
1144To draw a particular leaf, either double-click on the corresponding name in the \texttt{TBrowser} or use the \texttt{Draw} method of the corresponding \textit{tree}.
1145\begin{quote}
1146\begin{verbatim}
1147root [7] Trigger->Draw("TrigResult.Accepted");
1148\end{verbatim}
1149\end{quote}
1150Mathematical operations on several \textit{leaves} are possible within a given \textit{tree}:
1151\begin{quote}
1152\begin{verbatim}
1153root [8] Analysis->Draw("Muon.Px * Muon.Px");
1154root [9] Analysis->Draw("sqrt(pow(Muon.E,2) - pow(Muon.Pz,2) + pow(Muon.PT,2))");
1155\end{verbatim}
1156\end{quote}
1157Finally, to prepare an deeper analysis, the \texttt{MakeClass} method is useful:
1158\begin{quote}
1159\begin{verbatim}
1160root [10] Trigger->MakeClass()
1161Info in <TTreePlayer::MakeClass>: Files: Trigger.h and
1162 Trigger.C generated from TTree: Trigger
1163\end{verbatim}
1164\end{quote}
1165
1166To run the \texttt{Examples/Analysis\_Ex.cpp} code, the two following arguments are required: a text file containing the input \textsc{Delphes} \textsc{root} files to run, and the name of the output \textsc{root} file.
1167 \begin{quote}
1168\begin{verbatim}
1169me@mylaptop:~$ ./Analysis_Ex input_file.list output_file.root
1170\end{verbatim}
1171 \end{quote}
1172
1173\subsubsection{Adding the trigger information}
1174The \texttt{Examples/Trigger\_Only.cpp} code permits to run the trigger selection separately from the general detector simulation on output \textsc{Delphes} root files.
1175A \textsc{Delphes} root file is mandatory as an input argument for the \texttt{Trigger\_Only} routine.
1176The new \textit{tree} containing the trigger result data will be appended to this file.
1177The trigger datacard is also necessary. To run the code:
1178 \begin{quote}
1179\begin{verbatim}
1180me@mylaptop:~$ ./Trigger_Only input_file.root data/trigger.dat
1181\end{verbatim}
1182 \end{quote}
1183
1184\subsection{Running the FROG event display}
1185
1186\begin{itemize}
1187\item If the { \verb FLAG_frog } was switched on in the smearing card, two files have been created during the running of \textsc{Delphes}: {\verb DelphesToFrog.vis } and {\verb DelphesToFrog.geom }. They contain all the needed pieces of information to run \textsc{frog}.
1188\item To display the events and the geometry, you first need to compile \textsc{Frog}. Go to the {\verb Utilities/FROG } and type {\verb make }. This compilation is done once for all, with this geometry (i.e. as long as the \texttt{*vis} and \texttt{*geom} files do not change).
1189\item Go back into the main directory and type
1190\begin{quote}
1191\texttt{me@mylaptop:~\$ ./Utilities/FROG/frog}.
1192\end{quote}
1193\end{itemize}
1194
1195\end{document}
1196
1197%[25] ATLAS Collaboration, Detector and Physics Performance Technical Design
1198% Report, Vols. 1 and 2, CERN–LHCC–99–14 and CERN–LHCC–99–15.
1199%[26] CMS Collaboration, CMS Physics Technical Design Report, CERN/LHCC 2006–001.
1200%[27] A. Djouadi, J. Lykken, K. Monig, Y. Okada, M. J. Oreglia and S. Yamashita,
1201% International Linear Collider Reference Design Report Volume 2: PHYSICS
1202% AT THE ILC, arXiv:0709.1893 [hep-ph].
1203
1204% personnes qui pourraient être intéressées:
1205% Alice, Benjamin
1206% auteurs de arXiv:0801.3359
Note: See TracBrowser for help on using the repository browser.