Code Block

def procPeakDataArc(pk) :
    """ Process peak for ARC region; accumulate peak statistics in histogram arrays.
    """
    #===================
    # discard from all histograms except its own
    sp.lst_arc_atot.append(pk.atot)
    if pk.atot<2000 : return
    #===================
    sp.lst_arc_amax.append(pk.amax)
    sp.lst_arc_npix.append(pk.npix)
    sp.lst_arc_r   .append(pk.r)
...

EQU region

Code Block

def procPeakDataEqu(pk) :
    """ Process peak for EQU region; accumulate peak data
    """
    #===================
    # discard from all histograms except its own
    sp.lst_equ_atot.append(pk.atot)
    if pk.atot<2000 : return
    sp.lst_equ_r_raw.append(pk.r)
    if pk.r<100     : return
    #===================
    sp.lst_equ_r   .append(pk.r)
    sp.lst_equ_amax.append(pk.amax)
    sp.lst_equ_npix.append(pk.npix)
  ...

Peak selection for fit

...

ARC region

Code Block

def peakIsSelectedEqupeakIsSelectedArc(pk) :
    """Apply peak selection criteria to each peak from file
    """
    if pk.rms>60son<9     : return False
    if pk.atot<1800amax<150  : return False
    if pk.npix>500 atot<2000 : return False
    if pk.r<100   npix>500  : return False
    if pk.r>454r<435     : return False
    if math.fabs(pk.bkgd)>20r>443     : return False
    return True

ARC region

Code Block
def peakIsSelectedArc(pk) : if pk.rms>80 : return False if pk.bkgd<-20 : return False if pk.bkgd>50 : return False return True

Image AddedImage AddedImage Added

To fit peaks we use funcy_l1_v0(x, phi_deg, bet_deg, DoR=433/sp.DETD, sgnrt=-1.)

Image AddedImage AddedImage Added

EQU region

Code Block

def peakIsSelectedEqu(pk) :
    """Apply peak selection criteria to each peak from file
    """ file
    """
    if pk.son<9     : return False
    if pk.amax<150  : return False
    if pk.rms>60   atot<2000 : return False
    if pk.atot<1800npix>500  : return False
    if pk.npix>500r<100     : return False
    if pk.r<435r>454     : return False
    if pk.r>443 rms>80    : return False
    if math.fabs(pk.bkgd)>20 : return False
    return True

...

To fit peaks we use funcy_l0 which aotomatically select solution depending on sign of parameter B.

Image AddedImage AddedImage Added