Changeset 644 in ETALON for SPESO

Timestamp:

Apr 3, 2017, 9:54:44 AM (7 years ago)

Author:

hodnevuc

Message:

 

File:

• SPESO/ana2015/analyse_day.m (modified) (5 diffs)

SPESO/ana2015/analyse_day.m

@@ -52 +52 @@
 end % if ch(1) is on 
 
+if((theShot.shotInfo.ch(3).on)==1)
+        dayInfo.ch3signalValSPM(iShot)=dayInfo.shots(iShot).shotInfo.ch(3).values.SPM.signalVal;
+        dayInfo.ch3signalNegPeakSPM(iShot)=dayInfo.shots(iShot).shotInfo.ch(3).values.SPM.signalNegPeak;
+        dayInfo.ch3signalValLPM(iShot)=dayInfo.shots(iShot).shotInfo.ch(3).values.LPM.signalVal;
+        dayInfo.ch3signalNegPeakLPM(iShot)=dayInfo.shots(iShot).shotInfo.ch(3).values.LPM.signalNegPeak;
+end % if ch(3) is on 
         hourStr=line(16:21);
         htmltxt_table{str2num(hourStr(1:2))+1}=[ htmltxt_table{str2num(hourStr(1:2))+1} get_channels_html(theShot.shotInfo,0,line(7:21)) ];
@@ -104 +110 @@
 end % if ch(1) is on 
 
+if((theShot.shotInfo.ch(3).on)==1)
+  if  (~(isempty(theShot.shotInfo.ch(3).values)))
+        dayFullInfo.ch3signalValSPM(iFullShot)=dayFullInfo.shots(iFullShot).shotInfo.ch(3).values.SPM.signalVal;
+        dayFullInfo.ch3signalNegPeakSPM(iFullShot)=dayFullInfo.shots(iFullShot).shotInfo.ch(3).values.SPM.signalNegPeak;
+        dayFullInfo.ch3signalValLPM(iFullShot)=dayFullInfo.shots(iFullShot).shotInfo.ch(3).values.LPM.signalVal;
+        dayFullInfo.ch3signalNegPeakLPM(iFullShot)=dayFullInfo.shots(iFullShot).shotInfo.ch(3).values.LPM.signalNegPeak;
+        else % shot is empty
+                 error_sp(['Empty channel 3 shot on ' dateStr ' shot number ' num2str(iFullShot) ])
+        end % is shot empty 
+end % if ch(3) is on 
         dayFullInfo.beamInfo(iFullShot,:)=theFullShot.shotInfo.beamInfo(1,:);
     else
@@ -127 +143 @@
         thisHourShots=find(floor(dayInfo.hoursNum(:)/10000)==ihour);
         if (~isempty(thisHourShots))
+            figure(plotStr*10000+ihour)
+            clf;
+plot(floor(mod(dayInfo.hoursNum(thisHourShots),10000)/100),dayInfo.ch1signalValSPM(thisHourShots))
+            title(' Ch1 - signal val SPM') 
+            ylabel('Intensity [AU]')
+            xlabel('Minutes in the hour')
+            figure(plotStr*10000+100+ihour)
+            clf;
+plot(floor(mod(dayInfo.hoursNum(thisHourShots),10000)/100),dayInfo.ch1signalNegPeakSPM(thisHourShots))
+            title(' Ch1 - signal Neg Peak SPM') 
+            ylabel('Intensity [AU]')
+            xlabel('Minutes in the hour')
+            figure(plotStr*10000+200+ihour)
+            clf;
+plot(floor(mod(dayInfo.hoursNum(thisHourShots),10000)/100),dayInfo.ch1signalValLPM(thisHourShots))
+            title(' Ch1 - signal val LPM') 
+            ylabel('Intensity [AU]')
+            xlabel('Minutes in the hour')
+            figure(plotStr*10000+300+ihour)
+            clf;
+plot(floor(mod(dayInfo.hoursNum(thisHourShots),10000)/100),dayInfo.ch1signalNegPeakLPM(thisHourShots))
+            title(' Ch1 - signal Neg Peak LPM') 
+            ylabel('Intensity [AU]')
+            xlabel('Minutes in the hour')
             figure(plotStr*1000+ihour)
             clf;
@@ -151 +191 @@
             ylabel('Intensity [AU]')
             xlabel('Minutes in the hour')
+            for ifig=[ 0:3]
+  jfig=ifig*100+plotStr*10000+ihour;
+  figure(jfig)
+  print('-dpng', [ plots_basedir dateStr  '/'  sprintf('%02d',ihour) '/'  'hourly_figure_' num2str(ifig) '_ch1.png' ]);
+            end % for ifig
             for ifig=[ 0:3]
   jfig=ifig*100+plotStr*1000+ihour;
   figure(jfig)
-  print('-dpng', [ plots_basedir dateStr  '/'  sprintf('%02d',ihour) '/'  'hourly_figure_' num2str(ifig) '.png' ]);
+  print('-dpng', [ plots_basedir dateStr  '/'  sprintf('%02d',ihour) '/'  'hourly_figure_' num2str(ifig) '_ch2.png' ]);
             end % for ifig
         end % if length>0
@@ -197 +242 @@
         
         for ifig=[ 0:3]
-          imgname=[ 'hourly_figure_' num2str(ifig) '.png'       ];
+          imgname=[ 'hourly_figure_' num2str(ifig) '_ch1.png'   ];
+          fprintf(fid, [ '<A HREF="' imgname '" >'...
+                         '<IMG SRC="' imgname '" WIDTH=' num2str(imageSizes) ' ></A>\n' ] );
+        end % for ifig
+
+        for ifig=[ 0:3]
+          imgname=[ 'hourly_figure_' num2str(ifig) '_ch2.png'   ];
           fprintf(fid, [ '<A HREF="' imgname '" >'...
                          '<IMG SRC="' imgname '" WIDTH=' num2str(imageSizes) ' ></A>\n' ] );