getisotopephotopeak(isotope)

function ppeak = getisotopephotopeak(isotope)
% Note that the unit for the halflife is minute
% Reference:Physics of Medical Imaging J.T. Bushberg,et al, 2nd Edition
switch isotope
    case 'Tc-99m'
        ppeak = [140.5];
    case 'Co-57'
        ppeak = [122,136];
    case 'Co-68'
        ppeak = [122,136];
    case 'F-18'
        ppeak = 511;
    case 'N-13'
        ppeak = 511;
    case 'C-11'
        ppeak = 511;
    case 'O-15'
        ppeak = 511;
    case 'Ga-67'
        ppeak = [93,184,300];
    case 'Ge-68' % Germanium-68 EC 288 days
        ppeak = 511;
    case 'Ga-68'
        ppeak = 511;
    case 'Gd-153'
        ppeak = [97,103];
    case 'In-111'
        ppeak = [171,245];
    case 'I-123'
        ppeak = 159;
    case 'I-131'
        ppeak = [284,364,637];
    case 'I-124'
        ppeak = 511;
    case 'Tl-201'
        ppeak = 167;
    case 'Rb-82'
         ppeak = 511;
    case 'Cu-62'
         ppeak = 511;
    case 'Br-76'
        ppeak = 511;
    case 'Na-22'
        ppeak = 511;
    case 'Xe-133'
        ppeak = 81;
    case 'Cs-137'
        ppeak = 662;
    case 'In-113m'
        ppeak = 392;
    case 'Sn-113m'
        ppeak = 0;
    case 'Sr-82'
        ppeak = 0;
    case 'Kr-81m'
        ppeak = 0;
    case 'Zn-62'
        ppeak = 0;
    case 'P-32'
        ppeak = 1700;
    case 'Sm-153'
        ppeak = [70,103];
    case 'Y-90'
        ppeak = 1760;
    otherwise
        errordlg('unknown isotope');
        ppeak = [];
end