IslamModel.h

/********************************************************************************

    Copyright 2013 Jan Kašpar

    This file is part of Elegent (http://elegent.hepforge.org/).

    Elegent is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    Elegent is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with Elegent.  If not, see <http://www.gnu.org/licenses/>.
 
********************************************************************************/

#ifndef _elegent_islam_model_
#define _elegent_islam_model_

#include "Model.h"
#include "Math.h"

namespace Elegent
{


class IslamModel : public Model
{
    public:
        enum VariantType
        {
            vHP,    
            vLxG,   
            vLxG13, 
        } variant;

        enum ModeType
        {
            mDiff,      
            mCore,      
            mQuark,     
            mDiffCore,  
            mFull       
        } mode;

        IslamModel();
        ~IslamModel();

        void Configure(VariantType _v, ModeType _m);

        virtual void Init();

        static TComplex CEF(double a, double b, double c);

        void SetUnitarizationOrders(int qq, int cgc)
        {
            qqMaxOrder = qq;
            cgcMaxOrder = cgc;
        }

        virtual void Print() const;

        virtual TComplex Amp(double t) const;
        virtual TComplex Prf(double b) const;
        
    protected:
        TComplex R, a, Diff_fac_profile, Diff_fac;
        
        TComplex Abs_fac;
    
        double beta, m_omega_sq, Core_fac;
        bool multipleOmegaExchange;
        
        double m0sq;
    
        double r0, omega;
        TComplex Quark_fac;
        TComplex Quark_const;
        int qqMaxOrder;
            
        double lambda, m_c;
        TComplex cgc_fac;
        int cgcMaxOrder;
    
        double precision_b, precision_t, upper_bound_b, upper_bound_t;

        bool integ_workspace_initialized;
        unsigned long integ_workspace_size_b;
        gsl_integration_workspace *integ_workspace_b;
        unsigned long integ_workspace_size_t;
        gsl_integration_workspace *integ_workspace_t;

        TComplex T_diff(double t) const;
        TComplex GammaD(double b) const;
        static TComplex GammaD_J0(double b, double *par, const void *vobj);
        
        static TComplex T_core_integ(double b, double *par, const void *vobj);
        double F_sq(double t) const;
        TComplex T_core(double t) const;

        TComplex T_quark(double t) const;
        double I_integral(double qt, double al) const;
        static double F_cal_integ(double x, double *par, const void *vobj);
        double F_cal(int n, double qt, double om, double m0sq) const;
        
        static double T_hp_integ(double b, double *par, const void *vobj);
        TComplex T_hp_n(int n, double t) const;
        TComplex T_hp(double t) const;

        static double T_lxg_integ(double b, double *par, const void *vobj);
        TComplex T_lxg_n(int n, double t) const;
        TComplex T_lxg(double t) const;

        static TComplex Amp_J0(double t, double *par, const void *vobj);
};

} // namespace

#endif