#include "box.h"

void boxexample()
{
   TCanvas *c1 = new TCanvas("c1","Sin + Cos",1200,800);
   c1->Divide(2,2);

   TRandom *rand = new TRandom();
   TH2I *h1 = new TH2I("h1","Sin",36,0,360,300,-1.5,1.5);
   TH2I *h2 = new TH2I("h2","Cos",36,0,360,300,-1.5,1.5);
   h1->GetXaxis()->SetTitle("Deg");
   float myRand;
   for (int i = 0; i < 360; i+=10) {
      for (int j = 0; j < 100; j++) {
         myRand = rand->Gaus(sin(i*3.14/180),0.2);
         h1->Fill(i,myRand);
         myRand = rand->Gaus(cos(i*3.14/180),0.1);
         h2->Fill(i,myRand);
      }
   }
   c1->cd(1);
   h1->Draw();
   c1->cd(2);
   h2->Draw();

   c1->cd(3);
   BoxPlot *b = new BoxPlot("Sin + Cos",(TH2D*)h1,(TH2D*)h2);
   b->SetLineColor(kRed,0);
   b->SetLineColor(kBlue,1);
   b->Draw();
}