1 |
\section{Introduction}
|
2 |
The decay $B_s\rightarrow J/ \psi \phi\rightarrow \mu^+\mu^-K^+K^-$ is of
|
3 |
particular interest since it allows the study of many properties of the $B^0_s$
|
4 |
sytem, such as the differences between the widths and the masses of the two
|
5 |
weak eigenstates, $B_s^H$ and $B_s^L$.
|
6 |
The measurement of the interference between mixing and decay might reveal $CP$
|
7 |
violation other than predicted by the Standard Model of Particle Physics
|
8 |
indicating the presence of new generations of particles and forces between
|
9 |
them~\cite{physics}. Since $J/\psi \phi$ is a vector-vector final state with
|
10 |
an unknown mix of $CP$-even and $CP$-odd amplitudes the
|
11 |
extraction of $CP$ asymmetries requires an angular analysis~\cite{angular}.
|
12 |
|
13 |
The objective for the Summer'10 conferences is to observe the decay
|
14 |
statistically significant and if data permits to extract the $B_s$ lifetime.
|
15 |
These measurements do strongly depend on the secondary vertex reconstruction
|
16 |
capabilities of the CMS pixel detector and hence serve as benchmark
|
17 |
in direct comparison to measurements from the Tevatron~\cite{pdg}.
|
18 |
As in the $B^0_d$ system, the average lifetime of $B_S$ is
|
19 |
due to the $b$ quark; it is well measured by the CDF~($\tau=1.40\pm 0.14(stat)
|
20 |
\pm 0.02(syst)\times 10^{-12}$~s)
|
21 |
and D0 experiments~($\tau=1.444\pm 0.094(stat)\pm 0.020(syst)\times 10^{-12}$~s)
|
22 |
at the Tevatron of Fermilab.
|
23 |
|
24 |
\subsection{Strategy}
|
25 |
We attempt to measure the cross section for the channel, total and
|
26 |
as function of the transverse momentum, and the $B_s$ lifetime.
|
27 |
Previous studies~\cite{lotte} indicate that the potential background from
|
28 |
prompt $J/\psi$ and inclusive $b$ production can be
|
29 |
identified and separated from the signal using the invariant $J/\psi\phi$
|
30 |
mass and the proper decay length $ct$ of reconstructed $B_s$ candidates.
|
31 |
We study these and QCD backgrounds with full detector Monte Carlo simulations
|
32 |
corresponding to a luminosity at least comparable or exceeding the one
|
33 |
expected for real data to develop our strategy.
|
34 |
We propose two approaches: one based on optimized selection criteria
|
35 |
for event variables (''cut-and-count'') and a two-dimensional maximum likelihood
|
36 |
fit to the inavriant $J/\psi \phi$ mass and proper decay length $ct$ with relaxed
|
37 |
requirements on event variables.
|
38 |
|
39 |
We estimate of the number of $B_s\to J/\psi(\to\mu^+\mu^-) \phi(\to K^+ K^-))$
|
40 |
events produced in proton-proton collisions at a center-of-mass (CM) energy of 7~TeV
|
41 |
from the cross section for an integrated luminosity of 1~pb$^{-1}$ the following:
|
42 |
\begin{eqnarray}
|
43 |
& & \sigma(B_s\to J/\psi(\to \mu^+\mu^-) \phi(\to K^+ K^-)) = \nonumber \\
|
44 |
& & \hspace*{4mm} \sigma(pp\to \bar{b}b)\cdot
|
45 |
2 \cdot B(\bar{b}\to B_s)\cdot \Gamma(B_s\to J/\psi \phi)\cdot \Gamma(J/\psi\to\mu^+\mu^-)\cdot
|
46 |
\Gamma(\phi\to K^+ K^-)
|
47 |
\end{eqnarray}
|
48 |
The value of total cross section for the production of $\bar{b}b$ pairs in proton-proton
|
49 |
collisions at a CM energy of 7~TeV is expected to be approximately half of the value
|
50 |
at 14~TeV, namely 250~$\mu$b~\cite{nason,campbell}. With the branching fractions~\cite{pdg}
|
51 |
listed in Table~\ref{tab:pdgtab} we predict that about 2020 events are produced.
|
52 |
\begin{table}[ht]
|
53 |
\begin{center}
|
54 |
\begin{tabular}{ll}
|
55 |
\hline
|
56 |
$B(\bar{b}\to B_s)$ & $(10.7\pm 1.2)$\% \\
|
57 |
$\Gamma(B_s\to J/\psi \phi)$ & $(1.3\pm 0.4)\cdot 10^{-3}$ \\
|
58 |
$\Gamma(J/\psi\to \mu^+\mu^-)$ & $(5.93\pm 0.06)$\% \\
|
59 |
$\Gamma(\phi\to K^+ K^-)$ & $(48.9 \pm 0.5)$\% \\
|
60 |
\hline
|
61 |
\end{tabular}
|
62 |
\label{tab:pdgtab}
|
63 |
\caption{\sl Cross sections and branching fractions from PDG~\cite{pdg}.}
|
64 |
\end{center}
|
65 |
\end{table}
|
66 |
We establish the efficiency for the trigger and the two reconstruction chains
|
67 |
in this note. We find the product of both efficiencies to be about 2\% which
|
68 |
translates into about 40 observed events for each pb$^{-1}$ of integrated luminosity.
|
69 |
The luminosity is expected to be known to about 10\% which is comparable
|
70 |
to the statistical error for the first 3~pb$^{-1}$.
|
71 |
The largest uncertainties in the estimate above of up to 50\% are due to the unknown
|
72 |
total $pp\to \bar{b}b$ cross section at 7~TeV and the transverse momentum distribution
|
73 |
of the $b$ quarks~\cite{baines}.
|
74 |
|
75 |
The blind analysis method has been widely used by the BaBar collaboration and
|
76 |
experiments at the Tevatron to establish least biased searches, first or
|
77 |
repeated ones~\cite{aaron}. For an integrated luminosity of 1~pb$^{-1}$
|
78 |
we approach a possible signal significance of 5$\sigma$, which makes it
|
79 |
by definition an observation, a first one at CMS, that corroborates this
|
80 |
approach.
|
81 |
|
82 |
\subsection{Implementation}
|
83 |
We establish here the complete analysis with software releases and simulations
|
84 |
available. We expect that the Summer'10 analysis is likely an update with somewhat
|
85 |
improved code and newly generated events for signal and background as they
|
86 |
become available, but methods are the same and efficiencies will not change by much.
|
87 |
The studies presented here have been performed with analysis software release
|
88 |
CMSSW 3.3.6~\cite{cmssw}. Simulations of $pp$ hard scatter use PYTHIA
|
89 |
version 6.416~\cite{pythia}
|
90 |
|
91 |
|
92 |
|
93 |
\subsection{To Do}
|
94 |
For the Summer'10 analysis the following issues still need to be addressed:
|
95 |
|
96 |
\begin{itemize}
|
97 |
\item Reconstruction efficiencies with new software release CMSSW 3.6.x
|
98 |
\item QCD background contribution with newly generated samples exceeding the
|
99 |
expected luminosity in data
|
100 |
\item
|
101 |
\end{itemize}
|
102 |
|
103 |
|
104 |
|
105 |
|
106 |
|
107 |
|
108 |
|
109 |
|
110 |
|
111 |
|
112 |
|