GPC299 Code QA

0 : Basic Information
GPC299 Code in CVS : /offline/paper/psn0750/

This code QA work follows the procedure described in the README file.

1 : Code Compilation

- 2021/3/16
Executed 'cons' command following the README in the repository using the SL18f, SL19c, and SL20d as indicated in the xml files.
Failed cons as the following const functions push back some values to member variable vectors.

.sl73_gcc485/obj/StRoot/StSmallSystemTreeMaker/StSmallSystemTreeMaker.cxx: In member function 'void StSmallSystemTreeMaker::calQxQy(StMuTrack*) const':
.sl73_gcc485/obj/StRoot/StSmallSystemTreeMaker/StSmallSystemTreeMaker.cxx:991:37: error: passing 'const DoubleVec {aka const std::vector<double>}' as 'this'

.sl73_gcc485/obj/StRoot/StSmallSystemTreeMaker/StSmallSystemTreeMaker.cxx: In member function 'void StSmallSystemTreeMaker::calQxQy(StPicoTrack*, TVector3) const':
.sl73_gcc485/obj/StRoot/StSmallSystemTreeMaker/StSmallSystemTreeMaker.cxx:1019:37: error: passing 'const DoubleVec {aka const std::vector<double>}' as 'this'

See cons_SL18f_20210316.txt  cons_SL19c_20210316.txt  cons_SL20d_20210316.txt  for the entire error messages.

- 2021/3/18
Shengli removed unnecessary functions that caused compilation errors and some compilation warning.
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/StSmallSystemTreeMaker/StSmallSystemTreeMaker.cxx?r1=1.1&r2=1.2
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/StSmallSystemTreeMaker/StSmallSystemTreeMaker.h?r1=1.1&r2=1.2

Confirmed that cons command worked OK and did not show any error or warning messages. See cons_SL18f_20210318.txt  cons_SL19c_20210318.txt  cons_SL20d_20210318.txt.

2 : Mini-Tree Production and Analysis

- 2021/3/18
Spotted some inconsistencies in SL version between the PicoDst production and mini-tree production job in the xml files. Shengli fixed as follows:
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/xml/dAu200.xml?r1=1.1&r2=1.2
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/xml/HeAu200.xml?r1=1.1&r2=1.2
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/xml/pAussdmb.xml?r1=1.1&r2=1.2
https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/xml/ppssdmblong2_PicoDst.xml?r1=1.1&r2=1.2

Confirmed that cons command works fine with the updated SL versions.

- 2021/3/23
Given the limited institutional space on the pwg-disk, we decided to produce small number of mini-DSTs using exactly the same set of PicoDsts, and to compare those produced by Shengli and Takahito.

The data sets used for dAu and HeAu test miniDst productions:
dAu : /star/u/tdrk/GPC/299/dAu200new/17140003.list
HeAu : /star/u/tdrk/GPC/299/HeAu200new/15180040.list

Test miniDsts produced by Shengli:
dAu : /star/u/slhuang/pwg/small/paper_code/lib/dAu200/SL18f/test.root
HeAu : /star/u/slhuang/pwg/small/paper_code/lib/HeAu200/SL18h/test.root

Test miniDsts produced by Takahito:
dAu : /star/u/tdrk/GPC/299/dAu200new/test_SL18f.root
HeAu : /star/u/tdrk/GPC/299/HeAu200new/test_SL18h.root

Confirmed that the contents in Shengli's and Takahito's mini-trees are identical. See
dAuSL18f_Show.txt,
dAuSL18f_Show_Shengli.txt, HeAuSL18h_Show_Shengli.txt, HeAuSL18h_Show.txt for details.

We analyzed mini-trees following the procedure described in 2) in the README file. From this point forward we used Shengli's mini-trees for v2 and v3 measurements that can be found at:
HeAu200 : /star/u/slhuang/pwg/HeAu/picoDst2/
dAu200 : /star/scratch/slhuang/dAu200
pAussdmb : /star/scratch/slhuang/pAussdmb/
pAuhimult : /star/scratch/slhuang/pAuhimult/
ppssdmb : /star/scratch/slhuang/ppssdmblong2/

All tree analysis jobs successfully finished without error.

 3 : Vn Calculations and Systematic Uncertainty

We extracted v2 and v3, and their systematic uncertainties from the output of the mini-tree analysis, following the procedure described in 3) in the README file.
The results are consistent till the last a few digits between Shengli's analysis and this code QA work for both results using the standard TPC centrality and BBC centrality.

See the pages below for the direct comparison.

Shengli's results:
TPC centrality: https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/GPC/woc0c1/
BBC centrality: https://www.star.bnl.gov/cgi-bin/protected/viewvc.cgi/cvsroot/offline/paper/psn0750/GPC/BBC/woc0c1/

Code QA result:
TPC centrality: https://www.star.bnl.gov/protected/lfsupc/tdrk/GPC/299/CodeQA/datatable/final/
BBC centrality: https://www.star.bnl.gov/protected/lfsupc/tdrk/GPC/299/CodeQA/datatable/bbc/


Confirmed that PA's analysis results can be successfully reproduced using the codes, macros in the CVS area.

4 : PRL Figures

Following the guidance in this README file, we could successfully reproduce the plots in the PRL paper (Phys. Rev. Lett. 130, 242301 (2023)).






5 : PRC Figures

Following the guidance in this README file, we could also successfully reproduce the plots in the PRC paper draft.