Code Block

def procPeakDataArc(pk) :
    """ Process peak for ARC region; accumulate peak statistics in histogram arrays.
    """
    #===================
    sp.lst_arc_atot.append(pk.atot)
    if pk.atot<1500 : return
    #===================
    sp.lst_arc_amax.append(pk.amax)
    sp.lst_arc_npix.append(pk.npix)
    sp.lst_arc_r   .append(pk.r)
    sp.lst_arc_phi .append(pk.phi)
...

plots for pfv2/3/4r1

...

EQU

pre-selection for histograms

...

plots for pfv2/3/4r1

Fit results

ARC

pre-fit selection

Code Block

def peakIsSelectedArc(pk) :
    """Apply peak selection criteria to each peak from file
    """
    if sp.USE_PEAK_SELECTION :
        if pk.atot<1500 : return False
        if pk.npix>200  : return False
        if pk.r<435     : return False
        if pk.r>443     : return False
        if pk.rms<10    : return False
        if pk.rms>80    : return False
        if math.fabs(pk.bkgd)>25 : return False
    return True

plots for pfv2/3/4r1

EQU

pre-fit selection

Code Block

def peakIsSelectedEqu(pk) :
    """Apply peak selection criteria to each peak from file
    """
    if sp.USE_PEAK_SELECTION :
        if pk.atot<1500 : return False
        if pk.npix>200  : return False
        if pk.rms>80    : return False
        if math.fabs(pk.bkgd)>20 : return False
    if pk.r<100     : return False
    if pk.r>454     : return False
    return True