Changeset 621 in ETALON for SPESO


Ignore:
Timestamp:
Oct 11, 2016, 3:00:57 PM (8 years ago)
Author:
delerue
Message:

Add log var

Location:
SPESO/ana2015
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • SPESO/ana2015/all_data_summary.m

    r620 r621  
    213213
    214214header_fun(fid,[ 'All days summary - ' mask ' -  recent ref: ' datestr(recentDateNum) ' - run on ' datestr(now) ])
     215
     216
     217fprintf(fid,'<script language="javascript" type="text/javascript">\n');
     218fprintf(fid,'         <!-- \n');
     219fprintf(fid,['var hoursSinceLoaded = Math.round((((new Date).getTime()/1000)-' num2str((now-datenum(1970,1,1))*86400) '+7200)/(60*60));\n' ]);   
     220fprintf(fid,'if (hoursSinceLoaded<4){\n document.write("<font color=green size=2>");\n} else {\n document.write("<font color=red size=6>");\n}\n');
     221        fprintf(fid,'document.write("Page generated "+ hoursSinceLoaded+" h ago</font>")\n');
     222        fprintf(fid,' //-->\n');
     223        fprintf(fid,'</script>\n');
     224fprintf(fid,'<small><A HREF=latest_process.txt>Latest process start</A></small>\n');
     225fprintf(fid,'<BR/>');
     226
    215227
    216228fwarnname=[ plots_basedir 'warn_3days.html' ];

  • SPESO/ana2015/data_get_vals.m

    r607 r621  
    11function [valsused,vals]=data_get_vals(var1,var2)
     2
     3%%%Each variable constant is associated to a value through this function
     4%%% ND 2016
     5
    26
    37global allData;
    … …  
    913idx=find(allData.selected==1);
    1014
    11 varidx(1)=mod(var1,1000);
    12 vartype(1)=floor(var1/1000)*1000;
     15nloop=0;
     16for iloop=1:2       
     17  if (nargin>=iloop)
     18     if (iloop==1)
     19       varval(iloop)=var1;
     20     elseif (iloop==2)
     21     varval(iloop)=var2;
     22     else
     23       error('Error in for loop');
     24     end
     25varval
     26     if (varval(iloop)>dvar.LOG_VAR)
     27       log_requested(iloop)=1;
     28       varval(iloop)=varval(iloop)-dvar.LOG_VAR;
     29     else
     30       log_requested(iloop)=0;
     31     end 
     32varval
     33     varidx(iloop)=mod(varval(iloop),1000)
     34     vartype(iloop)=floor(varval(iloop)/1000)*1000
     35     nloop=iloop;
     36   end %if
     37end %for
    1338
    14 if (nargin==2)
    15 varidx(2)=mod(var2,1000);
    16 vartype(2)=floor(var2/1000)*1000;
    17 nloop=2;
    18  else
    19    nloop=1;
    20 end
     39
     40
     41
     42
    2143
    2244if (~isempty(allData.selection))
    … …  
    3456for iloop=1:nloop
    3557switch vartype(iloop)
    36         case dvar.BEAMINFO
    37                 vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))'; %'
    38         case dvar.DATAINFO
    39                 if (iloop==1)
    40                         loopvar=var1;
    41            else
    42               loopvar=var2;
    43       end
    44       switch loopvar
    45                         case dvar.DATA_VAL_LPM
    46                                  valsused(iloop,1)=1;
    47                                  valsused(iloop+2,1)=1;
    48                                  vals(iloop,:)=allData.ch1signalValLPM(idx);
    49                                  vals(iloop+2,:)=allData.ch2signalValLPM(idx);
    50                         case dvar.DATA_PEAK_LPM
    51                                  valsused(iloop,1)=1;
    52                                  valsused(iloop+2,1)=1;
    53                                  vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
    54                                  vals(iloop+2,:)=allData.ch2signalNegPeakLPM(idx);
    55                         case dvar.DATA_VAL_SPM
    56                                  valsused(iloop,1)=1;
    57                                  valsused(iloop+2,1)=1;
    58                                  vals(iloop,:)=allData.ch1signalValSPM(idx);
    59                                  vals(iloop+2,:)=allData.ch2signalValSPM(idx);
    60                         case dvar.DATA_PEAK_SPM
    61                                  valsused(iloop,1)=1;
    62                                  valsused(iloop+2,1)=1;
    63                                  vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
    64                                  vals(iloop+2,:)=allData.ch2signalNegPeakSPM(idx);
     58    case dvar.BEAMINFO
     59        vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))'; %'
     60    case dvar.DATAINFO       
     61            switch varval(iloop)
     62            case dvar.DATA_VAL_LPM
     63            warning_sp('Using dvar.DATA_VAL_LPM is obsolete; use dvar.DATA1_VAL_LPM or dvar.DATA2_VAL_LPM instead!!!')
     64            valsused(iloop,1)=1;
     65            valsused(iloop+2,1)=1;
     66                vals(iloop,:)=allData.ch1signalValLPM(idx);
     67                vals(iloop+2,:)=allData.ch2signalValLPM(idx);
     68            case dvar.DATA_PEAK_LPM
     69            warning_sp('Using dvar.DATA_PEAK_LPM is obsolete; use dvar.DATA1_PEAK_LPM or dvar.DATA2_PEAK_LPM instead!!!')
     70            valsused(iloop,1)=1;
     71            valsused(iloop+2,1)=1;
     72                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
     73                vals(iloop+2,:)=allData.ch2signalNegPeakLPM(idx);
     74            case dvar.DATA_VAL_SPM
     75            warning_sp('Using dvar.DATA_VAL_SPM is obsolete; use dvar.DATA1_VAL_SPM or dvar.DATA2_VAL_SPM instead!!!')
     76            valsused(iloop,1)=1;
     77            valsused(iloop+2,1)=1;
     78                vals(iloop,:)=allData.ch1signalValSPM(idx);
     79                vals(iloop+2,:)=allData.ch2signalValSPM(idx);
     80            case dvar.DATA_PEAK_SPM
     81            warning_sp('Using dvar.DATA_PEAK_SPM is obsolete; use dvar.DATA1_PEAK_SPM or dvar.DATA2_PEAK_SPM instead!!!')
     82            valsused(iloop,1)=1;
     83            valsused(iloop+2,1)=1;
     84                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
     85                vals(iloop+2,:)=allData.ch2signalNegPeakSPM(idx);
    6586
    66                         case dvar.DATA1_VAL_LPM
    67                                 valsused(iloop,1)=1;
    68                                 vals(iloop,:)=allData.ch1signalValLPM(idx);
    69                         case dvar.DATA1_PEAK_LPM
    70                                 valsused(iloop,1)=1;
    71                                 vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
    72                         case dvar.DATA1_VAL_SPM
    73                                 valsused(iloop,1)=1;
    74                                 vals(iloop,:)=allData.ch1signalValSPM(idx);
    75                         case dvar.DATA1_PEAK_SPM
    76                                 valsused(iloop,1)=1;
    77                                 vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
     87            case dvar.DATA1_VAL_LPM
     88                valsused(iloop,1)=1;
     89                vals(iloop,:)=allData.ch1signalValLPM(idx);
     90            case dvar.DATA1_PEAK_LPM
     91                valsused(iloop,1)=1;
     92                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
     93            case dvar.DATA1_VAL_SPM
     94                valsused(iloop,1)=1;
     95                vals(iloop,:)=allData.ch1signalValSPM(idx);
     96            case dvar.DATA1_PEAK_SPM
     97                valsused(iloop,1)=1;
     98                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
    7899
    79                         case dvar.DATA2_VAL_LPM
    80                                 valsused(iloop,1)=1;
    81                                 vals(iloop,:)=allData.ch2signalValLPM(idx);
    82                         case dvar.DATA2_PEAK_LPM
    83                                 valsused(iloop,1)=1;
    84                                 vals(iloop,:)=allData.ch2signalNegPeakLPM(idx);
    85                         case dvar.DATA2_VAL_SPM
    86                                 valsused(iloop,1)=1;
    87                                 vals(iloop,:)=allData.ch2signalValSPM(idx);
    88                         case dvar.DATA2_PEAK_SPM
    89                                 valsused(iloop,1)=1;
    90                                 vals(iloop,:)=allData.ch2signalNegPeakSPM(idx);
    91                         % Select time in minutes       
    92                         case dvar.DATA_TIME_MINUTES
    93                                 valsused(iloop,1)=1;
    94                                 %convert time HHMMSS to minutes
    95                                 vals(iloop,:)=( fix(allData.hoursNum(idx)/10000)*60+...
     100            case dvar.DATA2_VAL_LPM
     101                valsused(iloop,1)=1;
     102                vals(iloop,:)=allData.ch2signalValLPM(idx);
     103            case dvar.DATA2_PEAK_LPM
     104                valsused(iloop,1)=1;
     105                vals(iloop,:)=allData.ch2signalNegPeakLPM(idx);
     106            case dvar.DATA2_VAL_SPM
     107                valsused(iloop,1)=1;
     108                vals(iloop,:)=allData.ch2signalValSPM(idx);
     109            case dvar.DATA2_PEAK_SPM
     110                valsused(iloop,1)=1;
     111                vals(iloop,:)=allData.ch2signalNegPeakSPM(idx);
     112            % Select time in minutes   
     113            case dvar.DATA_TIME_MINUTES
     114                valsused(iloop,1)=1;
     115                %convert time HHMMSS to minutes
     116                vals(iloop,:)=( fix(allData.hoursNum(idx)/10000)*60+...
    96117                                fix(mod(allData.hoursNum(idx),10000)/100)+...
    97118                                mod(allData.hoursNum(idx),100)/60);
    98                         case dvar.DATA_TIME_30
    99             valsused(iloop,1)=1;
    100             %convert time HHMMSS to minutes
    101             vals(iloop,:)=mod(( fix(allData.hoursNum(idx)/10000)*60+...
     119            case dvar.DATA_TIME_30
     120                valsused(iloop,1)=1;
     121                %convert time HHMMSS to minutes
     122                vals(iloop,:)=mod(( fix(allData.hoursNum(idx)/10000)*60+...
    102123                                fix(mod(allData.hoursNum(idx),10000)/100)+...
    103124                                mod(allData.hoursNum(idx),100)/60),30);
    104125           
    105                         case dvar.DATA_RS1_OFFSET
    106                                 RS_params=get_RS_params(allData.timestamp(1,idx));
    107                                 vals(iloop,:)=allData.beamInfo(idx,5)'-RS_params.RS1_offset;
    108                         case dvar.DATA_RS2_OFFSET
    109                                 RS_params=get_RS_params(allData.timestamp(1,idx));
    110                                 vals(iloop,:)=allData.beamInfo(idx,7)'-RS_params.RS2_offset;
    111                         case dvar.DATA_RS3_OFFSET
    112                                 RS_params=get_RS_params(allData.timestamp(1,idx));
    113                                 vals(iloop,:)=allData.beamInfo(idx,9)'-RS_params.RS3_offset;
    114                         case dvar.DATA_RS4_OFFSET
    115                                 RS_params=get_RS_params(allData.timestamp(1,idx));
    116                                 vals(iloop,:)=allData.beamInfo(idx,11)'-RS_params.RS4_offset;
    117                        
    118                         case dvar.DATA_DET1_THETA
    119                                 SPOS_offset=-645000; % position of the grating in TS
    120                                 X0_pos=220000; % position of the grating in TX
    121                                 raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
    122                                                                 (allData.beamInfo(idx,13)-SPOS_offset)./...
    123                                                                 0.6.*0.5)*180/pi;
    124                                 % 0.6 - coeff between TS and mm
    125                                 % 0.5 - coeff between TX and mm
    126                                 vals(iloop,:)=(-raw_th+(1+sign(raw_th))/2*180);
    127                                  % conver from 0:-pi/2 - pi/2:0 to 0:pi
    128                         case dvar.DATA_DET2_THETA
    129                                 SPOS_offset=-470000; % position of the grating in TS
    130                                 X0_pos=220000; % position of the grating in TX
    131                                 raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
    132                                                                 (allData.beamInfo(idx,13)-SPOS_offset)./...
    133                                                                 0.6.*0.5)*180/pi;
    134                                 % 0.6 - coeff between TS and mm
    135                                 % 0.5 - coeff between TX and mm
    136                                 vals(iloop,:)=(-raw_th+(1+sign(raw_th))/2*180);
    137                                  % conver from 0:-pi/2 - pi/2:0 to 0:pi
    138                         case dvar.DATA_DET1_PHI
    139                                 ZPOS_offset=-66000; % position of the grating in TS
    140                                 X0_pos=220000; % position of the grating in TX
    141                                 raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
    142                                                                 (allData.beamInfo(idx,17)-ZPOS_offset)./...
    143                                                                 2.48.*0.5)*180/pi;
    144                                 % 2.48 - coeff between TZ and mm
    145                                 % 0.5 - coeff between TX and mm
    146                                 vals(iloop,:)=raw_th-sign(raw_th)*90;
    147                                  % conver from -0:pi/2 - pi/2:0 to -pi/2:pi/2
    148                         case dvar.DATA_DET2_PHI
    149                                 ZPOS_offset=-66000; % position of the grating in TS
    150                                 X0_pos=220000; % position of the grating in TX
    151                                 raw_th=atan((X0_pos-allData.beamInfo(idx,15))./...
    152                                                                 (allData.beamInfo(idx,17)-ZPOS_offset)./...
    153                                                                 2.48.*0.5)*180/pi;
    154                                 % 2.48 - coeff between TZ and mm
    155                                 % 0.5 - coeff between TX and mm
    156                                 vals(iloop,:)=raw_th-sign(raw_th)*90;
    157                                  % conver from 0:-pi/2 - pi/2:0 to -pi/2:pi/2
    158 
    159                         otherwise
    160                                 warning_sp([ 'The varidx ' varidx ' is unknown.' ]);
    161                                 error('end');
    162                         end %switch
     126            otherwise
     127            warning_sp([ 'The varidx ' num2str(varidx(iloop)) ' is unknown.' ]);
     128                error('end');
     129            end %switch
    163130    otherwise
    164                         warning_sp([ 'The vartype ' vartype ' is unknown.' ]);
    165                         error('end');
    166         end %switch
     131            warning_sp([ 'The vartype ' num2str(vartype) ' is unknown.' ]);
     132        error('end');
     133end %switch
     134if (log_requested(iloop)==1)
     135  vals(iloop,:)=log(abs(vals(iloop,:)));
     136 end
    167137end %for
    168138

  • SPESO/ana2015/data_histo.m

    r321 r621  
    66global allData;
    77
    8 %standart number of bins in histigram
     8%standart number of bins in histogram
    99binsNum=20;
    1010%load values

  • SPESO/ana2015/define_data_var.m

    r607 r621  
    66dvar.DATAINFO=2000;
    77dvar.RAWINFO=3000;
     8
     9dvar.LOG_VAR=10000;
    810
    911dvar.BTIME=dvar.BEAMINFO+1;

  • SPESO/ana2015/get_var_name.m

    r607 r621  
    33global dvar;
    44define_data_var();
     5
     6if (varidx>dvar.LOG_VAR)
     7  strlog='Log of ';
     8  varidx=varidx-dvar.LOG_VAR;
     9 else
     10   strlog='';
     11end
    512
    613switch(varidx)
    … …  
    95102end
    96103
     104str=[ strlog str];
     105
    97106end %function

  • SPESO/ana2015/header_fun.m

    r620 r621  
    99fprintf(fid,'</H2>\n');
    1010fprintf(fid,'</CENTER>\n');
    11 fprintf(fid,'<script language="javascript" type="text/javascript">\n');
    12 fprintf(fid,'         <!-- \n');
    13 fprintf(fid,['var hoursSinceLoaded = Math.round((((new Date).getTime()/1000)-' num2str((now-datenum(1970,1,1))*86400) '+7200)/(60*60));\n' ]);   
    14 fprintf(fid,'if (hoursSinceLoaded<4){\n document.write("<font color=green size=2>");\n} else {\n document.write("<font color=red size=6>");\n}\n');
    15         fprintf(fid,'document.write("Page generated "+ hoursSinceLoaded+" h ago</font>")\n');
    16         fprintf(fid,' //-->\n');
    17         fprintf(fid,'</script>\n');
    18 fprintf(fid,'<small><A HREF=latest_process.txt>Latest process start</A></small>\n');
    19 fprintf(fid,'<BR/>');
    2011
    2112

  • SPESO/ana2015/load_analysis_data.m

    r620 r621  
    3434    end
    3535      else
    36         disp('No lsmask');
     36        %%%disp('No lsmask');
    3737    end
    3838  end

  • SPESO/ana2015/makefile_day

    r620 r621  
    3939$(ROOTDAYFULLINFO): $(DAYFULLINFO) ../$(DATEDIR)_scope
    4040                  echo Datedir $(DATEDIR)
     41                  date
     42                  date > $(DATEFILE)
    4143                  ln -f $(DAYINFO) $(ROOTDAYINFO)           
    4244                  ln -f $(DAYFULLINFO) $(ROOTDAYFULLINFO)           
Note: See TracChangeset for help on using the changeset viewer.