tcv2r5x2

 This function builds R5X2 type camera data from TCV ap structure

0001 function hxrcam = tcv2r5x2(ap,icam)
0002 %
0003 % This function builds R5X2 type camera data from TCV ap structure
0004 %
0005 hxrcam.camnum = ap.cg.cams(icam);
0006 hxrcam.id = ['TCV HXR camera # ',num2str(hxrcam.camnum)];
0007 %
0008 hxrcam.chordsnum = ap.cg.chords{hxrcam.camnum};
0009 hxrcam.chordsind = ap.cg.chords_abs{hxrcam.camnum};
0010 %
0011 hxrcam.centralchord = ap.cg.cc{hxrcam.camnum};
0012 %
0013 nc = length(hxrcam.chordsind);
0014 %
0015 % focal point position
0016 hxrcam.R_hxr=ap.cg.rd(hxrcam.chordsind);
0017 hxrcam.Z_hxr=ap.cg.zd(hxrcam.chordsind);
0018 % chord angles
0019 hxrcam.alpha_hxr=ap.cg.alpha(hxrcam.chordsind);
0020 hxrcam.beta_hxr=ap.cg.beta(hxrcam.chordsind);
0021 % etendue data
0022 hxrcam.EG_hxr = ap.cg.etendue(hxrcam.chordsind);
0023 
0024 % keep data for HXRC detector for the moment
0025 hxrcam.kd_hxr = ones(nc,1)*[10,30,50,70,90,110,130,150,170,190,210,250];%Parameters approximately corresponding to  photon energy spectra of gamma sources 57Co and 241Am. Reference energies for the photofraction [m,p]
0026 hxrcam.phf_hxr = ones(nc,1)*[0.75,0.62,0.51,0.41,0.32,0.2,0.2,0.2,0.2,0.2,0.2,0.2];%Parameters approximately corresponding to  photon energy spectra of gamma sources 57Co and 241Am. Photofraction determined using the "MCDET.f" Monte-carlo hard x-ray absorption code at energies kd [m,p]
0027 hxrcam.res_hxr = ones(nc,1)*[6,-0.9]; % fit parameters of the energy resolution of the detectors [m,2]
0028 
0029 % filter, window and absorber
0030 hxrcam.ec_hxr = NaN(1,nc);          % thickness of the detector shield (mm) [1,m]
0031 hxrcam.tc_hxr = repmat('fm',nc,1);  % type of material for the detector shield (Al,C ou Be) [m,2]
0032 hxrcam.ef_hxr = NaN(1,nc);          % thickness of the vacuum window (mm) [1,m]
0033 hxrcam.tf_hxr = repmat('wm',nc,1);  % type of material for the vacuum window (Al ou Be) [m,2]
0034 hxrcam.ea_hxr = NaN(1,nc);          % thickness of the absorber (mm) [1,m] (WARNING: for 2 mm Cu: 2.83= 2*(29/26)^2*(dens_Cu/dens_Fe)  dens_Cu = 63.546, dens_Fe = 55.845)
0035 hxrcam.ta_hxr = repmat('am',nc,1);  % type of material for the absorber (Al,Fe,Ge,Be ou Pb) [m,2]
0036 hxrcam.ed_hxr = NaN(1,nc);          % thickness of the detector (mm) [1,m]
0037 hxrcam.td_hxr = repmat('detm',nc,1);% type of material for the detector (BGO , CsI , NaI , Ge  , CdTe) [m,4]
0038 for ic = 1:nc
0039     hxrcam.ec_hxr(ic) = ap.fch{hxrcam.chordsind(ic)}.filt{1}.d;
0040     hxrcam.tc_hxr(ic,:) = ap.fch{hxrcam.chordsind(ic)}.filt{1}.mt;
0041     hxrcam.ef_hxr(ic) = ap.fch{hxrcam.chordsind(ic)}.win{1}.d;
0042     hxrcam.tf_hxr(ic,:) = ap.fch{hxrcam.chordsind(ic)}.win{1}.mt;
0043     hxrcam.ea_hxr(ic) = ap.fch{hxrcam.chordsind(ic)}.abs{1}.d;
0044     hxrcam.ta_hxr(ic,:) = ap.fch{hxrcam.chordsind(ic)}.abs{1}.mt;
0045     hxrcam.ed_hxr(ic) = ap.det{hxrcam.chordsind(ic)}.d;
0046     hxrcam.td_hxr(ic,:) = ap.det{hxrcam.chordsind(ic)}.mt;
0047 end