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
Recommend
More recommend