Difference between revisions of "PAA Software"

From New IAC Wiki
Jump to navigation Jump to search
 
(One intermediate revision by the same user not shown)
Line 11: Line 11:
 
Christian is a good resource for questions on Gamma W
 
Christian is a good resource for questions on Gamma W
  
2.)
+
2.)Valeriia's matlab script used to estimate activation times
  
Valeriia's matlab script used to estimate activation times
+
[[Vstar_Matlab_ActTime]]
  
[[File:Vstar_Matlab_ActTime.txt]]
 
 
This example was done for Bi-206.
 
 
%%%%%%%%%%%%%%%%%%
 
clear all;
 
 
load bics.m % loads cross-section file: first column is energy, second column is c-s in b
 
load flpb.m %loads photon flux file: first column is energy, second column is flux in photons/cm^2 per source electron
 
energy=bics(:, 1);
 
flux=flpb(:, 2);
 
cs=bics(:,2);
 
 
 
 
%% Parameters (target parameters) (need to be entered every time program is compiled)
 
 
% Below enter target mass in grams;
 
M = 1;
 
% Below enter target average atomic mass in amu;
 
m = 208;
 
% Below enter natural abundance of an isotope of interest;
 
abundance = 0.55;
 
% Below enter a half-life of the isotope of interest in hours:
 
halflife = 6.24*24;
 
% %If you want to know activity after certain time of irradiation, then below enter irradiation time in hours:
 
T = 1;
 
% %If you want to see activity as a function of time, then below enter irradiation time in hours:
 
t = [0:0.1:100];
 
%% Activity calculation is conducted below (don't edit)
 
 
% Average number of electrons per second is calculated and stored in variable N (for 1 kW of
 
% beam power and 40 MeV electron beam):
 
N = (1 * 10^(-3)) / (40 * 1.6 * 10^(-19));
 
 
% Photon flux - cross section overlap integral I (in s^-1) is calculated and stored in variable I:
 
 
 
I=0;
 
for i=1:length(cs),
 
I=I+flux(i).*cs(i);
 
 
end
 
I = I*N * 10^(-27);
 
% Total number of the isotope of interest in the target is calculated and stored in variable Tot;
 
Tot = M * abundance / (m * 1.6605 * 10^(-24) );
 
 
% Production rate (in s^-1) is calculated and stored in variable R:
 
R = Tot * I;
 
 
% Activity obtained in uCi:
 
A0 = ( R * (1 - exp(-log(2)* T / halflife))) / (37000); %converted from Bq into microCi
 
 
sprintf('Activity obtained is %g uCi/(%g g kW %g hours)', A0, T)
 
 
A = ( R * (1 - exp(-log(2)* t / halflife))) / (37000); %converted from Bq into microCi
 
plot(t,A)
 
  
 
[[PAA_Research]]
 
[[PAA_Research]]

Latest revision as of 19:53, 3 April 2014

In Order of preference

1.) Gamma W


Need to check out a USB dive which contains the key needed to run the program.


Can read MPA files saved as *.MP not *.MPA

Christian is a good resource for questions on Gamma W

2.)Valeriia's matlab script used to estimate activation times

Vstar_Matlab_ActTime


PAA_Research