study of kinematic fit for higgs inv
play

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


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

  2. 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

  3. 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 |

  4. 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 |

  5. 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

  6. 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

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

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

Recommend


More recommend