#include <TMath.h>



Double_t Espo (Double_t x, Double_t csi)
{
  //    if(x==0) return -10.;
    return exp(-x/csi);
}


void Esponenziale()
{

  gROOT->Reset();



  TCanvas* c1 = new TCanvas("c1","Esponenziale csi=1",180,10,700,500);
  TCanvas* c2 = new TCanvas("c2","Esponenziale csi=2",200,30,700,500);
  TCanvas* c3 = new TCanvas("c3","Esponenziale csi=3",220,50,700,500);
  TCanvas* c4 = new TCanvas("c4","Esponenziale csi=4",240,70,700,500);
  
  // Intervalli di visualizzazione delle distribuzioni
  const Int_t xMin = 0;
  const Int_t xMax = 6;
  
  // Media della distribuzione esponenziale
  Double_t csi = 0.5;
  

  

  Double_t x[xMax-xMin+1], y[xMax-xMin+1];
  for (Int_t i = 0; i<xMax-xMin+1; i++) 
    {
      x[i] = xMin + i;
      y[i] = Espo(x[i],csi)/csi;
         }
  
  // Disegno la distribuzione  nella prima canvas
  c1->cd();
  TGraph* Expo1PDF = new TGraph (xMax-xMin+1,x,y);
  Expo1PDF->SetFillColor(35);
  Expo1PDF ->SetTitle("Esponenziale csi=0.5");
  Expo1PDF->SetMarkerStyle(21);
  Expo1PDF->SetMarkerColor(4);
  Expo1PDF->Draw("APC");
  

  // Seconda distribuzione
  const Int_t xMin = 0;
  const Int_t xMax = 14;

  // Media della distribuzione Esponenziale
  Double_t csi = 1.;

  // Calcolo i valori della seconda distribuzione
  
  Double_t x1[xMax-xMin+1], y1[xMax-xMin+1];
  for (Int_t i = 0; i<xMax-xMin+1; i++) 
    { 
      x1[i] = xMin + i;
      y1[i] = Espo (x1[i],csi);
    }

  // Disegno la seconda distribuzione nella seconda canvas
  c2->cd();
  TGraph* Expo2PDF = new TGraph (xMax-xMin+1,x1,y1);
  Expo2PDF->SetFillColor(35);
  Expo2PDF ->SetTitle("Esponenziale csi=1");
  Expo2PDF->SetMarkerStyle(21);
  Expo2PDF->SetMarkerColor(2);
  Expo2PDF->Draw("APC");

  // terza  distribuzione
  const Int_t xMin = 0;
  const Int_t xMax = 20;

  // Media della distribuzione di Poisson
  Double_t csi = 2;

  // Calcolo i valori della terza distribuzione
  
  Double_t x2[xMax-xMin+1], y2[xMax-xMin+1];
  for (Int_t i = 0; i<xMax-xMin+1; i++) 
    { 
      x2[i] = xMin + i;
      y2[i] = Espo (x2[i],csi);
    }

  // Disegno la terza  poisson nella terza canvas
  c3->cd();
  TGraph* Expo3PDF = new TGraph (xMax-xMin+1,x2,y2);
  Expo3PDF->SetFillColor(35);
  Expo3PDF ->SetTitle("Esponenziale csi=2");
  Expo3PDF->SetMarkerStyle(21);
  Expo3PDF->SetMarkerColor(3);
  Expo3PDF->Draw("APC");


  // quarta  distribuzione
  const Int_t xMin = 0;
  const Int_t xMax = 35;

  // Media della distribuzione
  Double_t csi = 4;

  // Calcolo i valori della quarta poissoniana

  Double_t x3[xMax-xMin+1], y3[xMax-xMin+1];
  for (Int_t i = 0; i<xMax-xMin+1; i++)
    {
      x3[i] = xMin + i;
      y3[i] = Espo (x3[i],csi);
    }

  // Disegno la quarta  poissoniana nella quarta canvas
  c4->cd();
  TGraph* Expo4PDF = new TGraph (xMax-xMin+1,x3,y3);
  Expo4PDF->SetFillColor(35);
  Expo4PDF ->SetTitle("Esponenziale csi=4");
  Expo4PDF->SetMarkerStyle(21);
  Expo4PDF->SetMarkerColor(7);
  Expo4PDF->Draw("APC");

  
  TCanvas* all = new TCanvas ("all","All in one");
  all->Divide(2,2);
   all->cd(1);
   Expo1PDF->Draw("APC");
  leg = new TLegend(0.45,0.75,0.89,0.89); 
  leg->SetHeader("Valore di aspettazione = 0.5");
  leg->SetTextSize(0.05);
  leg->Draw();
   all->cd(2);
   Expo2PDF->Draw("APC");
  leg1 = new TLegend(0.45,0.75,0.89,0.89); 
  leg1->SetHeader("Valore di aspettazione = 1");
  leg1->SetTextSize(0.05);
  leg1->Draw();

   all_3->cd();
   Expo3PDF->Draw("APC");
   leg2 = new TLegend(0.45,0.75,0.89,0.89); 
   leg2->SetHeader("Valore di aspettazione = 2");
   leg2->SetTextSize(0.05);
   leg2->Draw();

   all_4->cd();
  leg3 = new TLegend(0.55,0.80,0.98,0.98); 
   Expo4PDF->Draw("APC");
  leg3->SetHeader("Valore di aspettazione = 4");
   leg3->SetTextSize(0.05);
  leg3->Draw();

  
   // hist->Draw();
    //  Poisson1PDF->Draw("PC");
    //Poisson2PDF->Draw("PC");
    //Poisson3PDF->Draw("PC");
  //leg->Draw();
  
}