Changeset 634 in ETALON for SPESO


Ignore:
Timestamp:
Jan 25, 2017, 10:08:40 PM (7 years ago)
Author:
delerue
Message:

Updated SPESO analysis code

Location:
SPESO/ana2015
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • SPESO/ana2015/all_data_summary.m

    r621 r634  
    1010global dayData;
    1111
    12 mask='2016????_scope' ;
     12mask='201[6,7]????_scope' ;
     13%mask='201[6,7]10??_scope' ;
    1314%mask='201511??_scope';
    1415%recentDateNum=736248;
    … …  
    108109
    109110           daysHtmlTableNS= [ ' <TD> ' num2str(dayDataInfo.nShots(iday)) '</TD>'  daysHtmlTableNS ];
     111disp(['iday ' num2str(iday) ' -- ' thisDayStr ])
     112
    110113idayFull=find(dayData.dateNum==dayDataInfo.dateNum(iday));
     114
     115length(idayFull)
     116if (length(idayFull)>0)
    111117dayDataInfo.nFullShots(iday)=dayData.nFullShots(idayFull);
     118else
     119  dayDataInfo.nFullShots(iday)=0
     120end
    112121daysHtmlTableNFS= [  ' <TD> ' num2str(dayData.nFullShots(idayFull)) '</TD>' daysHtmlTableNFS ];
    113122end %for
    … …  
    232241fclose(fwarn);
    233242
    234 
    235243fprintf(fid,daysHtmlTable);
     244fprintf(fid,'<HR>');
     245fprintf(fid,'<H2>Computers status</H2><BR/>');
     246[ret,txt]=system('uname -n');
     247fprintf(fid,[ 'Machine: ' txt '<BR/>']);
     248[ret,txt]=system('uptime');
     249fprintf(fid,[ 'Uptime: ' txt  ' ---  '] );
     250[ret,txt]=system('date');
     251fprintf(fid,[ 'Date: ' txt  '<BR/>'] );
     252%[ret,txt]=system('top -b -n1 | head -15')
     253%  txt=strrep(txt',10, '<BR/>'); %'
     254%fprintf(fid,[ 'Uptime: ' txt  '<BR/><BR/>'] );
     255
     256fstatusname=[ basedir 'status.info' ];
     257fstatus=fopen(fstatusname,'r');
     258txt=fread(fstatus);
     259fprintf(fid, '<BR/><BR/>Remote status: <BR/>' );
     260poseol=strfind(txt',10); %'
     261            txt=char(txt(1:poseol(min(12,length(poseol))))'); %'
     262txt=strrep(txt,10, '<BR/>');
     263fwrite(fid,txt);
     264fprintf(fid, '<BR/>' );
     265fclose(fstatus);
     266
     267fstatusname=[ basedir 'warning_sp.txt' ];
     268fstatus=fopen(fstatusname,'r');
     269txt=fread(fstatus);
     270fprintf(fid, 'Latest warnings: <BR/>' );
     271poseol=strfind(txt',10); %'
     272               txt=char(txt(poseol(length(poseol)-20):end)'); %'
     273txt=strrep(txt,10, '<BR/>');
     274fwrite(fid,txt);
     275fprintf(fid, '<BR/>' );
     276fclose(fstatus);
     277
    236278
    237279fprintf(fid,'<HR>');

  • SPESO/ana2015/analyse_channel.m

    r312 r634  
    166166        subplot(3,1,1)
    167167        hold on
    168         plot(data_num(:,1), data_num(:,2), '-');
     168        plot(data_num(:,1), data_num(:,2), '-');     
    169169        xlabel('time[seconds]');
    170170        ylabel('voltage[V]');
    … …  
    209209        xlabel('time[seconds]');
    210210        ylabel('voltage[V]');
     211if (iPulseMode==1)
     212        xlabel('LPM (sgnl:red / bkgd:blck) - time[seconds]');
     213 else
     214        xlabel('SPM (sgnl:green / bkgd:cyan) - time[seconds]');
     215end
    211216       
    212217        print( '-dpng', [ imgName '.png' ]);

  • SPESO/ana2015/analyse_day.m

    r626 r634  
    167167end %if nFullShots>0
    168168
    169 %plot erroros and send alert if something wrong
     169%plot errors and send alert if something wrong
    170170analyse_dayInfo(dateStr, 200);
    171171
    … …  
    281281        %============================= if full data exists      =================================
    282282        if (nFullShots>0)
     283             fprintf(fid, [ ' nFullShots: ' nFullShots '</Br>' ]);
    283284          for ifig=[1 501:504 ]
    284285           imgname=[ 'analysis_figure_fullinfo_' num2str(ifig) '.png' ];
    … …  
    293294                                                                  '</A>\n' ] );
    294295         end % for each figure
     296           else
     297             fprintf(fid, [ ' <P> Full day info not yet available </P>' ]);
    295298        end %if nFullShots>0
    296299        % ----------------------------------- full data --------------------------------------

  • SPESO/ana2015/analyse_dayFullInfo.m

    r626 r634  
    2020        define_data_var();
    2121        define_alerts_var();
     22        cosmetics
    2223
    2324        % turn off figure windows
    2425        %set(0, 'DefaultFigureVisible', 'off');
     26        %set(0, 'DefaultFigureVisible', 'on');
    2527
    2628        data_types=[dvar.RESEAU_POS dvar.DETECT_SPOS dvar.DETECT_XPOS dvar.DETECT_ZPOS dvar.DETECT_RS1 dvar.DETECT_RS2 dvar.DETECT_RS3 dvar.DETECT_RS4 ];
    … …  
    3638        end                     
    3739
     40        figure(401)
     41        data_plot2d(dvar.TIMESTAMP,dvar.RESEAU_POS)
     42
    3843        for dataloop=1:length(data_types)
    3944                allData.selection=['Data for ' date_str(1:4) '/' date_str(5:6) '/' date_str(7:8) char(10)];
    4045       
    4146                %================ Plot vars vs time ====================
    42                 fih=figure;
     47          fih=figure(figure_idx+1);
     48                        clf
     49                        set(0,'DefaultAxesFontSize',14)
     50                        data_plot2d(dvar.TIMESTAMP,data_types(dataloop))
     51                        set(fih,'DefaultAxesFontSize',14)
     52                        close(fih)
     53          fih=figure(figure_idx);
    4354                        clf
    4455                        data_plot2d(dvar.TIMESTAMP,data_types(dataloop))
     56                        set(gcf,'AxesFontSize',14)
    4557                        set(gcf,'PaperUnits','inches','PaperPosition',[0 0 24 18])
    4658                        print(fih, '-dpng', '-r300',[ basedirfigures 'analysis_figure_fullinfo_' num2str(figure_idx+dataloop) '.png' ]);
     59pause(3)
    4760                close(fih);
    4861                %-----------------------------------------------------

  • SPESO/ana2015/analyse_info_one_day.m

    r443 r634  
    1717global dvar;
    1818define_data_var();
    19 % commented since here
    20 %{
    21 figure(1)
    22 data_histo(dvar.RESEAU_POS)
    23 figure(2)
    24 data_histo(dvar.DETECT_SPOS)
    25 figure(3)
    26 data_histo(dvar.DETECT_XPOS)
    27 figure(4)
    28 data_histo(dvar.DETECT_ZPOS)
    29 figure(5)
    30 data_histo(dvar.DETECT_R2)
    31 figure(6)
    32 data_histo(dvar.DETECT_R1)
    33 figure(7)
    34 data_histo(dvar.LPM)
    35 figure(8)
    36 data_histo(dvar.XBPM)
    37 figure(9)
    38 data_histo(dvar.QICT1)
    39 
    40 
    41 figure(101)
    42 data_plot2d(dvar.RESEAU_POS,dvar.DATA_PEAK_LPM)
    43 
    44 data_select(dvar.RESEAU_POS,'>',1.2e6)
    45 
    46 figure(111)
    47 data_plot2d(dvar.DETECT_SPOS,dvar.DATA_PEAK_LPM)
    48 figure(112)
    49 %data_plot2d(dvar.DETECT_XPOS,dvar.DATA_PEAK_LPM)
    50 %figure(113)
    51 data_plot2d(dvar.QICT1,dvar.DATA_PEAK_LPM)
    52 
    53 
    54 plotsdirday=[ plots_basedir dateStr  '/' ];
    55 
    56 for ifig=[1:9,101,111:112]
    57   figure(ifig)
    58   print('-dpng', [ plotsdirday 'analysis_figure_' num2str(ifig) '.png' ]);
    59 end % for each figure
    60 %}
    61 % end of comment
    6219data_select_reset();
    6320fih=figure(301);
    … …  
    6724close(fih);
    6825
     26
    6927end % function analyse_one_day

  • SPESO/ana2015/error_sp.m

    r316 r634  
    11function error_sp(txt)
    2     dbstack();
     2  st=dbstack();
     3txtst='';
     4  for iloop=2:length(st)
     5              txtst=[ txtst st(iloop).file ':' num2str(st(iloop).line) ' ' ];
     6  end
    37    fid=fopen('/sps/sera/delerue/smith-purcell/SPESO/data/warning_sp.txt','a+');
    4 fprintf(fid,'[Error] %s :\t %s \n', datestr(now), txt  );
     8fprintf(fid,'[Error] %s :\t %s \t { %s} \n', datestr(now), txt, txtst  );
    59    fclose(fid);
    610    error(txt);

  • SPESO/ana2015/get_parameters.m

    r619 r634  
    33  dateNum=str2double(dateStr);
    44  hourNum=str2double(hourStr);
    5 if (dateNum>=20170101)
     5if (dateNum>=20180101)
    66  error_sp([ 'NO parameters defined on ' dateStr ]);
     7elseif (dateNum>=20161122)
     8  myParameters.ch(1).on=1;
     9 
     10  myParameters.ch(1).SPM.sigStart=2.40e-8;
     11  myParameters.ch(1).SPM.sigEnd=2.8e-8;
     12  myParameters.ch(1).SPM.riseTime=1.0e-8;
     13
     14  myParameters.ch(1).LPM.sigStart=-1e-8;
     15  myParameters.ch(1).LPM.sigEnd=1e-8;
     16  myParameters.ch(1).LPM.riseTime=1e-8;
     17
     18  myParameters.ch(2).on=1;
     19 
     20  myParameters.ch(2).SPM.sigStart=1.45e-8;
     21  myParameters.ch(2).SPM.sigEnd=1.85e-8;
     22  myParameters.ch(2).SPM.riseTime=1.0e-8;
     23
     24  myParameters.ch(2).LPM.sigStart=-1e-8;
     25  myParameters.ch(2).LPM.sigEnd=1e-8;
     26  myParameters.ch(2).LPM.riseTime=1e-8;
     27
     28  myParameters.ch(3).on=0;
     29 
     30  myParameters.ch(3).SPM.sigStart=-0.73e-7;
     31  myParameters.ch(3).SPM.sigEnd=-0.61e-7;
     32  myParameters.ch(3).SPM.riseTime=2.0e-8;
     33
     34  myParameters.ch(3).LPM.sigStart=-1e-8;
     35  myParameters.ch(3).LPM.sigEnd=1e-8;
     36  myParameters.ch(3).LPM.riseTime=1e-8;
     37
     38  myParameters.ch(4).on=0;
    739elseif (dateNum>=20160801)
    840  myParameters.ch(1).on=1;
    … …  
    2759
    2860  myParameters.ch(3).on=0;
    29   myParameters.ch(4).on=0;
     61 
     62  myParameters.ch(3).SPM.sigStart=-0.73e-7;
     63  myParameters.ch(3).SPM.sigEnd=-0.61e-7;
     64  myParameters.ch(3).SPM.riseTime=2.0e-8;
     65
     66  myParameters.ch(3).LPM.sigStart=-1e-8;
     67  myParameters.ch(3).LPM.sigEnd=1e-8;
     68  myParameters.ch(3).LPM.riseTime=1e-8;
     69
     70  myParameters.ch(4).on=0;
     71
    3072elseif (dateNum>=20160910)
    3173  myParameters.ch(1).on=1;

  • SPESO/ana2015/load_analysis_data.m

    r621 r634  
    154154    dayData.nFullShots=[ dayData.nFullShots 0 ] ;
    155155end % load full info
     156disp( [ ' => '  num2str(length(allData.selected)) ' events loaded so far...'])
    156157end % for each line
    157158

  • SPESO/ana2015/makefile_all

    r435 r634  
    1111                  uptime
    1212                  cd $< && $(MAKE) -f '$(SVN_DIR)'/makefile_day target_no_info
     13                 tail -100 warning_sp.txt > warning_sp.txt
    1314
    1415%.dayFullInfo: %_scope

  • SPESO/ana2015/parameters_from_time_fun.m

    r259 r634  
     1error_sp('Replaced by get_parameter');
     2
    13%sigDuration of the LPM has to be pair!
    24function [myParameters] = parameters_from_time_fun(dateStr, hourStr)

  • SPESO/ana2015/warning_sp.m

    r345 r634  
    44
    55        define_dirs_analysis();
    6                
     6        st=dbstack();
     7txtst='';
     8  for iloop=2:length(st)
     9              txtst=[ txtst st(iloop).file ':' num2str(st(iloop).line) ' ' ];
     10  end
    711        fid=fopen([basedir 'warning_sp.txt'],'a+');
    8           fprintf(fid,'[Warning] %s :\t %s \n', datestr(now), txt  );
     12fprintf(fid,'[Warning] %s :\t %s \t { %s} \n', datestr(now), txt, txtst  );
    913        fclose(fid);
    1014        warning(txt);
Note: See TracChangeset for help on using the changeset viewer.