Content

Phi-beta fit function

pyimgalgos/src/FiberAngles.py

y = funcy_l1_v1(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

Arc region

wide range of x:

Image RemovedImage Removed

narrow range of x:

Image RemovedImage Removed

Equatorial region

Problem with two roots is resolved for EQU region; sign of the root is selected the same as sign of parameter B;

Image Removed

• For |beta|>~48 solution does not exist for entire range of x

Fit results

Distance between 2 peaks in Arc region

Image Removed

Image Removed

Radius of the median point of 2 peaks in Arc region

Image Removed

Fit results for for phi and beta in ARC region

Negative solution y = -B - sqrt(B*B-C)

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

Image Removed

Image RemovedImage Removed

Image RemovedImage Removed

Positive solution y = -B + sqrt(B*B-C)

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

Image Removed

Image RemovedImage Removed

Image RemovedImage Removed

Fit results for peak list from Meng

Equatorial region

Image RemovedImage RemovedImage RemovedImage Removed

Image Removed

Arc region (positive root)

Image RemovedImage RemovedImage RemovedImage RemovedImage RemovedImage RemovedImage Removed

Arc region (negative root)

Image RemovedImage RemovedImage RemovedImage RemovedImage Removed

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