Simulations of the pile-up effect on HXR photon energy spectrum as function of the count rate by Y.Peysson CEA-DRFC <yves.peysson@cea.fr> and Joan Decker CEA-DRFC (joan.decker@cea.fr)
0001 %Simulations of the pile-up effect on HXR photon energy spectrum as 0002 %function of the count rate 0003 % 0004 %by Y.Peysson CEA-DRFC <yves.peysson@cea.fr> and Joan Decker CEA-DRFC (joan.decker@cea.fr) 0005 % 0006 clear all 0007 pileup_id = 'TCV';%Set here the name (identifier) of the HXR diagnostic electronic system 0008 % 0009 % Pulse shape parameters 0010 % 0011 diag.tok = 'TCV'; 0012 diag.namediag = ''; 0013 diag.theta = 0.18;%Shaping time (microsecond) 0014 diag.tc0 = 0.0;%Coding time ADC (microsecond) 0015 diag.tpeak = 0.4; 0016 diag.tl0 = 7.0; 0017 diag.tp0 = 2.0; 0018 diag.beta = 5;%For pulse shape 0019 diag.kmin = 20;%Lower limit of the photon energy spectrum [kev] 0020 diag.kmin0 = diag.kmin; 0021 diag.kmax = 300;%Upper limit of the photon energy spectrum [kev] 0022 diag.nhisto = 14;%Number of channels 0023 diag.khistotphmin = 50;%Minimum photon energy of the interval for calculating the photon temperature [keV] 0024 diag.khistotphmax = 110;%Maximum photon energy of the interval for calculating the photon temperature [keV] 0025 % 0026 % Pile-up simulation parameters 0027 % 0028 pileupsim.pur = 1;%Pile-up rejection system On (1) or OFF (0) 0029 pileupsim.nreset = 3;%Random number generator reset number 0030 pileupsim.typspec = 2;%Line photon energy spectrum (1) or exponential (2) 0031 pileupsim.Ni0 = 50000;%Plasma photon count rate (photon/s) 0032 pileupsim.ech = 1;%Scaling for the Monte-Carlo calculations (for reducing the computational effort, increasing ech) 0033 pileupsim.k0 = 20;%Characteristic decay energy of the photon spectrum for exponential case or line energy for the line case (keV) 0034 pileupsim.modeMC = 2;%(1) number of input photon for the Monte-Carlo is the pile-up level from theoretical calculations, (0) Monte-Carlo use in addition NiMC 0035 pileupsim.NiMC = [1e5:0.5e5:1e6];%Plasma photon count rate for the Monte-Carlo (photon/s) 0036 pileupsim.seuil = 0.5;%Fraction of pile-up 0037 pileupsim.p_opt = 2;%Figure printing options 0038 pileupsim.filename = pileup_id ;%Figure printing options 0039 % 0040 eval(['save PILEUP_',pileup_id,'.mat diag pileupsim']); 0041