%% birth_pie
% Plots pie diagram for the cumulative energy investment at birth

%%
function [Hfig EMJHG ue0 info] = birth_pie(p, eb)
  % created at 2011/01/19 by Bas Kooijman; modified 2016/05/03, 2017/01/04
  
  %% Syntax
  % [Hfig EMJHG E0 info] = <../birth_pie.m *birth_pie*> (p, eb)
  
  %% Description
  % Shell around <get_EMJHG.html *get_EMJHG*> to draw pies.
  % Plots pie diagram for the cumulative energy investment at birth
  % See  <http://www.bio.vu.nl/thb/research/bib/Kooy2010.pdf *section 4.3.3 of the comments on the DEB-book*>. 
  % The pie pieces are also given numerically as fractions in the output, which sum to 1. 
  % Investments that are fully dissipated are drawn exploded. 
  % If the growth efficiency is specified in p(5), cumulative allocation to growth is splitted into that to growth overhead and structure.  
  %
  % Input
  %
  % * p: 4 or 5-vector with parameters: g, k, vHb, kap, kap_G
  % * eb: optional n-vector with scaled reserve densities at birth (default: 1)
  %
  % Output
  %
  % * Hfig: scalar with figure handle
  % * EMJHG: (n,5 or 6)-matrix with in the columns fractions of initial reserve at birth
  %
  %    - reserve left at birth, cumulatively allocated to som maint, mat maint, maturation, growth 
  %
  % * E0: scaler with initial amount of energy in reserve
  % * info: n-vector with 1's success and 0's otherwise for uE0 and tau_b-computations
  
  %% Remarks
  % See also <birth_pie_foetus.html *birth_pie_foetus*> for foetal development.
  
  %% Example of use
  % See <../mydata_birth_pie.m *mydata_birth_pie*>

  if exist('eb', 'var') == 0
    eb = 1; % maximum value as juvenile
  end
  n = length(eb);
  [EMJHG uE0 info] = get_EMJHG(p, eb);
  
  %close all
  if length(p) == 4
    pie_color = [1 1 .8; 1 0 0; 1 0 1; 0 0 1; 0 1 0];
    for i = 1:n
      pie_txt = { ...
        ['reserve ', num2str(EMJHG(i,1), '% 3.2f')], ['som maint ', num2str(EMJHG(i,2), '% 3.2f')], ['mat maint ', num2str(EMJHG(i,3), '% 3.2f')], ...
        ['maturity', num2str(EMJHG(i,4), '% 3.2f')] ,['growth ', num2str(EMJHG(i,5), '% 3.2f')]};
      Hfig = figure;
      colormap(pie_color);
      set(gca, 'FontSize', 15, 'Box', 'on')
      pie3s(EMJHG(i,:), 'Explode', [0 1 1 1 0], 'Labels', pie_txt, 'Bevel', 'Elliptical');
      title({['cum. investment at birth, e_b = ', num2str(eb(i))], ['E_0 = ', num2str(E0,'% 1.3e'), ' J'] }) 
    end
  else
    pie_color = [1 1 .8; 1 0 0; 1 0 1; 0 0 1; 0 1 0; .8 1 .8];
    for i = 1:n
      pie_txt = { ...
        ['reserve ', num2str(EMJHG(i,1), '% 3.2f')], ['som maint ', num2str(EMJHG(i,2), '% 3.2f')], ['mat maint ', num2str(EMJHG(i,3), '% 3.2f')], ...
        ['maturity ', num2str(EMJHG(i,4), '% 3.2f')],['growth overhead ', num2str(EMJHG(i,5), '% 3.2f')], ['structure ', num2str(EMJHG(i,6), '% 3.2f')]};
      Hfig = figure;
      colormap(pie_color);
      set(gca, 'FontSize', 15, 'Box', 'on')
      pie3s(EMJHG(i,:), 'Explode', [0 1 1 1 1 0], 'Labels', pie_txt, 'Bevel', 'Elliptical');
      title({['cum. investment at birth, e_b = ', num2str(eb(i))], ['E_0 = ', num2str(E0,'% 1.3e'), ' J'] })
    end
  end