300 likes | 438 Views
HBD performance study during run9 pp 200GeV. Katsuro Nakamura. HBD Sdphi , Sdz stability check during Run-9. Sdphi , Sdz … distance between track and cluster, normalized with RMS used whole electron tracks in EWG_DST on CCJ. electron track selection. abs(bbcz)<20cm quality==31||51||63
E N D
HBD performance study during run9 pp 200GeV Katsuro Nakamura
HBD Sdphi, Sdzstability check during Run-9 • Sdphi, Sdz … distance between track and cluster, normalized with RMS • used whole electron tracks in EWG_DST on CCJ
electron track selection • abs(bbcz)<20cm • quality==31||51||63 • pt>0.15GeV/c • n0>=2 • prob>0.01 • abs(emcsdphi_e)<4, abs(emcsdz_e)<4 • e/p>0.6 • hbdsize>=2
Sdphi fitting • fit with gaus + pol1 • check the run dependence of mean and sigma of gaus
Sdphi 0.0 < pt < 0.5 GeV/c 0.5 < pt < 1.0 GeV/c mean mean sigma sigma
Sdphi 1.0 < pt < 1.5 GeV/c 1.5 < pt < 2.0 GeV/c mean mean sigma sigma
Sdz fitting • fit with gaus + pol1 • check the run dependence of mean and sigma of gaus z
Sdz 0.0 < pt < 0.5 GeV/c 0.5 < pt < 1.0 GeV/c mean mean • displacement • of mean value • recalibrate this value • (but this displacement is enough small to be ignored.) sigma sigma dynamical moving during run-9 ?
Sdz after private recalibration 0.0 < pt < 0.5 GeV/c 0.5 < pt < 1.0 GeV/c mean mean new_sdz = (sdz - mean)/sigma sigma sigma
Sdz 1.0 < pt < 1.5 GeV/c 1.5 < pt < 2.0 GeV/c mean mean sigma sigma
Sdz after private recalibration 1.0 < pt < 1.5 GeV/c 1.5 < pt < 2.0 GeV/c mean mean sigma sigma
electron track selection • abs(bbcz)<20cm • quality==31||51||63 • pt>0.15GeV/c • n0>=2 • prob>0.01 • abs(emcsdphi_e)<4, abs(emcsdz_e)<4 • e/p>0.6 • hbdsize>=2 • abs(hbdsdphi)<4, abs(hbdsdz)<4
electron track distribution pad id (0-origin): 96 97 98 …………… Sect 9 South • electron track dist. on sect9 is almost uniform. fill hbdcharge to the nearest pad from p = {phbdx, phbdy, phbdz}. ///////////////////////////////// int padnum = NearestPad(phbdx,phbdy,phbdz); hist[padnum]->Fill(hbdcharge); /////////////////////////////////
I found the update of HbdWisClustrizer.C on CVS • HbdWisClustrizer.C • define clustering algorithm • HbdWisClustrizer.C –r1.12 (previously used) • cause nonuniform cluster peak distribution • detail is in slide sent to electron-l on 5/ 6/ 2010 • update • HbdWisClustrizer.C –r1.16 (latest version) • looks fine!
peak pad distribution after update pad id (0-origin): 96 97 98 …………… Sect 9 South • reasonable distribution fill hbdcharge to the peak pad in the cluster ////////////////////////// int max_pad = peak pad in the cluster; hist[max_pad]->Fill(hbdcharge); //////////////////////////
Dalitz Electron Pair Selection • Electron Selection • bbcz < 20 cm • quality==31||51||63 • pt>0.15GeV • n0>=2 • abs(emcsdphi)<4.0 , abs(emcsdz)<4.0 • prob>0.01 • e/p>0.6 • hbdsize>=2 • abs(hbdsdphi)<4.0 , abs(hbdsdz)<4.0 • Pair Selection • Mass<0.15GeV/c^2 select Dalitz region • PhiV cut reject conversion pairs • Unlike sign
Single e and Double e definition Single Double hbd_id(trk_i) == hbd_id(trk_j) trk_i trk_i • hbd_id(trk_i) != hbd_id(trk_j) trk_ j π π trk_ j notice: used HbdWisClusterizer.C –r1.12 in the following analysis
HBD charge distribution (East Arm)for Dalitz event (Raw dist.) Sect Sect blue: single red : double 0 3 1 4 5 2 Direction N N S S
HBD charge distribution (West Arm)for Dalitz event (Raw dist.) Sect Sect blue: single red : double 6 9 7 10 11 8 Direction S S N N
HBD charge distribution (East Arm)for Dalitz event (normalized) Sect Sect • normalized with integral of entry blue: single red : double slightly worse slightly worse 0 3 difficult to use 1 4 5 2 Direction N N S S
HBD charge distribution (West Arm)for Dalitz event (normalized) Sect Sect • normalized with integral of entry blue: single red : double 6 9 7 10 11 8 Direction S S N N
comparing with the distribution made by Ilia pair mass ~ 0 field ~ 0 • good agreement e HBD charge distribution e+ HBD charge distribution D e- π Dalitz event Dalitz event single double by I.Ravinovichi
For the next step • study Pt dependence of the ratio between single and double. • low pt 2 clusters tend to be separated • high pt 2 clusters tend to be merged • decide charge threshold sector-by-sector
HBD scintillation charge distribution • fit with a0*exp(-x/a1) • a0 = const, a1 = slope fit on 2.5 ~ 5 [p.e.]
run dependence of fitting slope • it’s stable • since it’s calibrated • no large difference on pad-by-pad • also “const” has no large difference HBD sect9 Yellow: pad107 Red: pad127 Blue: pad106 Green: pad105
central pad in cluster pad id (0-origin): 96 97 98 …………… Sect 9 South • you can see the pad dependence of yielded charge