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

#ifndef gamma_N_d13_vertex_h
#define gamma_N_d13_vertex_h

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


#define G1p_d13 5.33     /* 5.33 */    
#define G2p_d13 5.58     /* 5.58 */
#define G3p_d13 0.00     /* 0.0 */

#define G1n_d13 -1.508   /* -1.508 */ 
#define G2n_d13 -0.145   /* -0.145 */  
#define G3n_d13 0.00     /* 0.0 */

#define G1_d13_off 0.0  
#define G2_d13_off 0.0  
#define G3_d13_off 0.0  


spinor_mat gamma_N_d13_vertex(int iso, int mu, int beta, fvec q4, fvec p4_av);
spinor_mat gamma_d13_N_vertex(int iso, int mu, int beta, fvec q4, fvec p4_av);


#endif