#include <stdio.h>
#include "TH2F.h"
#include "TCanvas.h"
#include "TRint.h"
#include "TROOT.h"
#include "TRandom.h"
#include "TF1.h"
#include "TVirtualPad.h"
#include "TVirtualX.h"


void DynamicExec()
{
   // Example of function called when a mouse event occurs in a pad.
   // When moving the mouse in the canvas, a second canvas shows the
   // projection along X of the bin corresponding to the Y position
   // of the mouse. The resulting histogram is fitted with a gaussian.
   // A "dynamic" line shows the current bin position in Y.
   // This more elaborated example can be used as a starting point
   // to develop more powerful interactive applications exploiting Cling
   // as a development engine.
 
   static int pyold = 0;

	 if( gPad->GetEvent() != 11 ) return;    // Ignore everything except mouse button L release

	 float uxmin = gPad->GetUxmin();
   float uxmax = gPad->GetUxmax();
   int pxmin = gPad->XtoAbsPixel(uxmin);
   int pxmax = gPad->XtoAbsPixel(uxmax);
   int px = gPad->GetEventX();
   int py = gPad->GetEventY();
	 printf( "px, py = ( %d, %d )      uxmin, usmax = ( %f, %f )    ( %d , %d )\n", px, py, uxmin, uxmax, pxmin, pxmax );

	 TObject *select = gPad->GetSelected();
 
   gPad->GetCanvas()->FeedbackMode(kTRUE);
   if (pyold) {
      // erase line at old position
      gVirtualX->DrawLine(pxmin, pyold, pxmax, pyold);
      pyold = 0;
   }
 
   if(!select || !select->InheritsFrom(TH2::Class()))
      return;
 
   TH2 *h = (TH2*)select;
 
   // draw a line at current position
   gVirtualX->DrawLine(pxmin, py, pxmax, py);
   pyold = py;
 
   Float_t upy = gPad->AbsPixeltoY(py);
   Float_t y = gPad->PadtoY(upy);
 
   //create or set the new canvas c2
   TVirtualPad *padsav = gPad;
   TCanvas *c2 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("c2");
   if(c2) delete c2->GetPrimitive("Projection");
     else c2 = new TCanvas("c2","Projection Canvas",710,10,700,500);
   c2->SetGrid();
   c2->cd();
 
   //draw slice corresponding to mouse position
   Int_t biny = h->GetYaxis()->FindBin(y);
   TH1D *hp = h->ProjectionX("",biny,biny);
   hp->SetFillColor(38);
   hp->SetName("Projection");
   hp->SetTitle(TString::Format("Projection of biny=%d",biny));
   hp->Fit("gaus","ql");
   hp->GetFunction("gaus")->SetLineColor(kRed);
   hp->GetFunction("gaus")->SetLineWidth(6);
   c2->Update();
   padsav->cd();
}


#ifdef __CLING__

void Dynamic()
{

#else

int main( int argc, char *argv[] )
{

	TRint *app = new TRint( "Dynamic", &argc, argv);

#endif
	
   // Create a new canvas.
   TCanvas* c1 = new TCanvas("c1","Dynamic Slice Example", 10, 10, 700, 500);
 
   //create a 2-d histogram, fill and draw it
   TH2F *hpxpy  = new TH2F("hpxpy", "py vs px", 40,-4,4, 40,-4,4);
   hpxpy->SetStats(0);
   Double_t px,py;
   for (Int_t i = 0; i < 50000; i++) {
      gRandom->Rannor(px,py);
      hpxpy->Fill(px,py);
   }
   hpxpy->Draw("col");
   //Add a TExec object to the canvas
   c1->AddExec("dynamic", "DynamicExec()");

#ifndef __CLING__
	 
	 app->Run( kTRUE );

#endif
}