# Presentation: 320211

## Charmless Two Body b Decays from CDF

- Simone Donati, Michael J. Morello, Giovanni Punzi, Diego Tonelli
- S.DeCecco, M.R., S.Giagu, G.Salamanna, S.D'auria, D.Lucchesi,
- S.DaRonco, R.Carosi, P.Catastini, A.Ciocci, P.Squillacioti, S.Torre, M.Casarsa
- For the CDF collaboration
- Fermilab W&C 8/6/2004

## Outline

- A bit of history
- Hadronic b trigger (SVT)
- Motivation
- Method
- Results
- Prospects
- Summary

- W&C - Two body B decays with CDF data 6/8/04

## SUMMARY REPORT OF THE WORKING GROUP ON MEASUREMENT OF ANGLE ALPHA WORKSHOP ON B PHYSICS AT HADRON COLLIDERS, SNOWMASS, JUNE 1993

- SUMMARY REPORT OF THE WORKING GROUP ON MEASUREMENT OF ANGLE ALPHA WORKSHOP ON B PHYSICS AT HADRON COLLIDERS, SNOWMASS, JUNE 1993

- (2001-200?)

- W&C - Two body B decays with CDF data 6/8/04

## From design to reality

- W&C - Two body B decays with CDF data 6/8/04

## Silicon Vertex Tracker 

- Silicon Vertex Tracker

**Online**
**Track**
**Impact**
**Param**.

- 35m33m
-   resolbeam
-     47m

- W&C - Two body B decays with CDF data 6/8/04

## Motivations

## BdK vs Bd

- BdK vs Bd

- Original observation of Cleo:
- BR(K)~=4BR()
- made it clear that penguins are an important element in two body B decays.

- Tree is Cabibbo-suppressed Penguin dominant in BdK+ decays
- (Similar conclusion apply  to its SU(3) partner BsKK decay)

- Revised strategies to extract CKM matrix elements:
- use several Isospin related decays modes
- Measure  Bs and Bd SU(3) related modes (hadron colliders !)

- W&C - Two body B decays with CDF data 6/8/04

## Why Bhh’ at CDF (1) 

- Penguin-tree interference provide information on angle 
- A measure of the relative penguin to tree contribution needed
- Use of both Bshh’ processes and Bdhh’, relies on flavour SU(3) symmetry and provides a theoretically clean method
- CDF is the first experiment to be sensitive to Bshh decays
- Measure BR( BsK+K ) relative to BR(BdK+

- Significant direct CP violation observed in B decay in the channel BdK+ (BABAR hep-ex/0407057) !
- CDF measurement of ACP(BdK+ ) consistent with Babar & Belle.
- Will eventually be competitive with B-factories with more data
- CDF will eventually observe BsK- and measure CP asymmetry for this mode also
- Direct CP asymmetry fix theory parameters and allow more precise predictions for yet to be observed observables

- W&C - Two body B decays with CDF data 6/8/04

## Why Bhh’ at CDF (2)  

- Measurement of the time evolution of the untagged BsK+K- sample sensitive to s!
- If  s in BsK+K- turns out different from the time evolution of the CP even component in BsJ/ he may have hint for New Physics

- CDF is searching for charmless two body b decays where large CP violation is predicted by theory

- W&C - Two body B decays with CDF data 6/8/04

## Bd(s)hh’  penguin and tree

- Amplitudes related by U-spin symmetry of strong interactions (sd interchange) !

- W&C - Two body B decays with CDF data 6/8/04

## Bhh observables

- Many related observables determine the angle  using  BsKK/Bd CDF data alone
- Time dependent CP asymmetry requires b-flavor tagging and need more statistics
- (a major goal for CDF Run II)

- R.Fleisher hep-ph/0405091;
- D.London, J.Matias hep-ph/0404009

- Phase space factor = 0.92

- QCD sum rules: 1.76+0.15-0.17 (A.Khodyamirian et al., Phys.Rev D68 114007)

- 

- 

- 

- 

- 

- 

- 

- Branching Ratio measurements can constrain theory too!

- W&C - Two body B decays with CDF data 6/8/04

## BsKK vs Bd

- BsKK vs Bd

- Determine allowed region in the R vs ACPdir() plan from Babar & Belle measurements
- Check Theory (or claim New Physics...) by comparing the allowed range with CDF experiment data (Lepton Photon 03 prel. result)

- Inputs:
- 2 = 47o=65o±7o
- Current Babar/Belle measurement
- ±20% SU(3) breaking effect as additional theory error

- W&C - Two body B decays with CDF data 6/8/04

## Quadrant of CDF II Tracker

**Quadrant of CDF II Tracker**

**LAYER 00**: 1 layer of radiation-hard silicon at very small radius (1.5 cm)
-                   (expected 50 fs proper time resolution in Bs  Ds )

**COT:** large radius (1.4 m) Drift C.
- 96 layers, 200ns drift time
- Precise PT above 400 MeV/c
- Precise 3D tracking in |****|<1
- (1/PT) ~ 0.1%GeV –1; (hit)~150****m
- dE/dx info provides >1.3 sigma K/****  separation above 2 GeV

**SVX-II + ISL:**  6 (7) layers of double-side silicon (3cm < R < 30cm)
-  Standalone 3D tracking up to |****|= 2
-  Very good I.P. resolution: ~30****m (~20 ****m with Layer00)

**TOF:** 100ps resolution, 2 sigma K/**** separation for tracks below 1.6 GeV/c
- (significant improvement of Bs flavor tag effectiveness)

**TIME OF FLIGHT**

- W&C - Two body B decays with CDF data 6/8/04

## Measurement of the relative fractions and CP asymmetry in B0d,s  h+h’-  

- W&C - Two body B decays with CDF data 6/8/04

## B0  h+h’- mass plot – trigger like selection

- B0  h+h’- mass plot – trigger like selection

- Online trigger cuts are just a little looser!!!

- Two-Track data, 180 pb-1

- Selection optimized using sideband data and MC

- W&C - Two body B decays with CDF data 6/8/04

## B0  h+h’- mass plot – isolation cut added

- B0  h+h’- mass plot – isolation cut added

- Define B isolation as :

- Where the sum over all charged tracks within a cone of radius R=1 around candidate B meson direction

- W&C - Two body B decays with CDF data 6/8/04

## Strategy

- To measure Branching Franction and CP asymmetry we need to separate the 4 signals superimposed in the mass peak. Fit the composition of the B0  h+h’-  signal with a likelihood that combines the invariant mass (M ,the kinematics and PID information (dE/dx from drift chamber).

- Notice how the BsKK and Bd sit one on top of the other. dE/dx separation crucial for separate the two.

- W&C - Two body B decays with CDF data 6/8/04

## Handle 1: kinematics

- 

- Invariant mass ( hypothesis) vs signed momentum
- imbalance ******=[1–p1/p2] x q1, **discriminates among the
- four signals and the B flavour for flavour specif decays.

- W&C - Two body B decays with CDF data 6/8/04

## Handle 2: PID from dE/dx 

- Handle 2: PID from dE/dx

**1.39 ******** ********/K separation**
**For P>2 GeV/c**

** ****Note:**
**The combined fit with PID at his level provides an effective separation only 40% worse than a perfect PID would do...**

- W&C - Two body B decays with CDF data 6/8/04

## Handle 2 (aside): huge calibration sample from D*±D0± decays !

- CDF  is accumulating huge sample of clean D signals with the SVT
- Crucial point for accurate calibration and understanding of the dE/dx measurement
- Crucial for understanding trigger efficiencies and its dependence on particle type

- (_)

- World best BR and CP asymmetry in D0K+K- and  final state from CDF submitted for publication today!

- Purity 98 %

- 4 *105D0

- W&C - Two body B decays with CDF data 6/8/04

## The  Likelihood

- Fit the events falling in the range:
- (4.85 < M < 5.8) U ( -0.8 <  < 0.8) U (all dE/dx)

- Background fraction (float)

- Background Likelihood

- Signal Likel.

- W&C - Two body B decays with CDF data 6/8/04

## The Likelihood (cont’d)

- dE/dx term

- kinematics term

- fraction of events of the jth mode

- Fit the fraction fj of each of the 6 modes:
-           Bd + - Bd  K+ -	Bd  K- +                               	Bs  K+K-	Bs  K- +	Bs  K+ -
- With normalization condition

- Signal likelihood

- Background likelihood

- W&C - Two body B decays with CDF data 6/8/04

## The kinematics  term (signal)

- Expected value of the** ******invariant mass (function of** ******). Input the fit with its  analytic expression

- Pdf of ****** **variable. Template from MC

- Invariant mass width assuming tracks in each mode assigned with correct masses

- W&C - Two body B decays with CDF data 6/8/04

## The kinematics term (BCKG)

- Pdf of the **** variable.  Extract from separated fit  on data

- Mass shape of the background.  Parameters c_i floating in the fit

- W&C - Two body B decays with CDF data 6/8/04

## dE/dx Likelihood (signal)

- PID is the mean of the Gaussian:
- PIDK = 1, PID  = 0

*********  is the RMS of the Gaussian, it depends on the event and on the mass hypothesis

- For Bd,s  

- W&C - Two body B decays with CDF data 6/8/04

## dE/dx Likelihood (background)

- Assume the background made of pions and kaons only, the BCKG Likelihood is

- f  is the pion fraction in background
- Floating in the fit

- Correlations between dE/dx measurements of two tracks in the same event (a.k.a. common-mode fluctuations) potentially bias the result of the fit.
- The likelihood takes in to account this effect (although the formulas are considerably complicated and are not displayed here)

- W&C - Two body B decays with CDF data 6/8/04

## Fit Results

| Decay | # B |
| --- | --- |
| BdK+ | 509 |
| Bd+ | 134 |
| BsK+ | 232 |
| BsK- |  --- |

- Biggest sample of Bs decays!
- BdK+- 509/180 pb-1
- (compare Babar 1600/220 fb-1)

- W&C - Two body B decays with CDF data 6/8/04

## Extraction of ACP

- The RAW fit results need to becorrected for relative acceptance, trigger and selection efficiency:

- from Monte Carlo derive the ratio of efficiency for K+ and K- due to the different nuclear interaction rate for K+ and K- with detector material (1% correction)

- W&C - Two body B decays with CDF data 6/8/04

## Bd/BdK Branching Ratio Ratio

- The RAW fit results need to becorrected for relative acceptance, trigger and selection efficiency:

- Get from Monte Carlo simulations the ratio of efficiencies from Kinematics and Kaon vs pion decays in flight and interaction probability.

- Correct for specific ionization dependence of trigger efficiency in the level1 trigger (XFT). Use data from unbiased legs in D+K- control sample to derive correction.

- W&C - Two body B decays with CDF data 6/8/04

## BsKK/BdK Branching Ratio Ratio

- The RAW fit results need to becorrected for relative acceptance, trigger and selection efficiency:

- At colliders we can only measure the product of production fractions time BR. fs(d) is the probability that a b quark hadronize in a Bs(d) meson. Averages exists in the PDG (dominated by LEP measurements and b time integrated mixing). Interest in an indpendent CDF measurement!

- Fragmentation process might be different for
- Bs and Bd meson. Derive the efficiency of the Isolation cut from samples of fully reconstructed Bs/Bd mesons

- W&C - Two body B decays with CDF data 6/8/04

## Production Pt spectra

- B hh trigger accept very soft B  big samples available!
- Measurement of the production Pt spectrum from inclusive bJX in this region important for reliable MC simulation

- Overflow bin!

- W&C - Two body B decays with CDF data 6/8/04

## Efficiency corrections (3) 

- 

- W&C - Two body B decays with CDF data 6/8/04

## Results

- Annihilation dominated modes
- For Bs assume same time evolution as for the BsKK (see next)

- W&C - Two body B decays with CDF data 6/8/04

##  Sensitivity to a sizeable s

- CDF recently measured an anomalously large (even if not yet statistically compelling) value of the lifetime difference in the Bs system
- Use angular analysis to project out CP-odd and CP-even component of BsJ/
- What’s the implication on the two body decays measurement ?

- W&C - Two body B decays with CDF data 6/8/04

##  Sensitivity to a sizeable s

- = 0.50 ± 0.08 ± 0.09

- dN(BsKK)/dt  RLexp(-tL)+RHexp(-tH)
- The above central value assumes:
- RH=0  (no Heavy Bs decay to KK or, equivalently, no tree contribution)
-  L = 1/(s+s/2) = 1.45 ps
- (from SM s/s = 0.12 ± 0.06
- and s = d )
- Yellow band shows acceptance corrected BR for any assumed effective BsKK lifetime on the X axis

- Central value
- for L  ( CDF)

- W&C - Two body B decays with CDF data 6/8/04

## Summary of systematics 

- W&C - Two body B decays with CDF data 6/8/04

## Systematics  (1)

- a) Mass resolution: the mass resolution is input from MC. It is rescaled to match the D0 resolution on data.
- b) dE/dx - uncertainty on RMS(s): repeat the fit varying the correlation RMS from its minimum (0.24) to its maximum (0.52), quote the differences wrt to central fit.
- c) dE/dx - uncertainty on the shape p(s): repeat the fit at central value of RMS(s) = 0.38 assuming a Double Dirac delta for p(s) shape quote the difference wrt to central fit.
- d) dE/dx – tails: the central fit assumes Gaussian dEdx pdf. Repeat the fit with more accurate parameterization of the pdf; repeat the fit adding extra component.

- W&C - Two body B decays with CDF data 6/8/04

## Systematics  (2) 

- e) input masses: the fit is done on data in which the recipe used for mass measurement at CDF II was applied. Input masses in the kinematics pdf are those measured by CDF II. Repeat the fit varying M(Bd) and M(Bs) within their statistical uncertainties (0.92 and 1.29 MeV/c2).       Quote the differences wrt the central fit.
- f) Background model: the fit assumes mass spectrum of bckg = exp + C. Repeat the fit with p2,p3,p4 and quote the difference wrt central value.
- g) B lifetimes: relative kinematics efficiencies depend on the lifetime assumed in MC. Re-evaluate efficiencies after simultaneous shift of  Bs lifetime (+1) and  Bd (-1)  and viceversa.  is the PDG2004 uncertainty. Quote difference wrt central value.

- W&C - Two body B decays with CDF data 6/8/04

## Systematics  (3) 

- h) Isolation efficiency: has a ~ 10% from measurement on data. Reevaluate the efficiency at +/- 1 and quote difference wrt central value
- i) MC statistics: kinematics efficiencies have statistical error. Reevaluate them at +/- 1  and quote difference wrt the central fit.
- l) Trigger dE/dx correction:  the correction function have uncertainties. Stretch (push) K/ discrepancy shifting simultaneously the correction coefficients by 1, reevaluate the correction, and quote differences wrt the central fit

- W&C - Two body B decays with CDF data 6/8/04

##  Systematics  (4)

- m) ******s****/********s**   (Standard Model)
******s****/********s** : Standard Model predicts ~ 0.12±0.06 and
**B****s**** ********K****+****K****-** to be dominated by the short-lived component.
- We derive the systematic uncertainties from these assumptions by varying s/s from 0.06 to 0.18 , re-evaluating the relative efficiencies and quoting the differences wrt the central fit.

- W&C - Two body B decays with CDF data 6/8/04

## Experimental Comparison

|  | CDF | Babar | Belle |
| --- | --- | --- | --- |
| N(BdK+) | 509/180 pb-1 | 1600/200 fb-1 | 1030/140 fb-1 |
| ACP (BdK+) | -0.04±0.08±0.01 | -0.133±0.030±0.009 | -0.085±0.030±0.013 |
| BR()/BR(K) | 0.24±0.08 | 0.25± 0.04 | 0.24 ± 0.04 |

- ACP measurement with systematic uncertainty comparable to Babar/Belle
- Expect to reach comparable stat. precision with 0.8 fb-1 of data (2005?)
- Ratio of Bd Branching Ratio consistent with world average and provide valuable cross-check for the other Branching ratio measurements

- W&C - Two body B decays with CDF data 6/8/04

## What’s behind the corner?

- W&C - Two body B decays with CDF data 6/8/04

## Where are the b?

- Used the same data to look for evidence of   b’s
- Large direct CP asymmetries expected
- Theory predicts:
- BR(bpK), BR(bp) ~ 10-6 -2*10-6 (Mohanta, Phys. Rev. D63:074001, 2001)
- Current limits:
- BR(bpK) <50*10-6, BR(bp) <50*10-6 @90% C.L.
- Blind optimization to reduce background in the b mass region including the contribution from B0hh’
- Normalize to BR(Bd0->K
- Extract number of BR(Bd0->Kevents from a fit like the one described before.

-      Use BR(Bd0K)=(17.4±1.5) 10-5 and f/fd=0.25±0.04 Obtain :
-   BR(bp)+ BR(bp) < 22 *10-6

- Improved sensitivity in the future with combined TOF and dE/dx PID proton identification

- W&C - Two body B decays with CDF data 6/8/04

## Perspective for futher measurements

- CDF has twice the luminosity on tape wrt to this analysis, better tracking alignment and reconstruction make mass and vertex resolution even better:
- Bs-->Kshould be observable
- Direct CP asymmetry
- BsKK lifetime will be measured soon and will give interesting information on s/s
- Interesting resolution even with current statistics if s so large
- If s different from what is observed in J/   new physics
- Direct CP asymmetry will be competitive with current B-Factory with twice the data we have on tape by now
- Time dependent CP violation measurement interesting further down the road due to the low flavor tagging efficiency at hadron machines (needs full nominal Run II luminosity)
- standalone and clean measurment of CKM angle  is the ultimate goal

- Need to keep the trigger working fine even at high instantaneous luminosity
- Relentless day-to-day fight to save bandwidth!
- Important trigger/daq upgrades in the pipeline

- W&C - Two body B decays with CDF data 6/8/04

## Conclusion

- New exciting result from CDF:

- No evidence (yet) for the BsK decay
- No large annhilation in Bs
- Other charmless mode have been measured or searched for with available data, expect many new results with increasing data size and improved offline quality. Stay tuned...

- W&C - Two body B decays with CDF data 6/8/04

## Summary of CDF result for HFAG

- ...for new CDF result also on Heavy Flavour Averaging Group pages

- W&C - Two body B decays with CDF data 6/8/04

## Backup Slides

- W&C - Two body B decays with CDF data 6/8/04

## Slide 48

- W&C - Two body B decays with CDF data 6/8/04

## Raw fit results

- BCKG shape

- Pion fraction in BCKG

- W&C - Two body B decays with CDF data 6/8/04

## dE/dx in the Likelihood

- dE/dx information is included through the ID variable:

- ID can be written as a function of the dE/dx pulls

- W&C - Two body B decays with CDF data 6/8/04

## dE/dx in the Likelihood (cont’d)

- Write pull = pull (ID, *)

- We assume the pull distribution to be Gaussian. Actually pdf(pull) has a small tail at high values, consider this effect in the systematics. Therefore the pdf(ID)  becomes Gaussian with:
- Mean:  0 for pions and 1 for kaons
- RMS:  *  depends on the mass hypothesis and changes
-            event by event

- W&C - Two body B decays with CDF data 6/8/04

## Final results

- [LP03: 0.74 ± 0.2 ± 0.22]

- [LP03: 0.02 ± 0.15 ± 0.017]

- [LP03: 0.26 ± 0.11 ± 0.055]

- W&C - Two body B decays with CDF data 6/8/04

## Disentangling the  Bh+h- contributions (I)

**Disentangling the  B********h+h- contributions (I)**

**Must **_**disentangle **_**contributions from each mode **
**To do this we use:**
**-****Kinematical**** ****variable  separation**** ****M******** vs ********=(1-p****1****/p****2****)********q****1 **
**-dE/dx based K and ******** identification**** **

- W&C - Two body B decays with CDF data 6/8/04

## Disentangling the  Bh+h- contributions (II)

**M******** vs a for each B********h****+****h****-**** mode**

**dE/dx calibration using D*****± ******** ****D****0************, **
**D****0 ********K****±********±**** ****(******** from D*  unambiguously distinguishes K, ******** from D****0****)**** **

_**Sanity check: Measure Ratio of Branching Ratios**_
_**CDF :**_** ********(B****d ******** ****************-****)/********(B****d******** ****K****+********-****) = 0.26 ±0.11±0.055, ****PDG:**** **

_**Yield for each mode:**_
**B****d ******** **********************
**B****d ******** ****K****±**** ********±**********
**B****s ******** ****K****± ********±**********
**B****s ******** ****K****+****K****-  ********(stat) ********(stat)******
******First observation !**

_**Method works !**_**  Confirmed by Sanity **
**check against ratio of branching ratios**
_**Have first observation of B**__**s **_****** **_**K+K-**_
_**Its a CP Eigenstate: Can use this**_
_**To measure **__****__**s**__** as well !!**_

- W&C - Two body B decays with CDF data 6/8/04

## Where are the b (syst.)?

- Where are the b (syst.)?

- Largest source of syst is the uncertainty on the Pt spectra, and the production fraction.
- However, the impact on the limit is small.
- Limit most sensitive to background uncertainty.

- W&C - Two body B decays with CDF data 6/8/04