#ifndef ROOT_TGraphBentErrors
#define ROOT_TGraphBentErrors

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// TGraphBentErrors                                                     //
//                                                                      //
// a Graph with bent, asymmetric error bars                             //
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#ifndef ROOT_TGraph
#include "TGraph.h"
#endif

class TGraphBentErrors : public TGraph {

protected:
    Double_t    *fEXlow;        //[fNpoints] array of X low errors
    Double_t    *fEXhigh;       //[fNpoints] array of X high errors
    Double_t    *fEYlow;        //[fNpoints] array of Y low errors
    Double_t    *fEYhigh;       //[fNpoints] array of Y high errors

    Double_t    *fEXlowd;        //[fNpoints] array of X low displacements
    Double_t    *fEXhighd;       //[fNpoints] array of X high displacements
    Double_t    *fEYlowd;        //[fNpoints] array of Y low displacements
    Double_t    *fEYhighd;       //[fNpoints] array of Y high displacements

    virtual void     SwapPoints(Int_t pos1, Int_t pos2);

public:
        TGraphBentErrors();
        TGraphBentErrors(Int_t n);
        TGraphBentErrors(Int_t n, 
			 const Float_t *x, const Float_t *y, 
			 const Float_t *exl=0, const Float_t *exh=0, 
			 const Float_t *eyl=0, const Float_t *eyh=0,
			 const Float_t *exld=0, const Float_t *exhd=0, 
			 const Float_t *eyld=0, const Float_t *eyhd=0);
        TGraphBentErrors(Int_t n, 
			 const Double_t *x, const Double_t *y, 
			 const Double_t *exl=0, const Double_t *exh=0, 
			 const Double_t *eyl=0, const Double_t *eyh=0,
			 const Double_t *exld=0, const Double_t *exhd=0, 
			 const Double_t *eyld=0, const Double_t *eyhd=0);
        virtual ~TGraphBentErrors();
        virtual void    Apply(TF1 *f);
        virtual void    ComputeRange(Double_t &xmin, Double_t &ymin, 
				     Double_t &xmax, Double_t &ymax) const;
        Double_t        GetErrorX(Int_t bin) const;
        Double_t        GetErrorY(Int_t bin) const;
        Double_t       *GetEXlow()  const {return fEXlow;}
        Double_t       *GetEXhigh() const {return fEXhigh;}
        Double_t       *GetEYlow()  const {return fEYlow;}
        Double_t       *GetEYhigh() const {return fEYhigh;}
        virtual Int_t   InsertPoint(); // *MENU*
        virtual void    Paint(Option_t *chopt="");
        virtual void    Print(Option_t *chopt="") const;
        virtual Int_t   RemovePoint(); // *MENU*
        virtual Int_t   RemovePoint(Int_t ipoint);
        virtual void    SavePrimitive(ofstream &out, Option_t *option);
        virtual void    Set(Int_t n); 
        virtual void    SetPoint(Int_t i, Double_t x, Double_t y);
        virtual void    SetPointError(Double_t exl, Double_t exh, 
				      Double_t eyl, Double_t eyh); // *MENU*
        virtual void    SetPointError(Int_t i, 
				      Double_t exl, Double_t exh, 
				      Double_t eyl, Double_t eyh);

        ClassDef(TGraphBentErrors,1)  //A graph with bent, asymmetric error bars
};

#endif