TPC SCGL Calibration: Run 20 AuAu 9.2 GeV

Updated on Wed, 2021-09-08 13:54. Originally created by pokhrel on 2021-04-26 07:54.

This dataset has three periods. The plan is to do the calibration for the first period and test for the rest of the periods whether the parameter works.
--> Out of 21, 18 DAQ files were restored from the period 1 . The DAQ list is attached. (daq_p1.txt for the first pass and daq_p1_02.txt for the second pass!)
--> Working directory: /star/u/pokhrel/TPC_Cal/Run20/9p2GeV/Period1/
--> For the first pass, since the CINT files didn't have any suggested parameters, I started with the random scaler "vpdx" with "row.fullFieldB=0".
------------------------------------------------------------------------------------
First Pass Results:

USING STDDEV sc => 0.000193 :: gapf => 0.014826

SpaceCharge & GridLeak fit results {scaler: bbcw}:

g2 : -0.00291744 +1.43826/-1.50864 or +/- 0.272177

g2 : -0.00291744 +/- 1.47345

log(g5) : 4.7592 +0/-1.27487 or +/- 2.63544

log(g5) : 4.7592 +/- 2.63544

log(SC) : -16.4278 +0.155302/-0.183563 or +/- 0.167541

log(SC) : -16.4278 +/- 0.169433

SO : 6576.45 +3478.86/-3486.52 or +/- 3464.11

SO : 6576.45 +/- 3482.69

log(GL) : 2.13237 +1.25114/0 or +/- 2.48517

log(GL) : 2.13237 +/- 2.48517

*** FINAL CALIBRATION VALUES: ***

sc = (7.337e-08 +/- 1.243e-08) * ((bbcw) - ( 6576 +/- 3483))

with GL = 8.43 +/- 20.96

*** The following calibration values may not be trusted at this time... ***

*** Try the following calibration values: ***

sc = 3.229e-06 * ((bbcw) - ( 6575)) with GL = 1.00

sc = 3.229e-06 * ((bbcw) - ( 6575)) with GL = 0.50

sc = 3.229e-06 * ((bbcw) - ( 6575)) with GL = 0.00

--> The GL values looks awkward. Tried fixing GL =0.5.

USING STDDEV sc => 0.000193 :: gapf => 0.014825

SpaceCharge & GridLeak fit results {scaler: bbcw}:

g2 : -0.0628521 +22.6922/-23.7702 or +/- 8.36865

g2 : -0.0628521 +/- 23.2312

log(g5) : 10.581 +0/0 or +/- 167.963

log(g5) : 10.581 +/- 1.94591

log(SC) : -16.3617 +0.156544/-0.185687 or +/- 0.167538

log(SC) : -16.3617 +/- 0.171116

SO : 6565.19 +3506.12/-3506.15 or +/- 3466.23

SO : 6565.19 +/- 3506.13

*** FINAL CALIBRATION VALUES: ***

sc = (7.838e-08 +/- 1.341e-08) * ((bbcw) - ( 6565 +/- 3506))

with GL = 0.50 +/- 0.00

*** The following calibration values may not be trusted at this time... ***

*** Try the following calibration values: ***

sc = 3.229e-06 * ((bbcw) - ( 6577)) with GL = 1.00

sc = 3.229e-06 * ((bbcw) - ( 6577)) with GL = 0.50

sc = 3.229e-06 * ((bbcw) - ( 6577)) with GL = 0.00

--> The fit parameters looks ok. I got a message saying "parameters cannot be trusted".

-----------------------------------------------------------------------------------
Pass 02:
--> Added more daq (26 total) files for this pass to include more data. Ntuple for a couple of runs (3) were empty in the first pass.
--> Propagated first pass results, the result with a "bbcw" scaler (not the suggested results!).

1 Scaler Fit:

USING STDDEV sc => 0.000185 :: gapf => 0.013377

SpaceCharge & GridLeak fit results {scaler: vpdw}:

g2 : -0.00334557 +0/0 or +/- 0.103188

g2 : -0.00334557 +/- 3.04257e-06

log(g5) : -0.204077 +0/0 or +/- 0.00165836

log(g5) : -0.204077 +/- 6.70551e-06

log(SC) : -15.6968 +0.352571/-0.558781 or +/- 0.389601

log(SC) : -15.6968 +/- 0.455676

SO : -186.906 +2054.3/-1979.62 or +/- 2149.16

SO : -186.906 +/- 2016.96

log(GL) : -1.72859 +0/0 or +/- 0.00162189

log(GL) : -1.72859 +/- 2.82952e-05

*** FINAL CALIBRATION VALUES: ***

sc = (1.524e-07 +/- 6.944e-08) * ((vpdw) - (-186.9 +/- 2017))
with GL = 0.18 +/- 0.00

2Scaler:

USING STDDEV sc => 0.000120 :: gapf => 0.013315

SpaceCharge & GridLeak fit results {scaler: (9.40161e-08*(vpde-4620.44))+(1.24631e-07*(vpdw))}:

g2 : -1.815e-05 +0/0 or +/- 0.0729546

g2 : -1.815e-05 +/- 0.69129

log(g5) : 10.2905 +0/0 or +/- 153.43

log(g5) : 10.2905 +/- 1.62786

log(SC) : 0.0101588 +0.110423/-0.124272 or +/- 0.11593

log(SC) : 0.0101588 +/- 0.117348

SO : -3.04895e-08 +0.000169115/-0.00016912 or +/- 0.000167851

SO : -3.04895e-08 +/- 0.000169118

log(GL) : 0.857185 +0/0 or +/- 831.39

log(GL) : 0.857185 +/- 3.56602

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.1173)*(9.498e-08*(vpde-( 4620 +/- 1799))+(1.259e-07*(vpdw)))

with GL = 2.36 +/- 8.40

--> Large GL value.
--> Fixing GL=0.5 results in distorted fit.

{USING STDDEV sc => 0.044329 :: gapf => 0.013314

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 3.502)*(1.033e-07*(vpde-( 5533 +/- 6.084e+05))+(1.369e-07*(vpdw)))

with GL = 0.50 +/- 0.00
}

--> Though the 1-scaler fit parameter seems reasonable, the scaler is different than in the first pass.
--> Forcing the scaler as "bbcw" results in the following:

"bbcw" Fit:
--> Free GL fit:

USING STDDEV sc => 0.000251 :: gapf => 0.013268

*** FINAL CALIBRATION VALUES: ***

sc = (6.588e-08 +/- 1.648e-08) * ((bbcw) - ( 5953 +/- 5007))

with GL = 9.52 +/- 31.51

--> Large GL value.
--> Fixing GL = 0.5.

USING STDDEV sc => 0.000251 :: gapf => 0.013268

SpaceCharge & GridLeak fit results {scaler: bbcw}:

g2 : 1.17652e-06 +0/0 or +/- 0.0533819

g2 : 1.17652e-06 +/- 0.680944

log(g5) : 10.6174 +0/0 or +/- 230.353

log(g5) : 10.6174 +/- 1.18085

log(SC) : -16.4614 +0.220263/-0.282904 or +/- 0.244546

log(SC) : -16.4614 +/- 0.251583

SO : 5949.9 +5023.1/-5023.15 or +/- 4985.19

SO : 5949.9 +/- 5023.12

*** FINAL CALIBRATION VALUES: ***

sc = (7.094e-08 +/- 1.785e-08) * ((bbcw) - ( 5950 +/- 5023))

with GL = 0.50 +/- 0.00

--> Although the fit converged well, the three data points seems off (outliers!).

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

Final Pass:
--> "Bus Error" was occurred on most of the calibration jobs.
--> Gene looked at the log files and found that only about 1-2% events have passed the default 5 cm cut of VPD and primary vertex position.
--> I eliminate this selection assigning negative value to the vtxVpdAgree changing vtxVpdAgree(5.0) to vtxVpdAgree(-5.0) in the constructor of StSpaceChargeEbyEMaker.
--> This fixes the calibration jobs issues.

--> Result for "bbcw" scalar:

USING STDDEV sc => 0.000321 :: gapf => 0.009908

SpaceCharge & GridLeak fit results {scaler: bbcw}:

g2 : 8.08112e-05 +14.5512/-10.7044 or +/- 2.7401

g2 : 8.08112e-05 +/- 12.6278

log(g5) : 10.4937 +0/0 or +/- 215.604

log(g5) : 10.4937 +/- 1.80334

log(SC) : -16.3239 +0.184369/-0.226234 or +/- 0.201086

log(SC) : -16.3239 +/- 0.205302

SO : 7519.9 +5591.24/-5591.22 or +/- 5553.1

SO : 7519.9 +/- 5591.23

*** FINAL CALIBRATION VALUES: ***

sc = (8.14e-08 +/- 1.671e-08) * ((bbcw) - ( 7520 +/- 5591))

with GL = 0.50 +/- 0.00

--> For the "VPDW" scalar and fixing GL to 0.5. Free GL leads to very large GL.

USING STDDEV sc => 0.000187 :: gapf => 0.009763

SpaceCharge & GridLeak fit results {scaler: vpdw}:

g2 : 0.000722362 +0/0 or +/- 0.557945

g2 : 0.000722362 +/- 0.953358

log(g5) : 10.4076 +0/0 or +/- 176.977

log(g5) : 10.4076 +/- 1.89872

log(SC) : -15.642 +0.114934/-0.129916 or +/- 0.120982

log(SC) : -15.642 +/- 0.122425

SO : 1899.71 +1646.32/-1646.3 or +/- 1635

SO : 1899.71 +/- 1646.31

*** FINAL CALIBRATION VALUES: ***

sc = (1.61e-07 +/- 1.971e-08) * ((vpdw) - ( 1900 +/- 1646))

with GL = 0.50 +/- 0.00

--> all the fit parameters are reasonable. All data points are all aligned with low STD value of straight line fit.
--> Results for the "vpdw" are much better than "bbcw".

--> Final calibration parameters for Run 20, 9.2 GeV datasets?

sc = (1.61e-07 +/- 1.971e-08) * ((vpdw) - ( 1900 +/- 1646))

with GL = 0.50 +/- 0.00

Next Step:
--> Period 2 and 3 SC/GL calibration.
--> First to check whether Period 1 parameters works for the rest of the periods.
--> Is there a way that I can directly use these parameters to check for the rest of the periods?
------------------------------
Period 2
----------------------------
--> Final calibration parameters are propagated to the period 2 dataset. Since the experimental configuration is same as in the period 1, it is expected that the parameters should fit to this dataset as well.
--> Fit failed for the single scalar, not as expected.
--> Rechecked: Datasets are all fine, the produced histograms and the ntuples are all fine. Enough events have passed the vertex requirements.
--> The issue was tracked by Gene. His initial guess was poor SC GL guess and the poorly set offsets while fitting all three datasets.
--> Updated calibration macro can be found here (from Gene): ~genevb/public/FOR_BABU/Calib_SC_GL.C
--> Still with the updated calibration macro, the fit fails for the single rich scalar.
--> Try fitting with two rich scalar:
--> Fit converged, with the results below:

USING STDDEV sc => 0.000231 :: gapf => 0.011079

SpaceCharge & GridLeak fit results {scaler: (1.73619e-06*(bbcx-365.97))+(2.52128e-07*(bbcbb))}:

g2 : 1.40973e-05 +0/0 or +/- 0.117109

g2 : 1.40973e-05 +/- 0.363439

log(g5) : 10.6142 +0/0 or +/- 357.543

log(g5) : 10.6142 +/- 1.24997

log(SC) : 0.013768 +0.38013/-0.620822 or +/- 0.458982

log(SC) : 0.013768 +/- 0.500476

SO : -3.33476e-08 +0.000324055/-0.000324055 or +/- 0.0003216

SO : -3.33476e-08 +/- 0.000324055

log(GL) : -5.79076 +0/0 or +/- 253530

log(GL) : -5.79076 +/- 3.93304

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 0.5005)*(1.76e-06*(bbcx-( 366 +/- 186.6))+(2.556e-07*(bbcbb)))

with GL = 0.00 +/- 0.01

--> Not a well constrained fit, needs more iteration.
--> Calibration jobs submitted using the 2 scalers parameter for the third pass.

Pass03:
------------------

SpaceCharge fit results {scaler: (1.77167e-06*(bbcx-338.129))+(2.85315e-07*(bbcbb))}:

Finished outlier removal: 0 (of 72) data points removed (2 iterations)

SO : 0.000287816 +0.000602048/-0.00060204 or +/- 0.000593666

SO : 0.000287816 +/- 0.000602044

log(SCe) : 3.75915 +1.7865/0 or +/- 4.9765

log(SCe) : 3.75915 +/- 4.9765

log(g5) : 3.72052 +0/-1.81436 or +/- 5.11521

log(g5) : 3.72052 +/- 5.11521

VAR_all = 0.0003722 for (1.77167e-06*(bbcx-338.129))+(2.85315e-07*(bbcbb))

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

GridLeak fit results {scaler: (1.77167e-06*(bbcx-338.129))+(2.85315e-07*(bbcbb))}:

Finished outlier removal: 0 (of 72) data points removed (2 iterations)

g2 : 0.0990082 +0.470378/-0.488636 or +/- 0.556109

g2 : 0.0990082 +/- 0.479507

g1/g2 = SCxGL : 0.0518263 +0.242776/-0.253273 or +/- 0.286651

g1/g2 = SCxGL : 0.0518263 +/- 0.248025

GLO : 0.344483 +1.60818/-1.67862 or +/- 1.89931

GLO : 0.344483 +/- 1.6434

*** FIRST PASS CALIBRATION VALUES: ***

sc = (1.0 +/- 7.141)*(1.839e-06*(bbcx-( 500.6 +/- 339.8))+(2.962e-07*(bbcbb)))

with GL = 0.05 +/- 0.35

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

USING STDDEV sc => 0.003736 :: gapf => 0.009753

SpaceCharge & GridLeak fit results {scaler: (1.77167e-06*(bbcx-338.129))+(2.85315e-07*(bbcbb))}:

g2 : -1.90731 +25.4284/-29.9418 or +/- 35.3008

g2 : -1.90731 +/- 27.6851

log(g5) : 8.49207 +0/0 or +/- 1755.42

log(g5) : 8.49207 +/- 1.90827

log(SC) : 0.234951 +1.77453/0 or +/- 4.97462

log(SC) : 0.234951 +/- 4.97462

SO : -6.54275e-05 +0.00385229/-0.00388822 or +/- 0.00391231

SO : -6.54275e-05 +/- 0.00387026

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 4.975)*(2.241e-06*(bbcx-( 301.2 +/- 2185))+(3.609e-07*(bbcbb)))

with GL = 0.50 +/- 0.00

Fig: space charge fit
--> GL value is fixed to 0.5. (Without fixing it, GL turns out to be awkwardly large with very large error)
--> The space charge fits for all datasets shows the linear trend (brown points and the gray fitted line).
--> The combine fit, SC and GL, is now off than the space charge fit (black solid line). The overall combined fit results comes up with the higher STD value.

------------------------------
Period 3
------------------------------

Pass01:

USING STDDEV sc => 0.004047 :: gapf => 0.024774

SpaceCharge & GridLeak fit results {scaler: vpdw}:

g2 : -5.0968e-05 +0/0 or +/- 0.0861533

g2 : -5.0968e-05 +/- 1.22374

log(g5) : 10.3674 +0/0 or +/- 490.349

log(g5) : 10.3674 +/- 1.10985

log(SC) : -16.5241 +0.636063/-2.20426 or +/- 0.884598

log(SC) : -16.5241 +/- 1.42016

SO : 49132.2 +86169.9/-86141.2 or +/- 85676.1

SO : 49132.2 +/- 86155.6

*** FINAL CALIBRATION VALUES: ***

sc = (6.663e-08 +/- 9.462e-08) * ((vpdw) - (4.913e+04 +/- 8.616e+04))

with GL = 0.50 +/- 0.00

Pass02
----------------
Using initial value for SCe of 34599.3

SpaceCharge fit results {scaler: (9.00571e-08*(bbce-3841.43))+(-8.50798e-09*(vpdw))}:

Finished outlier removal: 0 (of 95) data points removed (2 iterations)

SO : 0.000252949 +0.00127277/-0.00127277 or +/- 0.00125933

SO : 0.000252949 +/- 0.00127277

log(SCe) : 3.76184 +0/0 or +/- 38.3662

log(SCe) : 3.76184 +/- 0.189351

log(g5) : 3.74888 +0/0 or +/- 41.137

log(g5) : 3.74888 +/- 0.254881

VAR_all = 0.000369 for (9.00571e-08*(bbce-3841.43))+(-8.50798e-09*(vpdw))

---------------
GridLeak fit results {scaler: (9.00571e-08*(bbce-3841.43))+(-8.50798e-09*(vpdw))}:

Finished outlier removal: 0 (of 95) data points removed (2 iterations)

g2 : 0.669425 +1.24441/-1.31292 or +/- 1.58089

g2 : 0.669425 +/- 1.27867

g1/g2 = SCxGL : 0.331559 +0.852084/-0.972241 or +/- 1.05577

g1/g2 = SCxGL : 0.331559 +/- 0.912163

GLO : 0.0197642 +0.050364/-0.0572968 or +/- 0.0226542

GLO : 0.0197642 +/- 0.0538304

*** FIRST PASS CALIBRATION VALUES: ***

sc = (1.0 +/- 0.3191)*(9.053e-08*(bbce-( 6650 +/- 1.413e+04))+(-8.553e-09*(vpdw)))

with GL = 0.33 +/- 0.92

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

USING STDDEV sc => 0.053298 :: gapf => 0.013825

SpaceCharge & GridLeak fit results {scaler: (9.00571e-08*(bbce-3841.43))+(-8.50798e-09*(vpdw))}:

g2 : 0.94704 +2.50455/-2.49695 or +/- 3.16688

g2 : 0.94704 +/- 2.50075

log(g5) : 10.5071 +0/0 or +/- 3597.84

log(g5) : 10.5071 +/- 1.01545

log(SC) : -1.42096 +1.72455/0 or +/- 5.37242

log(SC) : -1.42096 +/- 5.37242

SO : 0.0282284 +0.128904/-0.140524 or +/- 0.150068

SO : 0.0282284 +/- 0.134714

*** FINAL CALIBRATION VALUES: ***

sc = (1.0 +/- 5.372)*(2.175e-08*(bbce-(3.173e+05 +/- 1.496e+06))+(-2.055e-09*(vpdw)))

with GL = 0.50 +/- 0.00

Fig: space charge fit results

--> Pass03 results shows huge offset between data and the fit.
Pass04:

USING STDDEV sc => 0.000278 :: gapf => 0.011547

SpaceCharge & GridLeak fit results {scaler: bbce}:

g2 : -0.419745 +2.18697/-2.67082 or +/- 3.11152

g2 : -0.419745 +/- 2.4289

log(g5) : 10.5845 +0/0 or +/- 240.432

log(g5) : 10.5845 +/- 1.0362

log(SC) : -16.2642 +0.208769/-0.264185 or +/- 0.230887

log(SC) : -16.2642 +/- 0.236477

SO : -800.615 +4558.74/-4558.62 or +/- 4533.58

SO : -800.615 +/- 4558.68

*** FINAL CALIBRATION VALUES: ***

sc = (8.64e-08 +/- 2.043e-08) * ((bbce) - (-800.6 +/- 4559))

with GL = 0.50 +/- 0.00

Pass 4 fit result is in good agreement with the data.
Trying one more iteration with these fit results.

Pass05:

Period 1 and Period 2 simultaneous fit with VPDW:
Gene found that the period 1 and 2 fits well with the common rich scaler "VPDW". I checked this with the final pass from Period 1 (Pass03) and the first pass from period 2 (which has period 1 final parameters as an input). There is a good agreement between data and the fit. Below is the combined fit results:

USING STDDEV sc => 0.000266 :: gapf => 0.012187

SpaceCharge & GridLeak fit results {scaler: vpdw}:

g2 : -6.59881e-06 +0/0 or +/- 0.272729

g2 : -6.59881e-06 +/- 0.855691

log(g5) : 10.4329 +0/0 or +/- 225.03

log(g5) : 10.4329 +/- 1.24368

log(SC) : -15.6022 +0.200843/-0.251582 or +/- 0.221656

log(SC) : -15.6022 +/- 0.226213

SO : 3134.41 +2252.36/-2252.56 or +/- 2244.63

SO : 3134.41 +/- 2252.46

*** FINAL CALIBRATION VALUES: ***

sc = (1.675e-07 +/- 3.789e-08) * ((vpdw) - ( 3134 +/- 2252))

with GL = 0.50 +/- 0.00

pokhrel's blog
Login or register to post comments

The STAR experiment

pokhrel's blog

User login

Navigation

TPC SCGL Calibration: Run 20 AuAu 9.2 GeV