half EBT3 pseudo implement.

RCF_energyloss_GUI_2016.m

@@ -795,6 +795,11 @@ if handles.parallel==0
                  v(1,i) = E;
                  v(2,i) = E_stop;
                  v(4,i) = E_stop_corr;
+              case 'ES'
+                 [E_new, E_stop, E_stop_corr] = elo_EBT3half(E_new, loss.([ion '_in_Poly']), loss.([ion '_in_ActiveEBT3']), projmass);
+                 v(1,i) = E;
+                 v(2,i) = E_stop;
+                 v(4,i) = E_stop_corr;
               otherwise
                  set(handles.out_text,'String','Unknown layer.')
                  break
@@ -1017,8 +1022,11 @@ for i=1:length(handles.parameters.rcf_config)
         case 'HD'
             imgfile = ['icons/' 'hd810' '.png'];
             mName = 'HD';
-        case {'MD', 'M2'}
+        case 'MD'
             imgfile = ['icons/' 'md55' '.png'];
+            mName = 'MD';
+        case 'M2'
+            imgfile = ['icons/' 'md55v2' '.png'];
             mName = 'M2';
         case 'HS'
             imgfile = ['icons/' 'hs' '.png'];
@@ -1038,6 +1046,9 @@ for i=1:length(handles.parameters.rcf_config)
         case 'E3'
             imgfile = ['icons/' 'ebt3' '.png'];
             mName = 'E3';
+        case 'ES'
+            imgfile = ['icons/' 'ebt3half' '.png'];
+            mName = 'ES';
         otherwise
             set(handles.out_text,'String','Unknown layer.')
             break
@@ -1158,7 +1169,7 @@ for i=1:length(handles.parameters.rcf_config)
             save(filename,'data', '-ascii', '-double');
         end
          j=j+1;
-      case {'HD','MD','M2','HS','H2','H3','E3'}
+      case {'HD','MD','M2','HS','H2','H3','E3','ES'}
          corroffhdl=findobj('Tag', 'corroff');
          corronhdl=findobj('Tag', 'corron');
          if (get(corronhdl,'Value')==1)

elo_EBT3half.m

+function [E_new, E_stop, E_stop_corr] = elo_EBT3half(E_new, H_in_Poly, H_in_ActiveEBT3, projmass)
+    % Routine to calculate the energyloss and energydeposition in EBT3
+    
+    thickness_poly_ebt3 = 125; % microns
+    thickness_active_ebt3 = 28; % microns
+
+    %assume same correction like active layer in HD => needs to be checked
+    %experimentaly
+    tmp = elo(E_new, thickness_active_ebt3, H_in_ActiveEBT3.x, H_in_ActiveEBT3.c, H_in_ActiveEBT3.crange, H_in_ActiveEBT3.xrangeinv, H_in_ActiveEBT3.crangeinv, H_in_ActiveEBT3.k, H_in_ActiveEBT3.Nele, H_in_ActiveEBT3.Zeff, H_in_ActiveEBT3.Ieff, projmass, 1, H_in_ActiveEBT3.ccorr);
+    E_new = tmp(:,1);
+    E_stop = mean(tmp(:,2));
+    E_stop_corr = mean(tmp(:,3));
+    
+    tmp = elo(E_new, thickness_poly_ebt3, H_in_Poly.x, H_in_Poly.c, H_in_Poly.crange, H_in_Poly.xrangeinv, H_in_Poly.crangeinv, H_in_Poly.k, H_in_Poly.Nele, H_in_Poly.Zeff, H_in_Poly.Ieff, projmass, 0);
+    E_new = tmp(:,1);
+end