Skip to content
Snippets Groups Projects

Add tcv get ids soft x rays to tcv imas functions

Merged Add tcv get ids soft x rays to tcv imas functions
Add-tcv_get_ids_soft_x_rays-to-TCV_IMAS-functions into master
Merged Luke Simons requested to merge Add-tcv_get_ids_soft_x_rays-to-TCV_IMAS-functions into master
Compare
and
1 file
+ 139
− 0
Compare changes
  • Side-by-side
  • Inline
matlab/TCV_IMAS/tcv_get_ids_soft_x_rays.m 0 → 100644
+ 139
− 0
function [ids_soft_x_rays,ids_soft_x_rays_description] = ...
tcv_get_ids_soft_x_rays(shot, ids_soft_x_rays_empty, gdat_params,varargin)
%
% [ids_soft_x_rays,ids_soft_x_rays_description] = ...
% tcv_get_ids_soft_x_rays(shot, ids_soft_x_rays_empty,varargin);
%
% Get the RADCAM Soft X-Ray (SXR) diagnostics data
%
% ids_soft_x_rays_empty should at least be the empty Soft X-Ray ids
% structure in input
%
% gdat_params: gdat_data.gdat_params to get all params passed from original
% call, in particular error_bar options
%
% Date: 20/09/24
% Author: L. Simons
% Initialise structure with default ids_properties
% shot=70777; gdat_params=[]; varargin={};
% aa=gdat([],'ids','source','soft_x_rays');
% ids_soft_x_rays_empty=aa.soft_x_rays;
if exist('gdat_params')
[ids_soft_x_rays, ~] = tcv_ids_headpart(shot, ids_soft_x_rays_empty, ...
+1
'soft_x_rays','homogeneous_time',1,'gdat_params',gdat_params, ...
'comment','RADCAM Soft X-Ray 47um Be filered optodiode AXUV20ELGDS diodes ', ...
varargin{:});
else
[ids_soft_x_rays, ~] = tcv_ids_headpart(shot, ids_soft_x_rays_empty, ...
'soft_x_rays','homogeneous_time',0, ...
'comment','RADCAM Soft X-Ray 47um Be filered optodiode AXUV20ELGDS diodes', ...
varargin{:});
end
ids_soft_x_rays_description = struct();
% Load the soft_x_rays geometry
sxr_geom=sxr_load_geometry();
sxr_data=rc_load_diodes(shot,'diag_name','sxr');
status=numel(sxr_data.data)>0;
%get gdat version, FIXME: Need to write function
%radcam_repo_version=imas_get_code_version();
if status
nchannel = sum(rc_n_chn(shot,'sxr'));
ids_soft_x_rays.channel(1:nchannel) = ids_soft_x_rays.channel(1);
for ii = 1:nchannel
+6
ids_soft_x_rays.channel{ii}.name = sxr_geom.channel{ii}.name;
ids_soft_x_rays_description.channel{ii}.name = ...
['RADCAM Soft X-Ray (SXR) detector name, from camera labels and ' ...
'number VT: Vertical Top, UL: Upper Lateral, ML: Middle Lateral, ' ...
'DIV: Divertor, VB: Vertical bottom. Matches cabling.'];
ids_soft_x_rays.channel{ii}.identifier = ...
sxr_geom.channel{ii}.identifier;
ids_soft_x_rays_description.channel{ii}.identifier = ...
['RADCAM Soft X-Ray (SXR) system, global index from clock-wise ' ...
'ordering from inboard side'];
%ids_soft_x_rays_description.channel{ii}.description = ...
% [sprintf('SXR Diode detector %d, %s Line of Sight %d',ii,,)];
ids_soft_x_rays.channel{ii}.detector.geometry_type= sxr_geom.channel{ii}.detector.geometry_type;
ids_soft_x_rays.channel{ii}.detector.centre.phi= sxr_geom.channel{ii}.detector.centre.phi;
ids_soft_x_rays.channel{ii}.detector.centre.r= sxr_geom.channel{ii}.detector.centre.r;
ids_soft_x_rays.channel{ii}.detector.centre.z= sxr_geom.channel{ii}.detector.centre.z;
ids_soft_x_rays.channel{ii}.detector.x1_unit_vector.x= sxr_geom.channel{ii}.detector.x1_unit_vector.x;
ids_soft_x_rays.channel{ii}.detector.x1_unit_vector.y= sxr_geom.channel{ii}.detector.x1_unit_vector.y;
ids_soft_x_rays.channel{ii}.detector.x1_unit_vector.z= sxr_geom.channel{ii}.detector.x1_unit_vector.z;
ids_soft_x_rays.channel{ii}.detector.x2_unit_vector.x= sxr_geom.channel{ii}.detector.x2_unit_vector.x;
ids_soft_x_rays.channel{ii}.detector.x2_unit_vector.y= sxr_geom.channel{ii}.detector.x2_unit_vector.y;
ids_soft_x_rays.channel{ii}.detector.x2_unit_vector.z= sxr_geom.channel{ii}.detector.x2_unit_vector.z;
ids_soft_x_rays.channel{ii}.detector.x3_unit_vector.x= sxr_geom.channel{ii}.detector.x3_unit_vector.x;
ids_soft_x_rays.channel{ii}.detector.x3_unit_vector.y= sxr_geom.channel{ii}.detector.x3_unit_vector.y;
ids_soft_x_rays.channel{ii}.detector.x3_unit_vector.z= sxr_geom.channel{ii}.detector.x3_unit_vector.z;
ids_soft_x_rays.channel{ii}.detector.x1_width= sxr_geom.channel{ii}.detector.x1_width;
ids_soft_x_rays.channel{ii}.detector.x2_width= sxr_geom.channel{ii}.detector.x2_width;
ids_soft_x_rays.channel{ii}.detector.surface= sxr_geom.channel{ii}.detector.surface;
ids_soft_x_rays.channel{ii}.aperture{1}.geometry_type= sxr_geom.channel{ii}.aperture{1}.geometry_type;
ids_soft_x_rays.channel{ii}.aperture{1}.centre.r= sxr_geom.channel{ii}.aperture{1}.centre.r;
ids_soft_x_rays.channel{ii}.aperture{1}.centre.z= sxr_geom.channel{ii}.aperture{1}.centre.z;
ids_soft_x_rays.channel{ii}.aperture{1}.x1_unit_vector.x= sxr_geom.channel{ii}.aperture{1}.x1_unit_vector.x;
ids_soft_x_rays.channel{ii}.aperture{1}.x1_unit_vector.y= sxr_geom.channel{ii}.aperture{1}.x1_unit_vector.y;
ids_soft_x_rays.channel{ii}.aperture{1}.x1_unit_vector.z= sxr_geom.channel{ii}.aperture{1}.x1_unit_vector.z;
ids_soft_x_rays.channel{ii}.aperture{1}.x2_unit_vector.x= sxr_geom.channel{ii}.aperture{1}.x2_unit_vector.x;
ids_soft_x_rays.channel{ii}.aperture{1}.x2_unit_vector.y= sxr_geom.channel{ii}.aperture{1}.x2_unit_vector.y;
ids_soft_x_rays.channel{ii}.aperture{1}.x2_unit_vector.z= sxr_geom.channel{ii}.aperture{1}.x2_unit_vector.z;
ids_soft_x_rays.channel{ii}.aperture{1}.x3_unit_vector.x= sxr_geom.channel{ii}.aperture{1}.x3_unit_vector.x;
ids_soft_x_rays.channel{ii}.aperture{1}.x3_unit_vector.y= sxr_geom.channel{ii}.aperture{1}.x3_unit_vector.y;
ids_soft_x_rays.channel{ii}.aperture{1}.x3_unit_vector.z= sxr_geom.channel{ii}.aperture{1}.x3_unit_vector.z;
ids_soft_x_rays.channel{ii}.aperture{1}.x1_width= sxr_geom.channel{ii}.aperture{1}.x1_width;
ids_soft_x_rays.channel{ii}.aperture{1}.x2_width= sxr_geom.channel{ii}.aperture{1}.x2_width;
ids_soft_x_rays.channel{ii}.aperture{1}.surface= sxr_geom.channel{ii}.aperture{1}.surface;
ids_soft_x_rays.channel{ii}.etendue = sxr_geom.channel{ii}.etendue;
ids_soft_x_rays.channel{ii}.etendue_method.name = sxr_geom.channel{ii}.etendue_method.name;
ids_soft_x_rays.channel{ii}.etendue_method.description = sxr_geom.channel{ii}.etendue_method.description;
ids_soft_x_rays.channel{ii}.filter_window{1}.identifier = 1;
ids_soft_x_rays.channel{ii}.filter_window{1}.material.name = 'Be';
ids_soft_x_rays.channel{ii}.filter_window{1}.material.index = 10;
ids_soft_x_rays.channel{ii}.filter_window{1}.material.description = '';
ids_soft_x_rays.channel{ii}.filter_window{1}.thickness = '47e-6';
ids_soft_x_rays.channel{ii}.line_of_sight.first_point.phi = sxr_geom.channel{ii}.line_of_sight.first_point.phi;
ids_soft_x_rays.channel{ii}.line_of_sight.first_point.r = sxr_geom.channel{ii}.line_of_sight.first_point.r;
ids_soft_x_rays.channel{ii}.line_of_sight.first_point.z = sxr_geom.channel{ii}.line_of_sight.first_point.z;
ids_soft_x_rays.channel{ii}.line_of_sight.second_point.phi = sxr_geom.channel{ii}.line_of_sight.second_point.phi;
ids_soft_x_rays.channel{ii}.line_of_sight.second_point.r = sxr_geom.channel{ii}.line_of_sight.second_point.r;
ids_soft_x_rays.channel{ii}.line_of_sight.second_point.z = sxr_geom.channel{ii}.line_of_sight.second_point.z;
ids_soft_x_rays.channel{ii}.brightness.data = sxr_data.data(ii,:);
ids_soft_x_rays.channel{ii}.brightness.time = sxr_data.time;
ids_soft_x_rays_description.channel{ii}.brightness.data = ...
['From ATLAS::SYSTEM:SXR:CHANNEL_ data, Radiance measured by ' ...
'Each SXR diode in arbitrary units of [counts]'];
ids_soft_x_rays.channel{ii}.validity = 0; % FIXME: Always valid data
ids_soft_x_rays.channel{ii}.validity = 0*ones(1,nchannel);
ids_soft_x_rays.code.name = 'rc_load_diodes';
ids_soft_x_rays.code.description = ['rc_load_diodes, RADCAM gitlab repo, loads' ...
'the signal from the SXR diodes in arbitrary units.'];
ids_soft_x_rays.code.commit='';
ids_soft_x_rays.code.version='';
ids_soft_x_rays.code.repository= ...
'https://gitlab.epfl.ch/spc/tcv/diag/radcam/';
ids_soft_x_rays.code.parameters='''diag_name'',''sxr''';
end
+2
ids_soft_x_rays.time = sxr_data.time;
else
warning('Failed to load data for shot %d',shot);
end
end