TPC SC&GL Calibration: Run19 AuAu 200 GeV
This document summarize the space charge and grid leak calibration for Run 19 AuAu 200 GeV dataset.
Details on the parameters definition and calibration can be found in the link below (Gene's Blog):
https://drupal.star.bnl.gov/STAR/blog/genevb/2013/mar/29/unified-spacecharge-gridleak-calibration-winter-2013
1. Data set summary:
Run Period: Au+Au 14.5 GeV
Trigger Setup: production_AuAu200_2019
Run time: Jul, 8 - Jul, 14
Pass01:
USING STDDEV sc => 0.000554 :: gapf => 0.015006
SpaceCharge & GridLeak fit results {scaler: bbcbb}:
g2 : 0.101986 +/- 7.88255 (AT LIMIT!)
log(g5) : 2.84587 +0/-1.71368 or +/- 0.170455
log(g5) : 2.84587 +/- 0.170455
log(SC) : -14.1021 +0.150461/-0.156133 or +/- 0.0661674
log(SC) : -14.1021 +/- 0.153297
SO : -2296.3 +911.982/-1030.53 or +/- 793.112
SO : -2296.3 +/- 971.254
*** FINAL CALIBRATION VALUES: ***
sc = (7.509e-07 +/- 1.151e-07) * ((bbcbb) - ( -2296 +/- 971.3))
with GL = 0.50 +/- 0.00
*** The following calibration values may not be trusted at this time... ***
*** Try the following calibration values: ***
sc = 1.33e-05 * ((bbcbb) - ( -2297)) with GL = 1.00
sc = 1.33e-05 * ((bbcbb) - ( -2297)) with GL = 0.50
sc = 1.33e-05 * ((bbcbb) - ( -2297)) with GL = 0.00
Pass02:
USING STDDEV sc => 0.000541 :: gapf => 0.006226
SpaceCharge & GridLeak fit results {scaler: bbcbb}:
g2 : 1.23094 +0/0 or +/- 0.0136802
g2 : 1.23094 +/- 0.0226187
log(g5) : 2.98279 +0/0 or +/- 0.0883189
log(g5) : 2.98279 +/- 0.452396
log(SC) : -14.218 +0.164116/-0.173445 or +/- 0.0539804
log(SC) : -14.218 +/- 0.16878
SO : -2296.52 +1002.53/-1128.57 or +/- 547.339
SO : -2296.52 +/- 1065.55
*** FINAL CALIBRATION VALUES: ***
sc = (6.686e-07 +/- 1.129e-07) * ((bbcbb) - ( -2297 +/- 1066))
with GL = 0.50 +/- 0.00
Pass03:
--> Number of runs increased for this pass to 34.
--> Single rich scaler failed for the combined SC&Gl fit. The best scalers are the bbcbb and vpdx.
USING STDDEV sc => 0.000432 :: gapf => 0.005796
SpaceCharge & GridLeak fit results {scaler: (4.28059e-07*(bbcbb-109.105))+(8.04667e-07*(vpdx))}:
g2 : 1.29106 +0/0 or +/- 0.073646
g2 : 1.29106 +/- 0.100391
log(g5) : 2.94698 +0/0 or +/- 0.024355
log(g5) : 2.94698 +/- 0.698164
log(SC) : -0.0214875 +0.112261/-0.132581 or +/- 0.0368819
log(SC) : -0.0214875 +/- 0.122421
SO : -0.000103431 +0.000545474/-0.000523979 or +/- 0.000260594
SO : -0.000103431 +/- 0.000534726
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1224)*(4.19e-07*(bbcbb-(-132.5 +/- 1249))+(7.876e-07*(vpdx)))
with GL = 0.50 +/- 0.00
--> bbcbb(bbc blue background) should not be the dominant. We should find other rich scaler, in addition to bbcbb.
--> Calibration macro found vpdx with the bbcbb dependence.
USING STDDEV sc => 0.000431 :: gapf => 0.005797
SpaceCharge & GridLeak fit results {scaler: (8.04667e-07*(vpdx-58.0406))+(4.28059e-07*(bbcbb))}:
g2 : 1.29105 +0/0 or +/- 0.0734803
g2 : 1.29105 +/- 0.0888637
log(g5) : 2.94708 +0/-1.5892 or +/- 0.132465
log(g5) : 2.94708 +/- 0.132465
log(SC) : -0.0210332 +0.112235/-0.132584 or +/- 0.0296126
log(SC) : -0.0210332 +/- 0.122409
SO : -0.00010147 +0.000545167/-0.000523703 or +/- 0.000276431
SO : -0.00010147 +/- 0.000534435
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1224)*(7.879e-07*(vpdx-(-68.06 +/- 664.2))+(4.191e-07*(bbcbb)))
with GL = 0.50 +/- 0.00
Pass04:
1 scaler fit result:
USING STDDEV sc => 0.000501 :: gapf => 0.005777
SpaceCharge & GridLeak fit results {scaler: bbcbb}:
g2 : 1.41623 +0/0 or +/- 0.0979595
g2 : 1.41623 +/- 0.159779
log(g5) : 2.91896 +0/-1.68534 or +/- 5.8377
log(g5) : 2.91896 +/- 5.8377
log(SC) : -14.2378 +0/0 or +/- 0.00641418
log(SC) : -14.2378 +/- 0.112792
SO : -2533.8 +0/0 or +/- 20.8093
SO : -2533.8 +/- 498.321
*** FINAL CALIBRATION VALUES: ***
sc = (6.555e-07 +/- 7.394e-08) * ((bbcbb) - ( -2534 +/- 498.3))
with GL = 0.50 +/- 0.00
2 scaler fit result:
USING STDDEV sc => 0.000432 :: gapf => 0.005772
SpaceCharge & GridLeak fit results {scaler: (4.58495e-07*(bbcbb+233.028))+(7.16046e-07*(vpdx))}:
g2 : 1.41637 +2.55131/-2.61671 or +/- 0.251437
g2 : 1.41637 +/- 2.58401
log(g5) : 2.87697 +0/-1.53881 or +/- 0.238663
log(g5) : 2.87697 +/- 0.238663
log(SC) : -0.0258936 +0.112118/-0.132044 or +/- 0.132647
log(SC) : -0.0258936 +/- 0.122081
SO : -0.000118954 +0.000545006/-0.000523367 or +/- 0.000616458
SO : -0.000118954 +/- 0.000534186
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1221)*(4.468e-07*(bbcbb-(-492.5 +/- 1165))+(6.977e-07*(vpdx)))
with GL = 0.50 +/- 0.00
3 scaler fit:
USING STDDEV sc => 0.000400 :: gapf => 0.005773
SpaceCharge & GridLeak fit results {scaler: (-1.627e-07*(zdce-2461.34))+(3.49884e-06*(bbcx))+(8.01285e-07*(bbcbb))}:
g2 : 1.41635 +0/0 or +/- 0.0997554
g2 : 1.41635 +/- 0.404707
log(g5) : 2.91594 +0/-1.49877 or +/- 0.457368
log(g5) : 2.91594 +/- 0.457368
log(SC) : -0.0207615 +0.100825/-0.119729 or +/- 0.128375
log(SC) : -0.0207615 +/- 0.110277
SO : -0.00011094 +0.000495095/-0.000466171 or +/- 0.000568761
SO : -0.00011094 +/- 0.000480633
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1103)*(-1.594e-07*(zdce-( 3143 +/- -2954))+(3.427e-06*(bbcx))+(7.848e-07*(bbcbb)))
with GL = 0.50 +/- 0.00
-: Without fixing GL:
Fit failed for up to 2 rich scalers and converged for 3 scalers. The GL value that the calibration macro is small.
USING STDDEV sc => 0.000400 :: gapf => 0.005773
SpaceCharge & GridLeak fit results {scaler: (-1.6645e-07*(zdce-2463.3))+(3.57942e-06*(bbcx))+(8.19764e-07*(bbcbb))}:
g2 : 1.43365 +0/0 or +/- 0.0305587
g2 : 1.43365 +/- 0.0719667
log(g5) : 2.93655 +0/-1.40347 or +/- 6.59487
log(g5) : 2.93655 +/- 6.59487
log(SC) : -0.0194186 +0.100327/-0.12003 or +/- 0.128133
log(SC) : -0.0194186 +/- 0.110178
SO : -0.000119374 +0.000508244/-0.000475513 or +/- 0.000581458
SO : -0.000119374 +/- 0.000491879
log(GL) : -2.8572 +0/0 or +/- 19.9402
log(GL) : -2.8572 +/- 5.56547
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1102)*(-1.632e-07*(zdce-( 3180 +/- -2955))+(3.511e-06*(bbcx))+(8.04e-07*(bbcbb)))
with GL = 0.06 +/- 0.32
-: Trying to find scalers combination with "VPDX" , instead of "bbcbb".
USING STDDEV sc => 0.000432 :: gapf => 0.005767
SpaceCharge & GridLeak fit results {scaler: (4.67167e-07*(bbcbb+233.565))+(7.29541e-07*(vpdx))}:
g2 : 1.43365 +0/0 or +/- 0.0190483
g2 : 1.43365 +/- 0.0226262
log(g5) : 2.90868 +0/-1.53781 or +/- 0.67751
log(g5) : 2.90868 +/- 0.67751
log(SC) : -0.0250654 +0.112517/0 or +/- 0.00985544
log(SC) : -0.0250654 +/- 0.00985544
SO : -0.00012664 +0/0 or +/- 8.99147e-06
SO : -0.00012664 +/- 0.000721948
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.009855)*(4.556e-07*(bbcbb-(-504.6 +/- 1545))+(7.115e-07*(vpdx)))
with GL = 0.15 +/- 0.00
-:Best Alternative scalers combination:
USING STDDEV sc => 0.000475 :: gapf => 0.005766
SpaceCharge & GridLeak fit results {scaler: (-1.05823e-07*(zdcw+1169.82))+(3.8649e-07*(vpdw))}:
g2 : 1.43365 +0/0 or +/- 0.00496243
g2 : 1.43365 +/- 0.0316213
log(g5) : 3.00836 +0/-1.69777 or +/- 0.463134
log(g5) : 3.00836 +/- 0.463134
log(SC) : -0.0882245 +0.140457/-0.133848 or +/- 0.129494
log(SC) : -0.0882245 +/- 0.137152
SO : -0.000349999 +0.000599284/-0.00070045 or +/- 0.000639803
SO : -0.000349999 +/- 0.000649867
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1372)*(-9.689e-08*(zdcw-( 2138 +/- -6141))+(3.539e-07*(vpdw)))
with GL = 0.15 +/- 0.00
-: Checking for three scalers:
- Best way is to run the calibration macro with debug mode 2. This will print out all scalers combination and RMS value.
- Below is the last segment of the output of pca iteration 1, with the best alternative for bbcbb or bbcyb at the top and the best scalers combination at the end.
RMS = 0.000442148 :: EV = 3.12944e-15 ::
sc = (-2.89603e-07*(bbce+82.6381))+(7.71676e-07*(vpdw))+(-1.9005e-07*(zdcwnk))
RMS = 0.000441458 :: EV = 5.93969e-14 ::
sc = (9.57245e-07*(bbcx+246.991))+(1.5014e-07*(bbcyb))+(5.27054e-07*(bbcbb))
RMS = 0.000439105 :: EV = 7.05269e-14 ::
sc = (-7.42001e-07*(zdcx-186.82))+(2.36616e-06*(bbcx))+(5.37057e-07*(bbcbb))
RMS = 0.000437468 :: EV = 6.4837e-14 ::
sc = (2.43126e-06*(bbcx+49.1797))+(5.30576e-07*(bbcbb))+(-6.58482e-07*(zdcxnk))
RMS = 0.000436185 :: EV = 8.4599e-14 ::
sc = (4.66794e-07*(bbcbb+234.359))+(7.28887e-07*(vpdx))
RMS = 0.000436152 :: EV = 3.01495e-15 ::
sc = (1.23364e-09*(zdce+83997.9))+(4.63723e-07*(bbcbb))+(7.20786e-07*(vpdx))
RMS = 0.000436099 :: EV = 2.31759e-14 ::
sc = (4.68289e-07*(bbcbb+250.444))+(7.80083e-07*(vpdx))+(-1.21371e-08*(vpde))
RMS = 0.000436089 :: EV = 2.21524e-15 ::
sc = (4.61777e-07*(bbcbb+220.142))+(7.15448e-07*(vpdx))+(1.75066e-09*(zdcenk))
RMS = 0.000436048 :: EV = 4.79478e-15 ::
sc = (4.67149e-07*(bbcbb+252.9))+(7.88432e-07*(vpdx))+(-5.46125e-09*(zdcwnk))
RMS = 0.000435685 :: EV = 6.12704e-15 ::
sc = (-1.31704e-08*(zdcw-10545.8))+(4.66133e-07*(bbcbb))+(8.59077e-07*(vpdx))
RMS = 0.000434992 :: EV = 6.09512e-14 ::
sc = (4.60941e-07*(bbcbb+216.16))+(6.8601e-07*(vpdx))+(4.69222e-08*(zdcxnk))
RMS = 0.000434951 :: EV = 7.9492e-14 ::
sc = (1.61901e-07*(bbcx+596.969))+(4.65033e-07*(bbcbb))+(6.42247e-07*(vpdx))
RMS = 0.000434833 :: EV = 6.56142e-14 ::
sc = (5.83183e-08*(zdcx+1692.2))+(4.60029e-07*(bbcbb))+(6.83808e-07*(vpdx))
RMS = 0.000433605 :: EV = 5.50671e-14 ::
sc = (-4.75892e-08*(bbcyb-1578.33))+(4.48138e-07*(bbcbb))+(7.839e-07*(vpdx))
RMS = 0.000429181 :: EV = 1.10779e-14 ::
sc = (6.48728e-07*(bbcbb+390.148))+(1.57565e-06*(vpdx))+(-1.73514e-07*(vpdw))
RMS = 0.000427295 :: EV = 7.25114e-15 ::
sc = (-5.21289e-08*(bbcw-1594.84))+(4.96856e-07*(bbcbb))+(1.11242e-06*(vpdx))
RMS = 0.000423745 :: EV = 1.79717e-15 ::
sc = (-1.04716e-06*(bbcyb-52.1191))+(7.33083e-07*(vpdw))+(-2.32001e-07*(zdcenk))
RMS = 0.000420619 :: EV = 1.35274e-14 ::
sc = (-1.30242e-07*(bbce-908.052))+(6.07553e-07*(bbcbb))+(1.19709e-06*(vpdx))
RMS = 0.000415867 :: EV = 2.25514e-15 ::
sc = (-2.93367e-07*(zdce-1071.59))+(-1.09714e-06*(bbcyb))+(7.75948e-07*(vpdw))
RMS = 0.000407995 :: EV = 1.94289e-15 ::
sc = (3.39897e-06*(bbcx+70.4538))+(7.82787e-07*(bbcbb))+(-1.3115e-07*(zdcenk))
RMS = 0.000402972 :: EV = 2.583e-15 ::
sc = (-1.6564e-07*(zdce-2463.29))+(3.562e-06*(bbcx))+(8.15775e-07*(bbcbb))
The alternative combination is comparable to what we have with 2 scalers combination. Since the result is not improving with the number of scalers, we should use 2 scalers result as final parameters for this dataset.
-: With the PCA's best scaler combination above:
USING STDDEV sc => 0.000400 :: gapf => 0.005767
SpaceCharge & GridLeak fit results {scaler: (-1.6564e-07*(zdce-2463.29))+(3.562e-06*(bbcx))+(8.15775e-07*(bbcbb))}:
g2 : 1.43365 +0/0 or +/- 0.0187684
g2 : 1.43365 +/- 0.0240966
log(g5) : 2.95104 +0/-1.48717 or +/- 4.27221
log(g5) : 2.95104 +/- 4.27221
log(SC) : -0.0160389 +0.100438/-0.119953 or +/- 0.128138
log(SC) : -0.0160389 +/- 0.110196
SO : -0.000103521 +0.000503699/-0.000471902 or +/- 0.000576608
SO : -0.000103521 +/- 0.0004878
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1102)*(-1.63e-07*(zdce-( 3088 +/- -2945))+(3.505e-06*(bbcx))+(8.028e-07*(bbcbb)))
with GL = 0.15 +/- 0.00
There is no significant difference between 2 and 3 scalers that PCA finds.
------FINAL PARAMETERS---------
USING STDDEV sc => 0.000475 :: gapf => 0.005766
SpaceCharge & GridLeak fit results {scaler: (-1.05823e-07*(zdcw+1169.82))+(3.8649e-07*(vpdw))}:
g2 : 1.43365 +0/0 or +/- 0.00496243
g2 : 1.43365 +/- 0.0316213
log(g5) : 3.00836 +0/-1.69777 or +/- 0.463134
log(g5) : 3.00836 +/- 0.463134
log(SC) : -0.0882245 +0.140457/-0.133848 or +/- 0.129494
log(SC) : -0.0882245 +/- 0.137152
SO : -0.000349999 +0.000599284/-0.00070045 or +/- 0.000639803
SO : -0.000349999 +/- 0.000649867
*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.1372)*(-9.689e-08*(zdcw-( 2138 +/- -6141))+(3.539e-07*(vpdw)))
with GL = 0.15 +/- 0.00
---------VVVVVVVVVVVVV------------------------
- pokhrel's blog
- Login or register to post comments