- 44 Views
- Uploaded on
- Presentation posted in: General

Approval Presentation, 17.08.11

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

LHCb-ANA-2011-062

Approval Presentation, 17.08.11

- The final state Ds-K+ is accessible by both Bs and Bs:

- Both diagrams are similar in magnitude, hence large interference between them is possible.
- Using flavour tagging, we can measure four decay rates
- Bs or Bs to Ds+K- or Ds-K+

- From these rates, Î³ can be extracted in an unambiguous and theoretically clean way.
- A review of the Î³extraction can be found in e.g. LHCb note 2007-041.

- A reliable extraction of Î³ requires a considerable amount of data
- Aim is to present first Î³ measurement from Bsâ†’DsK at Moriond 2012.

- First step towards measuring Î³ is to observe Bsâ†’DsK in our data, and measure its branching ratio (BR) relative to Bsâ†’DsÏ€.
- Most systematics cancel in the ratio
- Main differences between the modes are the bachelor PID requirements and the smaller Bsâ†’DsK yield.

- Analysis strategy for the relative BR measurement follows that used for the hadronic fd/fs measurement from 2010 data.
- Independently, the hadronic and semileptonic fd/fs measurements from 2010 data can be combined to extract BR(Bsâ†’DsÏ€). This is then used to measure BR(Bsâ†’DsK) absolutely.

- PDG value is BR(Bsâ†’DsK) = (3.0 Â± 0.7)*10-4(23% relative error)
- Calculated by rescaling Belle(*) measurement: BR(Bsâ†’DsK)= (2.4 +1.2/-1.0(stat)Â± 0.3(syst) Â± 0.3(fs) )*10-4. This uses 7Â±3 signal events.
- In addition, CDF(**) measures BR(Bsâ†’DsK)/BR(Bsâ†’DsÏ€) = 0.097 Â± 0.018 Â± 0.009 with ~100 DsK candidates, using a combined mass-PID fit.

Bsâ†’DsK

Belle

Bsâ†’DsÏ€

CDF

(*) PRL 102 021801

(**) PRL 103 191802

- BR(Bsâ†’DsÏ€) = (3.2 Â± 0.5)*10-3(16% relative error), from combining
- Belle(*): (3.67 Â± 0.34(stat)Â± 0.43(syst) Â± 0.49(fs) )*10-3 with 160 events
- CDF(**) : (3.03 Â±0.21(stat)Â± 0.45(syst) Â± 0.46(fd/fs) )*10-4with 500 events

Belle

CDF

(*) PRL 102 021801

(**) PRL 98 061802, rescaled to new BR(Bdâ†’DÏ€)

- Today we are requesting approval of preliminary results for BR(Bsâ†’DsK)/BR(Bsâ†’DsÏ€), BR(Bsâ†’DsÏ€) and BR(Bsâ†’DsK).
- We plan to write a paper in the very near future.
- Plots for approval are marked with

For Approval

- The analysis uses ~336pb-1 of 2011 data.
- Trigger requirements are
- L0: Hadron TOS or Global TIS
- HLT: HLT1Track and HLT2 Topo BBDT TOS (2,3 or 4 body)

- Stripping lines from B2DX module (with D2hhh)
- No PID is used, to allow inclusive selection of all the relevant modes

- Relative efficiency of reconstruction, trigger and stripping is checked on MC that has been reprocessed with a 2011 TCK (0x006d0032)
- Results on later slides

- For the 2010 fd/fs analysis, TMVA was used to check the performance of different provide classifiers as an offline selection, using kinematic and geometrical variables
- Role of the MVA is to minimise combinatorics (not physics backgrounds)

- The best performing MVA was the Boosted Decision Tree with Gradient boosting (BDTG).

- In the current analysis, the 2010 BDTG is retained, but optimal cut is re-evaluated
- Optimal working point need not be the same, as the trigger has changed

- Re-optimisation for 2011 is performed using 10% of the Bsâ†’DsÏ€ data (uniformly distributed in time)

- Figure of merit is

- Here, B refers to combinatoric bkg only, and 1/14 is the expected Cabibbo suppression factor between Bsâ†’DsÏ€ and Bsâ†’DsK.

- Choose start of significance plateau (BDTG>0.1) as our working point.
- This cut loses 6% of signal, for a background reduction of 45%.

- Correctly calibrating the PID cut efficiencies is crucial for this analysis.
- Extensive use is made of the tools developed by the RICH group.
- The D* (for K and Ï€) and Î› (for proton) calibration samples are binned in momentum and pT, and the resulting efficiency map is used to weight signal events.

- Magnet Up and Magnet Down data are calibrated separately
- PID performance is not constant in time, as RICH calibration needs to be propagated to the more recent data

- Separation between K and Ï€is poor above 100GeV, hence apply a cut of p<100GeV on the bachelor.

K eff

Ï€ misID

Example: DLL(K- Ï€)>5 (1D binning for visualisation)

- For the moment, only the Dsâ†’KKÏ€ mode is considered
- Other modes could be added in the future

- To obtain clean samples of Bsâ†’Dsh, PID cuts need to be applied to the D daughters to suppress Bdâ†’D+h and Î›bâ†’ Î›ch.
- Hence on the Ds+â†’ K-K+Ï€+ candidate we require:
- DLL(K- Ï€) > 0 for the K- and DLL(K- Ï€) < 5 for the Ï€+ (to suppress combinatorics)
- DLL(K- Ï€) > 5 for the K+ (to suppress D+ â†’K-Ï€+Ï€+)
- DLL(p-K) < 0 for the K+ (to suppress Î›c+ â†’K- p Ï€+)
- K+ failing DLL(p-K) < 0 are retained if KpÏ€ mass is outside the Î›c mass window

- Applying these cuts and a mass window of [1944,1990]MeV gives:
- Efficiency of 78% for Bsâ†’DsÏ€ (using momentum distributions from MC),
- MisID of 1.2% for Bdâ†’D+Ï€(using momentum distributions from data),
- MisID of 1.7% for Î›bâ†’ Î›cÏ€ (using momentum distributions from MC).

- After these cuts, the Bsâ†’DsÏ€ peak is rather pure

- For the DsÏ€ fit, a cut of DLL(K-Ï€)<0 is applied, to eliminate any residual contamination from DsK.
- A hard cut of DLL(K-Ï€)>5 is applied before doing the DsK fit, to suppress the favoured DsÏ€ mode.
- As a cross-check, DsK fit is also done with a loose cut of DLL(K-Ï€)>0, and a very tight cut of DLL(K-Ï€)>10.
- The efficiencies of these cuts, applied after the p<100GeV cut, are:

- Ratio of generator level efficiencies is found to be 1.027Â±0.010. Until the reasons for this are understood, the 1.027 is used as a correction factor, and a systematic of 2.7% is applied.

- Ratio of efficiencies for reconstruction, trigger, BDT cut and upper momentum cut on the bachelor is 1.03Â±0.01. This correction factor is applied, and the associated systematic is conservatively set to 3%.

- The B mass uses the D(s) mass constraint (improves resolution).
- Different shapes are tested on the MC signal samples.
- Deafult shape is double Crystal Ball, with common mean & sigma
- Radiative tail is smaller for modes with bachelor K than bachelor Ï€.

Bsâ†’DsÏ€

Bsâ†’DsK

- The physics bkgs to the different modes often involve misidentified hadrons. So getting the misIDâ€™d mass shapes correct is important.
- Example: the shape for DsÏ€ bkg to DsK is obtained as follows:
- Firstly, a clean sample of DsÏ€ is extracted from the Dsh data by applying DLL(K-Ï€)<0 on the bachelor
- This cut biases the bachelor momentum, however the original momentum distribution can be recovered from the whole Dsh sample
- This works because the DsÏ€ and DsKbachelor momenta are very similar

- Then the mass is recomputed under the DsK hypothesis
- Next, the shape is weighted according to the momentum spectrum of the misidentified bachelors
- This from the original momentum distribution and the misID rate as a function of momentum

- The shape for DÏ€ bkg to DsÏ€ is obtained in a similar way, changing D daughter mass hypothesis instead of the bachelor.
- The shape for the DKbkg to DsK should be the same as the DÏ€ bkg to DsÏ€.

- The shapes for DsÏ€ and DÏ€ under the DsK are sufficiently similar that in the fit only the DsÏ€ shape is used

Under the DsK mass hypothesis

- A bump was seen in the DsKfit at around 5500MeV, that was not described by the misidentified DsÏ€ shape.
- The bump was investigated, and it turned out to be Î›bâ†’Dsp!

A peak is also seen at lower mass, compatible with Î›bâ†’Ds*p

- A peak is seen at the Î›b mass after switching to the Dsp mass hypothesis, applying extremely tight PID cuts (DLL(p-Ï€)>10 and DLL(p-K)>15) on the bachelor, and tightening the BDT cut.
- In the future a measurement will be made of the BR of this mode, but for nowâ€¦

- Cutting on DLL(p-K) would lose too much signal, so we must live with this background, and model its shape.
- The shape is taken from simulated events, after reweighting for the efficiency of the DLL(K-Ï€)>5 cut as a function of momentum.
- The Î›bâ†’Ds*p shape is obtained by shifting the Î›bâ†’Dsp shape down by 200MeV. As a baseline, the relative amount of Î›bâ†’Dsp and Î›bâ†’Ds*p is assumed to be the same.

- The amount of Î›bâ†’Dsp in the DsK fit is estimated by taking the 24 events from the previous slide, and correcting for the efficiency of the tight PID cuts and the BDTG cut.
- This gives an expectation of ~150 events (Î›bâ†’Dsp + Î›bâ†’Ds*p)

Î›bâ†’Dsp plus Î›bâ†’Ds*p

- For partially reconstructed physics bkgs, the shapes are taken from MC, with data-driven momentum reweighting applied where a misidentification is involved. PDFs are made using RooKeysPDF.

- One final type of physics background needs to be considered: charmless modes such as Bsâ†’K*KK
- These can appear if no cut is applied on the flight distance of the D from the B vertex
- They can peak under the signal

- To remove such backgrounds, a soft cut of FDÏ‡2(D from B) > 2 is applied. This will have the same efficiency for Bsâ†’DsÏ€ and Bsâ†’DsÏ€, so will not affect the ratio of BRâ€™s.

- The slope of the combinatoric background can floated in the DsÏ€ fit.
- However it must be fixed in the DsKfit, due to the low statistics and the presence of the misidentified Bsâ†’DsÏ€ in the right-hand sideband.
- Fitting to wrong-sign (same-sign D and bachelor) events passing the DsK selection and PID cuts, the slope is compatible with being flat.

- As a cross-check, the wrong-sign events passing the DsÏ€ selection are also fitted, and the slope agrees well with that found in the DsÏ€ signal fit.

DsK wrong-sign

- Since the PID efficiencies vary slightly between MagUp and MagDown, the misID background shapes change.
- In addition, the signal mean is found to shift by ~1MeV between MagUp and MagDown.
- So we split the data by polarity, and fit the two subsamples independently.
- About 55% (45%) of the data is MagDown (MagUp).
- In the following slides, the MagDown fit is on the left, and the MagUp on the right.

- This fit is needed to estimate the amount of background DÏ€ to DsÏ€, and to check the mean and sigma of the signal shape with high statistics.
- The tails of the signal mass shape are fixed from the MC fit, but the mean and sigma are floated
- Mean and sigma are allowed to be different for MagUp and MagDown

- The yields of all components are left free.
- The slope of the combinatoric background is also floated in the fit.

- The expected number of misID Bdâ†’DÏ€ is calculated using the fitted DÏ€ yield, a mass window factor (from MC), and misID from the PID calibration tools.
- It is constrained to be within 10% of this estimate

- The misID Bdâ†’DÏ€ shape is also reweighted to take misID curve vs momentum into account.
- The signal width is fixed to that found in the Bdâ†’DÏ€ fit, scaled by the ratio of widths for Bsâ†’DsÏ€ and Bdâ†’DÏ€ in the MC
- Signal mean and comb background slope are floating.
- A Î›bâ†’ Î›cÏ€component was allowed in the fit, but got fitted to zero.
- Some Bdâ†’DsÏ€ can also be seen
- re-use Bsâ†’DsÏ€ mass shape, and constrain yield to {known BR ratio*fd/fs} = 1/35 relative to Bsâ†’DsÏ€ yield.

For Approval

For Approval

- The amount of misID Bsâ†’DsÏ€ background is floated
- Provides x-check on misID rate estimate
- Any Bdâ†’DÏ€should be taken care of by the Bsâ†’DsÏ€ shape

- Treatment of signal shape is the same as for Bsâ†’DsÏ€
- Comb background slope is fixed to be flat (from wrong-sign)

- Amount of Bdâ†’DKisconstrained from the Bdâ†’DÏ€ under DsÏ€, using the Bdâ†’DK/Bdâ†’DÏ€BR ratio
- Relative yields of PartReco backgrounds are constrained using
- Relative reconstruction efficiencies (from MC) when e.g. a charged track or soft pion/photon is missed
- Bs branching ratios from Bdbranching ratios, using SU(3) symmetry
- The yields can ove by 33% from these estimates

- TheBdâ†’DsKyield is floated. The Bsâ†’DsKand Bdâ†’DsKshapes are the same.

- Amount of Î›bâ†’Ds(*)p is constrained as detailed earlier.

For Approval

For Approval

For Approval

For Approval

For Approval

For Approval

- While this component is clearly visible in the DsK fits, the amount of background underneath it makes a reliable fit to its yield very difficult, at least with the current dataset.
- Hence we cannot make a competitive measurement of its BR (error in PDG is ~13%).

- Ratio of trigger/stripping/(non-PID) selection efficiencies from MC
- Fit model systematics will be evaluated using a large number of toy fits (as was done for fd/fs analysis).
- But for the moment, we simply apply cross-checks on the data, and assign conservative systematics.
- For the PID, take systematic on efficiency curves quoted by the RICH group, evaluated at our cut values
- PID systematic can enter in three different ways:
- Final PID efficiency correction to obtain BR(DsK)/BR(DsÏ€)
- Shape of misID bkgs after reweighting
- Expected number of DÏ€/K under DsÏ€/K (constrained in the fit)

- The fit model systematic for DsKis the most involved part of the systematics calculation.
- The main contributions to this part are:
- The slope of the combinatoric is fixed to half of the DsÏ€ slope. This reduces the signal yield by 3%.
- The constraints on the partially reconstructed backgrounds are all varied by a factor of two. This changes the signal yield by Â±4%.
- Also, the ratio of the Î›bâ†’Ds*p component to the Î›bâ†’Dsp component was varied by a factor of two. The change to the signal yield is <0.5%.

- Averaging MagUp and MagDown, we get
- N(DsK) = 406Â±26, N(DsÏ€) = 6038Â±105
- ÎµPID(DsK) = 83.4Â±0.2%, ÎµPID(DsÏ€) = 85.0Â±0.2%
- Îµsel(DsÏ€)/ Îµsel(DsK) = 0.945Â± 0.014
- We obtain

- Basically we turn the 2010 fd/fs combination on its head, by combining the ratio of yields of Bsâ†’DsÏ€ and Bdâ†’DÏ€ from the hadronic analysis, and the fd/fs value from the semileptonic analysis

Input:

Output:

- Combining these two results, we obtain

- This agrees with the Belle result, but is significantly below the CDF result.

- Using 2010 data we measure
- With 336pb-1 of 2011 data we measure
- These are combined to yield
- These are all Worldâ€™s Best measurements.
- Last but not least, we would like to thank our referees, StefaniaVecchi and StephaneMonteil, for their quick work which has been very helpful in improving our analysis!

Strong phase difference

- The MVA was trained for the fd/fs analysis using a small (2pb-1) subsample of the 2010 data
- Several MVAs were tried, the Boosted Decision Tree with Gradient boosting (BDTG) was found to have the best performance

PID variable (uses de/dx)

These are for the bachelor momentum spectrum, after the p<100GeV has been applied.

Shape for Bdâ†’D*-Ï€+ from 2010 MC, under pion (left) and kaon (right) mass hypothesis for the bachelor

- As a side-product of the 2010 fd/fsanalysis, we measured:

- Whereas we now measure :
- Bdâ†’DKHas only one tree diagram, while the Bsâ†’DsKhas two
- So our result suggests that the two different tree diagrams contributing to the DsK final state interfere destructively