Content

Plots

Distance between 2 peaks in Arc region

Image Removed

Image Removed

Image Removed

Fit results for for phi and beta in Equ region

Image RemovedImage Removed

Fit results for for phi and beta in Arc region

y = funcy_l1(x, phi_deg, bet_deg, DoR=392/913.27, sgnrt=-1.)

Image Removed

Image RemovedImage Removed

Image RemovedImage Removed

y = funcy_l1(x, phi_deg, bet_deg, DoR=392/913.27, sgnrt=+1.)

Phi-beta function to fit

pyimgalgos/src/FiberAngles.py

y = funcy_l1(x, phi_deg, bet_deg, DoR=392/913.27, sgnrt=-1.)

Two solutions of quadratic equation

• POS: y = -B + sqrt(B*B-C)
• NEG: y = -B - sqrt(B*B-C)

Legend:

• dot curves for negative beta
• black bold curve for beta=0
• solid curves for positive beta

Image RemovedImage Removed

For beta>27 solution does not exist?

Image processing in cxif5315

The best options for image processing in cxif5315

• Water background is evaluated WITH COMMON MODE CORRECTION that enhance intensity of central 2x1... - needs to be re-calibrated

• DO NOT use common mode

• polarization factor map should be rotated by 90 degree because of requested by Meng fancy 90 degree rotation in the geometry file
• USE status mask
• USE radial background subtraction for single angular bin

from pyimgalgos.RadialBkgd  import RadialBkgd, polarization_factor
from pyimgalgos.GlobalUtils import subtract_bkgd

        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

        # Subtract background shape averaged for pure water
        nda = subtract_bkgd(nda, sp.nda_bkgd, mask=sp.nda_smask, winds=sp.winds_bkgd, pbits=0)

        # Subtract dynamically evaluater redial background
        #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)

Fit results

ARC region, fit converged to beta ~ -17°

Image AddedImage Added

ARC region, fit converged to beta ~ -9°

Image AddedImage Added

Distance between 2 peaks in Arc region

Image Added