#ifndef DIRE__Shower__Lorentz_FI_H
#define DIRE__Shower__Lorentz_FI_H

#include "DIRE/Shower/Lorentz.H"

namespace DIRE {

  class Lorentz_FI: public Lorentz {
  protected:

    double Jacobian(const Splitting &s) const;

    double PDFEstimate(const Splitting &s) const;

  public:

    Lorentz_FI(const Kernel_Key &k);

    int Construct(Splitting &s,const int mode) const;

    bool Cluster(Splitting &s,const int mode) const;

  };// end of class Lorentz_FI

  class Lorentz_FI_123: public Lorentz_FI {
  protected:

    double Jacobian(const Splitting &s) const;

  public:

    Lorentz_FI_123(const Kernel_Key &k);

    int Construct(Splitting &s,const int mode) const;

    bool Cluster(Splitting &s,const int mode) const;

    void SetMS(ATOOLS::Mass_Selector *const ms);

    bool Allowed(const Splitting &s) const;

  };// end of class Lorentz_FI_123

}// end of namepsace DIRE

#endif