O'Grady, Paul Christopher   
Fri 2/24/2017 12:16 PM
To:   Liang, Mengning;  
Cc:   Dubrovin, Mikhail;  
Suggested Meetings

Hi Meng,

I talked with Mikhail about the fiber-diffraction “sprint” to get 
results for Germany that you have strong confidence in.  We’ll come by 
your office around 4 today to launch.

We’ll also touch base Monday, Tuesday, Wednesday, Thursday next week at 
3pm to make sure we’re headed in the right direction.  Those are allowed
 to be brief (more important to do real work).

chris

    alg_arc = PyAlgos(windows=winds_arc, mask=mask_arc, pbits=2)
    alg_arc.set_peak_selection_pars(npix_min=20, npix_max=1000, amax_thr=0, atot_thr=2000, son_min=5)

    alg_equ = PyAlgos(windows=winds_equ, mask=mask_equ, pbits=0)
    alg_equ.set_peak_selection_pars(npix_min=5, npix_max=1000, amax_thr=0, atot_thr=600, son_min=5)

    for i, evt in enumerate(ds.events()) :
        nda_data = det.calib(evt, cmpars=(5,50))

        if nda_data is not None :
            nda =  np.array(nda_data, dtype=np.float32, copy=True)

            if BKG_MODE == BKG_MODE_NORMAL : 
                nda = subtract_bkgd(nda, nda_bkgd, mask=nda_smask, winds=winds_bkgd, pbits=0)

            peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=160, rank=8, r0=10, dr=2)
            peaks_equ = alg_equ.peak_finder_v4r2(nda, thr_low=60, thr_high=100, rank=4, r0=6, dr=2)

sit_setup ~yoon82/ana2 # some problem fixed in this local release
psocake -e cxif5315 -r 169 -d CxiDs2.0:Cspad.0 -n 468 -o psocake

python cxif5315/proc-cxif5315-r0169-data-pfvn-2017-02-27.py

Conclusion: within mask-permitted region analysis script finds the same peaks as psocake

r0=6  dr=2
 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 0
 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0
 0 0 1 1 1 0 0 0 0 0 0 0 1 1 1 0 0
 0 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 1 1 1 0 0 0 0 0 + 0 0 0 0 0 1 1 1
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 0
 0 0 1 1 1 0 0 0 0 0 0 0 1 1 1 0 0
 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0
 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 0
 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
Number of pixels to estimate background = 88

r0=10  dr=2
 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0
 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0
 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0
 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0
 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0
 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0
 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 1 1 1 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 1 1 1
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0
 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0
 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0
 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0
 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0
 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0
 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0
 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
Number of pixels to estimate background = 136

Liang, Mengning   
Tue 2/28/2017 4:09 PM
To:   Dubrovin, Mikhail;  

148 – double hit in equator – 3 peaks
162 – no peaks
525 – double hit in equator - 6 peaks
1456 – should find multiple arc peaks
1521 – double hit equator – 4 total peaks
2551 – no peaks
I could not find parameters to find the right peaks in these events – can see by eye where the peaks are. 
Thanks,
Meng

peaks_equ = alg_equ.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=4, r0=6, dr=2)
alg_equ.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)

        if pk.r<100     : return False
        if pk.r>454     : return False
        if pk.atot<600  : return False
        if pk.rms>25    : return False
        if pk.bkgd>15   : return False
        if pk.bkgd<-15  : return False

pk.son ?
sigma_hov ?

alg_arc.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)
peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=8, r0=10, dr=2)