Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

 

Content

Table of Contents

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 phi and beta in EQU region

Image RemovedImage 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

Note

The number of events with well defined angular errors in the fit is rughly half of all events...

beta distribution is shown for 0.1<beta_error<100

 

 

 

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

Arc region

Image RemovedImage Removed

  • For beta>27 solution does not exist...

 

Equatorial region

Image RemovedImage Removed

  • For |beta|>45 solution does not exist...

Fit results for peak list from Meng

Image RemovedImage RemovedImage Removed

 

Image RemovedImage Removed

Image 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
Code Block
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