#if !defined(B1PIANGAMP)
#define B1PIANGAMP

#include "IUAmpTools/Amplitude.h"
#include "IUAmpTools/AmpParameter.h"
#include "IUAmpTools/UserAmplitude.h"

#include <string>
#include <complex>
#include <vector>

#include "GPUManager/GPUCustomTypes.h"

using std::complex;
using namespace std;



#ifdef GPU_ACCELERATION
void
GPUb1piAngAmp_exec(dim3 dimGrid, dim3 dimBlock, GPU_AMP_PROTO,
		   int polBeam, GDouble polFrac,
		   int J_X, int Par_X, int L_X, int I_X, int epsilon_R, 
		   int Iz_b1, int Iz_pi,
		   GDouble u_rho_1, GDouble u_rho_3, 
		   GDouble u_omega_1, GDouble u_omega_3,
		   GDouble u_b1_0, GDouble u_b1_2, 
		   GDouble G0_omega, GDouble G0_b1, bool orthocheck);
#endif



class Kinematics;

class b1piAngAmp : public UserAmplitude< b1piAngAmp >
{

public:
  
  b1piAngAmp() : UserAmplitude< b1piAngAmp >() { }
  b1piAngAmp( const vector< string >& args );
  
  string name() const { return "b1piAngAmp"; }
  
  complex< GDouble > calcAmplitude( GDouble** pKin ) const;

  GDouble u_rho(int J_rho) const;
  GDouble u_omega(int L_omega) const;
  GDouble u_b1(int L_b1) const;

  inline complex<GDouble> BreitWigner(GDouble m0, GDouble Gamma0, int L,
			       HepLorentzVector &P1,HepLorentzVector &P2) const;
  inline GDouble CB(int j1, int j2, int m1, int m2, int J, int M) const;

  inline GDouble N(int J) const;

private:
  
  int mpolBeam;
  // GDouble mpolFrac;
  AmpParameter mpolFrac;
  int mJ_X, mPar_X, mL_X, mI_X, mepsilon_R;
  
  GDouble m_u_rho_1, m_u_rho_3;
  GDouble m_u_omega_1, m_u_omega_3;
  GDouble m_u_b1_0, m_u_b1_2, mG0_omega, mG0_b1;
  GDouble m_v_p, m_v_m;
  bool m_ORTHOCHECK, m_fastCalc;
  bool m_disableBW_omega, m_disableBW_b1;

  vector< int > mIz;

#ifdef GPU_ACCELERATION
  
  void launchGPUKernel( dim3 dimGrid, dim3 dimBlock, GPU_AMP_PROTO ) const;
  
  bool isGPUEnabled() const { return true; }
  
#endif // GPU_ACCELERATION
  

};

#endif