BSWModel.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_bsw_model_
#define _elegent_bsw_model_

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

#include <vector>

//#define DEBUG

namespace Elegent
{

class BSWModel : public Model 
{
    public:
        struct Trajectory
        {
            double C, b, a, ap, signature;

            double sign;

            void Init(double _C, double _b, double _a, double _ap, double _signature)
            {
                C=_C; b=_b; a=_a; ap=_ap; signature = _signature;
                sign = 0.;
            }
        };

        enum ModeType
        {
            mPomReg,    
            mPom,       
            mReg        
        } mode;

        BSWModel();
        
        ~BSWModel();

        void Configure(ModeType _mode = mPomReg, bool _presampled = true);
        
        virtual void Init();

        virtual void Print() const;

        virtual TComplex Amp(double t) const;   
        virtual TComplex Prf(double b) const;   

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

        double c, cp, a, f, m1, m2, asq, m1sq, m2sq;

        Trajectory A2, rho, omega;

        TComplex regge_fac;     
        signed int k_u;         
        signed int k_lnu;       

        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 Ft(double t) const;

        TComplex Rt(Trajectory tr, double t) const;

        TComplex R0t(double t) const;

        TComplex S0(double t) const;
        
        TComplex S00;

        TComplex Omega0t(double t) const;

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

        TComplex Omega0b(double b) const;

        TComplex prf0(double b) const;

        static TComplex prf0_J0(double b, double *par, const void *vobj);

        double data_db;

        unsigned int data_N;

        std::vector<double> data_re, data_im;

#ifdef DEBUG
        std::vector<double> data_b;
#endif

        void BuildSample(unsigned int samples);

        TComplex SampleEval(double b) const;
};

} // namespace

#endif