#include "fiTQun.h"

class PDK_MuGamma: public fiTQun {

 private:
  bool same3Vertex;
  static PDK_MuGamma* static_PDK_MuGamma;
  
 public:
  PDK_MuGamma();
  ~PDK_MuGamma();
  double DoFit(int ievent_, int itwnd_, int seedtype, int fittype);
  bool GetSame3Vertex() {return same3Vertex;}
  
 protected:
  double RunMinimizer(double *X, int& PCflg, int *flgparfix, double *arStps=NULL);
  double SetPDK_MuGammaSeed(double* fitpars, int& PCflg);
  static void   PDK_MuGammaWrapper(int& nDim, double* gout, double& result, double par[], int flg);
  void   SaveTestPDK_MuGammaFit(double* fitparams, int istage, double totmu, double nll, int PCflg); 
  double PDK_MuGammaGetTwoRngnglogL(double *X, int& PCflg);
  void   GetTruePDK_MuGammaParams(double* tParams);
  double MuGammaNglogL(double *X, int& PCflg);
};