/* This is the header file for functions that specify 
   the elementary pion photoproduction t matrices. */

#ifndef t_gaN_piN_h
#define t_gaN_piN_h


#include <default.h>
#include <complex.h>
#include <vector.h>
#include <t_piprod_born.h>
#include <t_piprod_delta.h>



dcomplex t_gaN_piN(int reaction, int operator, int model, 
                   int lam, double sn, double sh,
		   double ga_en, double pi_en,
                   vec k_ga, vec p_pi, vec p_nu);

double unitary_phase(double sqrt_s);


dcomplex t_gan_ppim(int operator, int model, int lam, double sn, double sh,
		    double ga_en, double pi_en,
		    vec k_ga, vec p_pi, vec p_nu);
dcomplex t_l_ppim(int operator, int model, int lam, vec p_pi, vec k_ga,
		  double po_s, vec pi_s, double E_s, 
                  double po_u, vec pi_u, double E_u);
dcomplex t_l_ppim_background(int lam, vec p_pi, vec k_ga, 
                             double po_s, double E_s, 
                             double po_u, vec pi_u, double E_u);
vec t_k_ppim(int operator, int model, int lam, 
             vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
	     double po_s, vec pi_s, double E_s, 
             double po_u, vec pi_u, double E_u);
vec t_k_ppim_background(int lam, vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
			double po_s, double E_s, 
                        double po_u, vec pi_u, double E_u);


dcomplex t_gap_npip(int operator, int model, int lam, double sn, double sh,
		    double ga_en, double pi_en,
		    vec k_ga, vec p_pi, vec p_nu);
dcomplex t_l_npip(int operator, int model, int lam, vec p_pi, vec k_ga,
		  double po_s, vec pi_s, double E_s, 
                  double po_u, vec pi_u, double E_u);
dcomplex t_l_npip_background(int lam, vec p_pi, vec k_ga, 
                             double po_s, double E_s, 
                             double po_u, vec pi_u, double E_u);
vec t_k_npip(int operator, int model, int lam, 
             vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
	     double po_s, vec pi_s, double E_s, 
             double po_u, vec pi_u, double E_u);
vec t_k_npip_background(int lam, vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
			double po_s, double E_s, 
                        double po_u, vec pi_u, double E_u);


dcomplex t_gap_ppio(int operator, int model, int lam, double sn, double sh,
		    double ga_en, double pi_en,
		    vec k_ga, vec p_pi, vec p_nu);
dcomplex t_l_ppio(int operator, int model, int lam, vec p_pi, vec k_ga,
		  double po_s, vec pi_s, double E_s, 
                  double po_u, vec pi_u, double E_u);
dcomplex t_l_ppio_background(int lam, vec p_pi, vec k_ga, 
                             double po_s, double E_s, 
                             double po_u, vec pi_u, double E_u);
vec t_k_ppio(int operator, int model, int lam, 
             vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
	     double po_s, vec pi_s, double E_s, 
             double po_u, vec pi_u, double E_u);
vec t_k_ppio_background(int lam, vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
			double po_s, double E_s, 
                        double po_u, vec pi_u, double E_u);


dcomplex t_gan_npio(int operator, int model, int lam, double sn, double sh,
		    double ga_en, double pi_en,
		    vec k_ga, vec p_pi, vec p_nu);
dcomplex t_l_npio(int operator, int model, int lam, vec p_pi, vec k_ga,
		  double po_s, vec pi_s, double E_s, 
                  double po_u, vec pi_u, double E_u);
dcomplex t_l_npio_background(int lam, vec p_pi, vec k_ga, 
                             double po_s, double E_s, 
                             double po_u, vec pi_u, double E_u);
vec t_k_npio(int operator, int model, int lam, 
             vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
	     double po_s, vec pi_s, double E_s, 
             double po_u, vec pi_u, double E_u);
vec t_k_npio_background(int lam, vec p_pi, vec k_ga, vec p_nu, double Q_sqr,
			double po_s, double E_s, 
                        double po_u, vec pi_u, double E_u);


#endif