NAME

SOOT::Examples::Graph - SOOT Examples for Graph

DESCRIPTION

This is a listing of all SOOT examples for Graph.

EXAMPLES

approx.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  sub DrawSmooth;
  sub _approx;
  
  my $vC1;
  my ($grxy, $grin, $grout);
  my @obj;
  
  # Macro to test interpolation function Approx()
  
  # test data (square)
  use constant n => 11;
  my $x = [1,2,3,4,5,6,6,6,8,9,10];
  my $y = [1,4,9,16,25,25,36,49,64,81,100];
  $grxy = TGraph->new(n, $x, $y);
  
  # x values, for which y values should be interpolated
  my $xout = [1.2,1.7,2.5,3.2,4.4,5.2,5.7,6.5,7.6,8.3,9.7,10.4,11.3,13];
  
  # create Canvas
  $vC1 = TCanvas->new("vC1","square",200,10,700,700);
  $vC1->Divide(2,2);
  
  # initialize graph with data
  $grin = TGraph->new(n,$x,$y);
  
  # interpolate at equidistant points (use mean for tied x-values)
  my $gs = TGraphSmooth->new("normal");
  $grout = $gs->Approx($grin,"linear");
  DrawSmooth(1,"Approx: ties = mean","X-axis","Y-axis");
  
  # re-initialize graph with data (since graph points were set to unique vales)
  $grin = TGraph->new(n,$x,$y);
  
  # interpolate at given points xout
  $grout = $gs->Approx($grin,"linear", scalar(@$xout), $xout, 0, 130);
  DrawSmooth(2,"Approx: ties = mean","","");
  
  # print output variables for given values xout
  my $vNout = $grout->GetN();
  
  #   Double_t vXout, vYout;
  #   for (Int_t k=0;k<vNout;k++) {
  #      $grout->GetPoint(k, vXout, vYout);
  #      cout << "k= " << k << "  vXout[k]= " << vXout
  #           << "  vYout[k]= " << vYout << endl;
  #   }
  
  # re-initialize graph with data
  $grin = TGraph->new(n,$x,$y);
  
  # interpolate at equidistant points (use min for tied x-values)
  #   grout = $gs->Approx(grin,"linear", 50, 0, 0, 0, 1, 0, "min");
  $grout = $gs->Approx($grin,"constant", 50, [0.0], 0.0, 0.0, 1, 0.5, "min");
  
  DrawSmooth(3,"Approx: ties = min","","");
  
  # re-initialize graph with data
  $grin = TGraph->new(n,$x,$y);
  
  # interpolate at equidistant points (use max for tied x-values)
  $grout = $gs->Approx($grin,"linear", scalar(@$xout), $xout, 0, 0, 2, 0, "max");
  DrawSmooth(4,"Approx: ties = max","","");
  
  $gApplication->Run;
  
  sub DrawSmooth {
     my ($pad, $title, $xt, $yt) = @_;
     $vC1->cd($pad);
     my $vFrame = $vC1->DrawFrame(0,0,15,150);
     $vFrame->SetTitle($title);
     $vFrame->SetTitleSize(0.2);
     $vFrame->SetXTitle($xt);
     $vFrame->SetYTitle($yt);
     $grxy->SetMarkerColor(kBlue);
     $grxy->SetMarkerStyle(21);
     $grxy->SetMarkerSize(0.5);
     $grxy->Draw("P");
     $grin->SetMarkerColor(kRed);
     $grin->SetMarkerStyle(5);
     $grin->SetMarkerSize(0.7);
     $grin->Draw("P");
     $grout->DrawClone("LP");
     push @obj, $vFrame; # hack to keep it alive until I fix the object ownership issue
  }
  

bent.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  $gROOT->Reset();
  
  my $n = 10;
  my $x    = [-0.22, 0.05, 0.25, 0.35, 0.5, 0.61,0.7,0.85,0.89,0.95];
  my $y    = [1.,2.9,5.6,7.4,9,9.6,8.7,6.3,4.5,1];
  my $exl  = [.05,.1,.07,.07,.04,.05,.06,.07,.08,.05];
  my $eyl  = [.8,.7,.6,.5,.4,.4,.5,.6,.7,.8];
  my $exh  = [.02,.08,.05,.05,.03,.03,.04,.05,.06,.03];
  my $eyh  = [.6,.5,.4,.3,.2,.2,.3,.4,.5,.6];
  my $exld = [.0,.0,.0,.0,.0,.0,.0,.0,.0,.0];
  my $eyld = [.0,.0,.05,.0,.0,.0,.0,.0,.0,.0];
  my $exhd = [.0,.0,.0,.0,.0,.0,.0,.0,.0,.0];
  my $eyhd = [.0,.0,.0,.0,.0,.0,.0,.0,.05,.0];
  my $gr = TGraphBentErrors->new($n,$x,$y,$exl,$exh,$eyl,$eyh,$exld,$exhd,$eyld,$eyhd);
  $gr->SetTitle("TGraphBentErrors Example");
  $gr->SetMarkerColor(4);
  $gr->SetMarkerStyle(21);
  $gr->Draw("ALP");
  $gApplication->Run;
  

graphApply.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  my $npoints = 3;
  my $xaxis  = [1.,2.,3.];
  my $yaxis  = [10.,20.,30.];
  my $errorx = [0.5,0.5,0.5];
  my $errory = [5.,5.,5.];
  my $exl    = [0.5,0.5,0.5];
  my $exh    = [0.5,0.5,0.5];
  my $eyl    = [5.,5.,5.];
  my $eyh    = [5.,5.,5.];
  
  my $gr1 = TGraph->new($npoints,$xaxis,$yaxis);
  my $gr2 = TGraphErrors->new($npoints,$xaxis,$yaxis,$errorx,$errory);
  my $gr3 = TGraphAsymmErrors->new($npoints,$xaxis,$yaxis,$exl,$exh,$eyl,$eyh);
  my $ff  = TF2->new("ff","-1./y");
  
  my $c1 = TCanvas->new("c1","c1");
  $c1->Divide(2,3);
  # TGraph
  $c1->cd(1);
  $gr1->DrawClone("A*");
  $c1->cd(2);
  $gr1->Apply($ff);
  $gr1->Draw("A*");
  
  # TGraphErrors
  $c1->cd(3);
  $gr2->DrawClone("A*");
  $c1->cd(4);
  $gr2->Apply($ff);
  $gr2->Draw("A*");
  
  # TGraphAsymmErrors
  $c1->cd(5);
  $gr3->DrawClone("A*");
  $c1->cd(6);
  $gr3->Apply($ff);
  $gr3->Draw("A*");
  
  $c1->Update;
  
  $gApplication->Run;
  

gtime.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  my @xa;
  my @ya;
  
  my $c1 = TCanvas->new("c1");
  
  my $ng = 100;
  my $kNMAX = 10000;
  
  my $cursor = $kNMAX;
  my $g = TGraph->new($ng);
  $g->SetMarkerStyle(21);
  $g->SetMarkerColor(4);
  my $x = 0.0;
  
  while (1) {
    $c1->Clear();
    if ($cursor > $kNMAX-$ng) {
       foreach (0..$ng-1) {
         push @xa, $x;
         push @ya, sin($x);
         $x   += 0.1;
       }
       $g->Draw("alp");
       $cursor = 0;
    } 
    else {
      $x += 0.1;
      push @xa, $x;
      push @ya, sin($x);
      $cursor++;
      
      my @nxa;
      my @nya;
      for my $i ($cursor..$cursor+$ng-1) {
        push @nxa, $xa[$i];
        push @nya, $ya[$i];
      }
      $g->DrawGraph($ng, \@nxa, \@nya, "alp");
    }
    $c1->Update();
    $gSystem->ProcessEvents();
    $gSystem->Sleep(10);
  }
  

markerwarning.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  # This script illustrates the danger of using asymmetric symbols
  # Non-symmetric symbols should be used carefully in plotting. 
  # These two graphs show how misleading a careless use of symbols can be. 
  # The two plots represent the same data sets but because of a bad symbol choice,
  # the two plots on the top appear further apart than for the bottom example.
  $gROOT->Reset;
  $gStyle->SetOptStat(0);
  
  my $Nph = 14;
  my $np_ph  = [353.4,300.2,254.3,215.2,181.0,151.3,125.2,102.7, 83.3, 66.7, 52.5, 40.2, 30.2, 22.0];
  my $nc_ph  = [3.890,3.734,3.592,3.453,3.342,3.247,3.151,3.047,2.965,2.858,2.701,2.599,2.486,2.328];
  my $npe_ph = [10.068,9.004,8.086,7.304,6.620,6.026,5.504,5.054,4.666,4.334,4.050,3.804,3.604,3.440];
  my $nce_ph = [0.235,0.217,0.210,0.206,0.213,0.223,0.239,0.260,0.283,0.318,0.356,0.405,0.465,0.545];
  
  my $Nbr = 6;
  my $np_br  = [357.0,306.0,239.0,168.0,114.0, 73.0];
  my $nc_br  = [3.501,3.275,3.155,3.060,3.053,3.014];
  my $npe_br = [8.000,11.000,10.000,9.000,9.000,8.000];
  my $nce_br = [0.318,0.311,0.306,0.319,0.370,0.429];
  
  my $phUP = TGraphErrors->new($Nph,$np_ph,$nc_ph,$npe_ph,$nce_ph);
  my $phDN = TGraphErrors->new($Nph,$np_ph,$nc_ph,$npe_ph,$nce_ph);
  my $brUP = TGraphErrors->new($Nbr,$np_br,$nc_br,$npe_br,$nce_br);
  my $brDN = TGraphErrors->new($Nbr,$np_br,$nc_br,$npe_br,$nce_br);
  
  my $Top_margin        = 0.;
  my $Left_margin       = 0.025;
  my $Right_margin      = 0.005;
  my $maxPlotPart       = 395;
  my $Marker_Size       = 1.3;
  my $Marker_Style      = 8;
  
  my $Et_200_Min        = 0.71;
  my $Et_200_Max        = 3.80;
  my $Et_130_Min        = 1.21;
  my $Et_130_Max        = 3.29;
  
  my $Nc_200_Min        = 1.31;
  my $Nc_200_Max        = 4.30;
  my $Nc_130_Min        = 1.51;
  my $Nc_130_Max        = 3.89;
  
  
  my $canvasNc = TCanvas->new("canvasNc", "Multiplicity",630,10,600,500);
  $canvasNc->SetFillColor(10);
  $canvasNc->SetBorderSize(0);
  $canvasNc->SetLeftMargin(0.02);
  $canvasNc->SetRightMargin(0.02);
  $canvasNc->SetTopMargin(0.02);
  $canvasNc->SetBottomMargin(0.02);
  
  
  # ------------>Primitives in Nc200 pad
  my $padNcUP = TPad->new("padNcUP","200 GeV",0.07,0.60,1.,1.00);
  $padNcUP->Draw();
  $padNcUP->cd();
  $padNcUP->SetFillColor(10);
  $padNcUP->SetFrameFillColor(10);
  $padNcUP->SetBorderSize(0);
  $padNcUP->SetLeftMargin($Left_margin);
  $padNcUP->SetRightMargin($Right_margin);
  $padNcUP->SetTopMargin($Top_margin+0.005);
  $padNcUP->SetBottomMargin(0.00);
  
  my $frameNcUP = TH1F->new("frameNcUP","",100,0,$maxPlotPart);
  $frameNcUP->GetYaxis()->SetLabelOffset(0.005);
  $frameNcUP->GetYaxis()->SetLabelSize(0.12);
  $frameNcUP->SetMinimum($Nc_200_Min);
  $frameNcUP->SetMaximum($Nc_200_Max);
  $frameNcUP->SetNdivisions(505,"Y");
  $frameNcUP->SetNdivisions(505,"X");
  $frameNcUP->Draw();
  
  $brUP->SetMarkerStyle(22);
  $brUP->SetMarkerSize (2.0);
  $brUP->Draw("P");                    
  
  $phDN->SetMarkerStyle(23);
  $phDN->SetMarkerSize (2);
  $phDN->Draw("P");
  
  $padNcUP->Modified();
  $canvasNc->cd();
  
  # ------------>Primitives in Nc130 pad
  my $padNcDN = TPad->new("padNcDN","130 GeV",0.07,0.02,1.,0.60);
  $padNcDN->Draw();
  $padNcDN->cd();
  $padNcDN->SetFillColor(10);
  $padNcDN->SetFrameFillColor(10);
  $padNcDN->SetBorderSize(0);
  $padNcDN->SetLeftMargin($Left_margin);
  $padNcDN->SetRightMargin($Right_margin);
  $padNcDN->SetTopMargin($Top_margin+0.005);
  $padNcDN->SetBottomMargin(0.30);
  
  my $frameNcDN = TH1F->new("frameNcDN","",100,0,$maxPlotPart);
  $frameNcDN->GetYaxis()->SetLabelOffset(0.005);
  $frameNcDN->GetYaxis()->SetLabelSize(0.07);
  $frameNcDN->GetXaxis()->SetLabelOffset(0.005);
  $frameNcDN->GetXaxis()->SetLabelSize(0.07);
  $frameNcDN->SetMinimum($Nc_200_Min);
  $frameNcDN->SetMaximum($Nc_200_Max);
  $frameNcDN->SetNdivisions(505,"Y");
  $frameNcDN->SetNdivisions(505,"X");
  $frameNcDN->Draw();
  
  $brDN->SetMarkerStyle(23);
  $brDN->SetMarkerSize (2.0);
  $brDN->Draw("P");                    
             
  $phUP->SetMarkerStyle(22);
  $phUP->SetMarkerSize (2);
  $phUP->Draw("P");
  
  my $t1 = TLatex->new();
  $t1->SetTextFont(12); 
  $t1->SetTextSize(0.0525);
  $t1->DrawLatex(-5,0.6,"Non-symmetric symbols should be used carefully in plotting. These two graphs show how misleading");
  $t1->DrawLatex(-5,0.4,"a careless use of symbols can be. The two plots represent the same data sets but because of a bad");
  $t1->DrawLatex(-5,0.2,"symbol choice, the two plots on the top appear further apart than for the bottom example.");
  
  $padNcDN->Modified();
  $canvasNc->cd();
  
  $canvasNc->Modified();
  $canvasNc->cd();
  
  $gApplication->Run;
  

motorcycle.pl

  use strict;
  use warnings;
  use SOOT ':all';
  use File::Spec;
  
  # Macro to test scatterplot smoothers: ksmooth, lowess, supsmu
  # as described in:
  #    Modern Applied Statistics with S-Plus, 3rd Edition
  #    W.N. Venables and B.D. Ripley
  #    Chapter 9: Smooth Regression, Figure 9.1 
  #
  # Example is a set of data on 133 observations of acceleration against time
  # for a simulated motorcycle accident, taken from Silverman (1985).
  
  # data taken from R library MASS: mcycle.txt
  my $inFile = shift
               || File::Spec->catfile($ENV{ROOTSYS},
                                      qw(share doc root tutorials graphs motorcycle.dat));
  
  # read file and add to fit object
  my $x = [];
  my $y = [];
  
  my ($vX, $vY);
  my $vNData = 0;
  
  open INPUT, "$inFile" or die "Could not open $inFile"; 
  while (<INPUT>) {
    ($vX,$vY) = split;
    push @$x, $vX*1.;
    push @$y, $vY*1.;
    $vNData++;
  }
  close INPUT;
  
  my $grin = TGraph->new($vNData,$x,$y);
  
  # draw graph
  my $can = TCanvas->new("can","Smooth Regression",200,10,900,700);
  $can->Divide(2,3);
  
  # Kernel Smoother
  # create new kernel smoother and smooth data with bandwidth = 2.0
  my $gs = TGraphSmooth->new("normal");
  my $grout = $gs->SmoothKern($grin,"normal",2.0);
  DrawSmooth($can, $grin, $grout, 1, "Kernel Smoother: bandwidth = 2.0", "times", "accel");
  
  # redraw ksmooth with bandwidth = 5.0
  $grout = $gs->SmoothKern($grin,"normal",5.0);
  DrawSmooth($can,$grin,$grout,2,"Kernel Smoother: bandwidth = 5.0","","");
  
  # Lowess Smoother
  # create new lowess smoother and smooth data with fraction f = 2/3
  $grout = $gs->SmoothLowess($grin,"",0.67);
  DrawSmooth($can,$grin,$grout,3,"Lowess: f = 2/3","","");
  
  # redraw lowess with fraction f = 0.2
  $grout = $gs->SmoothLowess($grin,"",0.2);
  DrawSmooth($can,$grin,$grout,4,"Lowess: f = 0.2","","");
  
  # Super Smoother
  # create new super smoother and smooth data with default bass = 0 and span = 0
  $grout = $gs->SmoothSuper($grin,"",0,0);
  DrawSmooth($can,$grin,$grout,5,"Super Smoother: bass = 0","","");
  
  # redraw supsmu with bass = 3 (smoother curve)
  $grout = $gs->SmoothSuper($grin,"",3);
  DrawSmooth($can,$grin,$grout,6,"Super Smoother: bass = 3","","");
  
  sub DrawSmooth {
     my ($can, $grin, $grout, $pad, $title, $xt, $yt) = @_;
     $can->cd($pad);
     my $vFrame = $can->DrawFrame(0,-130,60,70);
     $vFrame->SetTitle($title);
     $vFrame->SetTitleSize(0.2);
     $vFrame->SetXTitle($xt);
     $vFrame->SetYTitle($yt);
     $grin->Draw("P");
     $grout->SetMarkerColor(kRed);
     $grout->SetMarkerStyle(21);
     $grout->SetMarkerSize(0.5);
     $grout->DrawClone("P");
     $grout->DrawClone("LPX");
     $vFrame->keep;
  }
  
  $gApplication->Run;

multigraph.pl

  use strict;
  use warnings;
  use SOOT ':all';
  
  $gROOT->Reset;
  $gStyle->SetOptFit;
  my $c1 = TCanvas->new("c1","multigraph",200,10,700,500);
  $c1->SetGrid();
  
  # draw a frame to define the range
  my $mg = TMultiGraph->new();
  
  # create first graph
  my $n1 = 10;
  my $x1  = [-0.1, 0.05, 0.25, 0.35, 0.5, 0.61,0.7,0.85,0.89,0.95];
  my $y1  = [-1.,2.9,5.6,7.4,9,9.6,8.7,6.3,4.5,1];
  my $ex1 = [.05,.1,.07,.07,.04,.05,.06,.07,.08,.05];
  my $ey1 = [.8,.7,.6,.5,.4,.4,.5,.6,.7,.8];
  my $gr1 = TGraphErrors->new($n1,$x1,$y1,$ex1,$ey1);
  $gr1->SetMarkerColor(kBlue);
  $gr1->SetMarkerStyle(21);
  $gr1->Fit("pol6","q");
  $mg->Add($gr1);
  
  # create second graph
  my $n2 = 10;
  my $x2  = [-0.28, 0.005, 0.19, 0.29, 0.45, 0.56,0.65,0.80,0.90,1.01];
  my $y2  = [2.1,3.86,7,9,10,10.55,9.64,7.26,5.42,2];
  my $ex2 = [.04,.12,.08,.06,.05,.04,.07,.06,.08,.04];
  my $ey2 = [.6,.8,.7,.4,.3,.3,.4,.5,.6,.7];
  my $gr2 = TGraphErrors->new($n2,$x2,$y2,$ex2,$ey2);
  $gr2->SetMarkerColor(kRed);
  $gr2->SetMarkerStyle(20);
  $gr2->Fit("pol5","q");
  
  $mg->Add($gr2);
  
  $mg->Draw("ap");
  
   #force drawing of canvas to generate the fit TPaveStats
  $c1->Update();
  my $stats1 = $gr1->GetListOfFunctions()->FindObject("stats")->as('TPaveStats');
  my $stats2 = $gr2->GetListOfFunctions()->FindObject("stats")->as('TPaveStats');
  $stats1->SetTextColor(kBlue); 
  $stats2->SetTextColor(kRed); 
  $stats1->SetX1NDC(0.12); $stats1->SetX2NDC(0.32); $stats1->SetY1NDC(0.75);
  $stats2->SetX1NDC(0.72); $stats2->SetX2NDC(0.92); $stats2->SetY1NDC(0.78);
  $c1->Modified();
  
  $gApplication->Run;

SEE ALSO

SOOT

AUTHOR

Steffen Mueller, <smueller@cpan.org>

COPYRIGHT AND LICENSE

Copyright (C) 2010 by Steffen Mueller

SOOT, the Perl-ROOT wrapper, is free software; you can redistribute it and/or modify it under the same terms as ROOT itself, that is, the GNU Lesser General Public License. A copy of the full license text is available from the distribution as the LICENSE file.