make_pileup_TS

PURPOSE ^

Simulations of the pile-up effect on HXR photon energy spectrum as

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

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)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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 = 'TS-DSPX';%Set here the name (identifier) of the HXR diagnostic electronic system
0008 %
0009 % Pulse shape parameters
0010 %
0011 diag.tok = 'TS';
0012 diag.namediag = 'DSPX';
0013 diag.theta = 0.25;%Shaping time (microsecond)
0014 diag.tc0 = 3.0;%Coding time ADC (microsecond)
0015 diag.tpeak = 0.5; 
0016 diag.tl0 = 7.3;
0017 diag.tp0 = 2.35;
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 = 180;%Upper limit of the photon energy spectrum [kev]
0022 diag.nhisto = 8;%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 = 30;%Characteristic decay energy of the photon spectrum for exponential case or line energy for the line case  (keV)
0034 pileupsim.modeMC = 1;%(1) number of input photon for the Monte-Carlo is the pile-up level from theoretical calculations, (0) Monte-Carlo use NiMC
0035 pileupsim.NiMC = 100000;%Plasma photon count rate for the Monte-Carlo (photon/s)
0036 pileupsim.seuil = 0.6;%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

Community support and wiki are available on Redmine. Last update: 18-Apr-2019.