310 likes | 414 Views
Measurement of the top quark properties with jet flavour and charge tagging Christophe Clément (CERN/ Stockholm University). 13 th Nordic LHC Physics Workshop Helsinki October 26, 2006. V tb unconstrained without 3x3 unitarity constraint. B( t →Wb ). B( t → Wb )=1 usually assumed by
E N D
Measurement of the top quark properties with jet flavour and charge taggingChristophe Clément (CERN/ Stockholm University) 13th Nordic LHC Physics Workshop Helsinki October 26, 2006
Vtb unconstrained without 3x3 unitarity constraint B(t→Wb) • B(t→Wb)=1 usually assumed by CDF and D analyses • B(t→Wb) might deviate from unity: • Additional quark singlets or doublets • Non-SM processes in the production • Non-SM in the decay (H+,...) • Experimentally B(t→Wb) affects number of b-jets need to experimentally discriminate b/w t→Wb and t→Wqlight Vtb October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 2/19
Missing transverse energy One high pT isolated lepton W t t b-jets? light jets? W light, c-jets Deriving B(t→Wb) experimentally.... • Select a top-enriched sample • e+jets and +jets channels • larges statistics, good S/B • # events with 0, 1 and ≥2 b-jets • B(t→Wb) • b-tagging efficiency • Jet identification efficiency • Probability to tag background October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 3/19
Lxy= How to identify b-quark jets... • Efficiency to tag a jet is broken down into • xplicitely reconstruct displaced secondary vertices: Secondary Vertex Tagger • Taggability: • Efficiency to match a calorimeter jet • to a track jet*, sensitive to calorimeter noise, • tracking efficiency. • 2. Tagging efficiency • Probability for a taggable jet to contain • a SV with Lxy > n σLxy • B-tagging efficiency decoupled from • calorimeter noise, tracking efficiency • Similar algorithm used by CDF • * Track jet = ≥ 2 tracks in ΔR=0.5, Δz<2cm, ≥ 1 hit in the silicon detector, pT>0.5 GeV, ≥ 1 track with pT>1 GeV Taggability, tagger independent October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 4/19
b-tagging efficiency from data ”n-sample”: 1 jet with a Not -tagged n -tagged n Not -tagged, SVT tagged nSVT -tagged, SVT tagged n ,SVT ”p-sample”: 2 b-2-b jets 1 jet with a -tagged, SVT tagged p,SVT -tagged p Not -tagged, SVT tagged pSVT Not -tagged p October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 5/19
b-tagging efficiency From di-jet data: extract b-tagging efficiency for muonic b-jets We need the b-tagging efficiency for ”all kinds of b-jets” Light jet tagging rate derived from data as well... b→, data b Similar techniques described in Phys. Rev. D71, 052003 (2005) b,MC Pbtag(ET,) = Taggability CorrectionFlavor(b) --------b→, data b→,MC Transform semi-muonic b-tag efficiency into inclusive one October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 6/19
Lepton + jets sample composition • 2) Other small backgrounds • derived from MC • Lepton, jet efficiencies • calibrated on data • σfrom data or NLO • Nother True isolated lepton processes Ntrue • . • W+jets • Z+jets • WW, WZ, ZZ • single top • multijet fake isolated lepton processes Nfake • Fake isolated electron • Jets with leading /πo, convertions, γwith random tracks,... • Fake isolated • inside jets from heavy flavor or in flight decays • I) Determine Ntrue , Nfake on a statistical basis • from data by deriving the probability for jets to fake isolated leptons 3) Fit B(t→Wb), Ntt, Nwj to the Nn-tags Nbefore tag= Ntt + NWj + Nfake + Nother Nn-tags = Pntt (B(t→Wb)) Nntt + PnWj Nnwj + N’nfake + Pnother Nnother n-tags = 0, 1, 2 October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 7/19
≥4 jets 3 jets ≥4 jets, 0 tag Number of b-tagged jets • Fit B(t→Wb) and Ntt simultaneously to 12 bins: e/+ 3/≥4 jets, 0/1/≥2 tags • B(t→Wb) is constrained by relative 0/1/≥2 tags populations • Very poor S/B in 0-tag sample, Nttbar~√Nobs • Perform measurement of Ntt in 0-tag as additional constrain. Number of b-tagged jets Topological discriminant October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 8/19
Confidence contour plots in R, Ntt No 0-tag B(t→Wb) B(t→Wb) Phys. Lett. B 639 (2006) B(t→Wb) = 1.03+0.19-0.17 68% CL : B(t→Wb)>0.78 |Vtb|>0.88 95% CL : B(t→Wb)>0.61 |Vtb|>0.78 October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 9/19
OR ? Mtop~270GeV +2/3 -1/3 -4/3 t Q1 mixing b Q4 M(Q4)~175GeV Top Quark Charge Lift ambiguity present in all top analyses! t→W+b or ”t”→W-b Test exotic models... Phys.Rev. D65 (2002) 053002 October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 10/19
Double tagged events S/B~10 OR ? kinematic fit Qtop + = Ingredients Discriminate b/w |Qtop| = 2e/3 and |Q”top”|=4e/3 Two |Qtop| per event Compare to expected shape of observable Qtop for • SM top • ”top” with 4e/3 charge? October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 11/19
Algorithm: Jet Charge Algorithm pTi,qi Compute jet charge only for b-tagged jets • (2 per events) Jet charge = Optimizaton on MC gives a=0.6 Sum over tracks with • pT>0.5GeV, • ΔR(track, jet) <0.5 of the jet axis Derive expected shape of Qjet from data with minimal input from simulation October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 12/19
Tight di-jet sample Jet Charge Performance in Data • Tag and probe methodin ”pure” events In reality: • Is it pure ? ? flavor excitation? g→ ? • B→ • B →D → • → • B→ light hadrons → • Charge misidentification _ bb >3.0 Ideal case: sign of q = sign of qb _ bb _ cc _ bb ´ _ Bo→ Bo Charge flipping processes October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 13/19
- Is the triple tag sample pure ”bb”? • The fraction of c-jets in the triple tag sample is determined by pTrel fit of the order of a few percents, • Flavor excitation/ splitting? >3.0 Flavor creation dominate for large R.D. Field, Phys. Rev D 65, 094006 (2002) CDF, Phys. Rev. D 71, 092001 (2005) D.A. Wijngaarden D PhD thesis, unpublished Vary 2ndlargest systematic error... October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 14/19
Tight di-jet sample - p.d.f.’s Qb, Qb, Qc, Qc from data... - - - P+ (Qjet) = (1-xc) (1-xflip)Pb (Qjet)+ (1-xc)xflip Pb + xc Pc _ cc Fraction of charge flipping processes 30±1% from MC, Cross checked on data Fraction of derived from pTrel spectrum of (1+2-1%) p.d.f of Qjet in probe jet + Similar equation for P- (Qjet) 4 Unknown p.d.f’s Pb, Pb, Pc, Pc - - October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 15/19
- - P+ (Qjet) = 0.69 Pb (Qjet) + 0.30 Pb + 0.01 Pc P- (Qjet) = 0.30 Pb (Qjet) + 0.69 Pb + 0.01 Pc P´+(Qjet) = 0.567 Pb (Qjet) + 0.243 Pb + 0.19 Pc P´- (Qjet) = 0.243 Pb (Qjet) + 0.243 Pb + 0.19 Pc - - - - - p.d.f.’s Qb, Qb, Qc, Qcfrom data... - Tight di-jet sample loose di-jet sample • Correct for different • t→Wb and bb kinematics • (largest systematic error) October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 16/19
qB qB qB qB qb qb qb qb SM Top and Exotic Quark Charge Observables • Use kinematic fit to assign b-jets to correct W-bosons in MC qb = qjet of the b-jet associated to the leptonic side of the event qB = qjet of the b-jet asssociated to the hadronic side of the event October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 17/19
Admixture of Q4 and t-quark - not excluded by - fQ4 <0.52 at 68% C.L. σtt Result on ~0.4 fb-1 of D data • 21 l + jets double tagged events • 16 events with converged kinematic fit 32 measured top charges • Perform likelihood ratio test b/w 4e/3 and 2e/3 hypothesis Exclude 4e/3 at 92%C.L. (91 % expected) hep-ex/060844 October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 18/19
Phys. Lett. B 639 (2006) B(t→Wb) = 1.03+0.19-0.17 Exclude 4e/3 at 92%C.L. (91 % expected) hep-ex/060844 Conclusions • World best/only measurement • Acknowledge collaborators: J. Strandberg (Stockholm University, now at Michigan) and P. Hansson (Kungliga Tekniska Högskolan) • Calibration of b-tagging algorithm and jet charge algorithm on data was crucial to establish these measurements • Calibration of efficiencies, fake rates, Monte Carlo,... with DATA: often very importantat the heart of analyses at Tevatron and most tricky part. • Probably even more so at LHC experiments, large leap forward, need to understand data before we can trust our results October 26, 2006 C.Clément - Helsinki Nordic LHC Physics Workshop 19/19
All hadronic b t W q t W q q’ Lepton+jets q’ b b Dilepton t W q b t W l t q’ W l t W l v b v v b Top Pair Final States • How top quark was discovered! • Used for most measured top quark properties so far ”golden channel” e+jets & +jets ~32% 44% ee+e+ ~5%
In SM top decays via V-A charged current ⇒ Mostly left handed b-quarks in the decay • 3 quark generations + direct measurements of Vub and Vcb predict Vtb,~1 B(t→Wb)~1 • SM predicts FCNC decays are tiny, t→Wq is dominant Top Quark Decay
R=1.0, σ=7pb R=0.5, σ=7pb Illustration: How data constraints B(t→Wb) and σtt ... Single Tag l+jets Jet multiplicity Double Tag l+jets Jet multiplicity
ttbar → l+jets Take into account ”contamination” by → ll Pntt versus B(t→Wb) • Probability to see n-tags in events (Pntt ) depends on the number of b-jets • Pn tt = R2 Pn tag(tt→bb) + 2 R(1-R)Pn tag(tt→bql) + (1-R)2Pn tag(tt→qlql)
Prediction 7pb & B(t→Wb)=1 • Fit B(t→Wb) and Ntt simultaneously to 12 bins: e/+ 3/≥4 jets, 0/1/≥2 tags • B(t→Wb) is constrained by relative 0/1/≥2 tags populations • Very poor S/B in 0-tag sample, Nttbar~√Nobs • Perform measurement of Ntt in 0-tag as additional constrain.
Lower limit on B(t→Wb) and |Vtb| • Prior π(B(t→Wb))=0 outside [0,1] • Monte Carlo integration over 191 nuisance parameters associated to systematic errors • Provides a 2D p.d.f. for B(t→Wb) and Ntt • Limit on |Vtb| can be derived using |Vtb|=√B(t→Wb) (SM) 68% CL : B(t→Wb)>0.78 |Vtb|>0.88 95% CL : B(t→Wb)>0.61 |Vtb|>0.78
Systematics on template shapes Some systematic uncertainties can affect the template shapes... • Systematics on template • shapes on ttbar→ l+jet • JES • JetID • Jet energy resolution • W-modeling (W+jets) • Taggability • Tagging probabilities for b, c • and light jets October 16, 2006 C.Clément CERN EP seminar 21/39
Likelihood discriminant in 0-tag sample l + 4 jets before tagging 1. Sphericity S = 3(λ2+λ3)/2 , λ’s smallest eigenvalues of momentum tensor M (ttbar S~1) 2. K’Tmin K’Tmin = Rminjj/EWT with EWT = ElT + MET 3. Centrality C=HT/H , HT is scalar sum of jets ET and H is the sum of the jet energies. 4. H’T2 H’T2 = HT2/Hz, HT2 : pT sum of all jets but leading jet, Hz is the scalar sum of all jets |Ez| plus |Ez| of the neutrino (W-assumption) October 16, 2006 C.Clément CERN EP seminar 19/39
Discriminant Power Data Calibration Corrections MC truth on tag side Tag and probe Method: data Tag and probe Method: Z→bb October 16, 2006 C.Clément CERN EP seminar 31/39
Charge of highest pT hadron Charge of b-quark Why does it work? The charge of the quark is correlated with the charge of the highest pT hadron resulting of the hadronization Simple study carried out with Pythia: Generated QCD 2→2 process, pT>15GeV And look at hight pT b quarks produced in the process. Usually large number of hadrons Produced, most quite modest pT Thís is then smeared by detector effects... The original b-quark not always In highest pT hadron October 16, 2006 C.Clément CERN EP seminar 29/39
Systematic uncertainties Jet charge specific Kinematic fit October 16, 2006 C.Clément CERN EP seminar 38/39