/* This file contains the functions that specify the
   gamma_N_p33 electromagnetic vertex. */

#ifndef gamma_N_p33_vertex_h
#define gamma_N_p33_vertex_h

#include <default.h>
#include <complex.h>
#include <vector.h>
#include <four_vector.h>
#include <t_piprod_parameters.h>

#define G1_p33 4.745    /* 4.56 */ /* 4.661  */
#define G2_p33 5.70    /* 5.49 */ /* 7.178 */
#define G3_p33 0.00     /* 0.0 */

#define G1_p33_off 0.0  /* -2.89 */  /* 0.0 */
#define G2_p33_off 0.0  /* 0.03 */   /* 0.0 */
#define G3_p33_off 0.0   /* 0.0 */

#define GM_p33 2.80   /* 2.80 */
#define GE_p33 0.04   /* 0.04 */  

spinor_mat gamma_N_p33_vertex(int type,int mu,int beta, fvec q4, fvec p4_av);
spinor_mat gamma_p33_N_vertex(int type,int mu,int beta, fvec q4, fvec p4_av);

double form_factor_G1_p33(double Q_sqr);
double form_factor_G2_p33(double Q_sqr);
double form_factor_G3_p33(double Q_sqr);

double form_factor_GM_p33(double Q_sqr);
double form_factor_GE_p33(double Q_sqr);


#endif