To Make Likelihood Fitting

Next: To draw Feynman diagrams Up: Physical Process Previous: To use structure function Contents

To Make Likelihood Fitting

The follwoing three terms should be added to the file "process.def".

$\begin{picture}(490,80) \put(230,35){\oval(480,70)} \put(0,35){ \begin{tabular}{... ...er\_list:='( ..... ); & \%declare fited parameters \end{tabular}} \end{picture}$

Then to run the following command as usual.

$\begin{picture}(490,120) \put(230,50){\oval(480,110)} \put(0,50){ \begin{tabular... ... \\ $>$ make & \% compile and link Fortran source \end{tabular}} \end{picture}$

There are five Executable file generated after "make". The user can choosed some of them to run according to his need.

$\begin{picture}(110,40) \put(50,20){\oval(100,30)} \put(0,20){ \begin{tabular}{l} $>$ int \end{tabular}} \end{picture}$

To integrate cross section and generate distribution. the result of cross section is in the file "fresult.dat" and the convergence behaive is in the file "convergence.dat". The distribution data is in the file "plot.hbook" for pawX11 and the file "plot.dat for directly text view. The information for event generate is in the file "bases.data".

$\begin{picture}(110,40) \put(50,20){\oval(100,30)} \put(0,20){ \begin{tabular}{l} $>$ gevent \end{tabular}} \end{picture}$

To generate event by using the information in the file "bases.data" and the generated 4-momentum data of event is in the file "pdata1.dat". The file "pdata1.dat" is of the fortmat as:

$\begin{picture}(310,80) \put(150,35){\oval(300,70)} \put(0,35){ \begin{tabular}{l} Data of each Event \\ ..... ..... \end{tabular}} \end{picture}$

and each event is of the following format. There is no momentum of initial paticles, just all the momentum of final particles

$\begin{picture}(310,80) \put(150,35){\oval(300,70)} \put(0,35){ \begin{tabular}{... ..._{4,y}$ &$P_{4,z}$ &$E_4$\\ &... &... &... &... \end{tabular}} \end{picture}$

By performing the following program

$\begin{picture}(110,40) \put(50,20){\oval(100,30)} \put(0,20){ \begin{tabular}{l} $>$ kint \end{tabular}} \end{picture}$

to integrate cross section and generate distribution. the result of cross section is in the file "kfresult.dat" and the convergence behaive is in the file "kconvergence.dat". The distribution data is in the file "kplot.hbook" for pawX11 and the file "kplot.dat for directly text view. The information for event generate is in the file "kases.data".

$\begin{picture}(110,40) \put(50,20){\oval(100,30)} \put(0,20){ \begin{tabular}{l} $>$ gkevent \end{tabular}} \end{picture}$

To generate event by using the information in the file "kases.data" and the generated 4-momentum data of event is in the file "pdata1.mc".

$\begin{picture}(110,40) \put(50,20){\oval(100,30)} \put(0,20){ \begin{tabular}{l} $>$ fit \end{tabular}} \end{picture}$

To do Likelihood fitting by using 4-momentum data file "pdata1.dat" of Experiment events, phase space Moto Calo 4-momentum data file "pdata1.mc" and amplitude square Fortran source. The result of parameter fitting is in the file "pep.res". The file "pep.res" is of exactly same index as the input file "fpara.inp", but the "fpara.inp" is input and "pep.res" is the output. The file "flag.inp" is used to control the behaive of "fit".

Example Of Input Items Or Files

$resonace\_list$ is the list of resonace in this proecss, you can change the mass and width of resonaces in input data file "fort/reson.inp" without changing of fortran source.

$\begin{picture}(490,160) \put(230,75){\oval(480,150)} \put(0,75){ \begin{tabular... ...0 &0.1 &)\\ &(4 &pr4 &1.0 &0.1& )\\ &); & & & & \end{tabular}} \end{picture}$

$fit\_parameter\_list$ is the list of parameters which you want to change their value when you run program "fit" to fit with experiment data. The parameters are in the input file "fort/fpara.inp" and you can change them in the file.

$\begin{picture}(490,170) \put(230,80){\oval(480,160)} \put(0,80){ \begin{tabular... ...&(5 &f48 &0.1 &0.01 &-10 & 10 &) \\ &) & & & & &; \end{tabular}} \end{picture}$

The follow input files, generated by FDC as input data for generated Fortran source, are in the subdirectory "fort". User can change all the parameters in them to control and change parameters in the calculation.

flag.inp

$\begin{picture}(490,170) \put(230,80){\oval(480,160)} \put(0,80){ \begin{tabular... ...er of events in experiment data file ''pdata1.dat'' \end{tabular}} \end{picture}$

Where the detail of "do_mc" option is described in FDC resport.

fpara.inp

$\begin{picture}(490,130) \put(230,60){\oval(480,120)} \put(0,60){ \begin{tabular... ....1 &0.01 &-10 & 10 \\ &5 &0.1 &0.01 &-10 & 10 \\ \end{tabular}} \end{picture}$

reson.inp

$\begin{picture}(490,130) \put(230,60){\oval(480,120)} \put(0,60){ \begin{tabular... ...2 &1.0 &0.1 &\\ &3 &1.0 &0.1 &\\ &4 &1.0 &0.1 & \end{tabular}} \end{picture}$

Next: To draw Feynman diagrams Up: Physical Process Previous: To use structure function Contents

wang jian xiong 2003-01-10