|
|
| (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]] |
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