#include <iostream>
#include "TCanvas.h"
#include "TH1.h"
#include "TF1.h"
#include "TMath.h"

#include "read.h"

const int N_PAR_BG = 2;

void set_title_histo(TH1F *&histo, Int_t tvivo)
{
    Char_t titolo[80];
    sprintf(titolo,"Spettro fondo rivelatore. Tm %d min", (Double_t)tvivo/60.);  
    histo->SetTitle(titolo);
    histo->GetXaxis()->SetTitle("Canale");
    histo->GetYaxis()->SetTitle("-dN/dCanale"); 
}

Double_t rumore(Double_t *x, Double_t *par)
{
//    return par[0] + par[1] * x[0];  //linear noise
    return par[0]*TMath::Exp(par[1]*x[0]);  //exp noise
}

void gamma_R()
{    
    gROOT->SetStyle("Plain");
    
    char filename[]  = {"232Th03.h"};
    
    Int_t tvivo, treal, ch_lo, ch_hi;
    TH1F *h;
    
    get_data_from_file(filename, tvivo, treal, ch_lo, ch_hi, h);
    
    TCanvas *c1 = new TCanvas("c1","Istogramma",100,5,700,500); 
    h->Draw();

    TF1 *BackGround = new TF1("BackGround", rumore, ch_lo, ch_hi, N_PAR_BG);
    h->Fit("BackGround");
//    h->Fit("expo");
}