#include "rootlogon.h"

#include <TF1.h>
#include <TH1F.h>
#include <TCanvas.h>
#include <TRandom.h>
#include <TGaxis.h>

int main()
{
	rootlogon();

	TF1* g1 = new TF1("gaus1", "0.35/x^2-2.5/x+exp(x-2.5)", 0, 2);

	g1->GetXaxis()->SetTitle("a");
	g1->GetYaxis()->SetTitle("F#left(a#right)");
	g1->GetYaxis()->SetTitleOffset(0.65);

	const double x0 = 0.27;
	const double y0 = g1->Eval(x0);
	TCanvas* cnv = new TCanvas;
	g1->GetYaxis()->SetRangeUser(-5,0);
	g1->SetTitle("");
	g1->Draw();

	TLine* line = new TLine(0.21, -5, 0.21, y0+0.5);
	line->SetLineColor(kRed);
	line->Draw("same");

	line = new TLine(0.42, -5, 0.42, y0+0.5);
	line->SetLineColor(kRed);
	line->Draw("same");

	line = new TLine(0, y0+0.5, 2, y0+0.5);
	line->SetLineColor(kRed);
	line->Draw("same");


	line = new TLine(0.17, -5, 0.17, y0+2.0);
	line->SetLineColor(kBlue);
	line->Draw("same");

	line = new TLine(0.818, -5, 0.818, y0+2.0);
	line->SetLineColor(kBlue);
	line->Draw("same");

	line = new TLine(0, y0+2.0, 2, y0+2.0);
	line->SetLineColor(kBlue);
	line->Draw("same");

	cnv->SetMargin(0.15, 0.9, 0.20, 0.9);
	cnv->SaveAs("out/nonsym.eps");
}