Page History

Make tuples using the instructions here. You can stop after you create the tuple directly from the text file, cuts will be applied in the analysis code below.

In order to analyze, you need to clone the HPS-CODE repository.

git clone https://github.com/JeffersonLab/HPS-CODE

Run the following command on the tuple file.

python HPS-CODE/ANALYSIS/SVT/MakeVertexResPlots.py -c <output file base name> <input tuple file>

Options are as follows

    -t: tuple name (default "ntuple")

    -m: number of mass bins (default 10)

    -e: beam energy (default 1.056 GeV)

    -z: position of target in z (default 0.5 mm)

    -c: apply cuts (default false)

Cuts can manually be changed by editing the "cuts" variable directly in the script

Output:

<output file base name>.pdf - vertex resolution as a function of mass as well as the individual fits

<output file base name>.root - same as PDF but with root histograms and graphs

For SVT NIM Paper

File:

Cuts:

isPair1

eleClY*posClY<0

uncP>0.8*ebeam

max(eleTrkChisq,posTrkChisq)<60

eleP<0.75*ebeam

bscChisq<10

bscChisq-uncChisq<5

uncP<1.15*ebeam

max(eleMatchChisq,posMatchChisq)<10

abs(eleClT-eleTrkT-43)<4

abs(posClT-posTrkT-43)<4

abs(eleClT-posClT)<2

eleHasL2

posHasL2

eleHasL1

posHasL1

min(eleMinPositiveIso+0.5*(eleTrkZ0+{1}*elePY/eleP)*sign(elePY),posMinPositiveIso+0.5*(posTrkZ0+{1}*posPY/posP)*sign(posPY))>0

abs(eleP-posP)/(eleP+posP)<0.4

posTrkD0+ztarg*posPX/posP<1.5

Vertex Plots

Root File