//B. Pietras, GENP-USC 2014.

#include <iostream>
#include <string>
#include <stdlib.h> 
#include <fstream>
#include <sstream> 
#include <TFile.h>
#include <TH1I.h>
#include <TF1.h>
#include <TCanvas.h>
#include <TSpectrum.h>
#include <TVirtualFitter.h>
#include <TGraph.h>
#include <TStyle.h>
#include <TGraphErrors.h>
#include <TLegend.h>
#include <TMath.h>
#include <TFile.h>
#include <TTree.h>
#include <TApplication.h>

using namespace std;

// g++ -o ted ted.cpp `root-config --cflags --glibs` -lSpectrum  && ./ted 

    int main()
    {
    TApplication *dummy=new TApplication("dummy",0,0);
    const char *names[]={"350V0214","355V0213","360V0212","365V0211","370V0210","375V0209","380V0208","385V0207","390V0206","395V0205","400V0204"}; 
    string path="/home/genp/data/rootfiles/Co60_2014-06-20-2300_";
    string end=".root";
    TString inno,histi;
    TFile *m;
    TTree *T;
    char histy[40];

    int j,z,k,thresh,hith,npeaks,nfound,l;
    Float_t *xpeaks, xp1, xp2;
    double mean1,sig1,mean2,sig2;
    string oof,str,doof;
    TString sugar,sugarb,topstr;
    char X,tmpstr[1],hymn[30],funname[20];
    TH1I *h[64], *hb[64]; 
    TCanvas *can[11];
    npeaks=2;
    TF1 *gas[3];
    // DON'T change these anymore.
    gas[0]=new TF1("gas","gaus(0)+gaus(3)",5,5);

 //   cout <<"HERE1"<<endl;


    for(int i=0;i<11;i++)
	{

	inno=path;
	inno+=names[i];
	inno+=end;
	cout<<inno<<endl;
	TFile *m= new TFile(inno);
	TTree *T= (TTree*)m->Get("events");
        can[i] = new TCanvas("blah","blah",0,0,1700,950);
        can[i]->Divide(8,8);

        for(int j=0;j<64;j++)
	{
	can[i]->cd(j+1);
	sprintf(histy,"Energy[%d]>>h%d(500,20,60020)",j,j);
	histi=histy;
	T->Draw(histi);
	cout<<"drawing hist "<<j<<endl;
	}

	histi=names[i];
	histi+=".png";
	cout<<"HERE5"<<endl;
 	can[i]->SaveAs(histi);
	cout<<"HERE6"<<endl;
	can[i]->Close();
	T->Delete();
	m->Delete();
	}}