lamp_mac/BRISP/
sqw_pix.pro
top sqw_pix
result = sqw_pix(W_in, XY0=XY0, MASK=MASK, file_el=file_el, dE=dE, dQ=dQ, Qmax=Qmax, Ib=Ib, kikf=kikf, w_el=w_el)
******* ******* ** User callable ** Modifications (April 2011, FF): ** 1) time channels are not anymore discrete (1,2,...) but randomised ** If Ib used, ambient background is subtracted. ** 2) possibility of not making the ki/kf correction (default: ki/kf factor is included) ** 3) Elastic peak position is determined through an analisys of the channel intensity pixel by pixel, ** and not anymore after a radial integratrion procedure. Input data can be the same W_in, ** if statistic is sufficient, or another worksheet (eg vanadium or sample) ** 4) Data are better centred in energy ** ** OUTPUT: Sqw from raw data, converting I(theta, channels) of each pixels, without rebinning. ** If Ib used, ambient background is subtracted. With r without ki/kf factor. ** ** ***** ** INPUT: W_in = 3D raw data matrix ** ***** xy0 = [x0,y0] center determined by "brisp_center.pro" ** MASK = mask determined manually and by "brisp_center.pro" ** contains 1's for masked values otherwise 0's(default) ** ** **** OPTIONAL INPUTS ** Ib = 3D worksheet containing background data randomnly generated by raw_bck.pro ** w_el = 3D worksheet containing "good" statistic data. Needed to evaluate elastic channel position as a function of theta ** dE = energies step (default 0.5 meV) ** kikf = 'no': without ki/kf. Default is with ki/kf ** dQ = Q step (default 0.1 Ang^-1) ** Qmax = maximum Q at which Sqw is determined (default 2.0 Ang^-1) ** **** DIMENSIONS: Wout(dE,dQ) ** ** Call: w20=sqw_pix(w2,XY0=c,MASK=w40,kikf='no',w_el=w2,Ib=w4,dE=0.5,dQ=0.1,Qmax=2.0) ** ** (FF, April 2010) ** ** NUOVO RANDOM SUI CHANNELS **
Parameters
- W_in
Keywords
- XY0
- MASK
- file_el
- dE
- dQ
- Qmax
- Ib
- kikf
- w_el
File attributes
Modification date: | Tue Oct 1 14:45:32 2013 |
Lines: | 204 |