Stv vs. Sti with first look at Run 9 pp500 real data

Looking at hits on global tracks in Run 9 pp500 real data events (exact same sample of 40 events for both Sti and Stv).

Δrow is defined as the difference in row number between any hit and the previous hit on the track (the hits are arranged such that the first hit is the innermost, which itself does not have a previous hit).

Observations:
  • Standard tracks have a 10 hit minimum in Stv, but 11 in Sti
  • Special short tracks have a 6 hit minimum in Stv, but 5 in Sti
  • Stv struggles to span dead regions more than does Sti (best demonstrated in sectors 5 and 6)
  • Stv can have negative Δrow, which appears to occur on low pT loopers; this is not allowed in Sti
  • Stv has a minimum pT of ~40-50 MeV/c, but Sti goes down to ~10 MeV/c
    • This appears to be due to a maximum restriction on track curvature of 10-1.5 cm-1 in Stv
  • For adjacent hits on tracks which are within the same sector, Stv generally limits Δrow to <~8 where there are no dead regions, and <~15 where there are dead regions; Sti limits to <~13 nearly everywhere
  • For adjacent hits on tracks which are in different sectors (track crossed a sector boundary or the central membrane), Stv generally allows a slightly larger Δrow, but Sti allows a much larger Δrow, up to ~30
  • Signed DCA distributions are not as well peaked in Stv as in Sti
  • Very few tracks appear to cross the central membrane in Stv (first and last hits are on the same side of the TPC), but it may only be a problem for pile-up tracks (i.e. not a problem for triggered event tracks)
  • Sti assigns the primary vertex to be the "first hit" on primary tracks, but Stv uses the first TPC hit
-Gene

NB: additional columns have been added for Stv as of Sep. 10th and Sep. 23rd with some items addressed.

  Sti Stv Stv Sep. 10 Stv Sep. 23
# hits on tracks
row vs. sector of all hits on tracks
row vs. sector of all hits
on short (# hits <= 20) tracks 
Δrow
Δradius vs. Δrow

(multiple bands are due to inner
sector vs. outer sector row spacing)
Δrow vs. log10(pT)
Δrow vs. φ with both hits
in same sector
Δrow vs. φ with hits
in different sectors
log10(cuvature)  
log10(p) vs. track length  
signed DCA distribution    
zlast hit vs. zfirst hit    
Primary track
zlast hit vs. zfirst hit
 
Primary track
yfirst hit vs. xfirst hit
 


__________________

Other (older) Stv datasets' zlast hit vs. zfirst hit:
  • Parent vertex id == 1 are tracks from the "triggered" collision vertex; others are pile-up and/or secondaries
  • Steepness of slope in these plots is a clear correlator to the source vertex:
    • Steeper slopes, of ~190/60, are triggered collision primaries, as this is ~(radiusoutermost padrow / radiusinnermost padrow), shifted left and right by the z of the triggered collision primary vertex
    • Shallower slopes, randomly spread around 1, are the pile-up and/or secondaries
  • AuAu200 without pile-up shows no noticeable effects of the TPC central membrane
  • pp500 with pile-up shows some effects of the central membrane
    • The effect on primary tracks from the triggered collision is modest, but not absent (as in the AuAu200 without pile-up)
    • The effect on other tracks is perhaps just a little stronger, but not as strong as seen in the pp500 real data shown above

Stv simulations AuAu200
without pile-up
2013-08-05
50 events
pp500
with 1MHz pile-up
2013-06-06
5000 events
Globals
parent vtx id == 1
Globals
parent vtx id > 1
Primaries
parent vtx id == 1
Primaries
parent vtx id > 1


__________________

New Stv datasets' zlast hit vs. zfirst hit:
  • For the Stv datasets, all tracks had parent vertex id = 1
  • Without pile-up, things look as expected
  • Using Stv primaries as a proxy for the triggered collision tracks, it is not obvious that things are OK when there is pile-up
  • Sti globals with parent vertex id = 1 look fine, implying that the simulation of pile-up has no inherent problems

pp500 simulations
2013-08-29
500 events
without pile-up with 1MHz pile-up
Stv
Globals
parent vtx id == any
Sti
Globals
parent vtx id == any
Sti
Globals
parent vtx id == 1
Sti
Globals
parent vtx id > 1
Stv
Primaries
parent vtx id == any
Sti
Primaries
parent vtx id == any
Sti
Primaries
parent vtx id == 1
Sti
Primaries
parent vtx id > 1


__________________

Fraction of found hits placed on tracks:

Looking again at the same 40 events of real Run 9 pp500 data from the top of this blog page...
  • Sti and Stv start out with the exact same sample of TPC hits (a good check!)
  • Stv places zero "wrong side" hits on tracks, though the transition to zero happens over a length of a few centimeters, not exactly at z=0
    • "wrong side" hits are pile-up hits in the TPC recorded in timebins which come after the full drift length readout for the triggered event, and are thus reconstructed with a z position which is beyond the central membrane for the sector in which it was read out
  • Sti has a higher efficiency of placing outer sector (padrow > 13) hits on tracks, but Stv has a smoother transition in efficiency from outer to inner sectors
  • Both Sti and Stv have an inefficieny for placing hits at the innermost padrows on tracks (which may be OK: distortion corrections are incorrect for pile-up, and are largest at the innermost padrows, perhaps causing pile-up tracks to lose hits there), but Stv clearly shows a much stronger effect
  • Sti has a higher efficiency for using hits (real from the triggered event, real from pile-up, or noise) in the reconstructed z region of the gating grid not yet being open; for much of this region Stv's efficiency for putting a hit on a track is zero
fraction of hits
on tracks
Sti Stv
vs. sector and reconstructed z
near the central membrane
vs. row and sector
vs. row and sector
excluding "wrong side" hits
vs. row and reconstructed |z|
near the endcap



__________________

Sep. 7 patch for fraction of found hits placed on tracks:

First 20 events from st_W_10102066_raw_7170001:
  • Only sector 1 seems to really use the "wrong side" hits (though perhaps a few of these hits are used in sectors 4, 11, 12)
  • Sector 1 has also seen a big drop in efficiency (~50%) for placing any hits on tracks
fraction of hits
on tracks
eval13a eval13 + Sep 7 patch
vs. sector and reconstructed z
near the central membrane


__________________

Sep. 15 patch for fraction of found hits placed on tracks:

First 20 events from st_W_10102066_raw_7170001:
  • Looks good!
fraction of hits
on tracks
eval13a eval13 + Sep 15 patch
vs. sector and reconstructed z
near the central membrane


__________________

Signed DCAs (after the Sep. 15 patch):

All 1900 events from st_W_10102066_raw_7170001, all global tracks with flag>0 and with no wrong sign (post-membrane) hits (from the Offline QA histograms, using the DCA geometry's helix function geometricSignedDistance(xvtx,yvtx); units are [cm]):

NB: I had to skip the check on track topology map for having TPC hits, as Stv does not properly populate the topology map still.
  • Top left: log10(DCA)
  • Top right: sDCA of all charges
  • Bottom left: sDCA of neg. charge
  • Bottom right: sDCA of pos. charge
Observations:
  • Stv has a much larger spread in sDCA which is not attributable to charge splitting (it's there for both charges)
  • I added some simulations from earlier this year, only including sDCA of all charges
    • Stv is much better in these simulations, and actually better than Sti in the pp500 simulation with pile-up

Dataset: Sti Stv
Run 9 pp500
real data
Run 9 pp500
simulation
Run 10 AuAu200
simulation


________________
Sep 23 Stv update: Still a problem:



______________
Oct. 1 Stv update: resolved! The problem was an inversion of (x,y) of the position of the DCA point in Stv to (-x,-y).

In addition to the Run 9 pp500 shown below, we re-ran the Run 10 AuAu200 ideal simulation and I re-produced the plot for it but did not include it here because it looks nearly identical to the above-shown plot of sDCA for AuAu200 from August 5th. A closer look makes it a little clearer that the plot has flipped left-to-right, which tells me that while the sign flip was present in the old simulation, the fact that the simulation threw vertices at (0,0) meant that a sign inversion in the in the DCA position did not lead to a smearing of the sDCA values.