rwmichaels / trkcorrtest Goto Github PK

TrkCorrTest -- to test HrsTrkCorr July 2021 Robert Michaels

You'll need to get HrsTrkCorr from the separate repository.

/rootfiles directory with root files sieve slit runs 2239,2240,2241,21363,21364,21365 std Qsq runs 1983, 1984 The files you need are at JLab on the ifarm. ifarm:/w/halla-scifs17exp/parity/disk1/bob/trkcorr_rootfiles

AfileL.root AfileR.root usually a link

MoveHoles.C Test of HrsTrkCorr. Read a sieve slit ROOT file and move the holes to agree with the data.

MoveTracks.C Test of HrsTrkCorr. Read a sieve slit ROOT file and move the tracks to agree with the holes.

PlotResid.C Test of HrsTrkCorr. Plot the residuals found

QsqDistort.C Read the ROOT file of a Qsq run and plot the Qsq with and without track distortions from HrsTrkCorr

trk_scriptN.C N = 1, 2, ... Example of how to use HrsTrkCorr in a ROOT script.

holefiles: holes* = where the sieve holes are in tg_th vs tg_ph residuals* = residuals between track and holes

Question: How do we know that HrsTrkCorr is correct ? Answer: Use root macros MoveHoles.C and MoveTracks.C to verify that tracks lie on the holes for the sieve slit runs. Also, some vetting of the class was done with QsqDistort.C

Notes on the work flow -- example on ifarm

git clone https://[email protected]/rwmichaels/TrkCorrTest.git

delete contents of rootfiles, make links. rootfiles -> /w/halla-scifs17exp/parity/disk1/bob/trkcorr_rootfiles

For convenience, copy stuff here cp ../HrsTrkCorr/HrsTrkCorr.h . cp ../HrsTrkCorr/HrsTrkCorr.C .

AfileL.root or AfileR.root is a link AfileR.root -> ./rootfiles/prexRHRS_21364_siyu.root

Run MoveHoles.C, MoveTracks.C, or QsqDistort.C as per instructions at the top of those files.

control the variable Which_Spectrometer (=L or R) in the macros. control ixbeam (=0,1,2), same macros AfileL(R).root is a link; make it point to the file corresponding to ixbeam.