root [0] .x loadLibs_rodMeasurement.C
root [1] va_data va;
root [2] va.parse("/data0/studen/data/imageRec/studies/general/match_harmonized/Fischer_final.cfg")
root [] TH2D h2;
root [] getImage(h2,va,0);

root [6] va.img.set_value<int>("dirlistnData",100000,0)
root [8] va.img.set_value<double>("dirlistbeta",1e-8,0)
root [9] idList fid;iList fi
root [10] va.img.set_value<string>("dirlistFlag","[100:0.95][f0:0.05]",0)
#set Gate for simulation
root [11] va.img.set_value<string>("dataSource","um",i)
root [12] getData(fid,va,0);
root [14] makeList(fi,fid,va,0)
root [17] fi.front().Draw("COLZ")
#combine samples
root [] variationData vd;
root [] assemble(vd,fi)
root [] vd.GetMeanImage()->Draw("COLZ")
#get mean
root [] evalData d
# i is image id, j is rod id
root [] evalRod(d,fi,va,i,j)