/* This header file contains the prototypes of the functions  
   that are needed to describe the delta contribution to the 
   pion photoproduction t matrix. */
   
#ifndef t_piprod_delta_rel_h
#define t_piprod_delta_rel_h

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


double deltawidth(double sqrt_s);
double crossed_del_formf(double mom);

dcomplex t_l_ppim_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_s, vec pi_s);
dcomplex t_l_ppim_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_u, vec pi_u);

vec t_k_ppim_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_s, vec pi_s);
vec t_k_ppim_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_u, vec pi_u);


dcomplex t_l_npip_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_s, vec pi_s);
dcomplex t_l_npip_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_u, vec pi_u);

vec t_k_npip_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_s, vec pi_s);
vec t_k_npip_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_u, vec pi_u);


dcomplex t_l_ppio_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_s, vec pi_s);
dcomplex t_l_ppio_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_u, vec pi_u);

vec t_k_ppio_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_s, vec pi_s);
vec t_k_ppio_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_u, vec pi_u);


dcomplex t_l_npio_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_s, vec pi_s);
dcomplex t_l_npio_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                                double po_u, vec pi_u);

vec t_k_npio_sdelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_s, vec pi_s);
vec t_k_npio_udelta_M1_rel(int lam, vec p_pi, vec k_ga, 
                           double po_u, vec pi_u);


#endif