{
//=========Macro generated from canvas: CT1D/CT 1D canvas
//=========  (Tue Feb 24 16:29:38 2009) by ROOT version5.21/04
   TCanvas *CT1D = new TCanvas("CT1D", "CT 1D canvas",0,0,800,600);
   CT1D->SetHighLightColor(2);
   CT1D->Range(-37.5,-15.61875,37.5,140.5688);
   CT1D->SetBorderSize(2);
   
   TH1 *Test1D = new TH1F("Test1D","Test1D",128,-30,30);
   Test1D->SetBinContent(10,95);
   Test1D->SetBinContent(11,111);
   Test1D->SetBinContent(12,105);
   Test1D->SetBinContent(13,91);
   Test1D->SetBinContent(14,107);
   Test1D->SetBinContent(15,90);
   Test1D->SetBinContent(16,91);
   Test1D->SetBinContent(17,104);
   Test1D->SetBinContent(18,94);
   Test1D->SetBinContent(19,90);
   Test1D->SetBinContent(20,110);
   Test1D->SetBinContent(21,107);
   Test1D->SetBinContent(22,83);
   Test1D->SetBinContent(23,90);
   Test1D->SetBinContent(24,83);
   Test1D->SetBinContent(25,82);
   Test1D->SetBinContent(26,98);
   Test1D->SetBinContent(27,83);
   Test1D->SetBinContent(28,95);
   Test1D->SetBinContent(29,100);
   Test1D->SetBinContent(30,100);
   Test1D->SetBinContent(31,102);
   Test1D->SetBinContent(32,101);
   Test1D->SetBinContent(33,5);
   Test1D->SetBinContent(34,3);
   Test1D->SetBinContent(35,1);
   Test1D->SetBinContent(36,2);
   Test1D->SetBinContent(37,4);
   Test1D->SetBinContent(38,1);
   Test1D->SetBinContent(39,1);
   Test1D->SetBinContent(40,3);
   Test1D->SetBinContent(41,1);
   Test1D->SetBinContent(42,3);
   Test1D->SetBinContent(43,2);
   Test1D->SetBinContent(44,1);
   Test1D->SetBinContent(45,2);
   Test1D->SetBinContent(46,3);
   Test1D->SetBinContent(47,2);
   Test1D->SetBinContent(48,2);
   Test1D->SetBinContent(49,8);
   Test1D->SetBinContent(50,2);
   Test1D->SetBinContent(52,6);
   Test1D->SetBinContent(53,2);
   Test1D->SetBinContent(54,2);
   Test1D->SetBinContent(55,2);
   Test1D->SetBinContent(56,6);
   Test1D->SetBinContent(57,3);
   Test1D->SetBinContent(59,2);
   Test1D->SetBinContent(60,3);
   Test1D->SetBinContent(61,4);
   Test1D->SetBinContent(62,1);
   Test1D->SetBinContent(63,2);
   Test1D->SetBinContent(64,1);
   Test1D->SetBinContent(65,3);
   Test1D->SetBinContent(66,2);
   Test1D->SetBinContent(67,3);
   Test1D->SetBinContent(68,1);
   Test1D->SetBinContent(70,2);
   Test1D->SetBinContent(71,1);
   Test1D->SetBinContent(72,3);
   Test1D->SetBinContent(73,1);
   Test1D->SetBinContent(74,2);
   Test1D->SetBinContent(75,2);
   Test1D->SetBinContent(77,3);
   Test1D->SetBinContent(78,3);
   Test1D->SetBinContent(79,1);
   Test1D->SetBinContent(82,4);
   Test1D->SetBinContent(83,2);
   Test1D->SetBinContent(84,3);
   Test1D->SetBinContent(86,1);
   Test1D->SetBinContent(87,1);
   Test1D->SetBinContent(88,2);
   Test1D->SetBinContent(89,1);
   Test1D->SetBinContent(90,1);
   Test1D->SetBinContent(91,1);
   Test1D->SetBinContent(92,2);
   Test1D->SetBinContent(93,2);
   Test1D->SetBinContent(94,2);
   Test1D->SetBinContent(96,1);
   Test1D->SetBinContent(97,98);
   Test1D->SetBinContent(98,90);
   Test1D->SetBinContent(99,98);
   Test1D->SetBinContent(100,85);
   Test1D->SetBinContent(101,87);
   Test1D->SetBinContent(102,90);
   Test1D->SetBinContent(103,90);
   Test1D->SetBinContent(104,92);
   Test1D->SetBinContent(105,101);
   Test1D->SetBinContent(106,81);
   Test1D->SetBinContent(107,108);
   Test1D->SetBinContent(108,86);
   Test1D->SetBinContent(109,98);
   Test1D->SetBinContent(110,89);
   Test1D->SetBinContent(111,92);
   Test1D->SetBinContent(112,94);
   Test1D->SetBinContent(113,119);
   Test1D->SetBinContent(114,108);
   Test1D->SetBinContent(115,79);
   Test1D->SetBinContent(116,99);
   Test1D->SetBinContent(117,102);
   Test1D->SetBinContent(118,104);
   Test1D->SetBinContent(119,99);
   Test1D->SetEntries(4531);
   
   TPaveStats *ptstats = new TPaveStats(0.78,0.835,0.98,0.995,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text = ptstats->AddText("Test1D");
   text->SetTextSize(0.0368);
   text = ptstats->AddText("Entries = 4531   ");
   text = ptstats->AddText("Mean  = -0.1397");
   text = ptstats->AddText("RMS   =  20.46");
   ptstats->SetOptStat(1111);
   ptstats->SetOptFit(0);
   ptstats->Draw();
   Test1D->GetListOfFunctions()->Add(ptstats);
   ptstats->SetParent(Test1D->GetListOfFunctions());
   Test1D->Draw("CT1D");
   
   TPaveText *pt = new TPaveText(0.01,0.945,0.1392965,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(2);
   pt->SetFillColor(19);
   text = pt->AddText("Test1D");
   pt->Draw();
   CT1D->Modified();
   CT1D->cd();
   CT1D->SetSelected(CT1D);
   TH1 *hm =0;
   TVirtualFFT::SetTransform(0);
   hm = Test1D->FFT(hm, "MAG");
   TH1 *hp = 0;
   hp = Test1D->FFT(hp, "PH");
   hp->SetTitle("Phase of the 1st transform");
   int n=128;
   int holly=128;
   Double_t *re_full = new Double_t[n];
   Double_t *im_full = new Double_t[n];
   TVirtualFFT *fft = TVirtualFFT::GetCurrentTransform();
   fft->GetPointsComplex(re_full,im_full);
   
   //Now let's make a backward transform:
   TVirtualFFT *fft_back = TVirtualFFT::FFT(1, &holly, "C2R M K");
   fft_back->SetPointsComplex(re_full,im_full);
   fft_back->Transform();
   TH1 *hb = 0;
   //Let's look at the output
   hb = TH1::TransformHisto(fft_back,hb,"Re");
   hb->SetTitle("The backward transform result");
   hb->Draw();
}