Home » Blogs » pokhrel's blog

TPC Space Charge Calibration: 2019 AuAu 14.5 GeV

Updated on Mon, 2021-04-26 07:52. Originally created by pokhrel on 2021-03-26 08:46.


--> Working Directory: /star/u/pokhrel/TPC_Cal/Run19/AfterAlignFixed/14p5/Pass0*

First Pass :
Single Parameter Fit:  Didn't converge. 

Two Parameter Fit:

 USING STDDEV sc => 0.000252 :: gapf => 0.005187

SpaceCharge & GridLeak fit results {scaler: (-9.90553e-09*(bbce+168065))+(3.53425e-06*(zdcenk))}:
g2          : 0.945327 +6.55063/-6.79905 or +/- 0.889722
g2          : 0.945327 +/- 6.67484
log(g5)     : 3.84153 +0/-2.84413 or +/- 6.80541
log(g5)     : 3.84153 +/- 6.80541
log(SC)     : -0.208838 +0.332672/-0.503489 or +/- 0.399168
log(SC)     : -0.208838 +/- 0.418081
SO          : 0.000199262 +0.000426394/-0.000430237 or +/- 0.000432464
SO          : 0.000199262 +/- 0.000428315
log(GL)     : 0.203774 +2.00044/0 or +/- 6.45753
log(GL)     : 0.203774 +/- 6.45753

*** FINAL CALIBRATION VALUES: ***
sc = (1.0 +/- 0.4181)*(-8.039e-09*(bbce-(-1.882e+05 +/- -4.324e+04))+(2.868e-06*(zdcenk)))
  with GL =  1.23 +/-  7.92


------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------
 Second Pass Results:   

--> 2-scalers parameters from the first pass
--> 1-scaler fit failed.
  
--> 2-scalers fit converged with large GL value:

USING STDDEV sc => 0.000293 :: gapf => 0.005051

 

SpaceCharge & GridLeak fit results {scaler: (-1.40998e-07*(bbcbb+5440.09))+(6.11532e-07*(zdcwnk))}:

g2          : 0.0930733 +1.80635/-1.74629 or +/- 1.91167

g2          : 0.0930733 +/- 1.77632

log(g5)     : 3.12291 +0/-1.23776 or +/- 0.852025

log(g5)     : 3.12291 +/- 0.852025

log(SC)     : -0.0311908 +0.367599/-0.585343 or +/- 0.0801808

log(SC)     : -0.0311908 +/- 0.476471

SO          : -2.53613e-06 +0.000347722/-0.00034914 or +/- 0.000134658

SO          : -2.53613e-06 +/- 0.000348431

log(GL)     : 1.66928 +1.20192/0 or +/- 0.242101

log(GL)     : 1.66928 +/- 0.242101

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.4765)*(-1.367e-07*(bbcbb-( -5422 +/-  -2471))+(5.928e-07*(zdcwnk)))

  with GL =  5.31 +/-  1.29

 

WARNING! Symptoms of non-linearity in chosen luminosity scaler!

  Suggestion: use more PCA variables


---------------------------------------------------------------------------------------------------
-->Large GL value. Possibly due to the one outlier (Gene). Removing outlier (run =20099044) results following: 

USING STDDEV sc => 0.000282 :: gapf => 0.005084

 

SpaceCharge & GridLeak fit results {scaler: (-7.67518e-09*(bbce+111571))+(2.69729e-06*(zdcenk))}:

g2          : 0.765402 +0/0 or +/- 0.178189

g2          : 0.765402 +/- 0.595587

log(g5)     : 3.47941 +0/0 or +/- 13.4737

log(g5)     : 3.47941 +/- 1.06141

log(SC)     : 0.00716186 +0.609871/-1.85467 or +/- 0.855324

log(SC)     : 0.00716186 +/- 1.23227

SO          : -1.94927e-06 +0.00036133/-0.000365084 or +/- 0.000366299

SO          : -1.94927e-06 +/- 0.000363207

log(GL)     : 1.02164 +1.98045/0 or +/- 6.36327

log(GL)     : 1.02164 +/- 6.36327

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/-  1.232)*(-7.73e-09*(bbce-(-1.113e+05 +/- -4.732e+04))+(2.717e-06*(zdcenk)))

  with GL =  2.78 +/- 17.68

 

--> GL value is still large with large error. 
--> Decided to fix GL value to be same as 19.6 GeV GL value (GL = 0.12 for scalers)
--------------------------------------------------------------------------------------------
  2-scalers fit result for fix GL=0.12

USING STDDEV sc => 0.000293 :: gapf => 0.005050

 

SpaceCharge & GridLeak fit results {scaler: (-1.68037e-07*(bbcbb+5440.11))+(7.28787e-07*(zdcwnk))}:

g2          : 0.327587 +8.54787/0 or +/- 0.708427

g2          : 0.327587 +/- 0.708427

log(g5)     : 3.1224 +0/-3.39622 or +/- 4.33947

log(g5)     : 3.1224 +/- 4.33947

log(SC)     : -0.00173277 +0.367715/-0.585539 or +/- 0.182512

log(SC)     : -0.00173277 +/- 0.476627

SO          : -3.20851e-06 +0.00041442/-0.000416293 or +/- 0.000345677

SO          : -3.20851e-06 +/- 0.000415356

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.4766)*(-1.677e-07*(bbcbb-( -5421 +/-  -2472))+(7.275e-07*(zdcwnk)))

  with GL =  0.12 +/-  0.00

-------------------------------------------------------------------------------------------------- 

3-Scalers Fit:

USING STDDEV sc => 0.000199 :: gapf => 0.005029

 

SpaceCharge & GridLeak fit results {scaler: (-5.8129e-07*(bbcbb+2229.45))+(8.69927e-07*(zdcwnk))+(1.69922e-06*(zdcenk))}:

g2          : 0.388394 +4.32351/-4.37688 or +/- 1.39945

g2          : 0.388394 +/- 4.3502

log(g5)     : 3.17686 +0/-1.70986 or +/- 5.94774

log(g5)     : 3.17686 +/- 5.94774

log(SC)     : 0.00137251 +0.254412/-0.349918 or +/- 0.295718

log(SC)     : 0.00137251 +/- 0.302165

SO          : -3.09503e-07 +0.00026077/-0.000263376 or +/- 0.00026595

SO          : -3.09503e-07 +/- 0.000262073

log(GL)     : 0.565666 +1.60512/0 or +/- 4.05036

log(GL)     : 0.565666 +/- 4.05036

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.3022)*(-5.821e-07*(bbcbb-( -2229 +/- -450.8))+(8.711e-07*(zdcwnk))+(1.702e-06*(zdcenk)))

  with GL =  1.76 +/-  7.13

--> Different scalers set than first pass. 
--> Trying with the same scalers set as in the previous pass.

2-scalers fit with scalers: bbce, zdcenk

USING STDDEV sc => 0.000387 :: gapf => 0.005014

 

SpaceCharge & GridLeak fit results {scaler: (-1.16577e-08*(bbce+100493))+(4.07732e-06*(zdcenk))}:

g2          : 0.592923 +5.10785/-5.34954 or +/- 0.656412

g2          : 0.592923 +/- 5.2287

log(g5)     : 3.31996 +0/0 or +/- 0.123703

log(g5)     : 3.31996 +/- 1.474

log(SC)     : -0.000755798 +0.476578/-0.932444 or +/- 0.0842553

log(SC)     : -0.000755798 +/- 0.704511

SO          : -6.13492e-07 +0.000517458/-0.000521204 or +/- 0.000325331

SO          : -6.13492e-07 +/- 0.000519331

log(GL)     : 0.42656 +2.10181/0 or +/- 2.79657

log(GL)     : 0.42656 +/- 2.79657

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.7045)*(-1.165e-08*(bbce-(-1.004e+05 +/- -4.455e+04))+(4.074e-06*(zdcenk)))

  with GL =  1.53 +/-  4.28

--> Different Second pass scaler set results are better than the first pass scalers set, with larger GL value. 
--> Still the convergence is not satisfactory.
---------------------------------------------------------------------------------------------------
--------------------------------------------------------------------------------------------------- 
   Pass 03

--> 2-scalers fit with fix GL =0.12 result from second pass is propagated.   
sc = (1.0 +/- 0.4766)*(-1.677e-07*(bbcbb-( -5421 +/-  -2472))+(7.275e-07*(zdcwnk)))

  with GL =  0.12 +/-  0.00

--> 1- and 2-scalers fit didn't converge. 
     --Fit failed for gapf, err = 4
--> 3-scalers fit converged, the result is below for fixed Gl = 0.12. 

USING STDDEV sc => 0.000202 :: gapf => 0.005856

 

SpaceCharge & GridLeak fit results {scaler: (-6.65651e-07*(bbcbb+2233.2))+(9.03231e-07*(zdcwnk))+(2.07901e-06*(zdcenk))}:

g2          : -0.39951 +11.5439/-10.5267 or +/- 5.51963

g2          : -0.39951 +/- 11.0353

log(g5)     : 3.84785 +0/-3.41286 or +/- 1.76355

log(g5)     : 3.84785 +/- 1.76355

log(SC)     : 0.00555139 +0.263256/-0.366991 or +/- 0.307256

log(SC)     : 0.00555139 +/- 0.315123

SO          : -2.13474e-06 +0.000280723/-0.000280712 or +/- 0.000284353

SO          : -2.13474e-06 +/- 0.000280718

 

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.3151)*(-6.694e-07*(bbcbb-( -2230 +/- -421.7))+(9.083e-07*(zdcwnk))+(2.091e-06*(zdcenk)))

  with GL =  0.12 +/-  0.00

 

 
-------------------------------------------------------------------------------------------------
--> Macro fixed for failing gapf error (by Gene).
--> Gene's study: 

Gene tried running with 2 scalers on my data and found the same thing as above 3-scaler fit shows: A negative coefficient for bbcbb and a positive coefficient for zdcwnk. It's possible this makes physical sense, but he was not confident about it. So running with debug=2 to see how far off other possibilities were, and he found these for having 
 
root -l hist*/* 'Calib_SC_GL.C+(0,0,-2,2)'
...
RMS = 0.000336308 :: EV = 3.9635e-14 :: 
  sc = (5.17877e-07*(zdcwnk-2061.14))
...
RMS = 0.000330844 :: EV = 3.69704e-14 :: 
  sc = (4.73346e-07*(zdcwnk-2536.85))+(2.8871e-07*(zdcenk))
...
RMS = 0.00032088 :: EV = 1.56541e-14 :: 
  sc = (-1.46686e-07*(bbcbb+6462.78))+(6.73408e-07*(zdcwnk))
 
Finding:: the RMS isn't much worse for the 1-scaler zdcwnk, but it doesn't seem to fit happily. So the 2-scaler zdcwnk:zdcenk also is pretty close in RMS and has positive coefficients. His recommendation was running this fit and using the results:
root -l hist*/* 'Calib_SC_GL.C+(0,0,"zdcwnk:zdcenk")'






sc = (1.0 +/- 0.5841)*(4.769e-07*(zdcwnk-(  2530 +/-    979))+(2.907e-07*(zdcenk)))
  with GL =  0.12 +/-  0.00

With this, 14.5 GeV calibration work is wrapped up! 

»