#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", "-log(1/(2*pi*sqrt(1))*exp( -((x)/1)**2/2))",-2.5,2.5);

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

	const double y0 = g1->Eval(0);
	TCanvas* cnv = new TCanvas;
	g1->SetTitle("");
	g1->Draw();

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

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

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


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

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

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

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