#include "TRandom.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TMath.h"
#include <string>
using namespace std;
using namespace TMath;

void areacerchio1(){
	int n_iter=1000;
	double media_raggio=5.;
	double sigma_raggio=1.;
	double raggio_medio=0.;
	double circ_media=0.;
	double area_media=0.;
	double varianzar=0.;
	TRandom *rng = new TRandom(123456789);
	TH1D *hr = new TH1D("hraggio","Il mio raggio",50,-5.,15.);
	TH1D *ha = new TH1D("harea","La mia area",50,0.,250.);
	TH1D *hc = new TH1D("hcirconferenza","la mia circonferenza",50,0.,60.);
	for(int i=0;i<n_iter;i++){
		double raggio=rng->Gaus(media_raggio,sigma_raggio);
		double circ=2*Pi()*raggio;
		double area=Pi()*raggio*raggio;
		varianzar=varianzar+raggio*raggio;
		cout<<"Raggio generato all'iterazione "<<i<<" "<< raggio<<endl;
		raggio_medio=raggio_medio+raggio;
		circ_media=circ_media + circ;
		area_media=area_media + area;
		hr->Fill(raggio);
		hc->Fill(2*Pi()*raggio);
		ha->Fill(Pi()*raggio*raggio);
	}
	raggio_medio=raggio_medio/n_iter;
	circ_media=circ_media/n_iter;
	area_media=area_media/n_iter;
	varianzar=varianzar/n_iter;
	varianzar=varianzar-raggio_medio*raggio_medio;
	double raggiomh=hr->GetMean();
	double circmh=hc->GetMean();
	double areamh=ha->GetMean();
	TCanvas *cr = new TCanvas("cr","cr",800,600);
	hr->Draw();
	TCanvas *cc = new TCanvas("cc","cc",800,600);
	hc->Draw();
	TCanvas *ca = new TCanvas("ca","ca",800,600);
	ha->Draw();
	cout<<"raggio medio istogramma "<<raggiomh<<" media calcolata "<<raggio_medio<<" media popolazione "<< media_raggio<<endl;
	cout<<"circonferenza media istogramma "<<circmh<<" media calcolata "<<circ_media<<" media popolazione "<< 2*Pi()*media_raggio <<endl;
	cout<<"area media istogramma "<<areamh<<" media calcolata "<<area_media<<" media popolazione "<< Pi()*media_raggio*media_raggio <<endl;
	double correzione= Pi()*varianzar;
	double areaattesa = Pi()*raggio_medio*raggio_medio;
	double circattesa = 2*Pi()*raggio_medio;
	double areacorretta = areaattesa+correzione;
	cout<<" area attesa,correzione,areacorretta "<< areaattesa << " "<< correzione <<" "<< areacorretta << endl;
	cout<<" circonferenza attesa "<< circattesa << endl; 
}