EvtSLPoleFF.cpp

Go to the documentation of this file.

/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen 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.                                  *
*                                                                      *
* EvtGen 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 EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/

#include "EvtGenModels/EvtSLPoleFF.hh"

#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenBase/EvtReport.hh"

#include <math.h>
#include <stdlib.h>
#include <string>

EvtSLPoleFF::EvtSLPoleFF( int numarg, double* arglist )
{
    //arg - maybe ignore the last argument - if odd ... Sigh
    numSLPoleargs = numarg - ( numarg % 2 );
    for ( int i = 0; i < numSLPoleargs; i++ ) {
        SLPoleargs[i] = arglist[i];
    }

    return;
}

void EvtSLPoleFF::getscalarff( EvtId parent, EvtId, double t, double,
                               double* fpf, double* f0f )
{
    // Form factors have a general form, with parameters passed in
    // from the arguements.

    if ( numSLPoleargs != 8 ) {
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "Problem in EvtSLPoleFF::getscalarff\n";
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "wrong number of arguements!!!\n";
    }

    double mb = EvtPDL::getMeanMass( parent );
    double mb2 = mb * mb;

    double f0, af, bf, powf;

    f0 = SLPoleargs[0];
    af = SLPoleargs[1];
    bf = SLPoleargs[2];
    powf = SLPoleargs[3];
    *fpf = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );

    f0 = SLPoleargs[4];
    af = SLPoleargs[5];
    bf = SLPoleargs[6];
    powf = SLPoleargs[7];

    *f0f = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );

    return;
}

void EvtSLPoleFF::getvectorff( EvtId parent, EvtId, double t, double,
                               double* a1f, double* a2f, double* vf, double* a0f )
{
    if ( numSLPoleargs != 16 ) {
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "Problem in EvtSLPoleFF::getvectorff\n";
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "wrong number of arguements!!!\n";
        EvtGenReport( EVTGEN_ERROR, "EvtGen" ) << numSLPoleargs << "\n";
    }

    double mb = EvtPDL::getMeanMass( parent );
    double mb2 = mb * mb;

    double f0, af, bf, powf;

    f0 = SLPoleargs[0];
    af = SLPoleargs[1];
    bf = SLPoleargs[2];
    powf = SLPoleargs[3];
    *a1f = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );

    f0 = SLPoleargs[4];
    af = SLPoleargs[5];
    bf = SLPoleargs[6];
    powf = SLPoleargs[7];

    *a2f = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );

    f0 = SLPoleargs[8];
    af = SLPoleargs[9];
    bf = SLPoleargs[10];
    powf = SLPoleargs[11];

    *vf = f0 /
          ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                 powf ) );

    f0 = SLPoleargs[12];
    af = SLPoleargs[13];
    bf = SLPoleargs[14];
    powf = SLPoleargs[15];

    *a0f = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );
    return;
}

void EvtSLPoleFF::gettensorff( EvtId parent, EvtId, double t, double,
                               double* hf, double* kf, double* bpf, double* bmf )
{
    if ( numSLPoleargs != 16 ) {
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "Problem in EvtSLPoleFF::gettensorff\n";
        EvtGenReport( EVTGEN_ERROR, "EvtGen" )
            << "wrong number of arguements!!!\n";
    }

    double mb = EvtPDL::getMeanMass( parent );
    double mb2 = mb * mb;

    double f0, af, bf, powf;

    f0 = SLPoleargs[0];
    af = SLPoleargs[1];
    bf = SLPoleargs[2];
    powf = SLPoleargs[3];
    *hf = f0 /
          ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                 powf ) );

    f0 = SLPoleargs[4];
    af = SLPoleargs[5];
    bf = SLPoleargs[6];
    powf = SLPoleargs[7];

    *kf = f0 /
          ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                 powf ) );

    f0 = SLPoleargs[8];
    af = SLPoleargs[9];
    bf = SLPoleargs[10];
    powf = SLPoleargs[11];

    *bpf = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );

    f0 = SLPoleargs[12];
    af = SLPoleargs[13];
    bf = SLPoleargs[14];
    powf = SLPoleargs[15];

    *bmf = f0 /
           ( pow( 1.0 + ( af * t / mb2 ) + ( bf * ( ( t / mb2 ) * ( t / mb2 ) ) ),
                  powf ) );
    return;
}

void EvtSLPoleFF::getbaryonff( EvtId, EvtId, double, double, double*, double*,
                               double*, double* )
{
    EvtGenReport( EVTGEN_ERROR, "EvtGen" )
        << "Not implemented :getbaryonff in EvtSLPoleFF.\n";
    ::abort();
}

void EvtSLPoleFF::getdiracff( EvtId, EvtId, double, double, double*, double*,
                              double*, double*, double*, double* )
{
    EvtGenReport( EVTGEN_ERROR, "EvtGen" )
        << "Not implemented :getdiracff in EvtSLPoleFF.\n";
    ::abort();
}

void EvtSLPoleFF::getraritaff( EvtId, EvtId, double, double, double*, double*,
                               double*, double*, double*, double*, double*,
                               double* )
{
    EvtGenReport( EVTGEN_ERROR, "EvtGen" )
        << "Not implemented :getraritaff in EvtSLPoleFF.\n";
    ::abort();
}