SNR convergence

Go to coinc/root_scripts directory. Most of the libraries are loaded by default once root is run.

By hand, add:

root [] .L estimatePSF.C+
root [] .L scripts_buildTradeOffFromIterationData.C

Then, load the data headers.

root [] tagCollection tg
root [] tg.parseFile("/data0/studen/data/imageRec/layout/setup2/tagMapSi.txt")

I find it helpful to extract just the entry I want before loading the full data set.

root [] tagCollection tg1
root [] tg1.add(tg.at(0),tg.at(tg.at(0))

Replace 0 with a different number should you want a different entry. There are routines that extract according to dataSpec entry - just set the fields you want to filter on and do

root [] tg.select(tg1,ds)

which selects all appropriate to tg1. Then read the data

root [] iterationCollection ic
root [] ic.init(&tg1)

The SNR extracting routines are in a separate file, so:

root [] .L scripts_getSNR.C+

To get SNR as it develops through iteration at a fixed resolution of eta for at a source is=(0,1) you should do

root [] TGraph gr
root [] int is=0 //set the source id
root [] double eta=8 //set the resolution
root [] getSNR(gr,ic.at(0),ic.at(ic.at(0)),eta,is)

Wait for maybe a minute.