Changeset 512 in ETALON for SPESO


Ignore:
Timestamp:
Apr 21, 2016, 5:56:31 PM (8 years ago)
Author:
malovyts
Message:

Fixed polar plots in the offline analysis

Location:
SPESO/ana2015
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • SPESO/ana2015/data_get_vals.m

    r504 r512  
    107107                vals(iloop,:)=allData.beamInfo(idx,5)'-RS_params.RS1_offset;
    108108            case dvar.DATA_RS2_OFFSET
    109                 RS_params=get_RS_params(allData.timestamp);
     109                RS_params=get_RS_params(allData.timestamp(1,idx));
    110110                vals(iloop,:)=allData.beamInfo(idx,7)'-RS_params.RS2_offset;
    111111            case dvar.DATA_RS3_OFFSET
    112                 RS_params=get_RS_params(allData.timestamp);
     112                RS_params=get_RS_params(allData.timestamp(1,idx));
    113113                vals(iloop,:)=allData.beamInfo(idx,9)'-RS_params.RS3_offset;
    114114            case dvar.DATA_RS4_OFFSET
    115                 RS_params=get_RS_params(allData.timestamp);
     115                RS_params=get_RS_params(allData.timestamp(1,idx));
    116116                vals(iloop,:)=allData.beamInfo(idx,11)'-RS_params.RS4_offset;
    117117            otherwise

  • SPESO/ana2015/data_rose.m

    r507 r512  
    1616                {bins2(2,:)*pi/180, bins2(1,:)}},...
    1717                {{'LineWidth', 2, 'Color', 'r'},{'LineWidth', 2, 'Color', 'b'}},...
    18                 [40 140]*pi/180) % polar plot angle edeges
     18                [min(edges) max(edges)]*pi/180) % polar plot angle edeges
    1919
    2020legend(selection_text1, selection_text2, 'Location','southoutside');

  • SPESO/ana2015/offline_analysis/data_rose.m

    r498 r512  
    1313plot_title=plot_bins.title;
    1414
    15    custom_rose({{bins1(2,:)*pi/180, bins1(1,:)},...
    16                 {bins2(2,:)*pi/180, bins2(1,:)}},...
     15   custom_rose({{bins1(2,:)*pi/180, abs(bins1(1,:))},...
     16                {bins2(2,:)*pi/180, abs(bins2(1,:))}},...
    1717                {{'LineWidth', 2, 'Color', 'r'},{'LineWidth', 2, 'Color', 'b'}},...
    18                 [40 140]*pi/180) % polar plot angle edeges
     18                [min(edges) max(edges)]*pi/180) % polar plot angle edeges
    1919
    2020legend(selection_text1, selection_text2, 'Location','southoutside');

  • SPESO/ana2015/offline_analysis/offline_analysis.m

    r503 r512  
    3939%date_mask='20160{3[^1]?,31[^3],2??,1??}.dayFullInfo'; % months 01, 02, 03 except day 13th
    4040
    41 load_analysis_data(date_mask);
     41%load_analysis_data(date_mask);
    4242
    4343% statistics plots

  • SPESO/ana2015/offline_analysis/offline_histograms.m

    r503 r512  
    4141
    4242offline_histograms_html('main', selection_id); % create HTML page
     43
     44array_length=length(LPM_selection)*length(XPOS_selection)*length(QICT1_selection)*2;
     45
     46Q=cell(array_length,1);
     47Signal_PM=cell(array_length,1);
     48x_pos=cell(array_length,1);
     49theta=cell(array_length,1);
     50values=cell(array_length,1);
     51value_error=cell(array_length,1);
     52bin_elem_num=cell(array_length,1);
     53Detector=cell(array_length,1);
     54
     55idx=0;
    4356
    4457% make plots
    … …  
    8194                %============ comparison of two detectors LPM, SPM ==============
    8295                fih = figure(1);
    83                         bins=data_histc_multiDet(POS_var{POSloop},...
     96                        bins=data_plot2d_multiDet(POS_var{POSloop},...
    8497                                                 data_types((LPMloop-1)*2+1),...
    8598                                                 data_types((LPMloop-1)*2+2),...
    … …  
    135148                close(fih);
    136149                % -------------------------------------------------
     150       
     151                switch(POSloop+(DetLoop-1)*2)
     152                                case 1
     153                                        DATA_RS=dvar.DATA_RS3_OFFSET;                           
     154                                        RS_edges=[-50:1:50];
     155                                case 2
     156                                        DATA_RS=dvar.DATA_RS1_OFFSET;                           
     157                                        RS_edges=[40:2:140];
     158                                case 3
     159                                        DATA_RS=dvar.DATA_RS4_OFFSET;                           
     160                                        RS_edges=[-50:1:50];
     161                                case 4
     162                                        DATA_RS=dvar.DATA_RS2_OFFSET;                           
     163                                        RS_edges=[40:2:140];
     164                                otherwise
     165                                        display('ERROR!!!');
     166                                        return;
     167                                        DATA_RS=dvar.DATA_RS1_OFFSET;
     168               RS_edges=[0:180];
     169                end;
    137170
    138171                % ================ polar rose signal/bg ===============
    139172                fih = figure(1);
    140                         signal_bins=data_rose(dvar.DATA_RS1_OFFSET, data_types(1),...
     173                        signal_bins=data_rose(DATA_RS, data_types((LPMloop-1)*2+DetLoop),...
    141174                                   {...
    142175                                    {RESEAU_POS_selection{1}{1}},...
    143176                                    {RESEAU_POS_selection{2}{1}}...
    144177                                   },...
    145                                 [95:2:120]);
     178                                RS_edges);
    146179                        print_fig(fih, [ figures_path 'polar_Det' num2str(DetLoop) '_' pars_txt '.png' ]);
    147180                close(fih);
    148181                % -------------------------------------------------
     182                if POSloop==2
     183                        idx=idx+1;
     184                        Q{idx}=QICT1loop;
     185                        Signal_PM{idx}=LPMloop;
     186                        x_pos{idx}=XPOSloop;
     187                        theta{idx}=signal_bins.bins1(2,:);
     188                        values{idx}=signal_bins.bins1(1,:)-signal_bins.bins2(1,:);
     189                        bins_elem_num{idx}=signal_bins.bins1(4,:);
     190                        Detector{idx}=DetLoop;
     191                        value_error{idx}=signal_bins.bins1(3,:);
     192                        %signal_values=signal_bins.bins1(1,:);
     193                        %bg_values=signal_bins.bins2(1,:);
     194                        %positions=signal_bins.bins1(2,:);
     195                end;
     196
    149197          end; % DetLoop
    150198         end; % XPOSloop
    … …  
    153201 end; % POSloop
    154202
     203save('ExperimentSpectrum_20160421.mat',...
     204                'Q',...
     205                'Signal_PM',...
     206                'x_pos',...
     207                'theta',...
     208                'values',...
     209                'value_error',...
     210                'bins_elem_num',...
     211                'Detector',...
     212                '-mat');
     213
    155214data_select_reset();
    156215end % function
Note: See TracChangeset for help on using the changeset viewer.