// B.Pietras, GENP USC 2013

#include <iostream>
#include <string>
#include <fstream>
#include <sstream> 
#include <iomanip>
#include <math.h>
#include <cstdlib>
#include <TFile.h>
#include <TTree.h>
#include <TH2F.h>
#include <TAxis.h>
#include <TAttAxis.h>
#include <TCanvas.h>
#include <TSystem.h>
#include <TLatex.h>
#include <TStyle.h>

using namespace std;

/*
worx
g++ `root-config --cflags` -g -O2 paparazzi2.cpp -o paparazzi2 `root-config --libs`

less bitchin
g++ `root-config --cflags` -Wl,--no-as-needed -g -O2 paparazzi2.cpp -o paparazzi2 `root-config --libs`
*/

int main()
{ 
//gStyle->SetOptStat(0);
//gStyle->SetLabelSize(0.06, "XY"); 
TFile *d;
TTree *T;
TCanvas *can[29];
TH2F *h[32];
ifstream cct("input.txt");
int i,j,ii,nrg[29],run[29];
TLatex *t = new TLatex();  
int topdown[32]={30,28,26,24,22,20,18,16,31,29,27,25,23,21,19,17,14,12,10,8,6,4,2,0,15,13,11,9,7,5,3,1};

      if (!cct) 
           {
           cout <<" Could not open file!" << endl;
           exit(-1);
           }
          else { 

  while (!cct.eof()) 
  //    while (1)
    {
  cct >> nrg[ii] >> run[ii];
      ii++;
    }
  cct.close();
	}

  char filename[70], pic[70], bo[70], lo[70];

  TFile *a = new TFile("topdown.root","RECREATE");

for(i=0;i<29;i++)
{
sprintf(filename,"/media/cebolla/rootfiles/Krakow/Run%ical.root",run[i]);
d= new TFile(filename);
gSystem->Load("libTree.so");
T =(TTree*)d->Get("TreeHouse");
sprintf(lo,"can[%i]",i);
can[i]= new TCanvas(lo,lo,300,10,1680,960);
can[i]->Divide(8,4);
for(j=0;j<32;j++)
{
sprintf(bo,"cal[%i]:cal[%i]>>h[%i](240,0,120,240,0,120)",topdown[j],topdown[j+1],j);
can[i]->cd(j+1);
gPad->SetLogz();T->Draw(bo,"","colz");
//h[j]->SetMaximum(10000);
//h[j]->SetMinimum(1000);
//h[j]->GetZaxis()->SetRangeUser(1000,10000);
//h[j]->GetXaxis()->SetLabelSize(0.06);

if(j==31){
    can[i]->cd(j+1);
    TPad* pad1 = new TPad("pad1","The pad",0.0,0.5,1.0,0.9,21);
    pad1->Draw();
    t->SetNDC();
    t->SetTextAlign(22);
    t->SetTextFont(63);
    t->SetTextSizePixels(22);
    t->DrawLatex(0.5,0.8,"Proton energy");
    sprintf(bo,"%i MeV",nrg[i]);
    t->DrawLatex(0.5,0.7,bo);
    sprintf(bo,"Run %i",run[i]);
    t->DrawLatex(0.5,0.6,bo);
        }
}
sprintf(pic,"./pics/DE_DE_Energy_%i_run_%i.png",nrg[i],run[i]);
can[i]->SaveAs(pic);
a->cd();
can[i]->Write();
}
a->Write();
a->Close();
}