Data used for this analysis is 2006 p+p 200 GeV data taken with transverse polarisation, trigger setup "ppProductionTrans". This spanned days 97 (7th April) to 129 (9th May) inclusive. Trigger bemc-jp0-etot-mb-l2jet (ID 127622) is used. A file catalogue query with the following conditions gives a list of runs for which data is available:
trgsetupname=ppProductionTrans,tpc=1,year=2006,sanity=1,collision=pp200,
magscale=FullField,filename~physics,library=SL06e,production=P06ie
This generates a list of 549 runs. These runs are then compared against the spin PWG run QC (see http://www.star.bnl.gov/protected/spin/sowinski/runQC_2006) and are rejected if any of the following conditions are true:
This excludes 172 runs, leaving 377 runs to be analysed.
I use a Maker class to create TTrees of event objects with V0 and spin information for these runs. Code for the Maker and Event classes can be found at /star/u/tpb/StRoot/StTSAEventMaker/ and /star/u/tpb/StRoot/StV0NanoDst/ respectively. Events are accepted only if they fulfill the following criteria:
TTrees are produced for 358 runs (19 produce no/empty output), yielding 2,743,396 events.
The vertex distribution of events from each run are then checked by spin bits. A Kolmogorov test (using ROOT TH1::KolmogorovTest) is used to compare the vertex distributions for (4-bit) spin bits values 5, 6, 9 and 10. If any of the distributions are inconsistent, the run is rejected. Each run's mean event vertex z position is then plotted. Figure 1 shows the distribution, fitted with a Gaussian. A 3σ cut is applied and outlier runs rejected. 38 runs are rejected by these further cuts. The remaining 320 runs, spanning 33 RHIC fills and comprising 2,500,421 events, are used in the analysis.
Figure 1: Mean event vertex z for each run. The red lines indicate the 3σ cut. |