Xarr =  sp.det.coords_x(runnum)
        Yarr =  sp.det.coords_y(runnum)

        sp.rb = RadialBkgd(Xarr, Yarr, mask=sp.nda_smask, radedges=(5200, 80000), nradbins=200, nphibins=1)
        sp.pf = polarization_factor(sp.rb.pixel_rad(), sp.rb.pixel_phi()+90, sp.DETD_um)

    nda_raw = sp.det.raw(evt)

    if nda_raw is not None :

        nda =  np.array(nda_raw, dtype=np.float32, copy=True)
        nda -= sp.nda_peds
        #nda =  subtract_bkgd(nda, sp.nda_bkgd, mask=sp.nda_smask, winds=sp.winds_bkgd, pbits=0)

        nda = sp.rb.subtract_bkgd(nda.flatten() * sp.pf)
        nda.shape = sp.nda_peds.shape

        nda *= sp.nda_smask
        #sp.det.common_mode_apply(evt, nda)
        
        img = getImage(evt, nda)