//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtSLPoleFF.cc
//
// Description: Routine to implement semileptonic form factors
//              according to the model SLPoles
//
// Modification history:
//
//    DJL       April 17,1998       Module created
//
//------------------------------------------------------------------------
// 
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenModels/EvtSLPoleFF.hh"
#include <string>
#include "EvtGenBase/EvtPDL.hh"
#include <math.h>
#include <stdlib.h>

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 ) {
     report(ERROR,"EvtGen") << "Problem in EvtSLPoleFF::getscalarff\n";
     report(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 ) {
     report(ERROR,"EvtGen") << "Problem in EvtSLPoleFF::getvectorff\n";
     report(ERROR,"EvtGen") << "wrong number of arguements!!!\n";
     report(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 ) {
     report(ERROR,"EvtGen") << "Problem in EvtSLPoleFF::gettensorff\n";
     report(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*){
  
  report(ERROR,"EvtGen") << "Not implemented :getbaryonff in EvtSLPoleFF.\n";  
  ::abort();

}