Fork me on GitHub

Changeset 532 in svn


Ignore:
Timestamp:
Sep 18, 2009, 5:57:14 PM (15 years ago)
Author:
Xavier Rouby
Message:

Update commentaires VL

Location:
trunk/paper/CommPhysComp
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/paper/CommPhysComp/notes.tex

    r525 r532  
    3535\ead{severine.ovyn@uclouvain.be}
    3636
    37 \author{X. Rouby\fnref{freiburg}}
    38 \fntext[freiburg]{Now in Physikalisches Institut, Albert-Ludwigs-Universit\"at Freiburg}
     37\author{X. Rouby}
     38%\author{X. Rouby\fnref{freiburg}}
     39%\fntext[freiburg]{Now in Physikalisches Institut, Albert-Ludwigs-Universit\"at Freiburg}
    3940%\ead{xavier.rouby@cern.ch}
    4041
     
    5354
    5455\begin{abstract}
    55 It is always delicate to  know whether theoretical predictions are visible and measurable in a high energy collider experiment due to the complexity of the related detectors, data acquisition chain and software.
    56 We introduce here a new \texttt{C++}-based framework, \textit{Delphes}, for fast simulation of
    57 a general-purpose experiment. The simulation includes a tracking system, embedded into a magnetic field, calorimetry and a muon
    58 system, and possible very forward detectors arranged along the beamline.
    59 The framework is interfaced to standard file formats (e.g.\ Les Houches Event File or \texttt{HepMC}) and outputs observable objects for analysis, like missing transverse energy and collections of electrons or jets.
    60 The simulation of detector response takes into account the detector resolution, and usual reconstruction algorithms, such as 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 \textit{Hector} software. Finally, the \textsc{FROG} 2D/3D event display is used for visualisation of the collision final states.
    61 An overview of \textit{Delphes} is given as well as a few \textsc{LHC} use-cases for illustration.\\ \\
     56% It is always delicate to  know whether theoretical predictions are visible and measurable in a high energy collider experiment due to the complexity of the related detectors, data acquisition chain and software.
     57% We introduce here a new \texttt{C++}-based framework, \textit{Delphes}, for fast simulation of
     58% a general-purpose experiment. The simulation includes a tracking system, embedded into a magnetic field, calorimetry and a muon
     59% system, and possible very forward detectors arranged along the beamline.
     60% The framework is interfaced to standard file formats (e.g.\ Les Houches Event File or \texttt{HepMC}) and outputs observable objects for analysis, like missing transverse energy and collections of electrons or jets.
     61% The simulation of detector response takes into account the detector resolution, and usual reconstruction algorithms, such as 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 \textit{Hector} software. Finally, the \textsc{FROG} 2D/3D event display is used for visualisation of the collision final states.
     62% An overview of \textit{Delphes} is given as well as a few \textsc{LHC} use-cases for illustration.\\ \\
     63
     64It is sometimes difficult to know whether theoretical predictions can be observed in a high energy collider experiment, especially when expected experimental signature involve jets and missing transverse energy.
     65For this purpose, we have designed a new \texttt{C++}-based framework, \textit{Delphes}, performing a fast multipurpose detector response simulation.
     66The simulation includes a tracking system, embedded into a magnetic field, calorimeters and a muon system, and possible very forward detectors arranged along the beamline.
     67The framework is interfaced to standard file formats (e.g.\ Les Houches Event File or \texttt{HepMC}) and outputs observables such as isolated leptons, missing transverse energy and collection of jets which can be used for dedicated analyses.
     68The simulation of the detector response takes into account the effect of magnetic field, the granularity of the calorimeters and subdetector resolutions.
     69A simplified preselection can also be applied on processed events for trigger emulation. Detection of very forward scattered particles relies on the transport in beamlines with the \textit{Hector} software. Finally, the \textsc{FROG} 2D/3D event display is used for visualisation of the collision final states.
     70\\ \\
     71
     72
    6273\textit{Preprint:} \texttt{CP3-09-01}, \texttt{arXiv:0903.2225 [hep-ph]}\\ \\
    6374%\includegraphics[scale=0.8]{DELPHESLogoSml}\\
     
    117128\section{Introduction}
    118129
    119 Experiments 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. Such apparatus differ in their detection principles, technologies, geometrical acceptances, resolutions 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 of ``uninteresting'' events, 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 approaches.
    120 
    121 This 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 transport of the primary and secondary particles through the detector material accounting for the various detector inefficiencies, the dead material, the imperfections and the geometrical details. Moreover, control of the 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 expected signals and associated backgrounds.
     130% Experiments 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. Such apparatus differ in their detection principles, technologies, geometrical acceptances, resolutions 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 of ``uninteresting'' events, 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 approaches.
     131
     132Multipurpose detectors at high energy colliders are very complex systems. Their simulation is in general performed by means of the GEANT~\citep{bib:geant} package and final observables used for analyses usually require sophisticated reconstruction algorithms.
     133
     134
     135This complexity is handled by large collaborations, and data and the expertise on reconstruction and simulation software are only available to their members. Precise data analyses require a full detector simulation, including transport of the primary and secondary particles through the detector material accounting for the various detector inefficiencies, the dead material, the imperfections and the geometrical details.
     136%\textcolor{blue}{Moreover, control of the detector calibration and alignment are crucial}.
     137Such 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 expected signal signatures and their associated backgrounds.
    122138
    123139A new framework, called \textit{Delphes}~\citep{bib:delphes}, is introduced here, for the fast simulation of a general-purpose collider experiment.
    124140Using the framework, observables can be estimated for specific signal and background channels, as well as their production and measurement rates.
    125 Starting 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.
     141Starting 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 (i.e. those considered as stable by the event generator
     142\footnote{In the current \textit{Delphes} version, particles other than electrons ($e^\pm$), photons ($\gamma$), muons ($\mu^\pm$), neutrinos ($\nu_e$, $\nu_\mu$ and $\nu_\tau$) and neutralinos are simulated as hadrons for their interactions with the calorimeters. The simulation of stable particles beyond the Standard Model should therefore be handled with care~\citep{qr:invisibleparticles}.}). Tracks of charged particles and deposits of energy in calorimetric cells (or \textit{calotowers}) are then created. These two types of quantities are used for the reconstruction of jets and the isolation of leptons.
    126143
    127144\textit{Delphes} includes the most crucial experimental features, such as (Fig.~\ref{fig:FlowChart}):
    128145\begin{enumerate}
    129146\item the geometry of both central and forward detectors,
    130 \item magnetic field for tracks
     147\item magnetic field for tracks and energy flow
    131148\item reconstruction of photons, leptons, jets, $b$-jets, $\tau$-jets and missing transverse energy,
    132149\item lepton isolation,
     
    141158\caption{Flow chart describing the principles behind \textit{Delphes}. Event files coming from external Monte Carlo generators are read by a converter stage (top).
    142159The kinematics variables of the final-state particles are then smeared according to the tunable subdetector resolutions.
    143 Tracks are reconstructed in a simulated solenoidal magnetic field and calorimetric towers sample the energy deposits. Based on these low-level objects, dedicated algorithms are applied for particle identification, isolation and reconstruction.
     160Tracks are reconstructed in a simulated solenoidal magnetic field and calorimetric cells sample the energy deposits. Based on these low-level objects, dedicated algorithms are applied for particle identification, isolation and reconstruction.
    144161The transport of very forward particles to the near-beam detectors is also simulated.
    145162Finally, an output file is written, including generator-level and analysis-object data.
     
    152169Although 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.
    153170
    154 Four datafile formats can be used as input in \textit{Delphes}\footnote{\texttt{[code] }See the \texttt{HEPEVTConverter}, \texttt{HepMCConverter}, \texttt{LHEFConverter} and \texttt{STDHEPConverter} classes.}. In order to process events from many different generators, the standard Monte Carlo event structures \texttt{StdHEP}~\citep{bib:stdhep} and \texttt{HepMC}~\citep{bib:hepmc} can be used as an input. Besides, \textit{Delphes} can also provide detector response for events read in ``Les Houches Event Format'' (\textsc{LHEF}~\citep{bib:lhe}) and \texttt{*.root} files obtained from \texttt{*.hbook} using the \texttt{h2root} utility from the \textsc{ROOT} framework~\citep{bib:Root}.
     171Several datafile formats can be used as input in \textit{Delphes} \citep{qr:inputformat},
     172%\footnote{\texttt{[code] }See the \texttt{HEPEVTConverter}, \texttt{HepMCConverter}, \texttt{LHEFConverter} and \texttt{STDHEPConverter} classes.}.
     173in order to process events from many different generators. The standard Monte Carlo event structures \texttt{StdHEP}~\citep{bib:stdhep} and \texttt{HepMC}~\citep{bib:hepmc} can be used as an input. Besides, \textit{Delphes} can also provide detector response for events read in ``Les Houches Event Format'' (\textsc{LHEF}~\citep{bib:lhe}) and \texttt{*.root} files obtained from \texttt{*.hbook} using the \texttt{h2root} utility from the \textsc{ROOT} framework~\citep{bib:Root}.
    155174%Afterwards, \textit{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.
    156175
    157176\textit{Delphes} uses the \texttt{ExRootAnalysis} utility~\citep{bib:ExRootAnalysis} to create output data in a \texttt{*.root} ntuple.
    158177This output contains a copy of the generator-level data (\texttt{GEN} tree), the analysis data objects after reconstruction (\texttt{Analysis} tree), and possibly the results of the trigger emulation (\texttt{Trigger} tree).
    159 In option\footnote{\texttt{[code]} See the \texttt{FLAG\_LHCO} variable in the detector datacard. This text file format is shortly described in the user manual.}, \textit{Delphes} can produce a reduced output file in \texttt{*.LHCO} text format, which is limited to the list of the reconstructed high-level objects in the final states.
    160 
    161 The program is driven by input cards. The detector card (\texttt{data/DetectorCard.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/TriggerCard.dat}) lists the user algorithms for the simplified online preselection. Even if \textit{Delphes} has been developped for the simulation of general-purpose detectors at the \textsc{LHC} (namely, \textsc{CMS} and \textsc{ATLAS}), the input cards allow a flexible parametrisation for other cases, e.g.\ at future linear colliders.
    162 
    163 
    164 \section{Detector simulation}
    165 
    166 The overall layout of the general-purpose detector simulated by \textit{Delphes} is shown in Fig.~\ref{fig:GenDet3}.
    167 A central tracking system (\textsc{TRACKER}) is surrounded by an electromagnetic and a hadron calorimeters (\textsc{ECAL} and \textsc{HCAL}, resp., each with a central region and two endcaps). 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
    168 The fast simulation of the detector response takes into account geometrical acceptance of sub-detectors and their finite resolution, as defined in the detector data card\footnote{\texttt{[code] }See the \texttt{RESOLution} class.}.
    169 If no such file is provided, predefined values based on ``typical'' \textsc{CMS} acceptances and resolutions are used\footnote{\texttt{[code] }Detector and trigger cards for the \textsc{ATLAS} and \textsc{CMS} experiments are also provided in \texttt{data/} directory.}. The geometrical coverage of the various subsystems used in the default configuration are summarised in Tab.~\ref{tab:defEta}.
    170 
    171 \begin{table*}[t]
     178In option
     179%\footnote{\texttt{[code]} See the \texttt{FLAG\_LHCO} variable in the detector datacard. This text file format is shortly described in the user manual.},
     180\textit{Delphes} can produce a reduced output file in \texttt{*.lhco} text format, which is limited to the list of the reconstructed high-level objects in the final states~\citep{qr:lhco}.
     181
     182
     183
     184\section{Simulation of the detector response}
     185
     186The overall layout of the multipurpose detector simulated by \textit{Delphes} is shown in Fig.~\ref{fig:GenDet3}.
     187It consists in a central tracking system (\textsc{TRACKER}) surrounded by an electromagnetic and a hadron calorimeters (\textsc{ECAL} and \textsc{HCAL}, each with a central region and two endcaps). 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.
     188A detector card \citep{qr:detectorcard} allows a large spectrum of running conditions by modifying basic detector parameters, including calorimeter and tracking coverage and resolution, thresholds or jet algorithm parameters.
     189Even if \textit{Delphes} has been developped for the simulation of general-purpose detectors at the \textsc{LHC} (namely, \textsc{CMS} and \textsc{ATLAS}), this input parameter file interfaces a flexible parametrisation for other cases, e.g.\ at future linear colliders~\citep{qr:datacards}.
     190If no detector card is provided, predefined values based on ``typical'' \textsc{CMS} acceptances and resolutions are used.
     191%\footnote{\texttt{[code] }Detector and trigger cards for the \textsc{ATLAS} and \textsc{CMS} experiments are also provided in \texttt{data/} directory.}.
     192The geometrical coverage of the various subsystems used in the default configuration are summarised in Tab.~\ref{tab:defEta}.
     193
     194\begin{table}[t]
     195% \begin{table*}[t]
    172196\begin{center}
    173197\caption{Default extension in pseudorapidity $\eta$ of the different subdetectors.
    174198Full azimuthal ($\phi$) acceptance is assumed.
    175 The corresponding parameter name, in the detector card, is given. \vspace{0.5cm}}
    176 \begin{tabular}{llcc}
     199 \vspace{0.5cm}}
     200% \begin{tabular}{llcc}
     201% \hline
     202% Subdetector & & $\eta$ & $\phi$ \\
     203% \textsc{TRACKER}      & {\verb CEN_max_tracker }      & $[-2.5; 2.5]$         & $[-\pi ; \pi]$\\
     204% \textsc{ECAL}, \textsc{HCAL} & {\verb CEN_max_calo_cen }& $[-1.7 ; 1.7]$      & $[-\pi ; \pi]$\\
     205% \textsc{ECAL}, \textsc{HCAL} endcaps & {\verb CEN_max_calo_ec }& $[-3 ; -1.7] \& [1.7 ; 3]$   & $[-\pi ; \pi]$\\
     206% \textsc{FCAL}                 & {\verb CEN_max_calo_fwd }     & $[-5 ; -3]$ \& $[3 ;5]$     & $[-\pi ; \pi]$\\
     207% \textsc{MUON}                 & {\verb CEN_max_mu }           & $[-2.4 ; 2.4]$        & $[-\pi ; \pi]$\\ \hline
     208% \end{tabular}
     209\begin{tabular}{lcc}
    177210\hline
    178 Subdetector & & $\eta$ & $\phi$ \\
    179 \textsc{TRACKER}        & {\verb CEN_max_tracker }      & $[-2.5; 2.5]$         & $[-\pi ; \pi]$\\
    180 \textsc{ECAL}, \textsc{HCAL} & {\verb CEN_max_calo_cen }& $[-1.7 ; 1.7]$        & $[-\pi ; \pi]$\\
    181 \textsc{ECAL}, \textsc{HCAL} endcaps & {\verb CEN_max_calo_ec }& $[-3 ; -1.7] \& [1.7 ; 3]$     & $[-\pi ; \pi]$\\
    182 \textsc{FCAL}           & {\verb CEN_max_calo_fwd }     & $[-5 ; -3]$ \& $[3 ;5]$     & $[-\pi ; \pi]$\\
    183 \textsc{MUON}           & {\verb CEN_max_mu }           & $[-2.4 ; 2.4]$        & $[-\pi ; \pi]$\\ \hline
     211 & $\eta$ & $\phi$ \\ \hline
     212\textsc{TRACKER}        & $[-2.5; 2.5]$                 & $[-\pi ; \pi]$\\
     213\textsc{ECAL}, \textsc{HCAL} & $[-1.7 ; 1.7]$           & $[-\pi ; \pi]$\\
     214\textsc{ECAL}, \textsc{HCAL} endcaps & $[-3 ; -1.7]$ \& $[1.7 ; 3]$     & $[-\pi ; \pi]$\\
     215\textsc{FCAL}           & $[-5 ; -3]$ \& $[3 ;5]$       & $[-\pi ; \pi]$\\
     216\textsc{MUON}           & $[-2.4 ; 2.4]$                & $[-\pi ; \pi]$\\ \hline
    184217\end{tabular}
    185218\label{tab:defEta}
    186219\end{center}
    187 \end{table*}
     220% \end{table*}
     221\end{table}
    188222
    189223\begin{figure}[!ht]
     
    202236
    203237
    204 \subsubsection*{Magnetic field}
    205 In addition to the subdetectors, the effects of a solenoidal magnetic field are simulated for the charged particles\footnote{\texttt{[code] }See the \texttt{TrackPropagation} class.}. This affects the position at which charged particles enter the calorimeters and their corresponding tracks. The field extension is limited to the tracker volume and is in particular not applied for muon chambers. Howerver, this is not a limiting factor as the resolution applied for muon reconstruction is the one expected by the experiment, which consequently includes the effects of the magnetic field within the muon system.
    206 
     238\subsection{Magnetic field}
     239In addition to the subdetectors, the effects of a solenoidal magnetic field are simulated for the charged particles~\citep{qr:magneticfield}
     240%\footnote{\texttt{[code] }See the \texttt{TrackPropagation} class.}
     241. This affects the position at which charged particles enter the calorimeters and their corresponding tracks. The field extension is limited to the tracker volume and is in particular not applied for muon chambers. Howerver, this is not a limiting factor as the resolution applied for muon reconstruction is the one expected by the experiment, which consequently includes the effects of the magnetic field within the muon system.
    207242
    208243
    209244\subsection{Tracks reconstruction}
    210245Every stable charged particle with a transverse momentum above some threshold and lying inside the detector volume covered by the tracker provides a track.
    211 By default, a track is assumed to be reconstructed with $90\%$ probability\footnote{\texttt{[code]} The reconstruction efficiency is defined in the detector datacard by the \texttt{TRACKING\_EFF} term.} if its transverse momentum $p_T$ is higher than $0.9~\textrm{GeV}/c$ and if its pseudorapidity $|\eta| \leq 2.5$.
    212 
    213 
    214 \subsection{Simulation of central calorimeters}
    215 
    216 The 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.
    217 The response of each sub-calorimeter is parametrised as a function of the energy:
     246By default, a track is assumed to be reconstructed with $90\%$ probability
     247%\footnote{\texttt{[code]} The reconstruction efficiency is defined in the detector datacard by the \texttt{TRACKING\_EFF} term.}
     248if its transverse momentum $p_T$ is higher than $0.9~\textrm{GeV}/c$ and if its pseudorapidity
     249$|\eta| \leq 2.5$~\citep{qr:tracks}. No smearing is currently applied on tracks.
     250
     251
     252\subsection{Calorimetric cells}
     253
     254The response of the calorimeters to energy deposits of incoming particles depends on their segmentation and resolution. In CMS and ATLAS detectors, for instance, the calorimeter characteristics are not identical in every direction, with typically finer resolution and granularity in the central regions~\citep{bib:cmsjetresolution,bib:ATLASresolution}. It is thus very important to compute the exact coordinates of the entry point of the particles into the calorimeters, via the magnetic field calculations.
     255
     256The response of each sub-calorimeter is parametrised through a Gaussian smearing of the particle energy with a variance $\sigma$:
    218257\begin{equation}
    219258\frac{\sigma}{E} = \frac{S}{\sqrt{E}} \oplus \frac{N}{E} \oplus C,
    220259\label{eq:caloresolution}
    221260\end{equation}
    222 where $S$, $N$ and $C$ are the \textit{stochastic}, \textit{noise} and \textit{constant} terms, respectively, and $\oplus$ stands for quadratic additions.\\
    223 
    224 
    225 The particle four-momentum $p^\mu$ are smeared with a parametrisation directly derived from typical detector technical designs\footnote{\texttt{[code] } The response of the detector is applied to the electromagnetic and the hadronic particles through the \texttt{SmearElectron} and \texttt{SmearHadron} functions.} \citep{bib:cmsjetresolution,bib:ATLASresolution}.
     261where $S$, $N$ and $C$ are the \textit{stochastic}, \textit{noise} and \textit{constant} terms, respectively, and $\oplus$ stands for quadratic additions~\citep{qr:energysmearing}.\\
     262
     263%\footnote{\texttt{[code] } The response of the detector is applied to the electromagnetic and the hadronic particles through the \texttt{SmearElectron} and \texttt{SmearHadron} functions.}
    226264In the default parametrisation, the calorimeter is assumed to cover the pseudorapidity range $|\eta|<3$ and consists in an electromagnetic and hadronic parts. 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.
    227 Muons and neutrinos are assumed not to interact with the calorimeters\footnote{In the current \textit{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 therefore be handled with care.}.
     265Muons and neutrinos are assumed not to interact with the calorimeters~\citep{qr:invisibleparticles}.
     266%\footnote{In the current \textit{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 therefore be handled with care.}.
    228267The default values of the stochastic, noise and constant terms are given in Tab.~\ref{tab:defResol}.\\
    229268
    230269\begin{table}[!h]
    231270\begin{center}
    232 \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}).
    233 The corresponding parameter name, in the detector card, is given. \vspace{0.5cm}}
    234 \begin{tabular}[!h]{lllc}
     271\caption{Default values for the resolution of the central and forward calorimeters (for both electromagnetic and hadronic parts). Resolution is parametrised by the \textit{stochastic} ($S$), \textit{noise} ($N$) and \textit{constant} ($C$) terms (Eq.~\ref{eq:caloresolution})~\citep{qr:resolutionterms}.
     272%The corresponding parameter name, in the detector card, is given.
     273\vspace{0.5cm}}
     274\begin{tabular}[!h]{lccc}
    235275\hline
    236 \multicolumn{2}{c}{Resolution Term}   & Card flag & Value\\\hline
    237  \multicolumn{4}{l}{\textsc{ECAL}} \\
    238         & $S$ (GeV$^{1/2}$) & {\verb ELG_Scen }  & $0.05$ \\
    239         & $N$ (GeV)& {\verb ELG_Ncen }  & $0.25$ \\
    240         & $C$ & {\verb ELG_Ccen }  & $0.0055$ \\
    241  \multicolumn{4}{l}{\textsc{ECAL}, end caps} \\
    242         & $S$ (GeV$^{1/2}$) & {\verb ELG_Sec }  & $0.05$ \\
    243         & $N$ (GeV)& {\verb ELG_Nec }  & $0.25$ \\
    244         & $C$ & {\verb ELG_Cec }  & $0.0055$ \\
    245  \multicolumn{4}{l}{\textsc{FCAL}, electromagnetic part} \\
    246         & $S$ (GeV$^{1/2}$)& {\verb ELG_Sfwd }  & $2.084$ \\
    247         & $N$ (GeV)& {\verb ELG_Nfwd }  & $0$ \\
    248         & $C$ & {\verb ELG_Cfwd }  & $0.107$ \\
    249  \multicolumn{4}{l}{\textsc{HCAL}} \\
    250         & $S$ (GeV$^{1/2}$)& {\verb HAD_Scen } & $1.5$ \\
    251         & $N$ (GeV)& {\verb HAD_Ncen } & $0$\\
    252         & $C$ & {\verb HAD_Ccen } & $0.05$\\
    253  \multicolumn{4}{l}{\textsc{HCAL}, end caps} \\
    254         & $S$ (GeV$^{1/2}$)& {\verb HAD_Sec } & $1.5$ \\
    255         & $N$ (GeV)& {\verb HAD_Nec } & $0$\\
    256         & $C$ & {\verb HAD_Cec } & $0.05$\\
    257  \multicolumn{4}{l}{\textsc{FCAL}, hadronic part} \\
    258         & $S$ (GeV$^{1/2}$)& {\verb HAD_Sfwd }   & $2.7$\\
    259         & $N$ (GeV)& {\verb HAD_Nfwd }   & $0$ \\
    260         & $C$ & {\verb HAD_Cfwd }   & $0.13$\\
     276%\multicolumn{2}{c}{Resolution Term}   & Value\\\hline
     277        & $S$ (GeV$^{1/2}$) & $N$ (GeV) & $C$ \\\hline
     278 %\multicolumn{4}{l}{\textsc{ECAL}} \\
     279  ECAL            & $0.05$   & $0.25$ & $0.0055$ \\
     280 %\multicolumn{4}{l}{\textsc{ECAL}, end caps} \\
     281  ECAL, end caps  & $0.05$  & $0.25$ & $0.0055$ \\
     282 %\multicolumn{4}{l}{\textsc{FCAL}, electromagnetic part} \\
     283  FCAL, e.m. part & $2.084$ & $0$    & $0.107$ \\
     284 %\multicolumn{4}{l}{\textsc{HCAL}} \\
     285  HCAL            & $1.5$   & $0$    & $0.05$\\
     286 %\multicolumn{4}{l}{\textsc{HCAL}, end caps} \\
     287  HCAL, end caps  & $1.5$   & $0$    & $0.05$\\
     288 %\multicolumn{4}{l}{\textsc{FCAL}, hadronic part} \\
     289  FCAL, had. part & $2.7$   & $0$    & $0.13$\\
    261290\hline
    262291\end{tabular}
     
    265294\end{table}
    266295
    267 The 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}.
    268 Some long-living particles, such as the $K^0_s$ and $\Lambda$'s, with 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
     296% \begin{table}[!h]
     297% \begin{center}
     298% \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}).
     299% The corresponding parameter name, in the detector card, is given. \vspace{0.5cm}}
     300% \begin{tabular}[!h]{lllc}
     301% \hline
     302% \multicolumn{2}{c}{Resolution Term}   & Card flag & Value\\\hline
     303%  \multicolumn{4}{l}{\textsc{ECAL}} \\
     304%       & $S$ (GeV$^{1/2}$) & {\verb ELG_Scen }  & $0.05$ \\
     305%         & $N$ (GeV)& {\verb ELG_Ncen }  & $0.25$ \\
     306%         & $C$ & {\verb ELG_Ccen }  & $0.0055$ \\
     307%  \multicolumn{4}{l}{\textsc{ECAL}, end caps} \\
     308%         & $S$ (GeV$^{1/2}$) & {\verb ELG_Sec }  & $0.05$ \\
     309%         & $N$ (GeV)& {\verb ELG_Nec }  & $0.25$ \\
     310%         & $C$ & {\verb ELG_Cec }  & $0.0055$ \\
     311%  \multicolumn{4}{l}{\textsc{FCAL}, electromagnetic part} \\
     312%         & $S$ (GeV$^{1/2}$)& {\verb ELG_Sfwd }  & $2.084$ \\
     313%         & $N$ (GeV)& {\verb ELG_Nfwd }  & $0$ \\
     314%         & $C$ & {\verb ELG_Cfwd }  & $0.107$ \\
     315%  \multicolumn{4}{l}{\textsc{HCAL}} \\
     316%         & $S$ (GeV$^{1/2}$)& {\verb HAD_Scen } & $1.5$ \\
     317%         & $N$ (GeV)& {\verb HAD_Ncen } & $0$\\
     318%         & $C$ & {\verb HAD_Ccen } & $0.05$\\
     319%  \multicolumn{4}{l}{\textsc{HCAL}, end caps} \\
     320%         & $S$ (GeV$^{1/2}$)& {\verb HAD_Sec } & $1.5$ \\
     321%         & $N$ (GeV)& {\verb HAD_Nec } & $0$\\
     322%         & $C$ & {\verb HAD_Cec } & $0.05$\\
     323%  \multicolumn{4}{l}{\textsc{FCAL}, hadronic part} \\
     324%         & $S$ (GeV$^{1/2}$)& {\verb HAD_Sfwd }   & $2.7$\\
     325%         & $N$ (GeV)& {\verb HAD_Nfwd }   & $0$ \\
     326%         & $C$ & {\verb HAD_Cfwd }   & $0.13$\\
     327% \hline
     328% \end{tabular}
     329% \label{tab:defResol}
     330% \end{center}
     331% \end{table}
     332
     333
     334The energy of electrons and photons found in the particle list are smeared using only the \textsc{ECAL} resolution terms, while charged and neutral final-state hadrons interact with all calorimeters.
     335Some long-living particles, such as the $K^0_s$ and $\Lambda$'s, with lifetime $c\tau$ smaller than $10~\textrm{mm}$ are considered as stable particles by the generators 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
    269336\begin{equation}
    270337\left\{
     
    276343\end{equation}
    277344where $0 \leq F \leq 1$. The electromagnetic part is handled the same way for the electrons and photons.
    278 The resulting calorimetry 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\footnote{\texttt{[code]} To implement different ratios for other particles, see the \texttt{BlockClasses} class.}, the energy fraction is $F$ is assumed to be $0.7$.\\
    279 
    280 \subsection{Calorimetric towers}
    281 
    282 The 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 deposit energy in a calorimetric tower, either in \textsc{ECAL}, in \textsc{HCAL} or \textsc{FCAL}.
    283 As the detector is assumed to be cylindrical (e.g.\ symmetric in $\phi$ and with respect to the $\eta=0$ plane), the detector 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 calorimeter segmentation, which is common for the electromagnetic and hadronic sections at a given $(\eta,\phi)$.
     345The resulting calorimetry 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
     346%\footnote{\texttt{[code]} To implement different ratios for other particles, see the \texttt{BlockClasses} class.}
     347, the energy fraction is $F$ is assumed to be $0.7$~\citep{qr:emhadratios}.\\
     348
     349The smallest unit for geometrical sampling of the calorimeters is a \textit{cell}; it segments the $(\eta,\phi)$ plane for the energy measurement. No longitudinal segmentation is available in the simulated calorimeters. \textit{Delphes} assumes that ECAL and HCAL have the same segmentations and that the detector is symmetric in $\phi$ and with respect to the $\eta=0$ plane~\citep{qr:calorimetriccells}.
     350Fig.~\ref{fig:calosegmentation} illustrates the default calorimeter segmentation.
    284351
    285352\begin{figure}[!ht]
     
    292359\end{figure}
    293360
    294 The 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. Smearing is applied directly on the accumulated electromagnetic and hadronic energies of each calorimetric tower.
    295 
    296 \subsection{Very forward detector simulation}
    297 
    298 Most 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. In \textit{Delphes}, Zero Degree Calorimeters, roman pots and forward taggers have been implemented (Fig.~\ref{fig:fdets}).
    299 
    300 \begin{figure}[!ht]
    301 \begin{center}
    302 %\includegraphics[width=\columnwidth]{fdets}
    303 \includegraphics[width=\columnwidth]{fig4}
    304 \caption{Default location of the very forward detectors, including \textsc{ZDC}, \textsc{RP220} and \textsc{FP420} in the \textsc{LHC} beamline.
    305 Incoming (beam 1, red) and outgoing (beam 2, black) beams on one side of the fifth interaction point (\textsc{IP5}, $s=0~\textrm{m}$ on the plot).
    306 The Zero Degree Calorimeter is located in perfect alignment with the beamline axis at the interaction point, at $140~\textrm{m}$, the beam paths are separated. The forward taggers are near-beam detectors located at $220~\textrm{m}$ and $420~\textrm{m}$. Beamline simulation with \textit{Hector}~\citep{bib:hector}. All very forward detectors are located symmetrically around the interaction point. }
    307 \label{fig:fdets}
    308 \end{center}
    309 \end{figure}
    310 
    311 \begin{table*}[t]
    312 \begin{center}
    313 \caption{Default parameters for the forward detectors: distance from the interaction point and detector acceptance. The \textsc{LHC} beamline is assumed around the fifth \textsc{LHC} interaction point (\textsc{IP}). For the \textsc{ZDC}, the acceptance depends only on the pseudorapidity $\eta$ of the particle, which should be neutral and stable.
    314 The tagger acceptance is fully determined by the distance in the transverse plane of the detector to the real beam position~\citep{bib:hector}. It is expressed in terms of the particle energy ($E$).
    315 All detectors are located on both sides of the interaction point.
    316 \vspace{0.5cm}}
    317 \begin{tabular}{llcl}
    318 \hline
    319 Detector & Distance from \textsc{IP}& Acceptance & \\ \hline
    320 \textsc{ZDC}   & $\pm 140$ m & $|\eta|> 8.3$       & for $n$ and $\gamma$\\
    321 \textsc{RP220} & $\pm 220$ m & $E \in [6100 ; 6880]$ (GeV) & at $2~\textrm{mm}$\\
    322 \textsc{FP420} & $\pm 420$ m & $E \in [6880 ; 6980]$ (GeV) & at $4~\textrm{mm}$\\
    323 \hline
    324 \end{tabular}
    325 \label{tab:fdetacceptance}
    326 \end{center}
    327 \end{table*}
    328 
    329 
    330 \subsubsection*{Zero Degree Calorimeters}
    331 
    332 In direct sight of the interaction point, on both sides of the central detector, the Zero Degree Calorimeters (\textsc{ZDC}s) are located at zero angle, i.e.\ are aligned with the beamline axis at the interaction point. They are placed beyond the point where the paths of incoming and outgoing beams separate. 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{ATLAS} and \textsc{CMS}).
    333 
    334 The trajectory of the neutrals observed in the \textsc{ZDC}s is a straight line, while charged particles are deflected away from their acceptance window by the powerful magnets located in front of them. The fact that additional charged particles may enter the \textsc{ZDC} acceptance is neglected here.
    335 
    336 The \textsc{ZDC}s have the ability to measure the time-of-flight of the particle.
    337 This corresponds to the delay after which the particle is observed in the detector, with respect to the bunch crossing reference time at the interaction point ($t_0$). The measured time-of-flight $t$ is simply given by:
    338 \begin{equation}
    339  t = t_0 + \frac{1}{v} \times \Big( \frac{s-z}{\cos \theta}\Big),
    340 \end{equation}
    341 where $t_0$ is thus 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, $\theta$ is the particle emission angle. It is then assumed 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}$.
    342 The formula then reduces to
    343 \begin{equation}
    344  t = \frac{1}{c} \times (s-z).
    345 \end{equation}
    346 For example, a photon takes $0.47~\mu\textrm{s}$ to reach a \textsc{ZDC} located at $s=140~\textrm{m}$, neglecting $z$ and $\theta$. For the time-of-flight measurement, a Gaussian smearing can be applied according to the detector resolution (Tab.~\ref{tab:defResolZdc}). In the current version of \textit{Delphes}, only neutrons, antineutrons and photons are assumed to be able to reach the \textsc{ZDC}s, all other particles being neglected.
    347 
    348 The \textsc{ZDC}s are composed of an electromagnetic and a hadronic sections, for the measurement of photons and neutrons, respectively. The energy of the observed neutral is smeared according to Eq.~\ref{eq:caloresolution} and the corresponding section resolutions (Tab.~\ref{tab:defResolZdc}). The \textsc{ZDC} hits do not enter in the calorimeter tower list used for reconstruction of jets and missing transverse energy.
    349 
    350 \begin{table}[!h]
    351 \begin{center}
    352 \caption{Default values for the resolution of the zero degree calorimeters. Resolution on energy measurement is parametrised by the \textit{stochastic} ($S$), \textit{noise} ($N$) and \textit{constant} ($C$) terms (Eq.~\ref{eq:caloresolution}). The time-of-flight is smeared according to a Gaussian function.
    353 The corresponding parameter name, in the detector card, is given. \vspace{0.5cm}}
    354 \begin{tabular}[!h]{lllc}
    355 \hline
    356 \multicolumn{2}{c}{Resolution Term}   & Card flag & Value\\\hline
    357  \multicolumn{4}{l}{\textsc{ZDC}, electromagnetic part} \\
    358         & $S$ (GeV$^{1/2}$)& \texttt{ELG\_Szdc}  & $0.7$ \\
    359         & $N$ (GeV)& \texttt{ELG\_Nzdc}  & $0.0$ \\
    360         & $C$ & \texttt{ELG\_Czdc}  & $0.08$ \\
    361  \multicolumn{4}{l}{\textsc{ZDC}, hadronic part} \\
    362         & $S$ (GeV$^{1/2}$)& \texttt{HAD\_Szdc}   & $1.38$\\
    363         & $N$ (GeV)& \texttt{HAD\_Nzdc}   & $0$ \\
    364         & $C$ & \texttt{HAD\_Czdc}   & $0.13$\\
    365  \multicolumn{4}{l}{\textsc{ZDC}, timing resolution} \\
    366         & $\sigma_t$ (s) & \texttt{ZDC\_T\_resolution} & $0$ \\
    367 \hline
    368 \end{tabular}
    369 \label{tab:defResolZdc}
    370 \end{center}
    371 \end{table}
    372 
    373 \subsubsection*{Forward taggers}
    374 
    375 Forward taggers (called here \textsc{RP220}, for ``roman pots at $220~\textrm{m}$'' and \textsc{FP420} ``for forward proton taggers at $420~\textrm{m}$'', as at the \textsc{LHC}) are meant for the measurement of particles following very closely the beam path. Such devices, also used at \textsc{HERA} and Tevatron, are located very far away from the interaction point (further than $150$~m in the \textsc{LHC} case).
    376 
    377 To be able to reach these detectors, particles must have a charge identical to the beam particles, and a momentum very close to the nominal value of 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}).
    378 For instance, roman pots at $220~\textrm{m}$ from the  \textsc{IP} and $2~\textrm{mm}$ from the beam will detect all forward protons with an energy between $120$ and $900~\textrm{GeV}$~\citep{bib:hector}.
    379 In practice, in the \textsc{LHC}, only positively charged muons ($\mu^+$) and protons can reach the forward taggers as other particles with a single positive charge coming from the interaction points will decay before their possible tagging. In \textit{Delphes}, extra hits coming from the beam-gas events or secondary particles hitting the beampipe in front of the detectors are not taken into account.
    380 
    381 While 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 \textit{Hector} software~\citep{bib:hector}, which includes the chromaticity effects and the geometrical aperture of the beamline elements of any arbitrary collider.
    382 
    383 Forward taggers are able to measure the hit positions ($x,y$) and angles ($\theta_x,\theta_y$) in the transverse plane at the location of the detector ($s$ meters away from the \textsc{IP}), as well as the time-of-flight\footnote{It should be noted that for both \textsc{CMS} and \textsc{ATLAS} experiments, the taggers located at $220$~m are not able to measure the time-of-flight, contrary to \textsc{FP420} detectors.} ($t$). Out of these the particle energy ($E$) and the momentum transfer it underwent during the interaction ($q^2$) can be reconstructed\footnote{The reconstruction of $E$ and $q^2$ are not implemeted in \textit{Delphes} but can be performed at the analysis level.}. The time-of-flight measurement can be smeared with a Gaussian distribution (default value\footnote{\texttt{[code] } The resolution is defined by the \texttt{RP220\_T\_resolution} and \texttt{RP420\_T\_resolution} parameters in the detector card.} $\sigma_t = 0~\textrm{s}$).
     361No sharing between neighbouring cells is implemented when particles enter a cell very close to its geometrical edge. Due to the finite segmentation, the smearing, as defined in Eq.~\ref{eq:caloresolution}, is applied directly on the accumulated electromagnetic and hadronic energies of each calorimetric cell. The calorimetric cells directly enter in the calculation of the missing transverse energy (\textsc{MET}), and as input for the jet reconstruction algorithms.
     362
    384363
    385364
     
    387366\section{High-level object reconstruction}
    388367
    389 Analysis 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 \textit{Delphes}.
    390 In addition, some detector data are added: tracks, calorimetric towers and hits in \textsc{ZDC}, \textsc{RP220} and \textsc{FP420}.
    391 While electrons, muons and photons are easily identified, some other objects are more difficult to measure, like jets or missing energy due to invisible particles.
     368Analysis 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
     369%\footnote{\texttt{[code] }All these processed data are located under the \texttt{Analysis} tree.}
     370in the output file created by \textit{Delphes}~\citep{qr:analysistree}.
     371In addition, some detector data are added: tracks, calorimetric cells and hits in the very forward detectors (\textsc{ZDC}, \textsc{RP220} and \textsc{FP420}, Sec.~\ref{sec:vfd}). While electrons, muons and photons are easily identified, some other objects are more difficult to measure, like jets or missing energy due to invisible particles.
    392372
    393373For most of these objects, their four-momentum and related quantities are directly accessible in \textit{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).
    394  
    395 
    396374
    397375\subsection{Photon and charged lepton reconstruction}
    398 From 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.
     376From 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. The collections of electrons, photons and muons are filled in with candidates observing some fiducial and reconstruction cuts, and are based on the true particle ID provided by the generator. Consequently, no fake candidates enter these collections. However, when needed, fake candidates can be added into the collections at the analysis level, when processing \textit{Delphes} output data. As effects like bremsstrahlung are not taken into account along the lepton propagation in the tracker, no clustering is needed for the electron reconstruction in \textit{Delphes}.
    399377
    400378\subsubsection*{Electrons and photons}
    401 Electron ($e^\pm$) and photon 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}/c$). A calorimetric tower will be seen in the detector, as electrons will leave in addition a track. Subsequently, electrons and photons create a candidate in the jet collection.
    402 Assuming a good measurement of the track parameters in the real experiment, the electron energy can be reasonably recovered. In \textit{Delphes}, electron energy is smeared according to the resolution of the calorimetric tower where it points to, but independently from any other deposited energy is this tower. This approach is still conservative as the calorimeter resolution is worse than the tracker one.
     379Real electron ($e^\pm$) and photon candidates are identified if they fall into the acceptance of the tracking system and have a transverse momentum above a threshold (default $p_T > 10~\textrm{GeV}/c$). A calorimetric cell will be activated in the detector and electrons will leave in addition a track. Subsequently, electrons and photons create a candidate in the jet collection.
     380Assuming a good measurement of the track parameters in the real experiment, the electron energy can be reasonably recovered. In \textit{Delphes}, electron energy is smeared according to the resolution of the calorimetric cell where it points to, but independently from any other deposited energy is this cell. This approach is still conservative as the calorimeter resolution is worse than the tracker one.
    403381
    404382\subsubsection*{Muons}
    405383Generator-level muons entering the detector acceptance are considered as candidates for the analysis level.
    406384The 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}/c$) and of the pseudorapidity coverage of the muon system (default: $-2.4 \leq \eta \leq 2.4$).
    407 The 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. Multiple scattering is neglected. This implies that low energy muons have in \textit{Delphes} a better resolution than in a real detector. Furthermore, muons leave no deposit in calorimeters. At last, the particles which might leak out of the calorimeters into the muon systems (\textit{punch-through}) will not be see
    408 n as muon candidates in \textit{Delphes}.
     385The application of the detector resolution on the muon momentum depends on a Gaussian smearing of the $p_T$ variable~\citep{qr:muonsmearing}.
     386%\footnote{\texttt{[code]} See the \texttt{SmearMuon} method.}.
     387Neither $\eta$ nor $\phi$ variables are modified beyond the calorimeters: no additional magnetic field is applied. Multiple scattering is neglected. This implies that low energy muons have in \textit{Delphes} a better resolution than in a real detector. Furthermore, muons leave no deposit in calorimeters. At last, the particles which might leak out of the calorimeters into the muon systems (\textit{punch-through}) will not be seen as muon candidates in \textit{Delphes}.
    409388
    410389\subsubsection*{Charged lepton isolation}
     
    414393The result (i.e.\ \textit{isolated} or \textit{not}) is added to the charged lepton measured properties.
    415394In addition, the sum $P_T$ of the transverse momenta of all tracks but the lepton one within the isolation cone is
    416 provided\footnote{\texttt{[code] }See the \texttt{IsolFlag} and \texttt{IsolPt} values in the \texttt{Electron} or \texttt{Muon} collections in the \texttt{Analysis} tree, as well as the \texttt{ISOL\_PT} and \texttt{ISOL\_Cone} variables in the detector card.}:
     395provided~\citep{qr:isolflag}:
     396%\footnote{\texttt{[code] }See the \texttt{IsolFlag} and \texttt{IsolPt} values in the \texttt{Electron} or \texttt{Muon} collections in the \texttt{Analysis} tree, as well as the \texttt{ISOL\_PT} and \texttt{ISOL\_Cone} variables in the detector card.}
    417397$$ P_T = \sum_{i \neq \mu}^\textrm{tracks} p_T(i)$$
    418398
    419 No calorimetric isolation is applied, but the muon collection contains also the ratio $\rho_\mu$ between (1) the sum of the transverse energies in all calotowers in a $N \times N$ grid around the muon, and (2) the muon transverse
    420 momentum\footnote{\texttt{[code] }Calorimetric isolation parameters in the detector card are \texttt{ISOL\_Calo\_ET} and  \texttt{ISOL\_Calo\_Grid}.}:
     399No calorimetric isolation is applied, but the muon collection contains also the ratio $\rho_\mu$ between (1) the sum of the transverse energies in all calorimetric cells in a $N \times N$ grid around the muon, and (2) the muon transverse momentum~\citep{qr:caloisolation}:
     400%\footnote{\texttt{[code] }Calorimetric isolation parameters in the detector card are \texttt{ISOL\_Calo\_ET} and  \texttt{ISOL\_Calo\_Grid}.}:
    421401$$ \rho_\mu = \frac{\Sigma_i E_T(i)}{p_T(\mu)}~,~ i\textrm{ in }N \times N \textrm { grid centred on }\mu.$$
    422402
    423403
    424 \subsubsection*{Forward neutrals}
    425 
    426 The zero degree calorimeter hits correspond to neutral particles with a lifetime long enough to reach these detectors (default: $c \tau \geq 140~\textrm{m}$) and very large pseudorapidities (default: $|\eta|>8.3$). In current versions of \textit{Delphes}, only photons and neutrons are considered. Photons are identified thanks to the electromagnetic section of the calorimeter, and if their energy overpasses a given threshold (def. $20$~GeV). Similarly, neutrons are reconstructed according to the resolution of the hadronic section, if their energy exceeds a threshold\footnote{\texttt{[code]} These thresholds are defined by the \texttt{ZDC\_gamma\_E} and \texttt{ZDC\_n\_E} variables in the detector card.} (def. $50$~GeV).
     404% \subsubsection*{Forward neutrals}
     405%
     406% The zero degree calorimeter hits correspond to neutral particles with a lifetime long enough to reach these detectors (default: $c \tau \geq 140~\textrm{m}$) and very large pseudorapidities (default: $|\eta|>8.3$). In current versions of \textit{Delphes}, only photons and neutrons are considered. Photons are identified thanks to the electromagnetic section of the calorimeter, and if their energy overpasses a given threshold (def. $20$~GeV). Similarly, neutrons are reconstructed according to the resolution of the hadronic section, if their energy exceeds a threshold (def. $50$~GeV)~\citep{qr:fwdneutrals}.
     407% %\footnote{\texttt{[code]} These thresholds are defined by the \texttt{ZDC\_gamma\_E} and \texttt{ZDC\_n\_E} variables in the detector card.} (def. $50$~GeV).
    427408
    428409
     
    430411\subsection{Jet reconstruction}
    431412
    432 A realistic analysis requires a correct treatment of particles which have hadronised. Therefore, the most widely currently used jet algorithms have been integrated into the \textit{Delphes} framework using the FastJet tools~\citep{bib:FASTJET}.
    433 Six different jet reconstruction schemes are available\footnote{\texttt{[code] }The choice is done by allocating the \texttt{JET\_jetalgo } input parameter in the detector 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 differ in their sensitivity to soft particles or collinear splittings, and in their computing speed performances.
     413A realistic analysis requires a correct treatment of particles which have hadronised. Therefore, the most widely currently used jet algorithms have been integrated into the \textit{Delphes} framework using the FastJet tools\footnote{A more detailed description of the jet algorithms is given in the User Manual, in appendix.}.
     414Six different jet reconstruction schemes are available, with three cone algorithms and three recombination algorithms~\citep{bib:FASTJET,qr:jetalgo}.
     415%\footnote{\texttt{[code] }The choice is done by allocating the \texttt{JET\_jetalgo } input parameter in the detector card.}.
     416% The first three belong to the cone algorithm class while the last three are using a sequential recombination scheme.
     417For all of them, the calorimetric cells are used as inputs for the jet clustering. Jet algorithms differ in their sensitivity to soft particles or collinear splittings, and in their computing speed performances.
    434418By default, reconstruction uses a cone algorithm with $\Delta R=0.7$.
    435 Jets are stored if their transverse energy is higher\footnote{\texttt{[code] PTCUT\_jet }variable in the detector card.} than $20~\textrm{GeV}$.
     419Jets are stored if their transverse energy is higher
     420%\footnote{\texttt{[code] PTCUT\_jet }variable in the detector card.}
     421than $20~\textrm{GeV}$~\citep{qr:ptcutjet}.
    436422 
    437423\subsubsection*{Cone algorithms}
     
    439425\begin{enumerate}
    440426 
    441 \item {\it CDF Jet Clusters}~\citep{bib:jetclu}: Algorithm forming jets by associating together towers lying within a circle (default radius $\Delta R=0.7$) in the $(\eta$, $\phi)$ space.
    442 This so-called JetCLU cone jet algorithm is used by the \textsc{CDF} experiment in Run II.
    443 All 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.
    444 The existing FastJet code has been modified to allow easy modification of the tower pattern in $(\eta, \phi)$ space.
    445 In following versions of \textit{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 detector card.}.
    446  
    447 \item {\it CDF MidPoint}~\citep{bib:midpoint}: Algorithm developed for the \textsc{CDF} Run II to reduce infrared and collinear sensitivities compared to purely seed-based cone by adding `midpoints' (energy barycentres) in the list of cone seeds.
     427\item {\it CDF Jet Clusters}~\citep{bib:jetclu}: Cone algorithm forming jets by combining cells lying within a circle (default radius $\Delta R=0.7$) in the $(\eta$, $\phi)$ space. Jets are seeded by all cells with
     428 transverse energy $E_T$ overpassing a given threshold (default: $E_T > 1~\textrm{GeV}$)~\citep{qr:jetparams}.
     429 
     430\item {\it CDF MidPoint}~\citep{bib:midpoint}: Cone algorithm with additional ``midpoints'' (energy barycentres) in the list of seeds; this algorithm has reduced infrared and collinear sensitivities.
    448431 
    449432\item {\it Seedless Infrared Safe Cone}~\citep{bib:SIScone}: The \textsc{SISC}one algorithm is simultaneously insensitive to additional soft particles and collinear splittings, and fast enough to be used in experimental analysis.
    450  
    451433\end{enumerate}
    452  
     434
    453435\subsubsection*{Recombination algorithms}
    454  
    455 The three sequential recombination jet algorithms are safe with respect to 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$.
    456 
    457 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 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:
     436
     437The next three jet algorithms rely on recombination schemes where calorimeter cell pairs are successively merged (\textit{E-scheme recombination}):
     438
     439% Two such variables are defined: the distance $d_{ij}$ between each pair of cells $(i,j)$, and a variable $d_{iB}$ (\textit{beam distance}) depending on the transverse momentum of the cell $i$.
     440
     441% The jet reconstruction algorithm browses the calorimetric cell 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 cells $i$ and $j$ are merged into a single cell 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 cell is declared as a final jet and is removed from the input list. This procedure is repeated until no cells 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 calorimetric cell $i$ and $\Delta R_{ij}= \sqrt{(y_i-y_j)^2+(\phi_i-\phi_j)^2}$ as the jet-radius parameter:
    458442 
    459443\begin{enumerate}[start=4]
    460444 
    461 \item {\it Longitudinally invariant $k_t$ jet}~\citep{bib:ktjet}:
    462 \begin{equation}
    463 \begin{array}{l}
    464   d_{ij} = \min(k_{ti}^2,k_{tj}^2)\Delta R_{ij}^2/R^2 \\
    465   d_{iB}=k_{ti}^2 \\
    466 \end{array}
    467 \end{equation}
    468  
    469 \item {\it Cambridge/Aachen jet}~\citep{bib:aachen}:
    470 \begin{equation}
    471 \begin{array}{l}
    472 d_{ij} = \Delta R_{ij}^2/R^2\\
    473 d_{iB}=1 \\
    474 \end{array}
    475 \end{equation}
    476  
    477 \item {\it Anti $k_t$ jet}~\citep{bib:antikt}: where hard jets are exactly circular in the $(y,\phi)$ plane
    478 \begin{equation}
    479 \begin{array}{l}
    480 d_{ij} =  \min(1/k_{ti}^2,1/k_{tj}^2)\Delta R_{ij}^2/R^2 \\
    481 d_{iB}=1/k_{ti}^2 \\
    482 \end{array}
    483 \end{equation}
     445\item {\it Longitudinally invariant $k_t$ jet}~\citep{bib:ktjet},
     446% \begin{equation}
     447% \begin{array}{l}
     448%   d_{ij} = \min(k_{ti}^2,k_{tj}^2)\Delta R_{ij}^2/R^2 \\
     449%   d_{iB}=k_{ti}^2 \\
     450% \end{array}
     451% \end{equation}
     452 
     453\item {\it Cambridge/Aachen jet}~\citep{bib:aachen},
     454% \begin{equation}
     455% \begin{array}{l}
     456% d_{ij} = \Delta R_{ij}^2/R^2\\
     457% d_{iB}=1 \\
     458% \end{array}
     459% \end{equation}
     460 
     461\item {\it Anti $k_t$ jet}~\citep{bib:antikt}, where hard jets are exactly circular in the $(y,\phi)$ plane.
     462% \begin{equation}
     463% \begin{array}{l}
     464% d_{ij} =  \min(1/k_{ti}^2,1/k_{tj}^2)\Delta R_{ij}^2/R^2 \\
     465% d_{iB}=1/k_{ti}^2 \\
     466% \end{array}
     467% \end{equation}
    484468\end{enumerate}
     469
     470The recombination algorithms are safe with respect to soft radiations (\textit{infrared}) and collinear splittings. Their implementations are similar except for the definition of the \textit{distances} used during the merging procedure.
     471
    485472 
    486473\subsubsection*{Energy flow}
    487474
    488 In jets, several particle can leave their energy into a given calorimetric tower, which broadens the jet energy resolution. However, the energy of charged particles associated to jets can be deduced from their reconstructed track, thus providing a way to identify some of the components of towers with multiple hits. When the \textit{energy flow} is switched on in \textit{Delphes}\footnote{\texttt{[code]} Set \texttt{JET\_Eflow} to $1$ or $0$ in the detector card in order to switch on or off the energy flow for jet reconstruction.}, the energy of tracks pointing to calotowers is extracted and smeared separately, before running the chosen jet reconstruction algorithm. This option allows a better jet $E$ reconstruction.
     475In jets, several particle can leave their energy into a given calorimetric cell, which broadens the jet energy resolution. However, the energy of charged particles associated to jets can be deduced from their reconstructed track, thus providing a way to identify some of the components of cells with multiple hits. When the \textit{energy flow} is switched on in \textit{Delphes}
     476%\footnote{\texttt{[code]} Set \texttt{JET\_Eflow} to $1$ or $0$ in the detector card in order to switch on or off the energy flow for jet reconstruction.}
     477, the energy of tracks pointing to calorimetric cells is extracted and smeared separately, before running the chosen jet reconstruction algorithm. This option allows a better jet $E$ reconstruction~\citep{qr:energyflow}.
    489478 
    490479\subsection{$b$-tagging}
    491480\label{btagging}
    492481
    493 A 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. By default, 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{BTAG\_b}, \texttt{BTAG\_mistag\_c} and \texttt{BTAG\_mistag\_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.}.
     482A 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. By default, 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~\citep{qr:btag}.
     483%\footnote{\texttt{[code] }Corresponding to the \texttt{BTAG\_b}, \texttt{BTAG\_mistag\_c} and \texttt{BTAG\_mistag\_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.}.
    494484The (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$). In current version of \textit{Delphes}, the displacement of secondary vertices is not simulated.
    495485
     
    504494\caption{ Branching ratios for $\tau^-$ lepton~\citep{bib:pdg}. $h^\pm$ and $h^0$ refer to charged and neutral hadrons, respectively. $n \geq 0$ and $m \geq 0$ are integers.
    505495\vspace{0.5cm}  }
    506 \begin{tabular}[!h]{ll}
     496\begin{tabular}[!h]{lll}
    507497\hline
    508  \multicolumn{2}{l}{\textbf{Leptonic decays}}\\
    509  $ \tau^- \rightarrow e^- \ \bar \nu_e \ \nu_\tau$ & $17.9\% $ \\
    510  $ \tau^- \rightarrow \mu^- \ \bar \nu_\mu  \ \nu_\tau$ & $17.4\%$ \\
    511  \multicolumn{2}{l}{\textbf{Hadronic decays}}\\
    512  $ \tau^- \rightarrow h^-\ (n\times h^\pm) \ (m\times h^0) \  \nu_\tau$  & $64.7\%$ \\
    513  $ \tau^- \rightarrow h^-\ (m\times h^0) \ \nu_\tau$  & $50.1\%$ \\
    514  $ \tau^- \rightarrow h^-\ h^+ h^-  (m\times h^0) \ \nu_\tau$  & $14.6\%$ \\
     498 \multicolumn{3}{l}{\textbf{Leptonic decays}}\\
     499 & $ \tau^- \rightarrow e^- \ \bar \nu_e \ \nu_\tau$ & $17.9\% $ \\
     500 & $ \tau^- \rightarrow \mu^- \ \bar \nu_\mu  \ \nu_\tau$ & $17.4\%$ \\
     501 \multicolumn{3}{l}{\textbf{Hadronic decays}}\\
     502 & $ \tau^- \rightarrow h^-\ (n\times h^\pm) \ (m\times h^0) \  \nu_\tau$  & $64.7\%$ \\
     503 & $ \tau^- \rightarrow h^-\ (m\times h^0) \ \nu_\tau$  & $50.1\%$ \\
     504 & $ \tau^- \rightarrow h^-\ h^+ h^-  (m\times h^0) \ \nu_\tau$  & $14.6\%$ \\
    515505\hline
    516506\end{tabular}
     
    531521\begin{table}[!h]
    532522\begin{center}
    533 \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.
     523\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 calorimetric cells $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~\citep{qr:taujets}.
    534524\vspace{0.5cm}  }
     525% \begin{tabular}[!h]{lll}
     526% \hline
     527% Parameter  & Card flag & Value\\\hline
     528% \multicolumn{3}{l}{\textbf{Electromagnetic collimation}} \\
     529% $R^\textrm{em}$     & \texttt{TAU\_energy\_scone } & $0.15$\\
     530% min $E_{T}^\textrm{tower}$     & {\verb JET_M_seed }  & $1.0$~GeV\\
     531% $C_{\tau}$         & \texttt{TAU\_energy\_frac} & $0.95$\\
     532% \multicolumn{3}{l}{\textbf{Tracking isolation}} \\
     533% $R^\textrm{tracks}$ & \texttt{TAU\_track\_scone} & $0.4$\\
     534% min $p_T^\textrm{tracks}$      & \texttt{PTAU\_track\_pt } & $2$ GeV$/c$\\
     535% \multicolumn{3}{l}{\textbf{$\tau$-jet candidate}} \\
     536% $\min p_T$ & \texttt{TAUJET\_pt} & $10$ GeV$/c$\\
     537% \hline
     538% \end{tabular}
    535539\begin{tabular}[!h]{lll}
    536540\hline
    537 Parameter  & Card flag & Value\\\hline
    538541\multicolumn{3}{l}{\textbf{Electromagnetic collimation}} \\
    539 $R^\textrm{em}$     & \texttt{TAU\_energy\_scone } & $0.15$\\
    540 min $E_{T}^\textrm{tower}$     & {\verb JET_M_seed } & $1.0$~GeV\\
    541 $C_{\tau}$         & \texttt{TAU\_energy\_frac} & $0.95$\\
     542& $R^\textrm{em}$          & $0.15$\\
     543& min $E_{T}^\textrm{tower}$ & $1.0$~GeV\\
     544& $C_{\tau}$              & $0.95$\\
    542545\multicolumn{3}{l}{\textbf{Tracking isolation}} \\
    543 $R^\textrm{tracks}$ & \texttt{TAU\_track\_scone} & $0.4$\\
    544 min $p_T^\textrm{tracks}$      & \texttt{PTAU\_track\_pt } & $2$ GeV$/c$\\
     546& $R^\textrm{tracks}$      & $0.4$\\
     547& min $p_T^\textrm{tracks}$ & $2$ GeV$/c$\\
    545548\multicolumn{3}{l}{\textbf{$\tau$-jet candidate}} \\
    546 $\min p_T$ & \texttt{TAUJET\_pt} & $10$ GeV$/c$\\
     549& $\min p_T$              & $10$ GeV$/c$\\
    547550\hline
    548551\end{tabular}
     
    554557\subsubsection*{Electromagnetic collimation}
    555558
    556 To use the narrowness of the $\tau$-jet, the \textit{electromagnetic collimation} $C_{\tau}$ 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.
    557 To be taken into account, a calorimeter tower should have a transverse energy $E_T^\textrm{tower}$ above a given threshold.
     559To use the narrowness of the $\tau$-jet, the \textit{electromagnetic collimation} $C_{\tau}$ is defined as the sum of the energy of cells in a small cone of radius $R^\textrm{em}$ around the jet axis, divided by the energy of the reconstructed jet.
     560To be taken into account, a calorimeter cell should have a transverse energy $E_T^\textrm{tower}$ above a given threshold.
    558561A large fraction of the jet energy is expected in this small cone. This fraction, or \textit{collimation factor}, is represented in Fig.~\ref{fig:tau2} for the default values (see Tab.~\ref{tab:tauRef}).
    559562
     
    609612\end{equation}
    610613The \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.
    611 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\footnote{However, as tracks and calorimetric towers are available in the output file, the missing transverse energy can always be reprocessed a posteriori. }:
     614In a real experiment, calorimeters measure energy and not momentum. Any problem affecting the detector (dead channels, misalignment, noisy cells, cracks) worsens directly the measured missing transverse energy $\overrightarrow {E_T}^\textrm{miss}$. In this document, \textsc{MET} is based on the calorimetric cells and only muons and neutrinos are not taken into account for its evaluation:
    612615\begin{equation}
    613616\overrightarrow{E_T}^\textrm{miss} = - \sum^\textrm{towers}_i \overrightarrow{E_T}(i)
    614617\end{equation}
    615 
     618However, as muon candidates, tracks and calorimetric cells are available in the output file, the missing transverse energy can always be reprocessed a posteriori with more specialised algorithms.
    616619
    617620\section{Trigger emulation}
     
    621624%High statistics are required for data analyses, consequently imposing high luminosity, i.e.\ a high collision rate.
    622625As 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.
    623 This data selection is supposed to reject only well-known \textsc{SM} events\footnote{However, some bandwidth is allocated to minimum-bias and/or zero-bias (``random'') triggers that stores a small fraction of the events without any selection criteria.}.
     626This data selection is supposed to reject only well-known \textsc{SM} events\footnote{In real experiments, some bandwidth is allocated to minimum-bias and/or zero-bias (``random'') triggers that stores a small fraction of  random events without any selection criteria.}.
    624627Dedicated 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.
    625628
    626629Most 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}/c$ and $p_T^e > 50~\textrm{GeV}/c$.
    627630
    628 A trigger emulation is included in \textit{Delphes}, using a fully parametrisable \textit{trigger table}\footnote{\texttt{[code] }The trigger card is the \texttt{data/TriggerCard.dat} file.}. When enabled, this trigger is applied on analysis-object data.
     631A trigger emulation is included in \textit{Delphes}, using a fully parametrisable \textit{trigger table} \citep{qr:triggercard}
     632%\footnote{\texttt{[code] }The trigger card is the \texttt{data/TriggerCard.dat} file.}
     633. When enabled, this trigger is applied on analysis-object data.
    629634In a real experiment, the online selection is often divided into several steps (or \textit{levels}).
    630635This splits the overall reduction factor into a product of smaller factors, corresponding to the different trigger levels.
     
    635640Real triggers are thus intrinsically based on reconstructed data with a worse resolution than final analysis data.
    636641On the contrary, same data are used in \textit{Delphes} for trigger emulation and for final analyses.
     642
     643\section{\label{sec:vfd}Very forward detector simulation}
     644
     645Most 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. In \textit{Delphes}, Zero Degree Calorimeters, roman pots and forward taggers have been implemented (Fig.~\ref{fig:fdets}), similarly to the plans for CMS and ATLAS collaborations~\citep{bib:cmsjetresolution, bib:ATLASresolution}.
     646
     647\begin{figure}[!ht]
     648\begin{center}
     649%\includegraphics[width=\columnwidth]{fdets}
     650\includegraphics[width=\columnwidth]{fig4}
     651\caption{Default location of the very forward detectors, including \textsc{ZDC}, \textsc{RP220} and \textsc{FP420} in the \textsc{LHC} beamline.
     652Incoming (beam 1, red) and outgoing (beam 2, black) beams on one side of the fifth interaction point (\textsc{IP5}, $s=0~\textrm{m}$ on the plot).
     653The Zero Degree Calorimeter is located in perfect alignment with the beamline axis at the interaction point, at $140~\textrm{m}$, the beam paths are separated. The forward taggers are near-beam detectors located at $220~\textrm{m}$ and $420~\textrm{m}$. Beamline simulation with \textit{Hector}~\citep{bib:hector}. All very forward detectors are located symmetrically around the interaction point. }
     654\label{fig:fdets}
     655\end{center}
     656\end{figure}
     657
     658%\begin{table*}[t]  % the star (*) allows to arrange the table over the two columns
     659\begin{table}[t]
     660\begin{center}
     661\caption{Default parameters for the forward detectors: distance from the interaction point and detector acceptance. The \textsc{LHC} beamline is assumed around the fifth \textsc{LHC} interaction point (\textsc{IP}). For the \textsc{ZDC}, the acceptance depends only on the pseudorapidity $\eta$ of the particle, which should be neutral and stable.
     662The tagger acceptance is fully determined by the distance in the transverse plane of the detector to the real beam position~\citep{bib:hector}. It is expressed in terms of the particle energy ($E$).
     663All detectors are located on both sides of the interaction point.
     664\vspace{0.5cm}}
     665\begin{tabular}{llcl}
     666\hline
     667%Detector & Distance from \textsc{IP}& Acceptance & \\ \hline
     668Detector & Distance & Acceptance & \\ \hline
     669\textsc{ZDC}   & $\pm 140$ m & $|\eta|> 8.3$       & for $n$ and $\gamma$\\
     670\textsc{RP220} & $\pm 220$ m & $E \in [6100 ; 6880]$ (GeV) & at $2~\textrm{mm}$\\
     671\textsc{FP420} & $\pm 420$ m & $E \in [6880 ; 6980]$ (GeV) & at $4~\textrm{mm}$\\
     672\hline
     673\end{tabular}
     674\label{tab:fdetacceptance}
     675\end{center}
     676\end{table}
     677
     678
     679\subsection{Zero Degree Calorimeters}
     680
     681In direct sight of the interaction point, on both sides of the central detector, the Zero Degree Calorimeters (\textsc{ZDC}s) are located at zero angle, i.e.\ are aligned with the beamline axis at the interaction point. They are placed beyond the point where the paths of incoming and outgoing beams separate. 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{ATLAS} and \textsc{CMS}).
     682
     683The trajectory of the neutrals observed in the \textsc{ZDC}s is a straight line, while charged particles are deflected away from their acceptance window by the powerful magnets located in front of them. The fact that additional charged particles may enter the \textsc{ZDC} acceptance is neglected in the current versions of \textit{Delphes}.
     684
     685The \textsc{ZDC}s have the ability to measure the time-of-flight of the particle.
     686This corresponds to the delay $t$ after which the particle is observed in the detector, with respect to the bunch crossing reference time at the interaction point ($t_0$):
     687\begin{equation}
     688 t = t_0 + \frac{1}{v} \times \Big( \frac{s-z}{\cos \theta}\Big) \approx \frac{1}{c} \times (s-z),
     689\end{equation}
     690where $t_0$ is thus 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, $\theta$ is the particle emission angle. It is assumed that the neutral particle observed in the \textsc{ZDC} is highly relativistic and very forward.
     691% 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}$.
     692% The formula then reduces to
     693% \begin{equation}
     694%  t = \frac{1}{c} \times (s-z).
     695% \end{equation}
     696% For example, a photon takes $0.47~\mu\textrm{s}$ to reach a \textsc{ZDC} located at $s=140~\textrm{m}$, neglecting $z$ and $\theta$.
     697For the time-of-flight measurement, a Gaussian smearing can be applied according to the detector resolution (Tab.~\ref{tab:defResolZdc})~\citep{qr:resolutionterms}.
     698%In the current version of \textit{Delphes}, only neutrons, antineutrons and photons are assumed to be able to reach the \textsc{ZDC}s, all other particles being neglected.
     699
     700The \textsc{ZDC}s are composed of an electromagnetic and a hadronic sections, for the measurement of photons and neutrons, respectively. The energy of the observed neutral is smeared according to Eq.~\ref{eq:caloresolution} and the corresponding section resolutions (Tab.~\ref{tab:defResolZdc}). The \textsc{ZDC} hits do not enter in the calorimeter cell list used for reconstruction of jets and missing transverse energy.
     701
     702\begin{table}[!h]
     703\begin{center}
     704\caption{Default values for the resolution of the zero degree calorimeters. Resolution on energy measurement is parametrised by the \textit{stochastic} ($S$), \textit{noise} ($N$) and \textit{constant} ($C$) terms (Eq.~\ref{eq:caloresolution})~\citep{qr:resolutionterms}. The time-of-flight is smeared according to a Gaussian function.
     705\vspace{0.5cm}}
     706% \begin{tabular}[!h]{lllc}
     707% \hline
     708% \multicolumn{2}{c}{Resolution Term}   & Card flag & Value\\\hline
     709%  \multicolumn{4}{l}{\textsc{ZDC}, electromagnetic part} \\
     710%         & $S$ (GeV$^{1/2}$)& \texttt{ELG\_Szdc}  & $0.7$ \\
     711%         & $N$ (GeV)& \texttt{ELG\_Nzdc}  & $0.0$ \\
     712%         & $C$ & \texttt{ELG\_Czdc}  & $0.08$ \\
     713%  \multicolumn{4}{l}{\textsc{ZDC}, hadronic part} \\
     714%         & $S$ (GeV$^{1/2}$)& \texttt{HAD\_Szdc}   & $1.38$\\
     715%         & $N$ (GeV)& \texttt{HAD\_Nzdc}   & $0$ \\
     716%         & $C$ & \texttt{HAD\_Czdc}   & $0.13$\\
     717%  \multicolumn{4}{l}{\textsc{ZDC}, timing resolution} \\
     718%         & $\sigma_t$ (s) & \texttt{ZDC\_T\_resolution} & $0$ \\
     719% \hline
     720% \end{tabular}
     721\begin{tabular}[!h]{llcc}
     722\hline
     723 \multicolumn{3}{l}{\textsc{ZDC}, electromagnetic part} & hadronic part \\
     724        & $S$ (GeV$^{1/2}$)     & $0.7$ & $1.38$\\
     725        & $N$ (GeV)             & $0$   & $0$   \\
     726        & $C$                   & $0.08$& $0.13$ \\
     727  \multicolumn{4}{l}{\textsc{ZDC}, timing resolution} \\
     728        & $\sigma_t$ (s) & $0$ & \\
     729\hline
     730\end{tabular}
     731\label{tab:defResolZdc}
     732\end{center}
     733\end{table}
     734
     735% \subsubsection*{Forward neutrals}
     736
     737The reconstructed ZDC hits correspond to neutral particles with a lifetime long enough to reach these detectors (default: $c \tau \geq 140~\textrm{m}$) and very large pseudorapidities (default: $|\eta|>8.3$).
     738%In current versions of \textit{Delphes}, only photons and neutrons are considered.
     739Photons are identified thanks to the electromagnetic section of the calorimeter, and if their energy overpasses a given threshold (def. $20$~GeV). Similarly, neutrons are reconstructed according to the resolution of the hadronic section, if their energy exceeds a threshold (def. $50$~GeV)~\citep{qr:fwdneutrals}.
     740%\footnote{\texttt{[code]} These thresholds are defined by the \texttt{ZDC\_gamma\_E} and \texttt{ZDC\_n\_E} variables in the detector card.} (def. $50$~GeV).
     741
     742
     743\subsection{Forward taggers}
     744
     745Forward taggers (called here \textsc{RP220}, for ``roman pots at $220~\textrm{m}$'' and \textsc{FP420} for ``forward proton taggers at $420~\textrm{m}$'', as at the \textsc{LHC}) are meant for the measurement of particles following very closely the beam path. Such devices, also used at \textsc{HERA} and Tevatron, are located very far away from the interaction point (further than $150$~m in the \textsc{LHC} case).
     746
     747To be able to reach these detectors, particles must have a charge identical to the beam particles, and a momentum very close to the nominal value of 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}).
     748For instance, roman pots at $220~\textrm{m}$ from the  \textsc{IP} and $2~\textrm{mm}$ from the beam will detect all forward protons with an energy between $120$ and $900~\textrm{GeV}$~\citep{bib:hector}.
     749In practice, in the \textsc{LHC}, only positively charged muons ($\mu^+$) and protons can reach the forward taggers as other particles with a single positive charge coming from the interaction points will decay before their possible tagging. In \textit{Delphes}, extra hits coming from the beam-gas events or secondary particles hitting the beampipe in front of the detectors are not taken into account.
     750
     751While 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 \textit{Hector} software~\citep{bib:hector}, which includes the chromaticity effects and the geometrical aperture of the beamline elements of any arbitrary collider.
     752
     753Forward taggers are able to measure the hit positions ($x,y$) and angles ($\theta_x,\theta_y$) in the transverse plane at the location of the detector ($s$ meters away from the \textsc{IP}), as well as the time-of-flight\footnote{It is worth noting that for both \textsc{CMS} and \textsc{ATLAS} experiments, the taggers located at $220$~m are not able to measure the time-of-flight, contrary to \textsc{FP420} detectors.} ($t$). Out of these the particle energy ($E$) and the momentum transfer it underwent during the interaction ($q^2$) can be reconstructed at the analysis level (it is not implemented in the current versions of \textit{Delphes}. The time-of-flight measurement can be smeared with a Gaussian distribution (default value
     754%\footnote{\texttt{[code] } The resolution is defined by the \texttt{RP220\_T\_resolution} and \texttt{RP420\_T\_resolution} parameters in the detector card.}
     755$\sigma_t = 0~\textrm{s}$)~\citep{qr:protontaggers}.
     756
     757
    637758
    638759\section{Validation}
     
    710831The samples used to study the \textsc{MET} performance are identical to those used for the jet validation.
    711832It is worth noting that the contribution to $E_T^\textrm{miss}$ from muons is negligible in the studied sample.
    712 The 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.
     833The 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 cells.
    713834The quality of the \textsc{MET} reconstruction is checked via the resolution on its horizontal component $E_x^\textrm{miss}$.
    714835
     
    723844\includegraphics[width=\columnwidth]{fig10}
    724845\includegraphics[width=\columnwidth]{fig10b}
    725 \caption{$\sigma(E^\textrm{mis}_{x})$ as a function on the scalar sum of all towers ($\Sigma E_T$) for $pp \rightarrow gg$ events, for a \textsc{CMS}-like detector (top) and an \textsc{ATLAS}-like detector (bottom), for di-jet events produced with MadGraph/MadEvent and hadronised with \textit{Pythia}.}
     846\caption{$\sigma(E^\textrm{mis}_{x})$ as a function on the scalar sum of all cells ($\Sigma E_T$) for $pp \rightarrow gg$ events, for a \textsc{CMS}-like detector (top) and an \textsc{ATLAS}-like detector (bottom), for di-jet events produced with MadGraph/MadEvent and hadronised with \textit{Pythia}.}
    726847\label{fig:resolETmis}
    727848\end{center}
     
    734855where the $\alpha$ parameter depends on the resolution of the calorimeters.
    735856
    736 The \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 high-luminosity in the same bunch crossing.}~\citep{bib:cmsjetresolution}, which compares very well with the $\alpha = 0.63$ obtained with \textit{Delphes}. Similarly, for an \textsc{ATLAS}-like detector, a value of $0.53$ is obtained by \textit{Delphes} for the $\alpha$ parameter, while the experiment expects it in the range $[0.53~ ;~0.57]$~\citep{bib:ATLASresolution}.
     857The \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 (i.e. extra simultaneous $pp$ collision occurring at high-luminosity in the same bunch crossing)~\citep{bib:cmsjetresolution}, which compares very well with the $\alpha = 0.63$ obtained with \textit{Delphes}. Similarly, for an \textsc{ATLAS}-like detector, a value of $0.53$ is obtained by \textit{Delphes} for the $\alpha$ parameter, while the experiment expects it in the range $[0.53~ ;~0.57]$~\citep{bib:ATLASresolution}.
    737858
    738859\subsection{\texorpdfstring{$\tau$}{\texttau}-jet efficiency}
     
    770891\section{Visualisation}
    771892
    772 When 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}~\citep{bib:FROG} has been interfaced in \textit{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$.}.
     893When 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}~\citep{bib:FROG} has been interfaced in \textit{Delphes}, allowing an easy display of the defined detector configuration~\citep{qr:frog}.
     894%\footnote{\texttt{[code] } To prepare the visualisation, the \texttt{FLAG\_FROG} parameter should be equal to $1$.}.
    773895 
    774896% \begin{figure}[!ht]
     
    867989Moreover, the framework allows trigger emulation and 3D event visualisation.
    868990
    869 \textit{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.
    870 
    871 This framework has already been used for several analyses, in particular in photon-induced interactions at the \textsc{LHC}~\citep{bib:wtphotoproduction, bib:papierquisortirajamais, bib:papiersimon}.
     991\textit{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, a better $b$-tag description and possibly the implementation of an event mixing module for pile-up event simulation. This framework has already been used for several analyses~\citep{bib:wtphotoproduction, bib:papierquisortirajamais, bib:papiersimon}, in particular in photon-induced interactions at the \textsc{LHC}.
    872992
    873993
     
    8801000\begin{thebibliography}{99}
    8811001\addcontentsline{toc}{section}{References}
     1002
     1003\bibitem{bib:geant} J. Allison, et al., Nucl. Inst. \& Meth. in \textbf{Phys. Res. A} \href{http://dx.doi.org/10.1016/S0168-9002(03)01368-8}{506 (2003) 250-303}; \textbf{IEEE Trans. on Nucl. Sc.} \href{http://dx.doi.org/10.1109/TNS.2006.869826}{53:1 (2006) 270-278}.
    8821004 
    8831005\bibitem{bib:delphes} \textit{Delphes}, \href{http://www.fynu.ucl.ac.be/delphes.html}{www.fynu.ucl.ac.be/delphes.html}
    884 %hepforge:
    8851006\bibitem{bib:stdhep} L.A. Garren, M. Fischler, \href{http://cepa.fnal.gov/psm/stdhep/c++}{cepa.fnal.gov/psm/stdhep/c++}
    8861007\bibitem{bib:hepmc} M. Dobbs and J.B. Hansen, \textbf{Comput. Phys. Commun.} \href{http://dx.doi.org/10.1016/S0010-4655(00)00189-2}{134 (2001) 41}.
     
    9171038\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].
    9181039
    919 \bibitem{bib:papierquisortirajamais}J. de Favereau~et~al, \textbf{CP3-08-04} (2008), to be published in EPJ.
     1040\bibitem{bib:papierquisortirajamais}J. de Favereau~et~al, \href{http://arxiv.org/abs/0908.2020}{arXiv:0908.2020v1} [hep-ph] (2008), to be published in EPJ.
    9201041
    9211042%\bibitem{bib:papiersimon} ``Phenomenology of a twisted two-Higgs-doublet model'', Simon de Visscher, Jean-Marc Gerard, Michel Herquet, Vincent Lema\^itre, Fabio Maltoni, to be published.
     
    9231044
    9241045\bibitem{bib:mcfio} P. Lebrun, L. Garren, Copyright (c) 1994-1995 Universities Research Association, Inc.
    925 
    926 
    9271046\end{thebibliography}
     1047
     1048
     1049
     1050% references to code
     1051\renewcommand\refname{Internal code references}
     1052\begin{thebibliography}{2}
     1053\addcontentsline{toc}{section}{Internal code references}
     1054
     1055\bibitem[a]{qr:inputformat} See the following classes: \texttt{HEPEVTConverter}, \texttt{HepMCConverter}, \texttt{LHEFConverter}, \texttt{STDHEPConverter} and \texttt{DelphesRootConverter}.
     1056
     1057\bibitem[b]{qr:invisibleparticles} The list of particles considered as invisible is accessible in the \texttt{PdgParticle} class. This list currently contains the PIDs 12, 14, 16, 1000022,  1000023, 1000025, 1000035 and 1000045, in absolute values.
     1058
     1059\bibitem[c]{qr:lhco} Set the \texttt{FLAG\_LHCO} variable to $1$ or $0$ in the detector card to switch on/off the creation of \texttt{*.lhco} output file.
     1060
     1061\bibitem[d]{qr:detectorcard}The detector card is the \texttt{data/DetectorCard.dat} file. This file is parsed by the \texttt{SmearUtil} class.
     1062
     1063\bibitem[e]{qr:datacards} Detector and trigger cards for the \textsc{ATLAS} and \textsc{CMS} experiments are also provided in \texttt{data/} directory.
     1064
     1065\bibitem[f]{qr:resolutionterms}The resolution terms in the detector card are named \texttt{ELG\_Xyyy} or \texttt{HAD\_Xyyy},  refering to electromagnetic and hadronic terms (resp.); \texttt{X} is replaced by \texttt{S}, \texttt{N}, \texttt{C} for the stochastic, noise and constant terms; and finally \texttt{yyy} is \texttt{cen} for central part, \texttt{ec} for end-caps, \texttt{fwd} for the forward calorimeters and \texttt{zdc} for the zero-degree calorimeters.
     1066
     1067\bibitem[g]{qr:magneticfield} See the \texttt{TrackPropagation} class.
     1068
     1069\bibitem[h]{qr:tracks} See the \texttt{TRACK\_eff} and \texttt{TRACK\_ptmin} terms in the detector card.
     1070
     1071\bibitem[i]{qr:energysmearing} The response of the detector is applied to the electromagnetic and the hadronic particles through the \texttt{SmearElectron} and \texttt{SmearHadron} methods in the \texttt{SmearUtil} class.
     1072
     1073\bibitem[j]{qr:emhadratios} To implement different ratios for other particles, see the \texttt{BlockClasses} class.
     1074
     1075\bibitem[k]{qr:calorimetriccells} As the detector is assumed to be cylindrical (e.g.\ symmetric in $\phi$ and with respect to the $\eta=0$ plane), the detector card stores the number of calorimetric cells with $\phi=0$ and $\eta>0$ (default: $40$ cells). For a given $\eta$, the size of the $\phi$ segmentation is also specified. See the \texttt{TOWER\_number}, \texttt{TOWER\_eta\_edges} and \texttt{TOWER\_dphi} variables in the detector card.
     1076
     1077\bibitem[l]{qr:analysistree} All these processed data are located under the \texttt{Analysis} tree.
     1078
     1079\bibitem[m]{qr:muonsmearing} See the \texttt{SmearMuon} method in the \texttt{SmearUtil} class.
     1080
     1081\bibitem[n]{qr:isolflag} See the \texttt{IsolFlag} and \texttt{IsolPt} values in the \texttt{Electron} or \texttt{Muon} collections in the \texttt{Analysis} tree, as well as the \texttt{ISOL\_PT} and \texttt{ISOL\_Cone} variables in the detector card.
     1082
     1083\bibitem[o]{qr:caloisolation} Calorimetric isolation parameters in the detector card are \texttt{ISOL\_Calo\_ET} and  \texttt{ISOL\_Calo\_Grid} in the detector card.
     1084
     1085\bibitem[p]{qr:fwdneutrals} These thresholds are defined by the \texttt{ZDC\_gamma\_E} and \texttt{ZDC\_n\_E} variables in the detector card.
     1086
     1087\bibitem[q]{qr:jetalgo} The choice is done by allocating the \texttt{JET\_jetalgo } input parameter in the detector card.
     1088
     1089\bibitem[r]{qr:ptcutjet} See the \texttt{PTCUT\_jet }variable in the detector card.
     1090
     1091\bibitem[s]{qr:jetparams} See the \texttt{JET\_coneradius} and \texttt{JET\_seed} variables in the detector card. The existing FastJet code has been modified to allow easy modification of the cell pattern in $(\eta, \phi)$ space.
     1092In following versions of \textit{Delphes}, a new dedicated plug-in will be created on this purpose.
     1093
     1094\bibitem[t]{qr:energyflow} Set \texttt{JET\_Eflow} to $1$ or $0$ in the detector card in order to switch on or off the energy flow for jet reconstruction.
     1095
     1096\bibitem[u]{qr:btag} Corresponding to the \texttt{BTAG\_b}, \texttt{BTAG\_mistag\_c} and \texttt{BTAG\_mistag\_l} constants, for the efficiency of tagging of a $b$-jet, the efficiency of mistagging a $c$-jet as a $b$-jet, and the
     1097efficiency of mistagging a light jet ($u$,$d$,$s$,$g$) as a $b$-jet.
     1098
     1099\bibitem[v]{qr:taujets} See the following parameters in the detector card:\\
     1100\texttt{TAU\_energy\_scone } for $R^\textrm{em}$; \texttt{JET\_M\_seed }  for min $E_{T}^\textrm{tower}$;
     1101\texttt{TAU\_energy\_frac} for $C_{\tau}$; \texttt{TAU\_track\_scone} for $R^\textrm{tracks}$;
     1102 \texttt{PTAU\_track\_pt } for min $p_T^\textrm{tracks}$ and \texttt{TAUJET\_pt} for $\min p_T$.
     1103
     1104
     1105\bibitem[w]{qr:triggercard} The trigger card is the \texttt{data/TriggerCard.dat} file. Default trigger files are also available for CMS-like and ATLAS-like detectors
     1106
     1107\bibitem[x]{qr:protontaggers} The resolution is defined by the \texttt{RP220\_T\_resolution} and \texttt{RP420\_T\_resolution} parameters in the detector card.
     1108
     1109\bibitem[y]{qr:frog} To prepare the visualisation, the \texttt{FLAG\_FROG} parameter should be equal to $1$.
     1110
     1111\end{thebibliography}
     1112
     1113
     1114
    9281115
    9291116\onecolumn
    9301117\appendix
    931  
     1118
    9321119\section{User manual}
    9331120 
     
    9381125 
    9391126In order to run \textit{Delphes} on your system, first download its sources and compile them:\\
    940 \texttt{wget http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes\_V\_*.tar.gz}\\
    941 Replace the \texttt{*} symbol by the proper version number\footnote{Refer to the download page on the \textit{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}. Current version of Delphes for this manual is V 1.8 (July 2009)}.
     1127\begin{quote}\texttt{wget http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes\_V\_*.tar.gz}\end{quote}
     1128Replace the \texttt{*} symbol by the proper version number. Always refer to the download page on the \textit{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}. Current version of Delphes for this manual is V 1.8 (July 2009).
    9421129
    9431130\begin{quote}
     
    9791166 \end{itemize}
    9801167 
    981 If no datacard is provided by the user, the default smearing and running parameters are used:
    982 \begin{quote}
    983 \begin{verbatim}
    984 # Detector extension, in pseudorapidity units (|eta|)
     1168If no datacard is provided by the user, the default smearing and running parameters are used (corresponding to tables~\ref{tab:defEta},~\ref{tab:defResol}).\\
     1169Definition of the sub-detector extensions:
     1170\begin{quote}
     1171\begin{verbatim}
    9851172CEN_max_tracker    2.5     // Maximum tracker coverage
    9861173CEN_max_calo_cen   1.7     // central calorimeter coverage
     
    9881175CEN_max_calo_fwd   5.0     // forward calorimeter pseudorapidity coverage
    9891176CEN_max_mu         2.4     // muon chambers pseudorapidity coverage
    990  
     1177\end{verbatim}
     1178\end{quote}
     1179Definition of the sub-detector resolutions:
     1180\begin{quote}
     1181\begin{verbatim}
    9911182# Energy resolution for electron/photon in central/endcap/fwd/zdc calos
    9921183# \sigma/E = C + N/E + S/\sqrt{E}, E in GeV
     
    10321223\end{verbatim}
    10331224\end{quote}
    1034 
     1225Definitions related to the calorimetric cells:
    10351226\begin{quote}
    10361227\begin{verbatim}
    10371228# Calorimetric towers
    10381229TOWER_number         40
    1039 ### list of the edges of each tower in eta for eta>0 assuming
    1040 ###a symmetric detector in eta<0
    1041 ### the list starts with the lower edge of the most central tower
    1042 ### the list ends with the higher edged of the most forward tower
    1043 ### there should be NTOWER+1 values
    10441230TOWER_eta_edges 0. 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783
    10451231               0.870 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653
     
    10481234               5.000
    10491235 
    1050 ### list of the tower size in phi (in degrees), assuming that all
    1051 ### towers are similar in phi for a given eta value
    1052 ### the list starts with the phi-size of the most central tower (eta=0)
    1053 ### the list ends with the phi-size of the most forward tower
    1054 ### there should be NTOWER values
    10551236TOWER_dphi 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 10
    10561237           10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 20 20
    10571238\end{verbatim}
    10581239\end{quote}
    1059  
     1240\texttt{TOWER\_eta\_edges} is the list of the edges in $\eta$ of all cells, in the $\eta>0$ hemisphere (the detector is supposed to be symmetric with respect to the $\eta=0$ plane, as well as around the $z$-axis). Starts with the lower edge of the most central tower (default: $\eta = 0$) and ends with the higher edge of the most forward tower.
     1241\texttt{TOWER\_dphi} lists the tower size in $\phi$ (in degree), assuming that all cells are similar in $\phi$ for a given $\eta$.\\
     1242Thresholds applied for storing the reconstructed objects in the final collections:
    10601243\begin{quote}
    10611244\begin{verbatim}
     
    10701253ZDC_gamma_E      20
    10711254ZDC_n_E          50
    1072 
     1255\end{verbatim}
     1256\end{quote}
     1257Definitions of variables related to the charged lepton isolation:
     1258\begin{quote}
     1259\begin{verbatim}
    10731260# Charged lepton isolation. Pt and Et in GeV
    10741261ISOL_PT          2.0  //minimal pt of tracks for isolation criteria
     
    10771264ISOL_Calo_ET     2.0  //minimal tower E_T for isolation criteria. 1E99 means "off"
    10781265ISOL_Calo_Grid   3    //Grid size (N x N) for calorimetric isolation
    1079 
     1266\end{verbatim}
     1267\end{quote}
     1268Definitions of variables related to the jet reconstruction:
     1269\begin{quote}
     1270\begin{verbatim}
    10801271# General jet variable
    10811272JET_coneradius   0.7  // generic jet radius
     
    10891280JET_Eflow        1    // Energy flow: perfect energy assumed in the tracker coverage.
    10901281                      // 1 is 'on' ; 0 is 'off'
    1091 \end{verbatim}
    1092 \end{quote}
    1093 
    1094 \begin{quote}
    1095 \begin{verbatim}
     1282
    10961283# Tagging definition
    10971284BTAG_b           40    // b-tag efficiency (%)
    10981285BTAG_mistag_c    10    // mistagging (%)
    10991286BTAG_mistag_l    1     // mistagging (%)
    1100  
     1287\end{verbatim}
     1288\end{quote}
     1289Switches for options
     1290\begin{quote}
     1291\begin{verbatim}
    11011292# FLAGS
    11021293FLAG_bfield      1     //1 to run the bfield propagation else 0
     
    11061297FLAG_FROG        1     //1 to run the FROG event display
    11071298FLAG_LHCO        1     //1 to run the LHCO
    1108 
     1299\end{verbatim}
     1300\end{quote}
     1301Parameters for the magnetic field simulation:
     1302\begin{quote}
     1303\begin{verbatim}
    11091304# In case BField propagation allowed
    11101305TRACK_radius      129   // radius of the BField coverage, in cm
     
    11131308TRACK_bfield_y    0     // Y component of the BField, in T
    11141309TRACK_bfield_z    3.8   // Z component of the BField, in T
    1115 
     1310\end{verbatim}
     1311\end{quote}
     1312Parameters related to the very forward detectors
     1313\begin{quote}
     1314\begin{verbatim}
    11161315# Very forward detector extension, in pseudorapidity
    11171316# if allowed
    11181317VFD_min_zdc       8.3   // Zero-Degree neutral Calorimeter
    11191318VFD_s_zdc         140   // distance of the ZDC, from the IP, in [m]
    1120 \end{verbatim}
    1121 \end{quote}
    1122 
    1123 
    1124 \begin{quote}
    1125 \begin{verbatim}
     1319
    11261320#\textit{Hector} parameters
    11271321RP_220_s          220     // distance of the RP to the IP, in meters
     
    11391333RP_cross_ang_x    142.5   // half-crossing angle in horizontal plane, in microrad
    11401334RP_cross_ang_y    0       // half-crossing angle in vertical plane, in microrad
    1141 
    1142 
     1335\end{verbatim}
     1336\end{quote}
     1337Others parameters:
     1338\begin{quote}
     1339\begin{verbatim}
    11431340# In case FROG event display allowed
    11441341NEvents_FROG      100
    11451342# Number of events to process
    11461343NEvents           -1                    // -1 means 'all'
    1147 
    11481344
    11491345# input PDG tables
     
    12291425\textbf{Analysis \texttt{Tree}} & & \\
    12301426~~~Tracks     & Collection of tracks                       & {\verb TRootTracks }\\
    1231 ~~~CaloTower  & Calorimetric towers                        & {\verb TRootCalo }\\
     1427~~~CaloTower  & Calorimetric cells                         & {\verb TRootCalo }\\
    12321428~~~Electron   & Collection of electrons                    & {\verb TRootElectron }\\
    12331429~~~Photon     & Collection of photons                      & {\verb TRootPhoton }\\
     
    12841480   \texttt{~~~float PhiCalo } &\texttt{ // particle azimuthal angle in rad when entering the calo }\\
    12851481   \texttt{~~~float EHoverEE }&\texttt{ // hadronic energy over electromagnetic energy }\\
    1286    \texttt{~~~float EtRatio } &\texttt{ // calo Et in NxN-tower grid around the muon over the muon Et }\\
     1482   \texttt{~~~float EtRatio } &\texttt{ // calo Et in NxN-cell grid around the muon over the muon Et }\\
    12871483\end{tabular}
    12881484\end{quote}
     
    13131509\begin{tabular}{ll}
    13141510\multicolumn{2}{l}{\textbf{Leaves in the \texttt{CaloTower} branch (\texttt{Analysis} tree)}}\\
    1315     \texttt{~~~float Eta }     &\texttt{ // pseudorapidity of the tower }\\
    1316     \texttt{~~~float Phi }     &\texttt{ // azimuthal angle of the tower in rad }\\
    1317     \texttt{~~~float E }       &\texttt{ // tower energy in GeV }\\
    1318     \texttt{~~~float E\_em }   &\texttt{ // electromagnetic component of the tower energy in GeV}\\
    1319     \texttt{~~~float E\_had }  &\texttt{ // hadronic component of the tower energy in GeV}\\
    1320     \texttt{~~~float ET }      &\texttt{ // tower transverse energy in GeV }\\
     1511    \texttt{~~~float Eta }     &\texttt{ // pseudorapidity of the cell }\\
     1512    \texttt{~~~float Phi }     &\texttt{ // azimuthal angle of the cell in rad }\\
     1513    \texttt{~~~float E }       &\texttt{ // cell energy in GeV }\\
     1514    \texttt{~~~float E\_em }   &\texttt{ // electromagnetic component of the cell energy in GeV}\\
     1515    \texttt{~~~float E\_had }  &\texttt{ // hadronic component of the cell energy in GeV}\\
     1516    \texttt{~~~float ET }      &\texttt{ // cell transverse energy in GeV }\\
    13211517& \\
    13221518\multicolumn{2}{l}{\textbf{Leaves in the \texttt{ETmis} branch (\texttt{Analysis} tree)}}\\
     
    13671563The hit position is computed from the center of the beam position, not from the edge of the detector.
    13681564
    1369  
     1565\subsection{Deeper description of jet algorithms}
     1566
     1567In this section, we briefly describe the differences between the six jet algorithms interfaced in \textit{Delphes}, via the FastJet utiliy~\citep{bib:FASTJET}. Jet algorithms differ in their sensitivity to soft particles or collinear splittings, and in their computing speed performances. The first three belong to the cone algorithm class while the last three are using a sequential recombination scheme. For all of them, the calorimetric cells are used as inputs for the jet clustering.
     1568 
     1569\subsubsection*{Cone algorithms}
     1570 
     1571\begin{enumerate}
     1572 
     1573\item {\it CDF Jet Clusters}~\citep{bib:jetclu}: Basic cone reconstruction algorithm used by the \textsc{CDF} experiment in Run II). All cells lying in a circular cone around the jet axis with a transverse energy $E_T$ higher than a given threshold are used to seed the jet candidates. This algorithm is fast but sensitive to both soft particles and collinear splittings.
     1574
     1575\item {\it CDF MidPoint}~\citep{bib:midpoint}: Cone reconstruction algorithm developed for the \textsc{CDF} Run II to reduce infrared and collinear sensitivities compared to purely seed-based cone by adding `midpoints' (energy barycentres) in the list of cone seeds.
     1576 
     1577\item {\it Seedless Infrared Safe Cone}~\citep{bib:SIScone}: The \textsc{SISC}one algorithm is simultaneously insensitive to additional soft particles and collinear splittings, and fast enough to be used in experimental analysis.
     1578 
     1579\end{enumerate}
     1580
     1581
     1582\subsubsection*{Recombination algorithms}
     1583 
     1584The three sequential recombination jet algorithms are safe with respect to soft radiations (\textit{infrared}) and collinear splittings. They rely on recombination schemes where calorimeter cell pairs are successively merged.
     1585The 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 cells $(i,j)$, and a variable $d_{iB}$ (\textit{beam distance}) depending on the transverse momentum of the cell $i$.
     1586The jet reconstruction algorithm browses the calorimetric cell 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 cells $i$ and $j$ are merged into a single cell 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 cell is declared as a final jet and is removed from the input list. This procedure is repeated until no cells 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 calorimetric cell $i$ and $\Delta R_{ij}= \sqrt{(y_i-y_j)^2+(\phi_i-\phi_j)^2}$ as the jet-radius parameter:
     1587 
     1588\begin{enumerate}[start=4]
     1589 
     1590\item {\it Longitudinally invariant $k_t$ jet}~\citep{bib:ktjet}, with
     1591  $d_{ij} = \min(k_{ti}^2,k_{tj}^2) \times \frac{\Delta R_{ij}^2}{R^2}$ and $d_{iB}=k_{ti}^2$,
     1592\item {\it Cambridge/Aachen jet}~\citep{bib:aachen}, with $d_{ij} = \frac{\Delta R_{ij}^2}{R^2}$ and $d_{iB}=1$,
     1593\item {\it Anti $k_t$ jet}~\citep{bib:antikt}, where hard jets are exactly circular in the $(y,\phi)$ plane:
     1594$d_{ij} =  \min(1/k_{ti}^2,1/k_{tj}^2) \times \frac{\Delta R_{ij}^2}{R^2}$ and $d_{iB}=\frac{1}{k_{ti}^2}$.
     1595\end{enumerate}
     1596
     1597
    13701598\subsection{Running an analysis on your \textit{Delphes} events}
    13711599 
     
    15191747
    15201748\paragraph{9th column (\texttt{had/em})}
    1521 For jets, electrons and photons, the ninth column is the ration between hadronic and electromagnetic energies in the calorimetric towers associated to the object. This is always \texttt{0} for missing transverse energy.
     1749For jets, electrons and photons, the ninth column is the ration between hadronic and electromagnetic energies in the calorimetric cells associated to the object. This is always \texttt{0} for missing transverse energy.
    15221750For muons, this number (\texttt{aaa.bb}) reports two values related to the muon isolation (section \ref{sec:isolation}). The integer part (\texttt{aaa}) is transverse momentum sum $P_T$ (in GeV/$c$) and the fractional part (\texttt{bb}) is the energy ratio $\rho_\mu$.
    15231751
Note: See TracChangeset for help on using the changeset viewer.