#include <TH1.h>
#include <TH2.h>
#include <TCanvas.h>
#include <fstream>

void minLin(){
  //Define parameters common to all analysis configurations
	const Int_t points = 6;

	gStyle->SetOptStat(kFALSE);
	gStyle->SetOptFit(0000);
	gStyle->SetPadGridX(kTRUE);
	gStyle->SetPadGridY(kTRUE);
	TGaxis::SetMaxDigits(4);

	ifstream infile1("minLin_ZK4033.txt");
	ifstream infile2("minLin_ZK5363.txt");
	ifstream infile3("minLin_ZK5365.txt");
	ifstream infile4("minLin_ZK5369.txt");
	ifstream infile5("minLin_ZK5370.txt");
	ifstream infile6("minLin_ZK5401.txt");
	ifstream infile7("minLin_ZK5407.txt");
	ifstream infile8("minLin_ZK5553.txt");
	ifstream infile9("minLin_ZK5555.txt");


	float LL[9], HV[9], I_out[9];
	Double_t my_LL[9][points], my_HV[9][points], my_Iout[9][points];
	char header[256];
	for(int h = 0; h<3; h++) {//read the header info
		infile1 >> header;
		infile2 >> header;
		infile3 >> header;
		infile4 >> header;
		infile5 >> header;
		infile6 >> header;
		infile7 >> header;
		infile8 >> header;
		infile9 >> header;
	}

	for(int i=0; i<6; i++){
		infile1 >> LL[0] >> HV[0] >> I_out[0];
          my_LL[0][i] = LL[0];
		my_HV[0][i] = HV[0];
		my_Iout[0][i] = I_out[0];
	}
	for(int i=0; i<6; i++){
		infile2 >> LL[1] >> HV[1] >> I_out[1];
          my_LL[1][i] = LL[1];
		my_HV[1][i] = HV[1];
		my_Iout[1][i] = I_out[1];
	}
	for(int i=0; i<6; i++){
		infile3 >> LL[2] >> HV[2] >> I_out[2];
          my_LL[2][i] = LL[2];
		my_HV[2][i] = HV[2];
		my_Iout[2][i] = I_out[2];
	}
	for(int i=0; i<5; i++){
		infile4 >> LL[3] >> HV[3] >> I_out[3];
          my_LL[3][i] = LL[3];
		my_HV[3][i] = HV[3];
		my_Iout[3][i] = I_out[3];
	}
	for(int i=0; i<5; i++){
		infile5 >> LL[4] >> HV[4] >> I_out[4];
          my_LL[4][i] = LL[4];
		my_HV[4][i] = HV[4];
		my_Iout[4][i] = I_out[4];
	}
	for(int i=0; i<6; i++){
		infile6 >> LL[5] >> HV[5] >> I_out[5];
          my_LL[5][i] = LL[5];
		my_HV[5][i] = HV[5];
		my_Iout[5][i] = I_out[5];
	}
	for(int i=0; i<5; i++){
		infile7 >> LL[6] >> HV[6] >> I_out[6];
          my_LL[6][i] = LL[6];
		my_HV[6][i] = HV[6];
		my_Iout[6][i] = I_out[6];
	}
	for(int i=0; i<5; i++){
		infile8 >> LL[7] >> HV[7] >> I_out[7];
          my_LL[7][i] = LL[7];
		my_HV[7][i] = HV[7];
		my_Iout[7][i] = I_out[7];
	}
	for(int i=0; i<5; i++){
		infile9 >> LL[8] >> HV[8] >> I_out[8];
          my_LL[8][i] = LL[8];
		my_HV[8][i] = HV[8];
		my_Iout[8][i] = I_out[8];
	}

		infile1.close();
		infile2.close();
		infile3.close();
		infile4.close();
		infile5.close();
		infile6.close();
		infile7.close();
		infile8.close();
		infile9.close();

  TCanvas *c1 = new TCanvas("c1","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr1;
  hr1 = new TH2F("hr1","",100,-970,-500,100,1,35);
  hr1->SetTitle("Smallest non-Linearity for a given LL (ZK4033);HV (V);I_out(%)");
  hr1->GetXaxis()->CenterTitle();
  hr1->GetYaxis()->CenterTitle();
  hr1->Draw();

  TGraph *gr1 = new TGraph(6,my_HV[0],my_Iout[0]);
  gr1->SetLineColor(1);
  gr1->SetMarkerStyle(20);
  gr1->SetMarkerColor(1);
  gr1->SetMarkerSize(0.55);
  gr1->Draw("LP");
  c1->SaveAs("ZK4033.png");

  TCanvas *c2 = new TCanvas("c2","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr2;
  hr2 = new TH2F("hr2","",100,-900,-500,100,0,37);
  hr2->SetTitle("Smallest non-Linearity for a given LL (ZK5363);HV (V);I_out(%)");
  hr2->GetXaxis()->CenterTitle();
  hr2->GetYaxis()->CenterTitle();
  hr2->Draw();

  TGraph *gr2 = new TGraph(6,my_HV[1],my_Iout[1]);
  gr2->SetLineColor(2);
  gr2->SetMarkerStyle(20);
  gr2->SetMarkerColor(2);
  gr2->SetMarkerSize(0.55);
  gr2->Draw("LP");
  c2->SaveAs("ZK5363.png");

  TCanvas *c3 = new TCanvas("c3","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr3;
  hr3 = new TH2F("hr3","",100,-720,-550,100,0,29);
  hr3->SetTitle("Smallest non-Linearity for a given LL (ZK5365);HV (V);I_out(%)");
  hr3->GetXaxis()->CenterTitle();
  hr3->GetYaxis()->CenterTitle();
  hr3->Draw();

  TGraph *gr3 = new TGraph(6,my_HV[2],my_Iout[2]);
  gr3->SetLineColor(3);
  gr3->SetMarkerStyle(20);
  gr3->SetMarkerColor(3);
  gr3->SetMarkerSize(0.55);
  gr3->Draw("LP");
  c3->SaveAs("ZK5365.png");

  TCanvas *c4 = new TCanvas("c4","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr4;
  hr4 = new TH2F("hr4","",100,-700,-530,100,0,8);
  hr4->SetTitle("Smallest non-Linearity for a given LL (ZK5369);HV (V);I_out(%)");
  hr4->GetXaxis()->CenterTitle();
  hr4->GetYaxis()->CenterTitle();
  hr4->Draw();

  TGraph *gr4 = new TGraph(5,my_HV[3],my_Iout[3]);
  gr4->SetLineColor(4);
  gr4->SetMarkerStyle(20);
  gr4->SetMarkerColor(4);
  gr4->SetMarkerSize(0.55);
  gr4->Draw("LP");
  c4->SaveAs("ZK5369.png");

  TCanvas *c5 = new TCanvas("c5","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr5;
  hr5 = new TH2F("hr5","",100,-820,-550,100,0,37);
  hr5->SetTitle("Smallest non-Linearity for a given LL (ZK5370);HV (V);I_out(%)");
  hr5->GetXaxis()->CenterTitle();
  hr5->GetYaxis()->CenterTitle();
  hr5->Draw();

  TGraph *gr5 = new TGraph(5,my_HV[4],my_Iout[4]);
  gr5->SetLineColor(6);
  gr5->SetMarkerStyle(20);
  gr5->SetMarkerColor(6);
  gr5->SetMarkerSize(0.55);
  gr5->Draw("LP");
  c5->SaveAs("ZK5370.png");

  TCanvas *c6 = new TCanvas("c6","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr6;
  hr6 = new TH2F("hr6","",100,-800,-580,100,0,36);
  hr6->SetTitle("Smallest non-Linearity for a given LL (ZK5401);HV (V);I_out(%)");
  hr6->GetXaxis()->CenterTitle();
  hr6->GetYaxis()->CenterTitle();
  hr6->Draw();

  TGraph *gr6 = new TGraph(6,my_HV[5],my_Iout[5]);
  gr6->SetLineColor(7);
  gr6->SetMarkerStyle(20);
  gr6->SetMarkerColor(7);
  gr6->SetMarkerSize(0.55);
  gr6->Draw("LP");
  c6->SaveAs("ZK5401.png");

  TCanvas *c7 = new TCanvas("c7","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr7;
  hr7 = new TH2F("hr7","",100,-740,-630,100,0,28);
  hr7->SetTitle("Smallest non-Linearity for a given LL (ZK5407);HV (V);I_out(%)");
  hr7->GetXaxis()->CenterTitle();
  hr7->GetYaxis()->CenterTitle();
  hr7->Draw();

  TGraph *gr7 = new TGraph(5,my_HV[6],my_Iout[6]);
  gr7->SetLineColor(8);
  gr7->SetMarkerStyle(20);
  gr7->SetMarkerColor(8);
  gr7->SetMarkerSize(0.55);
  gr7->Draw("LP");
  c7->SaveAs("ZK5407.png");

  TCanvas *c8 = new TCanvas("c8","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr8;
  hr8 = new TH2F("hr8","",100,-960,-530,100,0,21);
  hr8->SetTitle("Smallest non-Linearity for a given LL (ZK5553);HV (V);I_out(%)");
  hr8->GetXaxis()->CenterTitle();
  hr8->GetYaxis()->CenterTitle();
  hr8->Draw();

  TGraph *gr8 = new TGraph(5,my_HV[7],my_Iout[7]);
  gr8->SetLineColor(9);
  gr8->SetMarkerStyle(20);
  gr8->SetMarkerColor(9);
  gr8->SetMarkerSize(0.55);
  gr8->Draw("LP");
  c8->SaveAs("ZK5553.png");

  TCanvas *c9 = new TCanvas("c9","nonLinearity vs HV",0,0,900,600);
  gPad->SetLogy(0);
  gPad->SetLogx(0);
  TH2F *hr9;
  hr9 = new TH2F("hr9","",100,-900,-500,100,0,6);
  hr9->SetTitle("Smallest non-Linearity for a given LL (ZK5555);HV (V);I_out(%)");
  hr9->GetXaxis()->CenterTitle();
  hr9->GetYaxis()->CenterTitle();
  hr9->Draw();

  TGraph *gr9 = new TGraph(5,my_HV[8],my_Iout[8]);
  gr9->SetLineColor(32);
  gr9->SetMarkerStyle(20);
  gr9->SetMarkerColor(32);
  gr9->SetMarkerSize(0.55);
  gr9->Draw("LP");

  c9->SaveAs("ZK5555.png");

  return ;
	}