//----------------------------------------------------------------------------------
//
// Module: EvtLb2Lll.hh
//
// Desription: Routine to implement Lambda_b0 -> Lambda_0 l+ l- decays accroding to
//             several models: Chen. Geng.
//                             Aliev. Ozpineci. Savci.
//
// Modification history:
//
//  15/09/2004  PR   Module created according to PHSP model
//  20/02/2005  PR   Added parameters, created matrix element (without polarization)
//  04/03/2005  PR   LD contrib., corrected WC eff. according to Chen. Geng.
//
//----------------------------------------------------------------------------------

#ifndef EVTLB2LLL_HH
#define EVTLB2LLL_HH

#include "EvtGenBase/EvtDecayAmp.hh"
#include "EvtGenBase/EvtTensor4C.hh"
#include "EvtGenBase/EvtDiracSpinor.hh"
#include "EvtGenModels/EvtWilsonCoeficients.hh"

class EvtLb2Lll:public EvtDecayAmp {

public:

  EvtLb2Lll() {}
  virtual ~EvtLb2Lll();

  std::string getName();
  EvtDecayBase* clone();

  void decay(EvtParticle *p);
  void init();
  void initProbMax();
  void calcAmp(EvtAmp *amp,EvtParticle *parent);

  EvtTensor4C EvtLeptonTG5Current(const EvtDiracSpinor &d,const EvtDiracSpinor &dp);

private:

  double m_polarizationLambdab0;
  double m_maxProbability;
  double m_poleSize;
  long   m_noTries;
  double m_omega;

  std::string m_decayName;
  std::string m_polarizationIntroduction;
  std::string m_HEPmodel;
  std::string m_FFtype;
  std::string m_effectContribution;
  
  EvtWilsonCoeficients m_WC;

};

#endif