Fast Evolution of Parton Distributions: Gauss

00001 
00008 /*
00009  *  Copyright (C) 2002 Stefan Weinzierl
00010  *
00011  *  This program is free software; you can redistribute it and/or modify
00012  *  it under the terms of the GNU General Public License as published by
00013  *  the Free Software Foundation; either version 2 of the License, or
00014  *  (at your option) any later version.
00015  *
00016  *  This program is distributed in the hope that it will be useful,
00017  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  *  GNU General Public License for more details.
00020  *
00021  *  You should have received a copy of the GNU General Public License
00022  *  along with this program; if not, write to the Free Software
00023  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
00024  */
00025 
00026 #ifndef __PDF_GAUSS_LAGUERRE_H__
00027 #define __PDF_GAUSS_LAGUERRE_H__
00028 
00029 #include "matrix_d.h"
00030 
00031 #include "partondistribution.h"
00032 
00033 namespace pdf {
00034 
00035   double evolve_nonsinglet(const partondistribution & f, double Q, double x, int number_of_points = 5);
00036   double evolve_nonsinglet_3(const partondistribution & f, double Q, double x);
00037   double evolve_nonsinglet_5(const partondistribution & f, double Q, double x);
00038   double evolve_nonsinglet_10(const partondistribution & f, double Q, double x);
00039   double evolve_nonsinglet_20(const partondistribution & f, double Q, double x);
00040   double evolve_nonsinglet_30(const partondistribution & f, double Q, double x);
00041 
00042   vector_d evolve_singlet(const partondistribution & s, const partondistribution & g, double Q, double x, int number_of_points = 5);
00043   vector_d evolve_singlet_3(const partondistribution & s, const partondistribution & g, double Q, double x);
00044   vector_d evolve_singlet_5(const partondistribution & s, const partondistribution & g, double Q, double x);
00045   vector_d evolve_singlet_10(const partondistribution & s, const partondistribution & g, double Q, double x);
00046   vector_d evolve_singlet_20(const partondistribution & s, const partondistribution & g, double Q, double x);
00047   vector_d evolve_singlet_30(const partondistribution & s, const partondistribution & g, double Q, double x);
00048 
00049   vector_d evolve_singlet_with_qed(const partondistribution & Delta, 
00050                                    const partondistribution & s, const partondistribution & g, 
00051                                    const partondistribution & phot, const partondistribution & lept, 
00052                                    double Q, double x, int number_of_points = 5);
00053   vector_d evolve_singlet_with_qed_3(const partondistribution & Delta, 
00054                                      const partondistribution & s, const partondistribution & g, 
00055                                      const partondistribution & phot, const partondistribution & lept, 
00056                                      double Q, double x);
00057   vector_d evolve_singlet_with_qed_5(const partondistribution & Delta, 
00058                                      const partondistribution & s, const partondistribution & g, 
00059                                      const partondistribution & phot, const partondistribution & lept, 
00060                                      double Q, double x);
00061   vector_d evolve_singlet_with_qed_10(const partondistribution & Delta, 
00062                                       const partondistribution & s, const partondistribution & g, 
00063                                       const partondistribution & phot, const partondistribution & lept, 
00064                                       double Q, double x);
00065   vector_d evolve_singlet_with_qed_20(const partondistribution & Delta, 
00066                                       const partondistribution & s, const partondistribution & g, 
00067                                       const partondistribution & phot, const partondistribution & lept, 
00068                                       double Q, double x);
00069   vector_d evolve_singlet_with_qed_30(const partondistribution & Delta, 
00070                                       const partondistribution & s, const partondistribution & g, 
00071                                       const partondistribution & phot, const partondistribution & lept, 
00072                                       double Q, double x);
00073 
00074 } // namespace pdf
00075 
00076 #endif // ndef __PDF_GAUSS_LAGUERRE_H__
00077
Gauss_Laguerre.h
