iluke_hxrcamimport

 script iluke_hxrcamimport

0001 % script iluke_hxrcamimport
0002 %
0003 hxr.cam = {};
0004 %
0005 if isfield(external.hxr,'ap') && isfield(external.hxr,'hxrparam'),% TCV data loaded by luke_get_hxr_diag_tcv
0006     %
0007 %     hxr.param = external.hxr.hxrparam;
0008     %
0009     ncams = length(external.hxr.ap.cams);
0010     %
0011     hxr.cam = cell(1,ncams);
0012     %
0013     for icam = 1:ncams
0014         %
0015         hxr.cam{icam} = tcv2r5x2(external.hxr.ap,icam);
0016         %
0017         hxr.cam{icam}.selectedchord = find(hxr.cam{icam}.chordsnum == hxr.cam{icam}.centralchord);
0018         %
0019     end
0020     %
0021     kphot_diag_S = [external.hxr.hxr.kdiag_hxr(1,1:(end-1)),external.hxr.hxr.kdiag_hxr(2,end-1)];% photon energies representing energy bin limits
0022     kphot_diag = (kphot_diag_S(1:end-1) + kphot_diag_S(2:end))/2;% photon energies representing energy bin centers
0023     %
0024     hxr.param.kphot = unique([hxr.param.kphot,kphot_diag,kphot_diag_S]);% energies for bremsstrahlung calculations
0025     %
0026 %     hxr.param.kmin = external.hxr.hxr.kdiag_hxr(1,end);
0027 %     hxr.param.kmax = external.hxr.hxr.kdiag_hxr(2,end);
0028     %
0029     clear ncams icam kphot_diag
0030     %
0031 end
0032     
0033 
0034