Page History
Content
Table of Contents |
---|
Fit results
Distance between 2 peaks in Arc region
Radius of the median point of 2 peaks in Arc region
Fit results for phi and beta in EQU region
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.)
Positive solution y = -B + sqrt(B*B-C)
y = funcy_l1(x, phi_deg, bet_deg, DoR=392/913.27, sgnrt=+1.)
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
- For beta>27 solution does not exist...
Equatorial region
- For |beta|>45 solution does not exist..
Fit results for peak list from Meng
Equatorial region (negative root)
Arc region (positive root)
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°
ARC region, fit converged to beta ~ -9°
Distance between 2 peaks in Arc region
Overview
Content Tools