Study of kinematic fit for Higgs->inv. Yu Kato The Univ. of - PowerPoint PPT Presentation

Study of kinematic fit for Higgs->inv. Yu Kato The Univ. of Tokyo Asian Physics and Software Meeting Nov. 28, 2017

2 fitter engine : NewtonFitter ( use Newton-Raphson iterative method ) Study of Kinematic Fit for Higgs->inv. 2017/12/1 kinematic fit (make my own kinematic fitter based on Kurata-san’s kinematic fitter) � 3. consider ISR effect ( e.g. include aa_lowpt overlay) 2. apply BW mass constraint & Z natural width 1. samples : ffH_ZZ_4n { ILD_l4_v02, ILD_s4_v02, (ILD_o1_v05) } p We want to improve recoil mass resolution of Higgs->inv. analysis. 𝑓 " 𝑓 # → 𝑎𝐼, 𝑎 → 𝑟𝑟, 𝐼 → 𝑗𝑜𝑤𝑗𝑡𝑗𝑐𝑚𝑓 − 𝐹 8 3 − 𝑞 3 ⃗ 8 3 𝑁 012 = 𝑡 p Apply kinematic fit minimize χ 2 under constraints => method of Lagrange multipliers 𝜓 3 = 𝜃 C 𝑊 "E 𝜃 ⃗ C I 𝑔 ⃗ ⃗ =>? − 𝜃 ⃗ @1AB ⃗ =>? − 𝜃 ⃗ @1AB + 2𝜇 • measured parameters 𝜃 ⃗ : 𝐹 KE , 𝜄 ME , 𝜚 KE , 𝐹 K3 , 𝜄 M3 , 𝜚 K3 ⃗ : 𝑛 KE,K3 = 𝑛 8 = 91.2 GeV , 𝛾 KE = |𝑞 KE |/𝐹 • constraint function 𝑔 KE , 𝛾 K3 = |𝑞 K3 |/𝐹 K3 • covariance matrix 𝑊 : 𝜏 [ (𝐹, 𝑑𝑝𝑡𝜄), 𝜏 ` (𝐹, 𝑑𝑝𝑡𝜄), 𝜏 a (𝐹, 𝑑𝑝𝑡𝜄) p Use MarlinKinfit p Status • include jet resolution for each event ( from Z->uds 2jet sample data ) • recoil mass sigma becomes smaller a little p Plan

Jet resolution study ILD_l4_v02 2017/12/1 Study of Kinematic Fit for Higgs->inv. 3 sv01-19-04_lcgeo.mILD_l4_v02 sv01-19-04_lcgeo.mILD_l4_v02 15 15GeV 15GeV ) [%] MC 20GeV 20GeV θ 0.06 - 30GeV 30GeV j (E REC 45.5GeV 45.5GeV 90 10 60GeV 60GeV θ ) / Mean 0.04 80GeV 80GeV = 100GeV 100GeV θ δ 120GeV 120GeV 5 150GeV 150GeV 0.02 j (E 175GeV 175GeV 90 200GeV 200GeV RMS 250GeV 250GeV 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 θ θ |cos | |cos | sv01-19-04_lcgeo.mILD_l4_v02 sv01-19-04_lcgeo.mILD_l4_v02 0.3 0.06 15GeV 15GeV MC MC 20GeV 20GeV θ φ 0.05 0.25 - cos - 30GeV 30GeV REC 45.5GeV 45.5GeV 0.2 0.04 60GeV 60GeV φ REC 80GeV 80GeV = 0.15 0.03 100GeV 100GeV θ φ = cos δ 120GeV 120GeV 0.1 0.02 150GeV 150GeV 175GeV 175GeV 0.05 θ 0.01 200GeV 200GeV cos 250GeV 250GeV 0 0 δ 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 θ θ |cos | |cos |

Jet resolution study ILD_s4_v02 Study of Kinematic Fit for Higgs->inv. 4 2017/12/1 sv01-19-04_lcgeo.mILD_s4_v02 sv01-19-04_lcgeo.mILD_s4_v02 15 15GeV 15GeV ) [%] MC 20GeV 20GeV θ 0.06 30GeV - 30GeV j (E 45.5GeV REC 45.5GeV 90 10 60GeV 60GeV ) / Mean θ 80GeV 0.04 80GeV = 100GeV 100GeV θ δ 120GeV 120GeV 5 150GeV 150GeV 0.02 j (E 175GeV 175GeV 90 200GeV 200GeV RMS 250GeV 250GeV 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 θ θ |cos | |cos | sv01-19-04_lcgeo.mILD_s4_v02 sv01-19-04_lcgeo.mILD_s4_v02 0.3 0.06 15GeV 15GeV MC MC 20GeV 20GeV φ θ 0.25 0.05 - - cos 30GeV 30GeV REC 45.5GeV 45.5GeV 0.2 0.04 60GeV 60GeV φ REC 80GeV 80GeV = 0.15 0.03 100GeV 100GeV θ φ = cos δ 120GeV 120GeV 0.1 0.02 150GeV 150GeV 175GeV 175GeV θ 0.05 0.01 200GeV 200GeV cos 250GeV 250GeV 0 0 δ 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 θ θ |cos | |cos |

MarlinKinfit result : jet beta constraint ILD_l4_v02 Study of Kinematic Fit for Higgs->inv. 5 right pol. 2017/12/1 left pol. ILD_o1_v05 ILD_s4_v02 sv01-19-04_lcgeo.mILD_l4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_l4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) Events / 2.00 GeV Events / 2.00 GeV Events / 2.00 GeV 2000 before fit: before fit: 2500 before fit: 2000 mean = 130.4 mean = 130.8 mean = 130.3 sigma = 12.022 sigma = 12.062 sigma = 11.669 2000 1500 1500 after fit: after fit: after fit: mean = 129.6 mean = 129.6 mean = 129.2 1500 sigma = 10.051 1000 sigma = 10.110 sigma = 10.099 1000 1000 500 500 500 0 0 0 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 Recoil Mass [GeV] Recoil Mass [GeV] Recoil Mass [GeV] sv01-19-04_lcgeo.mILD_l4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_l4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) Events / 2.00 GeV Events / 2.00 GeV Events / 2.00 GeV 2000 before fit: before fit: before fit: 2000 2500 mean = 130.8 mean = 131.1 mean = 130.2 sigma = 12.223 sigma = 12.368 sigma = 11.424 2000 1500 1500 after fit: after fit: after fit: mean = 129.5 mean = 129.4 mean = 129.1 1500 sigma = 10.051 sigma = 10.432 1000 sigma = 9.978 1000 1000 500 500 500 0 0 0 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 Recoil Mass [GeV] Recoil Mass [GeV] Recoil Mass [GeV] Cut : no error fit && 𝜓 3 < 500

MarlinKinfit result : jet massless constraint Study of Kinematic Fit for Higgs->inv. 6 ILD_l4_v02 ILD_s4_v02 ILD_o1_v05 right pol. left pol. 2017/12/1 sv01-19-04_lcgeo.mILD_l4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_l4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eL.pR ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) 2000 Events / 2.00 GeV 2500 Events / 2.00 GeV Events / 2.00 GeV before fit: before fit: before fit: mean = 130.4 mean = 130.8 mean = 130.3 2000 1500 sigma = 12.021 1500 sigma = 12.062 sigma = 11.667 after fit: after fit: after fit: 1500 mean = 127.9 mean = 128.0 mean = 127.6 1000 1000 sigma = 10.726 sigma = 10.827 sigma = 10.822 1000 500 500 500 0 0 0 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 Recoil Mass [GeV] Recoil Mass [GeV] Recoil Mass [GeV] sv01-19-04_lcgeo.mILD_l4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_l4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_s4_v02.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) sv01-19-04_lcgeo.mILD_o1_v05.eR.pL ((mkfprob > -0.5)&&(mkfchi2 < 500))&&(mkfmrec > 0) Events / 2.00 GeV Events / 2.00 GeV Events / 2.00 GeV 2500 2000 before fit: before fit: before fit: 2000 mean = 130.8 mean = 131.1 mean = 130.2 2000 sigma = 12.224 sigma = 12.368 sigma = 11.424 1500 1500 after fit: after fit: after fit: 1500 mean = 127.9 mean = 127.8 mean = 127.5 1000 sigma = 10.701 sigma = 11.160 sigma = 10.616 1000 1000 500 500 500 0 0 0 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 Recoil Mass [GeV] Recoil Mass [GeV] Recoil Mass [GeV] Cut : no error fit && 𝜓 3 < 500

backup 2017/12/1 Study of Kinematic Fit for Higgs->inv. 7

2017/12/1 Study of Kinematic Fit for Higgs->inv. 8