Analyzing Day 74 ZDC data
Updated on Mon, 2009-03-16 17:03. Originally created by aliceb on 2009-03-16 14:18.Runs 10074004,5
New PID Asymmetries
Update: I coded up the statistical uncertainty calculation incorrectly in this post. I forgot to divide by the purity! The optimal cuts changed by a lot when I fixed that error. I also generalized the formula to account for the presence of signal in the sidebands, and I included the lowest p_{T} bin in the analysis. For more details see New PID Asymmetries II. Bottom Line — the results on this page are wrong!
In an earlier post I explained the new method for calculating identified particle yields that I’m using in my A_{LL} analysis. I began that study because I planned to calculate A_{LL} differently than I had been in the past. Specifically, I wanted to incorporate the proton/kaon/electron backgrounds into the statistical uncertainty instead of assigning a separate (statistics-limited) systematic uncertainty to account for their presence. I’m using the following formulas for A_{LL} and its statistical uncertainty:
where the p_{T}-dependent background fractions are defined as:
I wrote a small function to estimate the statistical precision on A_{LL} given the p_{T} bin, pion acceptance window, and sideband acceptance windows. I didn’t care about the absolute statistical precision, so I just used 1/sqrt(N) for the uncertainty on each A_{LL}. I used Minuit2 to minimize this function and extract the optimal acceptance windows, with the constraint that the purity in each sideband is never below 90%. In principle, this approach would yield four momentum-dependent cuts. In practice, the p+K sideband cut and the left side of the pion acceptance window only had a small momentum dependence, so for the sake of simplicity I keep them fixed. I also choose to fix the other two cuts in each p_{T} bin instead of letting them vary with momentum. In the end I employ the following cuts
pT bin | π window | max p+K | min electron |
---|---|---|---|
3.18 - 4.56 | (-1.90, 2.40) | -1.90 | 2.40 |
4.56 - 6.32 | (-1.90, 2.25) | -1.90 | 2.50 |
6.32 - 8.80 | (-1.90, 2.00) | -1.90 | 2.60 |
8.80 - 12.84 | (-1.90, 1.50) | -1.90 | 2.60 |
These cuts are significantly wider than the (-1.0, 2.0) acceptance window I had been using in the past. Apparently the reduction in purity is more than offset by the extra efficiency.
The electron side of the acceptance window is interesting. As momentum increases the pion band moves closer to the electron band. As a result, we need to move the electron sideband cut further out to maintain the 90% purity. This cuts down on the electron background A_{LL} statistics. The minimizer compensates for that uncertainty by restricting the right side of the pion acceptance window and thus reducing the electron background fraction.
I compared the uncertainties obtained by the minimizer with the uncertainties from my old method (a flat (-1.0, 2.0) cut that does not subtract out the background asymmetries). It turns out that the uncertainties from the new method are actually smaller in every p_{T} bin. I haven’t calculated a systematic uncertainty for this method, but if there is one it will be far smaller than the systematic from the old method (~ background fraction * sigma of background A_{LL}). In other words, using the new method is a no-brainer.
Oh, and one plot, just because I think it’s pretty:
New PID Cuts
Updated on Sun, 2009-03-22 10:56. Originally created by kocolosk on 2009-03-16 09:09.Update:I determined the electron peak positions in the p < 3.00 bins by letting that parameter float instead of fixing it to the Bichsel value. The difference was about 0.75σ in the lowest bin. The χ2 for that bin is dramatically better. Some of the other higher momentum slices in the first bin still have large χ2 values, but for the moment I’m working under the assumption that the background fractions are basically correct. So that bin’s back in. I’ve updated the PDF with the new fits.
Historically, I’ve used separate triple-Gaussian fits for positively- and negatively-charged tracks (18 free parameters in total) to extract the pion, proton/kaon, and electron yields for my A_{LL} analyses. There are several problems with that approach
- To first order, dE/dx resolution is independent of particle species.
- The peak positions for pions, kaons, protons, and electrons are not charge-dependent.
- Those peak positions (or at least the separations between them) have been determined with a great deal of accuracy in other analyses.
- dE/dx scales with momentum, not p_{T}. At y = 1.0, p ~ 1.5*p_{T}. This difference was not being taken into account.
I set out to address those three points by redoing the fits in the manner employed by the lfspectra working group. I fill a 3D histogram of p_{T} vs. p vs. (nσ(π) + 6*track.charge()), and then fit individual xy slices of this histogram with a function comprised of 8 Gaussians. That function is subject to the following constraints:
- all Gaussian widths are identical
- π+ mean == π- mean (after subtracting imposed offset of 12)
- K+ mean == K- mean (ditto)
- proton mean == pbar mean (ditto)
- positron mean == electron mean (ditto)
- π mean - K mean is a fixed function of momentum
- π mean - p mean is a fixed function of momentum
- π mean - e mean is a fixed function of momentum
That leaves 24-14 = 10 free parameters for each fit. I received the momentum-dependent particle separations from Yichun Xu. Her work is documented in this NIM draft, and the actual values for the separation are posted at (columns are p, p/mass, and separation):
π - electron separations
π - kaon separations
π - proton separations
The results of the fits are attached as a PDF. The fits in the first p_{T} bin (2.00 - 3.18) are not good. Part of the reason is that Yichun’s analysis doesn’t go below 3.0 GeV/c, so I used a Bichsel parameterization for the particle separations instead of taking them from data. This bin is complicated because at the low momenta the protons are entering the 1/β^2 region. I think the sensible choice may be to drop the bin from the A_{LL} analysis.
The fits for the higher p_{T} bins are much better, particularly at mid-rapidity. I don’t really believe in the p/π and K/π ratios that come out, but in the end I think the only important thing for my analysis is that the (p+K)/π ratio is correct. I’ll explore the effect of these new fits on A_{LL} in a separate post.
Observation of hypertriton and anti-hypertriton at RHIC
Updated on Sat, 2009-03-28 09:18 by jhchen. Originally created by bedanga on 2009-03-16 07:18.Reference
Speaker : Jinhui Chen ( Kent State University )
Talk time : 11:45, Duration : 00:30
Draft of Junior's Talk 2009
Updated on Fri, 2009-03-20 16:00. Originally created by pagebs on 2009-03-15 02:40.Here is the draft of my talk for the Juiors meeting at the 2009 collaboration meeting.
04 Spectra from Problem PMT Channels
Updated on Sat, 2009-03-14 22:55. Originally created by trent on 2009-03-14 22:55. Under:We set the HV = 800 Volts for all channels which have been masked out, repaired, or otherwise had problems in 2008. I have examined the spectra for these channels and give pdfs with each of t
Analyzing Day 73 ZDC Runs
Updated on Sat, 2009-03-14 15:47. Originally created by aliceb on 2009-03-14 15:19.Run 10073029
BSMD Performance Table
Updated on Mon, 2009-04-13 22:51. Originally created by wleight on 2009-03-14 11:38.For plots, see here.
STAR Management Meeting
Updated on Thu, 2009-03-19 11:01. Originally created by nxu on 2009-03-14 09:50.BNL, 50-I-189, at 15:00 (GMT), duration : 01:00
EVO: Title: STAR Management Meeting // Community: STAR // Password: star
Phone Bridge: +1 626 395 2112 (CA) / +1 631 344 6100 (NY) // ID: 826668 // Password: 0952
Time:
* --===================================--*
Agenda:
1) Run9 Status - Bill / Bernd
2) QM09 STAR Talk & Posters status - Bedanga
3) TPC review committee charges
* --===================================--*
0) Attending: Bedanga, Jerome, Bill, Nu, Bernd, Jamie, Zhangbu, Olga,
1) Run Status: Prepare few slides for next Tuesday CAD planning meeting, Bernd present.
2) QM09 status: (I) Parallel Talks: UPC: ok // LFSpectra: ok // BulkCorr: ok // HF: NPE problematic; Upsilon ok, J/Psi dAu result problem; e-h correlation: statistics; FMS result // JetCorr: ridge ok; jet recostruction talkS, not clear; Maties talk not posted; FMS: correlations from both pp and dAu, needs details, simulation results // (II) Plenary talks: S Voloshin; G Wang; J Putschke //
Notes:
1) STAR Collaboration meeting: March 23-27, 2009, BNL
2) China-STAR ToF meeting: April 26-29, Hangzhou, China
3) STAR Regional meeting: May 6-9, 2009, Pusan, Korea
Updated QM talk
Updated on Fri, 2009-03-13 17:53. Originally created by bbiritz on 2009-03-13 17:53.Hi All,
Here is an updated version of my QM '09 talk – thanks for all the feedback.
To answer a question brought up during the HFT phone conference:
The STAR Forward GEM Tracker
Updated on Fri, 2009-03-13 16:12. Originally created by surrow on 2009-03-13 16:12.Presentatation at Frontier Detectors for Frontier Physics, May 24-30, 2009
http://agenda.infn.it/conferenceDisplay.py?confId=759
Rapidity distribution and centrality in dAu 2008 vs dAu 2003
Updated on Fri, 2009-03-13 15:28. Originally created by canson on 2009-03-13 12:39.The rapidity distribution of highest pt particles (pt > 1.5 GeV/c) for dAu2008 and dAu2003 appears to be
D0 efficiency with TOF PID.
Updated on Fri, 2009-03-13 12:14. Originally created by bouchet on 2009-03-13 12:14.mini production
Updated on Fri, 2009-03-13 12:15. Originally created by bouchet on 2009-03-13 11:57.Reference
Speaker : jonathan bouchet (ksu)
Talk time : 13:50, Duration : 00:01
Mini production
I ran AuAu Hijing events with new geometry (details from Gerrit are in the readme file) for:
TPC hits distortion
Updated on Fri, 2009-03-13 11:57. Originally created by bouchet on 2009-03-13 11:57.Reference
Speaker : jonathan bouchet (ksu)
Talk time : 13:30, Duration : 00:20
BSMD and BPRS runs 10066028-40
Updated on Mon, 2009-03-16 13:25. Originally created by aliceb on 2009-03-13 11:14.TCDPhase[ns] BSMD Delay4 Run
70 97 10066028
70 98 29
70 99 30
35 100 31
70 100 32
0 101 33
35 101 34
70 101 35
0 102 36
35 102 37
Some slides on ZDC SMD
Updated on Fri, 2009-03-13 11:30. Originally created by aliceb on 2009-03-13 10:28.Draft for RSC meeting
Plots for RSC meeting
Updated on Fri, 2009-03-13 03:18. Originally created by ogrebeny on 2009-03-13 03:18.HFT Software
Updated on Thu, 2009-03-12 21:58. Originally created by margetis on 2009-03-12 21:58.510 4867333, at 17:30 (GMT), duration : 00:00
Here
Time | Talk | Presenter |
---|---|---|
13:30 | TPC hits distortion ( 00:20 ) 1 file | jonathan bouchet (ksu) |
13:50 | mini production ( 00:01 ) 2 files | jonathan bouchet (ksu) |
13:51 | D0 efficiency with TOF PID. ( 00:20 ) 1 file | Yifei Zhang |