| 1 | #Script to perform data reduction for MAPS |
|---|
| 2 | |
|---|
| 3 | from qtiGenie import * |
|---|
| 4 | import time |
|---|
| 5 | |
|---|
| 6 | #instrument name: |
|---|
| 7 | inst='mer' |
|---|
| 8 | iliad_setup(inst) |
|---|
| 9 | ext='.raw' |
|---|
| 10 | |
|---|
| 11 | #map file |
|---|
| 12 | mapfile='/usr/local/mprogs/Libisis/InstrumentFiles/merlin/one2one_084' #single crystal mapping file |
|---|
| 13 | mv_mapfile='/usr/local/mprogs/Libisis/InstrumentFiles/merlin/mono_van_map' |
|---|
| 14 | |
|---|
| 15 | # latest white beam vanadium file for bad detector diagnosis |
|---|
| 16 | wbvan=3763 |
|---|
| 17 | |
|---|
| 18 | |
|---|
| 19 | #Run numbers can be specified as a list: |
|---|
| 20 | #runno=[17422,17423, etc] |
|---|
| 21 | runno=[4466] |
|---|
| 22 | |
|---|
| 23 | monovan=[4581] |
|---|
| 24 | |
|---|
| 25 | sam_rmm=5.433 |
|---|
| 26 | sam_mass=349.15 |
|---|
| 27 | |
|---|
| 28 | #Incident energy list e.g. ei=[20,30,40] |
|---|
| 29 | ei=[300] |
|---|
| 30 | |
|---|
| 31 | #Lo,step,hi for energy bins of output. Needs to be an array of arrays if more than one Ei used. |
|---|
| 32 | rebin_pars=[[-30,3,279]] |
|---|
| 33 | |
|---|
| 34 | |
|---|
| 35 | for i in range(len(runno)): |
|---|
| 36 | #w1=iliad_abs(wbvan,runno[i],monovan[i],wbvan,sam_rmm,sam_mass,ei[i],str(rebin_pars[i]).strip('[]'),mapfile,mv_mapfile,bkgd_range=[14000,19000],diag_remove_zero=False,save_format='None') |
|---|
| 37 | w1=iliad(wbvan,runno[i],ei[i],str(rebin_pars[i]).strip('[]'),mapfile,bkgd_range=[14000,19000],diag_remove_zero=False,save_format='None') |
|---|
| 38 | #save_file=inst+str(runno[i])+'_ei'+str(ei[i])+'.spe' |
|---|
| 39 | #SaveSPE('w1',save_file) |
|---|
| 40 | #save_file=inst+str(runno[i])+'_ei'+str(ei[i])+'.nxspe' |
|---|
| 41 | #SaveNXSPE('w1',save_file) |
|---|
| 42 | |
|---|
| 43 | |
|---|
| 44 | |
|---|
| 45 | |
|---|
