R&D on GridLeaks
productionHalfHigh: sector PHI =
productionHigh: sector PHI =
productionHalfHigh: |z| =
productionHigh: |z| =
productionHalfHigh:
We don't learn a lot more from these, but for documentation-sake:
Here are the same plots limited to |z|=25-50cm:
Here are the same plots limited to 0-1 EIGHTH (described earlier, all z):
productionHigh (all hits for each sector):
productionHalfHigh, 1 o'clock: sector PHI =
productionHalfHigh, 4 o'clock: sector PHI =
productionHalfHigh, 7 o'clock: sector PHI =
productionHalfHigh, 10 o'clock: sector PHI =
[Left plots are azimuthal dependence for |z|<25, where horizontal axis is "o'clock" as viewed from the west side looking east. So 9 o'clock is north, towards the inside of the RHIC ring. Blue is for sectors 1-12 (west), and red is for sectors 13-24 (east). Middle plots are z-dependence for all sectors. Right plots are azimuthal and z-dependence in one plot (its middle 2 z bins are the data in the left plot).]
productionHalfHigh: (BEFORE CORRECTION)
Beautifully linear z-dependence!
productionHalfHigh: (AFTER CORRECTION)
Remaining gap seems to be prominent in the 9-12 o'clock region.
productionHigh: (AFTER CORRECTION)
Reversed field flips the sign of the Gap. In the 3D plot, I have flipped the sign of the Gap so that any similarities to the forward field plots can be more easily seen. There is perhaps a similar effect of notable remaining Gaps at 9-12 o'clock and no Gaps at 4-7 o'clock, but the 1-3 o'clock region is markedly different. Also, there are notable issues at large z values on both sides of the TPC (the worst at sector 24, east side, negative z, 12 o'clock).
productionHigh: Day 46 (AFTER CORRECTION)
Forward field, so no sign flip. Looks like some really odd problem at large negative Z, but the 3D plot shows that most of this problem is isolated to one sector at 2 o'clock (sector 22). Otherwise similar to the productionHalfHigh taken on day 58.
Run 6035071: test (high gated grid rate, ZDCx = approx. 20k)
A little worse than last year's AuAu, but not an order of magnitude (compare 2mm gap at the central membrane with the un-GridLeak-corrected productionHalfHigh AuAu value of about 1.2mm). ZDCx rates are significantly higher for CuCu, though event sizes are smaller.
Run 6035027: control (normal rate, almost same ZDCx, taken several hours earlier during previous fill)
Not much different from the test file!!! So the gap appears to be GG rate independent.
Run 6035042: slightly lower ZDCx, taken between above two runs
Just making sure we still see a luminosity dependence, which we do.
Now applying SpaceCharge to run 6035071:
I expected there to be no effect on the observed GridLeak gap from applying SpaceCharge, but I was wrong. This may be due to how I measure the gap (difference of two residuals at padrows 12 and 14 instead of closer to the sector boundaries). [NOTE: SpaceCharge was determined by a prepass on the first few events of the file, but event-by-event SpaceCharge was not operational for as yet unknown reasons.]
Now applying SpaceCharge + GridLeak to run 6035071:
We can remove the gap using the same GridLeakStrength (15.0) that we found for the AuAu data...on average. The azimuthal dependence we have always seen is still there and is actually worse than ever. Runs 6035027 (my control) and 6035071 (the high GG rate test) share similar phi dependence for the gap, so the phi dependence appears to be GG rate independent.
Small note: ioratio = 0.75 worked best.
Now applying SpaceCharge to run 6035027:
Unfinished section...Not getting SpaceCharge prepass to work: no SpaceCharge means no GridLeak correction...
Now applying SpaceCharge + GridLeak to run 6035027:
I again applied the GridLeakStrength of 15.0, and this time EventByEvent SpaceCharge worked for some reason. It looks like we actually on average overcorrect (the V-shape in the Z dependence plot points upwards, though it points downwards before corrections). The strong azimuthal dependence remains.
GridLeak Distortion
I recommend using tabbed browsing to view these plots, whereby you open a series of plots in a series of tabs, and then move through the tabs to see the variations/dependencies.
Yuri's discussion of the problem and explanation of these plots: here
EbyE volume space charge is running in the following plots: Yuri's plot showing luminosity and B-field dependence
One of my first plots: Half field
The proposed location for GridLeaks can be seen by the gaps in the gating grid seen on the schematic diagrams found in this postscript file.
Other distortion corrections in place (no GridLeak, no volume SC):
Attempting to remove residual gap by varying GLS on half field run:
GLS (GridLeakageStrength) is a multiplier which defines how strong the grid leak charge is relative to the volume SpaceCharge. Tying the two together like this allows the charges/distortions/corrections to fluctuate along with the luminosity (as defined via the volume SpaceCharge correction).
The EbyE SpaceCharge correction is left running to bring the mean physical signed DCAs of global tracks to zero. The prepass must also be run because the GridLeak distortion changes DCAs, so a different guess of the necessary SpaceCharge magnitude is necessary (NOTE: StMagUtilities::PredictSpaceCharge() must call the function to undo the GridLeak distortion as well as the SpaceCharge).
During this process, I learned that I needed to reduce the InnerOuterRatio in StMagUtilities::PredictSpaceCharge() down to 0.5 (from 1.3) in order to get the prepass to work.
Application of GLS from half field solution to other luminosities and fields using the azimuthally symmetric GridLeak correction. There appears to be some luminosity-dependent issues still:
First attempt at moving the radius of the azimuthally symmetric GridLeak to higher radii using high luminosity run. This moved significant distortions into the outer sector. Using high luminosity run with GLS=65.
Now using sector-faceted GridLeak (no longer azimuthally symmetric, and GLS magnitude scale changed by a factor of ~4-5 to a value of ~15 to remove the gap) and allowing for inner/outer leaks. iGLS is the fraction of GLS for the strength of the inner sector leak, and oGLS for the outer. Inner sheet is at radius 55 +/- 3, and outer is at 195 +/- 3. Varying iGLS and oGLS together using high luminosity run:
Since the outer leak did not seem to help, now setting oGLS=0 and varying only iGLS:
SC = SpaceCharge. Had used only the Wieman (~1/r^2) shape until this point.
High luminosity, RFF run, GLS=15, no inner or outer grid leaks:
Change GLS to compensate for different volume SC necessary to bring DCAs to zero:
SpaceCharge shape = 1/r^3, GLS=15, no inner or outer grid leaks:
Looks pretty good! So the charge distributions must be something similar to this (1/r^3 moves more charge to small radii).
Following the hypothesis that an inner leak plus Wieman volume SC might be able to create similar distortions to the 1/r^3 volume space charge field, we started looking more closely at really large inner leaks (large iGLS). In the inner/outer leaks section above, iGLS of 25.0 and 50.0 changed the DCAs so much that the volume SC magnitude changed, thereby affecting the inner/outer gap because GLS (essentially the ratio of gap charge density to volume charge density) stayed constant at 15.0. Here, we try changing GLS while holding iGLS constant (and large) to see what happens.
iGLS=50:
iGLS=25:
From that data, it looks like the optimum relation for keeping the gap zero may be approximately:
GLS = 15.0 - iGLS/10.0:
Not perfect, so the above relation is likely inexact. And some structure between padrows 2 and 7 remains. A little tweaking of parameters:
Looking at those last two vesus luminosity: iGLS=35, GLS=11.0:
iGLS=40, GLS=10.5: