//*-- AUTHOR : B. Spruck
//*-- Modified : 26/10/2005 by B. Spruck
//*-- Modified : 15/nov/2005 by B. Spruck

#ifndef HHYP_PP_EP_EM_GAMMA_KINE_H
#define HHYP_PP_EP_EM_GAMMA_KINE_H

#include <TObject.h>
#include <TArrayI.h>
#include "hhypbasealgorithm.h"
#include "hhyplistfiller.h"

#include "hhypKineFitPPEpEmGamma.h"
#include <TNtuple.h>

// ***************************************************************************

#define	HHYPEPEM_MINUIT_SIZE	16
// 12 should be enough already... stay on the save side

class HHypPPEpEmGammaKine:public HHypBaseAlgorithm
{

  // ***************************************************************************
private:
	TNtuple       * qa;// QA ntuple (pulls, chi2, ..)
	TNtuple       *kinetest;

	bool use_err_fixed;// use fixed lin error per sector
	
	Float_t input[HHYPEPEM_MINUIT_SIZE];
	Double_t output[HHYPEPEM_MINUIT_SIZE];
	Double_t output2[HHYPEPEM_MINUIT_SIZE];
	Float_t errInput[HHYPEPEM_MINUIT_SIZE];
	Double_t errOutput[HHYPEPEM_MINUIT_SIZE];
	Float_t derOut[HHYPEPEM_MINUIT_SIZE];
	Double_t gInput[HHYPEPEM_MINUIT_SIZE];

//	Float_t mom_err_sec_0,mom_err_sec_1,mom_err_sec_2,
//			mom_err_sec_3,mom_err_sec_4,mom_err_sec_5;
	
	HHypKineFitPPEpEmGamma *fit;

	Double_t getMomErr(Int_t sec,Double_t P,Int_t pid);
	Bool_t DoTheRefit(TVector3 *momentum,Int_t idx_p1,Int_t idx_p2,Int_t idx_em,Int_t idx_ep,Float_t &chi2,Float_t &chi24,Float_t *pulls);

public:

	HHypPPEpEmGammaKine(char *name_i = NULL, Option_t par[] = NULL);
	
	// virtual 
	~ HHypPPEpEmGammaKine(void);
	
	Bool_t          execute();
	Bool_t          init();
	Bool_t          reinit();
	Bool_t          finalize();
	
	ClassDef(HHypPPEpEmGammaKine, 0) //Kinematic refit front-end alg
};

//============================================================================

#endif