260 likes | 390 Views
HKS collaboration meeting @JLab 2008/5/16. Akihiko&Yoko @ Wedding Party 2008/5/10. Topics. Bug in REPLAY Blind analysis by Lulin. Bug in REPLAY. Multiplicity problem. Run# 56089, eventID 1954 3 events @ JLab ENGINE 2events @ Tohoku ENGINE. Ntuples are inconsistent
E N D
HKS collaboration meeting @JLab2008/5/16 Akihiko&Yoko @ Wedding Party 2008/5/10
Topics • Bug in REPLAY • Blind analysis by Lulin
Multiplicity problem • Run# 56089, eventID 1954 • 3 events @ JLab ENGINE • 2events @ Tohoku ENGINE • Ntuples are inconsistent parameter files and source codes should be checked
Bug in analyzer • Cut condition for ENGE momentum (depends on “default” matrix) • Different “default” matrix between USA and Japan • Bug depends on multiplicity of ENGE Event#1 Event#2 Accept Event#3 (out of acceptance) Event#4 Reject …
Old ntuple vs. New ntuple • Fixed bug and re-analyze +12.1% +14.7% +17.5% +16.6%
Old RAP(12Jan08) vs. New RAP(04Apr08) • Number of events …. • Coin ntuple Multiplicity (ENGE) • RAP Coincidence between HKS and ENGE
Simulation data for blind analysis • Focal plane resolution in sigma 86 mm (exf), 0.7 mrad (expfp), 210mm(eyf), 2.8 mrad (eypf) 110 mm (hxf), 0.5 mrad (hxpfp), 110mm(hyf), 0.7 mrad (hypf) • Fixed beam energy (1853.1 MeV) • NO energy loss @ target • # of event & S/N (Updated) • moderately (L : ~2150 counts, S : ~400 counts, 12LB g.s. :~600 counts) • Sieve Slit data is ready • Raster (0.25 cm × 0.25 cm) for CH2
Lulin’s blind analysis • Tune CH2 & 12C • Add peaks of 12C one by one (maximum 6 peaks) • Raster correction • only momentum of kaon • Weight in the final step of tuning • 1.0(L), 1.0(S), 5.0 (B.E.=11.43), 0.0(B.E.=13.63), 6.0(B.E.=16.70), 3.0(B.E.=20.35), 4.5(B.E.=23.04), 6.0(B.E.=23.70) Detailed will be shown in Lulin’s talk * In the following pictures …… GOD : answer of simulation Lulin-tune : result from Lulin’s tuning
Lulin-tune CH2(tuned event) Top : GOD, Bottom : Lulin-tune for All event Top : GOD. Bottom Lulin-tune RED : coin. , BLACK : B.G.
Lulin-tune CH2(NOT tuned background event) TOP : GOD MM BOTTOM : LulinMM TOP : Lulin MM vs. GOD MM BOTTOM : wide region
Each component (xpt,ypt,dp)Lulin vs. GOD hxpt expt eypt hypt edp hdp offset due to different central momentum
Lulin-tune C12(tuned event) Top : GOD, Bottom : Lulin-tune for All event Top : GOD. Bottom Lulin-tune RED : coin. , BLACK : B.G.
Oka-tune C12(NOT tuned background event) RED : coin. , BLACK : B.G.
How to estimate contamination • Fit each peak by gaussian + offset • Fitted counts : Integration of gaussian • Real counts : # of real events in the peak(3s) (simulation data can be identified) • Contamination [%] = ((Fitted counts) – (Real counts)) / (Fitted counts)×100
Step by step missing mass • Step 1 Initial matrix ( *.1110 ) • Step 2 After involving L and S (*.02121) • Step 3 After involving B.E.=11.43, 23.04, 23.70 (*.03307) • Step4 After involving B.E.=20.35 (*.03308) • Step5 After involving B.E.=16.70 (*.04022)
Missing mass of step1(Initial) Top : CH2 for All event Bottom : Coin.(RED), B.G.(BLACK) Top : C12 for All event Bottom : Coin.(RED), B.G.(BLACK)
Missing mass of step2(tune L,S) Top : CH2 for All event Bottom : Coin.(RED), B.G.(BLACK) Top : C12 for All event Bottom : Coin.(RED), B.G.(BLACK)
Missing mass of step3(tune L,S,2 major peaks) Top : CH2 for All event Bottom : Coin.(RED), B.G.(BLACK) Top : C12 for All event Bottom : Coin.(RED), B.G.(BLACK)
Missing mass of step4(tune L,S,2 major peaks, core#1) Top : CH2 for All event Bottom : Coin.(RED), B.G.(BLACK) Top : C12 for All event Bottom : Coin.(RED), B.G.(BLACK)
Missing mass of step5(tune L,S,2 major peaks, core#1,core#2) Top : CH2 for All event Bottom : Coin.(RED), B.G.(BLACK) Top : C12 for All event Bottom : Coin.(RED), B.G.(BLACK)
Summary • The trend is the same as Okayasu’s result • The linearity between reference masses (L and S) seems OK but the outside is deformed • Background contamination in 12C peaks • The background events that not used for tuning do not make peak • Accuracy of binding energy : < 0.4 MeV
Comments on 7LHe Lulin Spectrum Lulin Mass : 6716.71 MeV/c2 (atomic masses are used) Lulin Binging energy : -5.53 MeV/c2 Lulin 7He + L : 6722.2413 MeV/c2 Okayasu Spectrum Okayasu Mass : 6716.15 MeV/c2 (nuclear masses are used) Okayasu Binging energy : -5.07 MeV/c2 7He + L : 6721.20 MeV/c2 D = 0.46 MeV