Sort by:[Date]

First look at some Simulations

At the moment my main focus is the off line calibration of the barrel EMC using neutral pions from 2006.  My (very) rough plan of how to do this is as follows:

Some Luminosity Statistics

===========================================================================================
trigId L_int_mb ε_vtx ε_mb vz vz_mb L_samp_mb < ps > L_sampled
-------------------------------------------------------------------------------------------
117001 4.794 pb^-1 0.505 0.505 0.678 0.678 6.191 μb^-1 1.0 6.191 μb^-1
137221 111.932 nb^-1 0.876 0.475 0.622 0.679 360.330 mb^-1 24894.4 8.970 nb^-1
137222 4.682 pb^-1 0.936 0.507 0.616 0.678 5.830 μb^-1 35544.9 207.237 nb^-1
137611 3.606 pb^-1 0.896 0.500 0.641 0.675 2.620 μb^-1 39894.6 104.544 nb^-1
137622 4.682 pb^-1 0.977 0.507 0.637 0.678 5.830 μb^-1 24114.8 140.597 nb^-1
===========================================================================================

where

  • L_int_mb is the integrated luminosity seen by the minbias trigger for runs in which the specified trigger was active
  • eps_vtx is the vertex finding efficiency for this trigger
  • eps_mb is the vertex finding efficiency for mb-triggered events in runs in which the specified trigger was active
  • vz is the fraction of events with a reco vertex that had fabs(vz) < 60 cm
  • vz_mb is the same quantity for mb triggered events
  • L_samp_mb is N_good_vertex_mb / (eps_mb * sigmaBBC))
  • < ps > is sum_{runs} (ps_mb * n_mb) / sum_{runs} (ps_trig * n_trig)
  • L_sampled is L_samp_mb * < ps >

Data were obtained from the 289 spinTree runs without any additional QA selection.  If I restrict to the 188 runs in the latest jet runlist I get

===========================================================================================
trigId L_int_mb ε_vtx ε_mb vz vz_mb L_samp_mb < ps > L_sampled
-------------------------------------------------------------------------------------------
117001 003.286 pb^-1 0.513 0.513 0.676 0.676 4.757 μb^-1 1.0 4.757 μb^-1
137221 065.594 nb^-1 0.957 0.518 0.625 0.683 208.913 mb^-1 24838.9 5.189 nb^-1
137222 003.220 pb^-1 0.942 0.512 0.613 0.676 4.548 μb^-1 36362.2 165.371 nb^-1
137611 002.253 pb^-1 0.905 0.510 0.636 0.669 1.649 μb^-1 40234.9 066.358 nb^-1
137622 003.220 pb^-1 0.981 0.512 0.637 0.676 4.548 μb^-1 24376.0 110.859 nb^-1
===========================================================================================

The vertex finding efficiency for minbias triggers climbs a little bit, but it's still only 51%.  After doing a little digging it seems this is consistent with Jan's findings in his evaluation of PPV for 2006:

http://www.star.bnl.gov/HyperNews-star/protected/get/starspin/2820.html

BTOW pedestal comparison -- emcOnline vs. L2ped

I wrote a script that compares each Run 7 BTOW pedestal table with a corresponding L2ped printout. The script queries the RunLog_onl database for the runnumber immediately preceding the table timestamp, then it looks for a log file corresponding to that runnumber in the L2ped directory on online.star.bnl.gov. If it can’t find the precise runnumber, it compares the online pedestals to the next closest L2ped runnumbers both before and after, and then prints the comparison for the better match.

Here’s a summary log printing the number of channels per table where

  1. the difference in pedestal value is greater than 4, or
  2. the difference in pedestal sigma is greater than 1.5

There are a handful of channels per run with drastic differences. For example, emcOnline will declare a zero pedestal, but l2ped will say there’s a peak at 3500 ADCs. Channels like these are obviously bad. Of greater concern are the runnumbers where several hundred channels do not match. Most of the time it’s because the widths are different, but occasionally there are runs where a few hundred towers will have pedestals that differ by more than 4 ADCs. The full verbose output (listing every tower that fails the cut) can be found here.

.

.

Kaon-pion correlations in d+Au: varying the momentum cut + historic summary

In
this update
a set of results with Adam's suggested momentum cut is added to the pool of analysis results for d+Au.

First Look at Electron Jets

I've been working with Priscilla on adding her electron candidates to the spin trees.  spinTrees with electrons are now available, and we went ahead and made these plots investigating the properties of jets with identified electrons.  We analyzed BJP1 and L2Gamma triggers and plotted the p_T spectra, track multiplicity, and neutral energy fraction for all jets and for jets associated with an electron candidate (deltaR < 0.7).  Full-size plots available by clicking on each image.

Electron Cuts

  • Global dE/dx cut changing with momentum
  • nFitPoints >= 15
  • nDedxPoints >= 10
  • nHits / nPoss >= 0.52
  • track Chi2 < 4
  • DCAGlobal < 2
  • NEtaStrips > 1 && NPhiStrips > 1
  • Primary dE/dx cut changing with momentum
  • 0.3 < P/E < 1.5
  • -0.01287 < PhiDist < 0.01345
  • ZDist in [-5.47,1.796] (West) or [-2.706,5.322] (East)

Jet Cuts

  • 0 < R_T < 0.99
  • -0.7 < detEta < 0.9
  • Jet points at fired jet patch (BJP1 only)

l2gamma_rt bjp1_rt l2gamma_track_mult bjp1_track_mult l2gamma_ptComp bjp1_ptComp

Conclusions

  • The neutral energy fraction for L2gamma electron jets peaks at the same point as the ridge in the inclusive jet plot.
  • There are a significant number of 1-track jets in the L2gamma sample
  • In the plot of jet p_T vs. electron p_T we see the maximum of the distribution occurring at a electron pT of ~5.5 GeV and a jet pT of ~11 GeV.  This is in stark contrast to the BJP1 distribution. 

Average Partonic Pt Carried by Charged Pions

Here's a fragmentation study looking at the ratio of reconstructed charged pion p_{T} and the event partonic p_{T} in PYTHIA.  Cuts are

  • fabs(mcVertexZ()) < 60
  • fabs(vertexZ()) < 60
  • geantId() == 8 or 9 (charged pions)
  • fabs(etaPr()) < 1
  • fabs(dcaGl()) < 1
  • fitPts() > 25

fragmentation_parton_reco

Error bars are just the errors on the mean partonic p_{T} in each reconstructed pion p_{T} bin.  Next step is to look at the jet simulations to come up with a plot that is more directly comparable to a real data measurement.

Effect of Triggers on Relative Subprocess Contributions

These histograms plot the fraction of reconstructed charged pions in each pion pT bin arising from gg, qg, and qq scattering.  I use the following cuts:

  • fabs(mcVertexZ()) < 60
  • fabs(vertexZ()) < 60
  • geantId() == 8 or 9 (charged pions)
  • fabs(etaPr()) < 1
  • fabs(dcaGl()) < 1
  • fitPts() > 25

I analyzed Pythia samples from the P05ih production in partonic pT bins 3_4 through 55_65 (excluded minbias and 2_3).  The samples were weighted according to partonic x-sections and numbers of events seen and then combined.  StEmcTriggerMaker provided simulations of the HT1 (96201), HT2 (96211), JP1 (96221), and JP2 (96233) triggers.  Here are the results.  The solid lines are MB and are identical in each plot, while the dashed lines are the yields from events passing a particular software trigger.  Each image is linked to a full-resolution copy:

HT1HT2
JP1JP2

Conclusions

  1. Imposing an EMC trigger suppresses gg events and enhances qq, particularly for transverse momenta < 6 GeV/c.  The effect on qg events changes with pT.  The explanation is that the ratio (pion pT / partonic pT) is lower for EMC triggered events than for minimum bias.
  2. High threshold triggers change the subprocess composition more than low-threshold triggers.
  3. JP1 is the least-biased trigger according to this metric.  There aren't many JP1 triggers in the real data, though, as it was typically prescaled by ~30 during the 2005 pp run.  Most of the stats in the real data are in JP2.

Checklist for Run 6 preliminary result

  1. QA for BJP0 trigger and runlist generation
  2. Confirm vertex and nSigma cuts with systematic studies
  3. Run over new PYTHIA simulation, make basic data/simu comparisons
  4. Trigger Bias studies
  5. Extract A_TT from transverse data
  6. Randomized spin pattern studies
  7. PID background estimation

Notes

  • Adding the L2gamma trigger would improve statistics but we don't have an officially sanctioned simulation of that trigger available yet.
  • This checklist assumes the preliminary result is an inclusive measurement using only Run 6 data; other possibilities are to combine Runs 5 and 6 or to pursue the jet+pion measurement.