GodizovModel.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_godizov_model_
#define _elegent_godizov_model_

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

namespace Elegent
{

class GodizovModel : public Model
{
    public:
        GodizovModel();
        ~GodizovModel();
        
        void Configure(bool _presampled = true);

        virtual void Init();

        virtual void Print() const;

        virtual TComplex Amp(double t) const;

        virtual TComplex Prf(double b) const;

    protected:
        double De;      
        double ta_a;
        double Ga_P0;
        double ta_g;
        double s0;

        bool presampled;
        
        virtual void ForcePresampling(bool value)
        {
            presampled = value;
        }

        TComplex delta_t(double t) const;

        static TComplex delta_t_J0(double t, double *par, const void *obj);
        
        TComplex delta_b(double b) const;

        TComplex prf0(double b) const;

        static TComplex prf_J0(double b, double *par, const void *obj);

        double upper_bound_t, precision_t;
        double upper_bound_b, precision_b;

        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;

        double prf0_sample_db;

        unsigned int prf0_sample_N;

        std::vector<double> prf0_sample_re, prf0_sample_im;

        void Prf0SampleBuild(unsigned int samples);

        TComplex Prf0SampleEval(double b) const;
};

} // namespace

#endif