Peak/event selection parameters

Peak/event selection parameters 
def peak_is_selected_arc(pk) :
    if pk.reg!='ARC': return False 
    if pk.r<433     : return False
    if pk.r>445     : return False
    if pk.rms>30    : return False
    if pk.bkgd>15   : return False
    if pk.bkgd<-15  : return False
    if pk.son<12    : return False
    if pk.npix<10   : return False
    return True
#------------------------------
def peak_is_selected_equ(pk) :
    if pk.reg!='EQU': return False
    if pk.r<105     : return False
    if pk.r>450     : return False
    if pk.rms>25    : return False
    if pk.bkgd>10   : return False
    if pk.bkgd<-10  : return False
    return True
#------------------------------
def select_equ() :
    sel_arc = len(sp.peaks_arc_sel)
    sel_equ = len(sp.peaks_equ_sel)
    # special case of 2 peaks in arc
    if sel_arc==2 and sel_equ<5 : return True
    if sel_equ>4 : return False
    if sel_equ<2 : return False
    return True
#------------------------------
def select_arc() :
    sel_arc = len(sp.peaks_arc_sel)
    if sel_arc>2 : return False
    return True

Processing of RAW peaks after pfv4r2

cxif5315/filter-peaks-2017-03-07.py -i work/xpfv4r2-cxif5315-r0169-2017-03-03-e085409.txt -p
----
Summary
Input file       : work/xpfv4r2-cxif5315-r0169-2017-03-03-e085409.txt
Events processed : 83727
Events selected  : 9735
Consumed time    : 152.409 sec
----
cp results/peaks-filtred.txt work/xpfv4r2-cxif5315-r0169-2017-03-03-e085409-filtered.txt
mv results 2017-03-10-figs-raw

Processing of FILTERED peaks for cross-check

cxif5315/filter-peaks-2017-03-07.py -i work/xpfv4r2-cxif5315-r0169-2017-03-03-e085409-filtered.txt  -o results/peaks-filtred-2 -p
----
Summary
Input file       : work/xpfv4r2-cxif5315-r0169-2017-03-03-e085409-filtered.txt
Events processed : 9735
Events selected  : 9735
----
mv results 2017-03-10-figs-filtered

Histograms for RAW and FILTERED peaks processing

ARC

 

EQU

 

 

 

 

 

  • No labels