brisp_pix2sqw.prox

; Example of RAW DATA DIRECTLY TRANSFORMED IN S(Q,w)
;(without any rebin or radial integration)
; 1) Get beam center and mask from transmission measurement
w1=rdrun(24420) & w1=total(w1,3)
c=brisp_center(W1,Rb_stop=8,mask_in='msk128',mask_out=W60)
; 2) Sum data and set right lambda and distance
;(lambda/Ang= PG002: 0.988, PG004: 1.977, Cu111=1.28)
w2=rdsum(24421,24433) & w2=brisp_par(w2,lam=0.988,dst=4.)
; 3) look at 2D projection
w52=total(w2,3)*(1-w60);
; 4) individuate bckg channel range
w12=brisp_rad(W2,XY0=C,INCR=6,MASK=w60)
; 5) Calculate random raw data with the expected bckg counting rate
; averaged in a pixels window (kernel)
w4=raw_bck(w2,bck_ch=[1500,1600],mask=w60,kernel=[9,6])
; 6) Calculate S(Q,w) (default: ki/kf is taken into account)
; default: dE=0.5, dQ=0.1,Qmax=2.0: W9 (optional) may contain elastic data
w22=sqw_pix(w2,XY0=c,MASK=w60,Ib=w4,w_el=w9) & w22=brisp_norm(w22)
; 7) Evaluate S(Q)
w32=total(w12,1)
; export w22 in ASCII "Q files"
out_Sqw, w22, file='namefile'