//makro do narysowania 1-wymiarowego histogramu
//dane wczytywane z pliku
//dodatkowo dopasowanie Gaussa do histogramu (gotowa funkcja 'gaus')

{
gROOT->Reset(); 

ifstream my_input;    
my_input.open("mean_pt_pbpb.text");   

TH1D *pt_pbpb = new TH1D("pt_pbpb", "", 100, 0.22, 0.4);


double pt;  
for(int i=1; i<=150000; i=i+1)
{
my_input>>pt;        
pt_pbpb->Fill(pt);   
}

pt_pbpb->SetFillColor(33);
pt_pbpb->GetXaxis()->SetTitle("M(p_{T}) (GeV/c)");
pt_pbpb->GetYaxis()->SetTitleOffset(1.2); //odleglosc opisu od osi
pt_pbpb->GetXaxis()->SetTitleOffset(1.1);
pt_pbpb->GetYaxis()->SetTitle("entries (arb. units)");
pt_pbpb->Draw();


//deklarujemy 1-wym. funkcje postaci Gaussa (typ "gaus")
//funkcja ma nazwe gauss_dop

TF1 *gauss_dop = new TF1("gauss_dop", "gaus"); //inne gotowe funkcje do obejrzenia w klasie TMath
pt_pbpb->Fit("gauss_dop");               //ta komenda dopasowujemy funkcje typu gaus (o nazwie gauss_dop)
                                         //do histogramu pt_pbpb
                                         //parametry dopasowania wypisane na ekranie


}