File:Photon Mass Attenuation Coefficients.svg

From Wikimedia Commons, the free media repository
Jump to navigation Jump to search

Original file(SVG file, nominally 600 × 400 pixels, file size: 64 KB)

Captions

Captions

Add a one-line explanation of what this file represents

Summary[edit]

Description Photon Mass Attenuation Coefficient for photons in energy range from 1 keV to 20 MeV for Elements Z = 1 to 100. Based on [1]. Also shown are locations of photons absorption edges (and not Compton edges as mentioned in the comments from the Matlab code).
Date
Source Own work
Author Jarekt
Other versions PNG version of the same image can be found at Image:Photon Mass Attenuation Coefficients.png.
SVG development
InfoField
 
The SVG code is valid.
 
This diagram was created with MATLAB.
 
The file size of this SVG image may be irrationally large because its text has been converted to paths inhibiting translations.
Source code
InfoField

MATLAB code

figure
Z = 1:100;  % elements with Z in 1-100 range - number of columns
nr = 500;   % number of rows to use in the plot
E = logspace(log10(0.001), log10(20), 500);  % define energy grid
[mac, CEdge] = PhotonAttenuationQ(Z, E);
colormap(jet(128)) % use hi-res color palette
imagesc(log10(mac)); 
grid on;
axis xy; % put small numbers on y axis on the bottom
title('Photon Mass Attenuation Coefficients (in cm^2/g) and Compton Edges');
xlabel('Atomic Number of Elements');
ylabel('Energy in MeV');
% Add X-Axis 
EPos = [6 16 26 35 47 55 65 74 82 94];   % define array to store label location
ELab = { '6-C','16-S','26-Fe','35-Br','47-Ag','55-Cs','65-Tb','74-W','82-Pb','94-Pu'}; %Define Energy labels for y-axis
set(gca,'XTick'     ,EPos);
set(gca,'XTickLabel',ELab);
% Add Y-Axis
ELab = [0.001 0.003 0.01 0.03 0.1 0.3 1 3 10]; %Define Energy labels for y-axis
EPos = size(ELab);   % define array to store label location
for i=1:length(ELab), [tmp EPos(i)]=min(abs(E-ELab(i))); end
set(gca,'YTick'     ,EPos);
set(gca,'YTickLabel',ELab);
% add Colorbar
cbar_axes = colorbar;
set(cbar_axes,'YTick'     ,     -1:4 ); % The image is a log10 of the MAC ...
set(cbar_axes,'YTickLabel',10.^(-1:4)); % ... so add proper labels
hold on
% Add Conpton Edges to the plot
ed = accumarray([CEdge(:,1),CEdge(:,2)],CEdge(:,3)); % get per element energies of 14 compton edges 
ed = 500*(log(ed')-log(0.001))/(log(20)-log(0.001)); % convert energy to row numbers of the image
K=plot(ed(:,1)    ,'m','LineWidth',3); %Plot K Compton edge
L=plot(ed(:, 2: 4),'b','LineWidth',2); %Plot 3 L Compton edges
M=plot(ed(:, 5: 9),'g','LineWidth',2); %Plot 5 M Compton edges
N=plot(ed(:,10:14),'Y','LineWidth',2); %Plot first 5 N Compton edges
legend([K(1),L(1),M(1),N(1)], {'K','L','M','N'}, 'Location', 'NorthWest');
plot2svg('Photon_Mass_Attenuation_Coefficients_.svg')

With help of two external libraries:

  1. PhotonAtenuattion2
  2. Plot2SVG archive copy at the Wayback Machine

Licensing[edit]

I, the copyright holder of this work, hereby publish it under the following licenses:
GNU head Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts. A copy of the license is included in the section entitled GNU Free Documentation License.
w:en:Creative Commons
attribution share alike
This file is licensed under the Creative Commons Attribution-Share Alike 3.0 Unported, 2.5 Generic, 2.0 Generic and 1.0 Generic license.
You are free:
  • to share – to copy, distribute and transmit the work
  • to remix – to adapt the work
Under the following conditions:
  • attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
  • share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
You may select the license of your choice.

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current02:10, 1 October 2007Thumbnail for version as of 02:10, 1 October 2007600 × 400 (64 KB)Jarekt (talk | contribs)
17:07, 26 September 2007Thumbnail for version as of 17:07, 26 September 2007512 × 330 (64 KB)Jarekt (talk | contribs){{Information |Description={{en|Photon '''Mass Attenuation Coefficient''' for photons in energy range from 1 keV to 20 MeV for Elements Z = 1 to 100. Based on http://physics.nist.gov/PhysRefData/XrayNoteB.html}} |Source=self-made |Date=9/27/2007 |Auth