//
// GaussianResonance.C

#include <stdio.h>
#include "TGraphErrors.h"
#include "TCanvas.h"
#include "TFile.h"
#include <iostream>
#include "TFile.h"
#include "TGraph.h"
#include <fstream>

using namespace std;
void GaussianResonance(void) {

TFile *CrossSection = new TFile("CrossSection2.root", "RECREATE");
TCanvas *c2 = new TCanvas();
c2->SetLogy();

int i, j=0, k=0;
ifstream infile;
double m[1200], pt07[1200],pt10[1200],pt15[1200],res[1200];

infile.open("Gaussian_Resonance.txt");

for(i=0;i<1200;i++) { m[i]=0; res[i]=0, pt07[i]=0; pt10[i]=0; pt15[i]=0;}

 for (std::string line; std::getline(infile, line); ++k) {
    std::istringstream strstrm(line);
    strstrm >> m[k] >> res[k] >> pt07[k] >> pt10[k] >> pt15[k];
// cout<<"Mass: "<<m[k]<<" "<<c[k]<<" "<<br[k]<<" "<<a[k]<<endl;
//cout <<"k: "<<k;
}

//k--;
cout<<"Test\n";



TGraph *gr3= new TGraph(k-1,m, pt07);
gr3->SetNameTitle("gr3","Cross Section*Branching Ratio*Acceptance vs Mass jj; Mass jj (GeV);Cross Section*branching ratio*acceptance(pb)");
// gr3->SetName("gr3");

gr3->SetMaximum(10e+5);
gr3->SetMinimum(10e-10);
gr3->SetLineColor(kGreen);
gr3->SetLineWidth(4);
gr3->SetMarkerColor(kBlack);
gr3->SetMarkerStyle(21);
// gr3->Write();

TGraph *gr4= new TGraph(k-1,m,pt10);
gr4->SetNameTitle("gr4","Cross Section*Branching Ratio*Acceptance vs Mass jj; Mass jj (GeV);Cross Section*branching ratio*acceptance(pb)");
gr4->SetLineColor(kViolet+6);
gr4->SetLineWidth(4);
gr4->SetMarkerColor(kBlack);


TGraph *gr5= new TGraph(k-1,m,pt15);
gr4->SetNameTitle("gr4","Cross Section*Branching Ratio*Acceptance vs Mass jj; Mass jj (GeV);Cross Section*branching ratio*acceptance(pb)");
gr4->SetLineColor(kBlue);
gr4->SetLineWidth(4);
gr4->SetMarkerColor(kBlack);

c2->cd();
gr3->Draw("ALP");
gr4->Draw("LP");
gr5->Draw("LP");

c2->Modified();
c2->Update();

gDirectory->Append(gr3);
CrossSection->Write();

}