| 1 |
fisk |
1.1 |
\section {Analysis Demonstrations}
|
| 2 |
|
|
|
| 3 |
acosta |
1.15 |
A wide variety of physics analysis demonstrations were prepared by the
|
| 4 |
|
|
physics groups for CSA06 in order to test the analysis workflow for
|
| 5 |
|
|
CSA06. These demonstrations also proved to be useful training
|
| 6 |
|
|
exercises for collaborators in the new software and computing tools.
|
| 7 |
|
|
The list of specialized calibration and
|
| 8 |
|
|
alignment streams produced at the Tier-0 and analyzed is discussed in
|
| 9 |
|
|
Section~\ref{sec:offlineswalca}. The list of general physics analysis
|
| 10 |
|
|
skims produced at the Tier-1 centres is covered in
|
| 11 |
|
|
Table~\ref{tab:tier1skim}. Table~\ref{tab:analyses} lists the physics
|
| 12 |
|
|
analyses that were conducted as part of CSA06. Details on these
|
| 13 |
|
|
analyses are contained in the following subsections. Unless noted,
|
| 14 |
|
|
CRAB was used to submit the analysis jobs.
|
| 15 |
|
|
|
| 16 |
|
|
|
| 17 |
|
|
\begin{table}[phtb]
|
| 18 |
|
|
\centering
|
| 19 |
|
|
\caption{List of CSA06 analyses.}
|
| 20 |
|
|
\label{tab:analyses}
|
| 21 |
|
|
\begin{tabular}{|l|l|l|l|l|}
|
| 22 |
|
|
\hline
|
| 23 |
|
|
Group & Analysis & People & Dataset & Comments \\ \hline
|
| 24 |
|
|
hg & Selection of Z$\rightarrow 2\tau$ &
|
| 25 |
|
|
K.Petridis (IC), A.Kalinowski (Warsaw) & EWK skims & \\
|
| 26 |
|
|
hg & jet$\rightarrow \tau$ mis id &
|
| 27 |
|
|
F.Blekman (IC), C.Siamitros (Brunel) & Z+jet& \\
|
| 28 |
|
|
hg & tau tagging efficiency from Z$\rightarrow 2\tau$ &
|
| 29 |
|
|
S.Gennai and G.Bagliesi (Pisa), & EWK skims & \\
|
| 30 |
|
|
& & A.Goussiou and R.Vasques Sierra (FNAL) & & \\
|
| 31 |
|
|
hg & tau HLT with pixel trigger using Z$\rightarrow 2\tau$ &
|
| 32 |
|
|
D.Kotlinski and P.Trueb (PSI) & EWK skims & \\
|
| 33 |
|
|
hg & H$\rightarrow WW\rightarrow \ell\ell$ &
|
| 34 |
|
|
F. Stoeckli (ETH) & EWK & CERN batch\\
|
| 35 |
|
|
hg & Z+2Jet background to qqH, H$\rightarrow$inv &
|
| 36 |
|
|
J.Brooke and S.Metson (Bristol), & EWK skims & \\
|
| 37 |
|
|
& & K. Mazumdar, S.Bansal (TIFR) & & \\
|
| 38 |
|
|
\hline
|
| 39 |
|
|
\end{tabular}
|
| 40 |
|
|
\end{table}
|
| 41 |
|
|
|
| 42 |
|
|
|
| 43 |
ndefilip |
1.3 |
|
| 44 |
fisk |
1.1 |
\subsection {Calibration}
|
| 45 |
|
|
|
| 46 |
acosta |
1.14 |
\label{sec:calib}
|
| 47 |
|
|
|
| 48 |
malgeri |
1.6 |
\input{calib_ana.tex}
|
| 49 |
|
|
|
| 50 |
fpschill |
1.5 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
| 51 |
|
|
|
| 52 |
|
|
\subsection {Alignment}
|
| 53 |
|
|
|
| 54 |
acosta |
1.14 |
\label{sec:align}
|
| 55 |
|
|
|
| 56 |
fpschill |
1.5 |
A tracker alignment CSA06 exercise was carried out with the goal to
|
| 57 |
|
|
demonstrate the full work- and dataflow of the alignment process.
|
| 58 |
|
|
The exercise followed closely the ideas and concepts
|
| 59 |
|
|
developed during the T0-RTAG~\cite{tier0rtag}.
|
| 60 |
|
|
The exercise comprised the following steps:
|
| 61 |
|
|
\begin{itemize}
|
| 62 |
|
|
\item Reading of alignment constants from the offline database during prompt
|
| 63 |
|
|
reconstruction;
|
| 64 |
|
|
\item Writing dedicated AlCaReco streams for alignment;
|
| 65 |
|
|
\item Defining a misalignment scenario and insertion of the corresponding
|
| 66 |
|
|
object into the offline database;
|
| 67 |
|
|
\item Running an alignment algorithm at the Tier-0;
|
| 68 |
|
|
\item Inserting the resulting alignment corrections into the database;
|
| 69 |
|
|
\item Running re-reconstruction at a Tier-1 centre reading this alignment
|
| 70 |
|
|
object;
|
| 71 |
|
|
\item Running analysis jobs in which ideal, misaligned
|
| 72 |
|
|
and aligned distributions are compared.
|
| 73 |
|
|
\end{itemize}
|
| 74 |
|
|
|
| 75 |
|
|
The steps regarding production of the dedicated AlCaReco streams
|
| 76 |
|
|
as well as enabling the prompt reconstruction to read alignment
|
| 77 |
|
|
constants from the offline database were already described in
|
| 78 |
|
|
section~\ref{sec:offlineswalca}. In the following, a short
|
| 79 |
|
|
summary of the exercise is given. For more details,
|
| 80 |
|
|
see~\cite{alcacsa06note}.
|
| 81 |
|
|
|
| 82 |
|
|
The following misalignment scenario was defined for this exercise: The
|
| 83 |
|
|
sensors of the Tracker Inner Barrel TIB as well as the Rods of the
|
| 84 |
|
|
Tracker Outer Barrel TOB were misaligned. Random shifts drawing from
|
| 85 |
|
|
a flat distribution in the range $\pm 100 \rm\ \mu m$ were applied in
|
| 86 |
|
|
$u$, $v$, $w$ (local coordinate system) for layers built from double
|
| 87 |
|
|
sided modules, and in $u$ (precise coordinate) only for layers built
|
| 88 |
|
|
from single sided modules. In addition, random rotations around all
|
| 89 |
|
|
three local coordinate axes of size $\pm 10 \rm\ mrad$ were applied to
|
| 90 |
|
|
the modules/rods of both single and double sided TIB and TOB layers.
|
| 91 |
|
|
The pixel detector was kept fixed
|
| 92 |
|
|
in order to define a reference system. In addition, the outermost TOB
|
| 93 |
|
|
layer was also kept fixed in order to improve the convergence. This
|
| 94 |
|
|
misalignment scenario was inserted as an alignment object into the
|
| 95 |
|
|
offline database. In addition, another object corresponding to the
|
| 96 |
|
|
ideal tracker geometry was inserted, to be used during prompt
|
| 97 |
|
|
reconstruction.
|
| 98 |
|
|
|
| 99 |
|
|
\begin{figure}
|
| 100 |
|
|
\centering
|
| 101 |
malgeri |
1.6 |
\includegraphics[angle=270,width=0.8\linewidth]{figs/csa06_convergence}
|
| 102 |
fpschill |
1.5 |
\caption{Alignment exercise: The convergence of the alignment
|
| 103 |
|
|
algorithm in $\Delta u, \Delta v, \Delta w$ is shown for double sided
|
| 104 |
|
|
TIB modules as a function of the iteration number (left), as well as
|
| 105 |
|
|
projected for initial misalignment (labeled ``after 0 Iterations'')
|
| 106 |
|
|
and after 4, 7 and 10 iterations (right). }
|
| 107 |
|
|
\label{fig:alignment_convergence}
|
| 108 |
|
|
\end{figure}
|
| 109 |
|
|
|
| 110 |
|
|
The alignment was performed running the HIP alignment
|
| 111 |
|
|
algorithm~\cite{hipnote} as implemented in CMSSW\_1\_0\_6 over approx.
|
| 112 |
|
|
$1 \rm\ M$ AlCaReco $Z^0\rightarrow\mu^+\mu^-$ events produced during
|
| 113 |
|
|
prompt reconstruction at the Tier-0, reading the above mentioned
|
| 114 |
|
|
misalignment object from the offline database. The algorithm was run
|
| 115 |
|
|
on 20 dedicated CPUs in parallel at CERN, iterating 10 times over the
|
| 116 |
|
|
data sample. The result of the alignment was obtained in less than $5
|
| 117 |
|
|
\rm\ h$, and the corresponding tracker alignment object
|
| 118 |
|
|
was inserted into the database. Figure~\ref{fig:alignment_convergence}
|
| 119 |
|
|
illustrates the convergence of the alignment for the double sided TIB
|
| 120 |
|
|
sensors.
|
| 121 |
|
|
|
| 122 |
|
|
\begin{figure}
|
| 123 |
|
|
\centering
|
| 124 |
malgeri |
1.6 |
\includegraphics[angle=270,width=0.8\linewidth]{figs/csa06_aliexercise}
|
| 125 |
fpschill |
1.5 |
\caption{Invariant mass distribution from $Z^0\to \mu^+ \mu^-$ events,
|
| 126 |
|
|
obtained from events produced by the prompt reconstruction at
|
| 127 |
|
|
the Tier-0 (``ideal''), from events processed with misalignment
|
| 128 |
|
|
as used as input for the alignment algorithm (``misaligned'') and
|
| 129 |
|
|
from events re-reconstructed at a Tier-1 centre (PIC) using the alignment
|
| 130 |
|
|
constants derived from the alignment algorithm (``realigned'').
|
| 131 |
|
|
}
|
| 132 |
|
|
\label{fig:aliexercise}
|
| 133 |
|
|
\end{figure}
|
| 134 |
|
|
|
| 135 |
|
|
Once the alignment and calibration constants were inserted in the
|
| 136 |
|
|
database, they were deployed to the Tier-1/2 centres via Frontier.
|
| 137 |
|
|
Subsequently, re-reconstruction of some of the CSA06 datasets was
|
| 138 |
|
|
launched at various Tier-1 centres. For instance, the
|
| 139 |
|
|
$Z^0\rightarrow\mu^+\mu^-$ data set was re-reconstructed at PIC
|
| 140 |
|
|
(Barcelona) using the new alignment object. In order to demonstrate
|
| 141 |
|
|
the final missing piece of the workflow, grid analysis jobs were
|
| 142 |
|
|
submitted to PIC to process these re-reconstructed
|
| 143 |
|
|
$Z^0\rightarrow\mu^+\mu^-$ events. The reconstructed invariant di-muon
|
| 144 |
|
|
mass is presented in Figure~\ref{fig:aliexercise} for three cases:
|
| 145 |
|
|
\begin{itemize}
|
| 146 |
|
|
\item Using the ideal geometry, reading the RECO produced during prompt
|
| 147 |
|
|
reconstruction;
|
| 148 |
|
|
\item Using the misaligned geometry, reading the AlCaReco and the misalignment
|
| 149 |
|
|
scenario database object;
|
| 150 |
|
|
\item Using the realigned geometry, reading the events re-reconstructed
|
| 151 |
|
|
with the alignment database object.
|
| 152 |
|
|
\end{itemize}
|
| 153 |
|
|
As can be seen, the invariant mass resolution is degraded in the case
|
| 154 |
|
|
of misalignment. After the alignment algorithm has corrected the
|
| 155 |
|
|
tracker geometry, the resolution is recovered close to the original
|
| 156 |
|
|
value. This demonstrates that the work- and dataflow of the full
|
| 157 |
|
|
alignment was successfully carried out.
|
| 158 |
|
|
|
| 159 |
|
|
|
| 160 |
acosta |
1.16 |
\input{muonalignment.tex}
|
| 161 |
|
|
|
| 162 |
fpschill |
1.5 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
| 163 |
fisk |
1.1 |
|
| 164 |
|
|
\subsection {Physics Analysis Exercises}
|
| 165 |
|
|
|
| 166 |
ndefilip |
1.2 |
\subsubsection {Effect of tracker misalignment on track reconstruction performances}
|
| 167 |
|
|
|
| 168 |
ndefilip |
1.3 |
The alignment uncertainties of the CMS Tracker detector, made of a huge amount of independent silicon sensors with an excellent position resolution, affect the performances of the track reconstruction and track parameters measurement.
|
| 169 |
|
|
The analysis exercise performed by the team at Bari during the CSA06 had the purpose of study the effect of the CMS tracker misalignment on the performances of the track reconstruction \cite{misalignment}. Realistic estimates for the expected displacements of the tracking systems were supplied in different scenarios as specified in the following:
|
| 170 |
|
|
|
| 171 |
|
|
\begin{itemize}
|
| 172 |
|
|
\item the ideal scenario with a perfect tracker geometry;
|
| 173 |
|
|
\item the short term misalignment scenario supposed to reproduce the mis-alignment conditions during the first
|
| 174 |
|
|
data taking when the uncertainties on the position of the sub-structures of the CMS
|
| 175 |
|
|
tracker will be between $10 \, \mu$ for pixel detectors and $400 \, \mu$ for microstrip silicon detectors in
|
| 176 |
|
|
the endcaps. Detector position and errors are read
|
| 177 |
|
|
from the offline database at CERN by caching the needed information locally via frontier/squid
|
| 178 |
ndefilip |
1.7 |
software.
|
| 179 |
ndefilip |
1.3 |
|
| 180 |
|
|
\item the long term scenario when the alignment uncertainties are supposed to be a factor 10 smaller because of the
|
| 181 |
|
|
improvement obtained by using aligmnent algorithms with a high statistics of tracks.
|
| 182 |
|
|
|
| 183 |
|
|
\item the CSA06 aligned scenario by using the tracker module position and errors as obtained by the output of the
|
| 184 |
|
|
alignment procedure that was run at CERN Tier-0 to verify the efficiency of the alignment procedure on the track
|
| 185 |
|
|
reconstruction. The refit of tracks is performed also in this case.
|
| 186 |
|
|
|
| 187 |
|
|
\end{itemize}
|
| 188 |
|
|
|
| 189 |
ndefilip |
1.4 |
Track reconstruction is based on the Kalman Filter formalism \cite{Kalman} for trajectory building, cleaning and
|
| 190 |
|
|
smoothing steps and uses hits from pixel detector as seeds to provide initial trajectory candidates.
|
| 191 |
|
|
Because of the misalignment the analysis requires to refit tracks with a misaligned tracker geometry.
|
| 192 |
|
|
Global efficiency of track recostruction and track parameter resolutions for muons were
|
| 193 |
|
|
compared in all the cases. The association between simulated track and reconstruted tracks is performed
|
| 194 |
|
|
by comparing the corresponding track parameters at the closest approach point and choosing the pair which gives
|
| 195 |
|
|
the minimum $\chi^2$ from the best fit procedure.
|
| 196 |
ndefilip |
1.3 |
|
| 197 |
ndefilip |
1.4 |
Events from CSA06 $Z\rightarrow \mu \mu $ sample were firstly skimmed by selecting events with
|
| 198 |
|
|
Hep MC muons from Z decay with pseudorapidity, $\eta$, in the tracker acceptance, $|\eta| < 2.55$,
|
| 199 |
|
|
with transverse momentum larger than $5 \, \mathrm{GeV}/c^2$ and di-muons invariant mass in the following
|
| 200 |
|
|
range aroung the Z peak: $50 < m_{\mu\mu}(\mathrm{GeV}/c^2) < 130$; the efficiency of the previous selection is
|
| 201 |
|
|
between 50 and 60 \% mainly due to the cut on the acceptance, for a final statistics of 1 million events.
|
| 202 |
|
|
The output files in RECOSIM format were needed for the subsequent analysis.
|
| 203 |
ndefilip |
1.3 |
|
| 204 |
|
|
|
| 205 |
ndefilip |
1.4 |
Jobs executing the misalignment analysis were submitted at Bari with CRAB\_1\_4\_0 in the LCG infrastructure.
|
| 206 |
|
|
A total of about 2.5 thousands jobs (45 at most in parallel) ran with a grid efficiency of 90 \% and an
|
| 207 |
|
|
application efficiency of 80\%, by accessing detector position and errors from the offline database via
|
| 208 |
|
|
frontier.
|
| 209 |
ndefilip |
1.3 |
|
| 210 |
ndefilip |
1.4 |
|
| 211 |
|
|
Some results of the misalignment analysis were summarized below. The global efficiency of track reconstruction of muons coming from
|
| 212 |
|
|
Z decay is shown in Fig.~\ref{eff} as a function of the pseudorapidity, $\eta$, in the tracker acceptance.
|
| 213 |
|
|
In the case of a perfect geometry the global track reconstruction was not fully efficient over all the $\eta$
|
| 214 |
ndefilip |
1.12 |
range because of the track associator algorithm itself which discards tracks with $\chi^2$ of the fit larger
|
| 215 |
ndefilip |
1.4 |
than 25. The effect of misalignment is relevant in the short term scenario and causes a partial inefficiency of the
|
| 216 |
|
|
track reconstruction; that can be recovered if the intrinsic position resolution of the tracker
|
| 217 |
|
|
detector is combined with the alignment uncertainties to make larger the error on the position of the
|
| 218 |
|
|
reconstructed hit (called alignment position error, APE) so improving the track fit at the
|
| 219 |
|
|
expence of a larger rate of fake tracks.
|
| 220 |
|
|
|
| 221 |
|
|
|
| 222 |
|
|
The tranverse momentum resolution as a function of the transverse momentum is reported in Fig.~\ref{respt};
|
| 223 |
|
|
the degradation of the tranverse momentum resolution at large $p_{T}$ because of the misalignment is in a factor
|
| 224 |
|
|
between 2 and 3 with respect to the perfect geometry case. At low transverse momentum (less than few $\mathrm{GeV}/c$
|
| 225 |
|
|
the multiple scattering is the
|
| 226 |
|
|
most important contribution to the resolution so the effect of misalignment is overwhelmed at all.
|
| 227 |
|
|
|
| 228 |
|
|
The residual of Z mass obtained as the invariant mass of muons coming from Z decay in the case of perfect
|
| 229 |
|
|
tracker geometry and in short-term and long term misalignment scenarios is shown in Fig.~\ref{mz}; the $\sigma$ of
|
| 230 |
|
|
the Gaussian fit of the residual ditribution can be quoted as the Z mass resolution which is degradated of a factor 2
|
| 231 |
|
|
because of the tracker misalignment in the short term scenario.
|
| 232 |
|
|
|
| 233 |
|
|
\begin{2figures}{hbt}
|
| 234 |
malgeri |
1.6 |
\resizebox{\linewidth}{!}{\includegraphics{figs/Eff_eta}} &
|
| 235 |
|
|
\resizebox{\linewidth}{!}{\includegraphics{figs/SigmapT_pT}} \\
|
| 236 |
ndefilip |
1.3 |
\caption{Global track reconstrution efficiency vs pseudorapidity for muons coming from Z decay in the case of perfect
|
| 237 |
|
|
tracker geometry and in short-term and long term misalignment scenarios when the APE is not used.}
|
| 238 |
|
|
\label{eff} &
|
| 239 |
ndefilip |
1.4 |
\caption{$P_{T}$ resolution vs $p_{T}$ in the case of perfect
|
| 240 |
|
|
tracker geometry and in short-term and long term misalignment scenarios.}
|
| 241 |
|
|
\label{respt} \\
|
| 242 |
|
|
\end{2figures}
|
| 243 |
|
|
|
| 244 |
|
|
|
| 245 |
|
|
\begin{figure}[htb]
|
| 246 |
|
|
\begin{center}
|
| 247 |
malgeri |
1.6 |
\resizebox{0.7\linewidth}{!}{\includegraphics{figs/Residual_mZ_mu}}
|
| 248 |
ndefilip |
1.4 |
\end{center}
|
| 249 |
ndefilip |
1.3 |
\caption{Residual of Z mass obtained as the invariant mass of muons coming from Z decay in the case of perfect
|
| 250 |
|
|
tracker geometry and in short-term and long term misalignment scenarios.}
|
| 251 |
ndefilip |
1.4 |
\label{mz}
|
| 252 |
|
|
\end{figure}
|
| 253 |
|
|
|
| 254 |
acosta |
1.10 |
% min bias studies
|
| 255 |
|
|
\input{mbue.tex}
|
| 256 |
acosta |
1.8 |
|
| 257 |
|
|
% tau analyses:
|
| 258 |
ndefilip |
1.4 |
|
| 259 |
acosta |
1.8 |
\input{z_2tau.tex}
|
| 260 |
|
|
|
| 261 |
|
|
\input{taumisid.tex}
|
| 262 |
|
|
|
| 263 |
|
|
\input{tau_validation.tex}
|
| 264 |
|
|
|
| 265 |
acosta |
1.9 |
% muon analyses
|
| 266 |
|
|
|
| 267 |
|
|
\input{wmunu.tex}
|
| 268 |
|
|
|
| 269 |
acosta |
1.13 |
\input{dimuon.tex}
|
| 270 |
|
|
|
| 271 |
meridian |
1.11 |
% electron analyses
|
| 272 |
|
|
\input{zee.tex}
|
| 273 |
|
|
|
| 274 |
ndefilip |
1.12 |
\input{wmass.tex}
|
