Changeset 430 in ETALON for SPESO


Ignore:
Timestamp:
Mar 30, 2016, 8:54:21 PM (8 years ago)
Author:
malovyts
Message:

Added comparison plot for different TX

File:
1 edited

Legend:

Unmodified
Added
Removed
  • SPESO/ana2015/offline_analysis/offline_analyse_colormaps_2016JanFeb.m

    r428 r430  
    1414
    1515analyse_date=datestr(now, 'yyyymmdd');
    16 figures_path=[plots_basedir 'offline_analysis' '/' '2016JanFeb' '/' 'analysis_for' analyse_date '/'];
     16figures_path=[plots_basedir 'offline_analysis' '/' '2016JanFeb' '/' 'analysis_for' analyse_date '_02/'];
    1717
    1818mkdir(figures_path);
     
    109109fprintf(fid, '\t%s\n', '<head></head>');
    110110fprintf(fid, '\t%s\n', '<body>');
     111fprintf(fid, '\t\t%s\n', ['<H1>Simulation time:' datestr(now) '</H1>']);
    111112
    112113% all selection
     
    207208
    208209                        % for different XPOS signal/bg
    209                         for XPOSloop=1:length(XPOS_selection)
     210                        for RESEAUloop=1:length(XPOS_selection)
    210211                                fprintf(fid, '\t\t\t%s\n', '<tr>');
    211212                                fprintf(fid, '\t\t\t\t%s\n','<td>');
     
    217218                                for QICT1loop=1:length(QICT1_selection)
    218219                                        img_name=[data_names{dataloop} '_monthly_'...
    219                                                                          sprintf('%07d',(0+QICT1loop*10+LPMloop*100+POSloop*1000+XPOSloop*10000))...
     220                                                                         sprintf('%07d',(0+QICT1loop*10+LPMloop*100+POSloop*1000+RESEAUloop*10000))...
    220221                                                                         '.png'];
    221222                                        fprintf(fid, '\t\t\t\t%s\n', '<td>');
     
    237238                                for QICT1loop=1:length(QICT1_selection)
    238239                                        img_name=[data_names{dataloop} 'histc_monthly_'...
    239                                                                          sprintf('%07d',(1+QICT1loop*10+LPMloop*100+POSloop*1000+XPOSloop*10000))...
     240                                                                         sprintf('%07d',(1+QICT1loop*10+LPMloop*100+POSloop*1000+RESEAUloop*10000))...
    240241                                                                         '.png'];
    241242                                        fprintf(fid, '\t\t\t\t%s\n', '<td>');
     
    249250                                % each 2nd datatype - comparison of two detectors
    250251                                if mod(dataloop,2)==0
     252                                for XPOSloop=1:length(XPOS_selection)
    251253                                        fprintf(fid, '\t\t\t%s\n', '<tr>');
    252254                                        fprintf(fid, '\t\t\t\t%s\n','<td>');
     
    258260                                        for QICT1loop=1:length(QICT1_selection)
    259261                                                img_name=['histc_monthly_'...
    260                                                                                  sprintf('%07d',(3+QICT1loop*10+LPMloop*100+POSloop*1000+XPOSloop*10000))...
     262                                                                                 sprintf('%07d',(3+QICT1loop*10+LPMloop*100+POSloop*1000+RESEAUloop*10000+100000*XPOSloop))...
    261263                                                                                 '.png'];
    262264                                                fprintf(fid, '\t\t\t\t%s\n', '<td>');
     
    267269                                        end; % QICTloop
    268270                                        fprintf(fid, '\t\t\t%s\n', '</tr>');
     271                                end; %XPOSloop
    269272                                end; % mod(dataloop,2)==0
    270                         end;%XPOSloop
     273                        end;%RESEAUloop
    271274
    272275                end; % LPMloop
     
    305308                        % base selection
    306309                        base_selected=allData.selected;
    307                         base_selected_text=allData.selected;
     310                        base_selected_text=allData.selection;
    308311
    309312        if dataloop==1
     313                for XPOSloop=1:length(XPOS_selection)
    310314                %============ comparison of two detectors LPM, SPM ==============
     315                data_select(XPOS_selection{XPOSloop}{1}{:});   
     316                data_select(XPOS_selection{XPOSloop}{2}{:});   
    311317                fih = figure(1);
    312318                        bins=data_histc_multiDet(POS_var{POSloop}, data_types((RESEAUloop-1)*2+1), data_types((RESEAUloop-1)*2+2),...
     
    318324                        set(gcf,'PaperUnits','inches','PaperPosition',[0 0 12 9]);
    319325                        print(fih, '-dpng', '-r300',[ figures_path 'histc_monthly_'...
    320                                                               sprintf('%07d',(3+QICT1loop*10+LPMloop*100+POSloop*1000+10000*RESEAUloop))...
     326                                                              sprintf('%07d',(3+QICT1loop*10+LPMloop*100+POSloop*1000+10000*RESEAUloop+100000*XPOSloop))...
    321327                                                              '.png' ]);
    322328                close(fih);
    323329                % --------------------------------------------
     330                allData.selected=base_selected;
     331                allData.selection=base_selected_text;
     332                end; % XPOSloop
    324333        end; % if dataloop
    325 
    326334
    327335
     
    356364                %---------------------------------------------------------
    357365                allData.selected=base_selected;
    358                 allData.selected=base_selected_text;
     366                allData.selection=base_selected_text;
    359367
    360368                data_select(XPOS_selection{RESEAUloop}{1}{:}); 
Note: See TracChangeset for help on using the changeset viewer.