1 / 19

Study of B s  mmf fK + K -

Study of B s  mmf fK + K -. Supreet Pal Singh Prof. Kajari Mazumdar Prof. J.B.Singh. Introduction to Bs  mmf. The branching ratio is predicted to be 1.6 x 10 -6 in Standard Model. source : C. Q. Geng and C. C. Liu, J. Phys. G 29, 1103 (2003).

lharms
Download Presentation

Study of B s  mmf fK + K -

An Image/Link below is provided (as is) to download presentation 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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Study of BsmmffK+ K- Supreet Pal Singh Prof. KajariMazumdar Prof. J.B.Singh India LHC weekely meeting (6th July 2009)

  2. Introduction to Bs mmf • The branching ratio is predicted to be 1.6 x 10-6 in Standard Model. source : C. Q. Geng and C. C. Liu, J. Phys. G 29, 1103 (2003). • The latest update on the branching ratio is an upper limit from D0 experiment: B( Bs mmf) < 3.2 x 10-6 • s(ppbb) ~ 440 mb at 10 TeV. • Prob(b Bs) = 10% So ,for 100 pb-1 luminosity we expect about 4400 events of the type Bs mmfto be produced at 10 TeV pp collision. • A fraction of these will be within the detector acceptance. • But this number will further reduce due to non-zero trigger thresholds on final state particles, ie, muons. India LHC weekely meeting (6th July 2009)

  3. General strategy Bs mass = 5.369 GeV • Trigger events with 2 muons only. m+ m+ • Bs decays  Secondary vertex properties utilized during offline analysis in this study. • (displaced vertex trigger may be implemented after initial phase) • Combine 2 muons + 2 energetic tracks from f (fK+ K- , BR ~ 50 %) • Since CMS does not have particle identification assign kaon mass to hard tracks other than muons and then try to form the f resonance. • Most of the time the other b-quark in the event is expected to produce jet in the opposite hemisphere. • B-tagging is not used in the present study though. . B0s . SV K+ PV K- India LHC weekely meeting (6th July 2009)

  4. CMS software, data samples etc. for this work • The analysis has been done using PAT in CMSSW_2_2_6 via CRAB. • Signal events have been generated privately through full reconstruction chain in CMSSW_2_2_3 on the line of official data production. • (Tier2 facility at Uni. Wisconsin has been used mainly.) • Pythia for generic QCD, (MSEL =1, no thresholds on hard scattered partons)  bb pair produced due to gluon fusion, flavour excitation and gluon splitting processes. • EvtGenLHC generator has been used to generate the decay mode with standard format for user_decay_file : “GeneratorInterface/EvtGenInterface/data/Bs_phimumu_KK.dec” • Requirement at gen. level: • Presence of Bs meson • 2 muons in the event with PT > 2.5 GeV and |h| < 2.5 •  Generator efficiency = 1.73 x 10 -4 No. of signal events analyzed = 8665  Integrated Luminosity = 1.0 fb-1 India LHC weekely meeting (6th July 2009)

  5. More on Signal and Background events No. of signal events analyzed = 8665 • Integrated Luminosity = 1.0 fb-1 • Average momentum of Bs ~ 10 GeV • The daughter particles in 3-body decay are quite soft! • Not a convenient situation for analysis! PT vs. h plot for Bs mesons Main source of background QCD Background samples analyzed : InclusivePPmuX/Summer08_IDEAL_V9_v4/GEN-SIM-RECO For this dataset: MSEL=1 At least one m in the event with PT > 2.5 GeV and |h| < 2.5 Integrated Luminosity = 0.04 pb-1  Large scale factor needed to extrapolate results to 10 pb-1. India LHC weekely meeting (6th July 2009)

  6. Kinematic variables from signal sample at generator level Distributions are normalized to unity m m f f K K India LHC weekely meeting (6th July 2009)

  7. PT and h of the leading reconstructed muons Events are normalized to unity Leading Muon 2nd Leading muon Events lost in the central region India LHC weekely meeting (6th July 2009)

  8. Analysis Criteria and Candidate Selection • Two opposite charged leading pair of muons are selected. • Trigger condition not applied, but same cuts at offline on recoedmuons : PTm > 3.0 GeV & |hm|< 2.4 • All track pairs of opposite charges with PTtk > 0.4 GeV and lying in the hemisphere opposite to the leading jet (|htk|. |hjet|< 0) are chosen to form (K+K- system ) invariant mass by assigning kaon mass. • There are many combinations possible, (for background events)! The pair closest to f mass (1.005 < Mf <1.035 GeV) are taken as candidate kaons from Bs decay. • Events are studied in sideband region of the signal mass window (5.245 GeV – 5.491 GeV) to avoid biases due to peaking backgrounds. B-Physics Meeting 2nd June 2009

  9. formation of Bs in signal sample Invariant Mass Distributions formation of f Huge dominance of QCD over signal India LHC weekely meeting (6th July 2009)

  10. Analysis Strategy • From invariant mass distribution of dimuons, peaking backgrounds of J/ and Ψ(2S) are vetoed out. • The events for the background are then studied in the sideband region of signal mass region (5.369 ± 2.5 s) so as to reduce bias and double counting of peaking backgrounds. • The four tracks (m±,K±) are expected from one mother Bs  kinematically fitted to a common vertex using Kalman Vertex Fitter technique . • Corresponding decay length of the secondary vertex from the primary vertex in the XY-plane along with χ2 are required to satisfy quality criteria. • The PT of the reconstructed shd. not be too soft  Pt(Bs ) > 11 GeV • muons and f should be close by  ΔR between mm and f < 0.4 • Around Bs candidate do not expect too much extra activity isolation parameter using tracker information: India LHC weekely meeting (6th July 2009)

  11. Preliminary Cut 1.4 GeV < Mmm< 4.4 GeV (excl. J/ and Ψ(2S) resonances) India LHC weekely meeting (6th July 2009)

  12. QCD events: Invariant mass distribution of (m m k k) system in the relevant mass range Signal mass window: 5.245 to 5.491 GeV Integral fraction: Nev in signal window/Nev in full range = 0.132 We proceed to estimate the number of background events within signal mass window by first estimating from events within the whole mass range (excluding the Bs mass region), the efficiencies for subsequent selections.  Side band estimation. India LHC weekely meeting (6th July 2009)

  13. PT distribution of Bs candidate & ΔR distribution between mm and f system PT of the Bs candidate > 11 GeV ΔRmmf < 0.4 India LHC weekely meeting (6th July 2009)

  14. Displaced Vertex Tracker Isolation χ2 < 2.0 I > 0.78 lxy < 0.1 Decay Length (mmKK) < 0.1 India LHC weekely meeting (6th July 2009)

  15. Un-correlation between PT(Bs), I, ΔRmmf and displaced vertex parameters These are uncorrelated to one another  the cuts can be applied Factorially. Total efficiency is the Product of individual Efficiencies. India LHC weekely meeting (6th July 2009)

  16. The relative efficiency at each step of event selection for Signal and QCD background We got various efficiencies, Number of signal and background events expected to survive in the signal mass window (5.245-5.391 GeV) after all the cuts for a given luminosity can be calculated now. India LHC weekely meeting (6th July 2009)

  17. Number of Signal and Background events for 100 pb-1 • No. of signal events left after surviving all the cuts: Nsig = Ns . ЄA . ЄB . ЄC . ЄD where Ns = the number of signal events surviving in the signal mass region (5.245-5.492 GeV) after preliminary cuts = 391 Nsig = 7.95 for 1 fb-1 Total cumulative signal efficiency, εs = 7.95/8665  9.2 x 10-4 . For 100 pb-1 number of signal events expected to survive: ns = 0.795 • Total number of background events expected after applying all the selections for 100 pb-1 luminosity : nb = 4.75 India LHC weekely meeting (6th July 2009)

  18. Calculations for 100 pb-1 luminosity • The upper limit on branching fraction of the process Bsmmf BUL = NUL(90%) / (εs .egen . NBs) εs = signal efficiency = 9.2 x 10-4 εgen = Kinematic acceptance that the produced Bs meson decays into dimuon and f according to generator level cuts = 0.13 NBs = No. of Bs mesons expected to be produced for 100 pb-1 luminosity. NUL (90%) = upper limit on mean of Poisson variable n at 90% confidence level. • NUL is obtained using TRolke method in root package taking into account the number of background events (nb) and total number of observed events (ntot) using relation ntot = nb + ns . NUL = 4.94 Now, NBs = sbb . L . fs = 8.8 x 1011 where sbb = b production cross section = 440 mb L = integrated luminosity = 100 pb-1 fs = 0.1 (percentage of b-bbar to form Bs meson) India LHC weekely meeting (6th July 2009)

  19. Conclusion • Our analysis derived a signal selection efficiency of 9.2 x 10-4 • Since we considered data of accumulated lumi only 100 pb-1 and the SM expectation of Br is about 10-6 , again it is clear that we cannot observe a significant number of signal events with early data in CMS. • The upper limit on the branching fraction : B(Bsmmf ) <6.7 x 10-6 . • However the analysis method is established to search for this rare decay and with more data one expects to get better result. • Also at higher energy signal selection will be more efficient. India LHC weekely meeting (6th July 2009)

More Related