// @(#)root/hist:$Name:  $:$Id: TF12.cxx,v 1.4 2005/04/29 20:34:51 brun Exp $
// Author: Rene Brun   05/04/2003

/*************************************************************************
 * Copyright (C) 1995-2003, Rene Brun and Fons Rademakers.               *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/

#include "TF12.h"
#include "TH1.h"
#include "TVirtualPad.h"

ClassImp(TF12)

//______________________________________________________________________________
//
// a projection of a TF2 along X or Y
// has the same behaviour as a TF1
//
//      Example of a function
//
//   TF2 *f2 = new TF2("f2","sin(x)*sin(y)/(x*y)",0,5,0,5);
//   TF12 *f12 = new TF12("f12",f2,0.1,"y");
//   f12->Draw();

//______________________________________________________________________________
 TF12::TF12(): TF1()
{
// TF12 default constructor

   fCase = 0;
   fF2   = 0;
   fXY   = 0;
}


//______________________________________________________________________________
 TF12::TF12(const char *name, TF2 *f2, Double_t xy, Option_t *option)
      :TF1(name,"x",0,0)
{
// TF12 normal constructor
//   Create a TF12 (special TF1) from a projection of a TF2
//   for a fix value of Y if option="X" or X if option="Y"
// This value may be changed at any time via TF12::SetXY(xy)
   
   SetName(name);
   fF2 = f2;
   TString opt=option;
   opt.ToLower();
   if (!f2) {
      Error("TF12","Pointer to TF2 is null");
      return;
   }
   SetXY(xy);
   if (opt.Contains("y")) {
      fXmin = f2->GetYmin();
      fXmax = f2->GetYmax();
      fCase = 1;
   } else {
      fXmin = f2->GetXmin();
      fXmax = f2->GetXmax();
      fCase = 0;
   }
}

//______________________________________________________________________________
 TF12::~TF12()
{
//*-*-*-*-*-*-*-*-*-*-*F2 default destructor*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
//*-*                  =====================

}

//______________________________________________________________________________
 TF12::TF12(const TF12 &f12) : TF1(f12)
{
   ((TF12&)f12).Copy(*this);
}

//______________________________________________________________________________
 void TF12::Copy(TObject &obj) const
{
//*-*-*-*-*-*-*-*-*-*-*Copy this F2 to a new F2*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
//*-*                  ========================

   TF1::Copy(obj);
   ((TF12&)obj).fXY      = fXY;
   ((TF12&)obj).fCase    = fCase;
   ((TF12&)obj).fF2      = fF2;
}

//______________________________________________________________________________
 TF1 *TF12::DrawCopy(Option_t *option) const
{
//*-*-*-*-*-*-*-*Draw a copy of this function with its current attributes*-*-*
//*-*            ========================================================
//*-*
//*-*  This function MUST be used instead of Draw when you want to draw
//*-*  the same function with different parameters settings in the same canvas.
//*-*
//*-* Possible option values are:
//*-*   "SAME"  superimpose on top of existing picture
//*-*   "L"     connect all computed points with a straight line
//*-*   "C"     connect all computed points with a smooth curve.
//*-*
//*-* Note that the default value is "F". Therefore to draw on top
//*-* of an existing picture, specify option "SL"
//*-*
//*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

   TF12 *newf2 = new TF12();
   Copy(*newf2);
   newf2->AppendPad(option);
   newf2->SetBit(kCanDelete);
   return newf2;
}

//______________________________________________________________________________
 Double_t TF12::Eval(Double_t x, Double_t /*y*/, Double_t /*z*/, Double_t /*t*/) const
{
//*-*-*-*-*-*-*-*-*-*-*Evaluate this formula*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
//*-*                  =====================
//*-*
//*-*   Computes the value of the referenced TF2 for a fix value of X or Y
//*-*
//*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

  if (!fF2) return 0;
  if (fCase == 0) {
     return fF2->Eval(x,fXY,0);
  } else {
     return fF2->Eval(fXY,x,0);
  }
}

//______________________________________________________________________________
 Double_t TF12::EvalPar(const Double_t *x, const Double_t *params)
{
// Evaluate this function at point x[0]
//   x[0] is the value along X if fCase =0, the value along Y if fCase=1
// if params is non null, the array will be used instead of the internal TF2
// parameters
   
  if (!fF2) return 0;
  Double_t xx[2];
  if (fCase == 0) {
     xx[0] = x[0];
     xx[1] = fXY;
  } else {
     xx[0] = fXY;
     xx[1] = x[0];
  }
  fF2->InitArgs(xx,params);
  return fF2->EvalPar(xx,params);
}


//______________________________________________________________________________
 void TF12::SavePrimitive(ofstream & /*out*/, Option_t * /*option*/)
{
    // Save primitive as a C++ statement(s) on output stream out

   Error("SavePrimitive","Function not yet implemented");
}


//______________________________________________________________________________
 void TF12::SetXY(Double_t xy) 
{
   // set the value of the constant for the TF2
   //   constant in X when projecting along Y
   //   constant in Y when projecting along X
   //  The function title is set to include the value of the constant
   //  The current pad is updated
   
   fXY = xy;
   if (!fF2) return;
   if (fCase == 0) SetTitle(Form("%s (y=%g)",fF2->GetTitle(),xy));
   else            SetTitle(Form("%s (x=%g)",fF2->GetTitle(),xy));
   if (fHistogram) fHistogram->SetTitle(GetTitle());
   if (gPad) gPad->Modified();
}