#include "TFile.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TGaxis.h"
#include "TRandom.h"
#include "TLegend.h"
#include "TPaveStats.h"
#include "TGraph.h"

void hnenmumed()
{
   
   char datafilegamma[200]= "G:/cors/gamma/gammatot/distribuzioni-emmu.txt";
   char datafileprot[200]= "G:/cors/prot/prottot/distribuzioni-emmu.txt";
   char dataout[200]="G:/cors/gamma/gammaprot_hfirstint-med.png";
   
   gStyle->SetOptFit();
   TCanvas *c36 = new TCanvas("c36","multigraph",1280,1024);
   TPad *pad1  = new TPad("pad1", "pad1", 0, 0.3, 1, 1.0);
   float offx=1.1;
   float offy=1;
   float margr=0.18;
   float margl=0.12;
   //gPad->SetLeftMargin(margl);
   pad1->SetGridx();         // Vertical grid
   pad1->Draw();             // Draw the upper pad: pad1
   pad1->cd();    
   //gPad->SetRightMargin(margr);
   //c36->SetGrid();
   //TMultiGraph *mg = new TMultiGraph();
    //gamma secondari
   TGraphErrors *gamma = new TGraphErrors(datafilegamma,"%lg %lg %lg %lg");
   gamma->SetMarkerColor(kBlue);
   gamma->SetLineColor(kBlue);
   gamma->SetMarkerStyle(23);
   gamma->SetMarkerSize(2);
   gamma->GetXaxis()->SetTitle("N_{e.m.}");
   gamma->GetYaxis()->SetTitle("N_{#mu^{#pm}}");
   gamma->GetHistogram()->SetTitle("Numero di #gamma e e^{#pm} vs numero di #mu^{#pm} in EAS indotti da #gamma e #font[12]{p}");
   //gr1->Fit("pol1","q");
   gamma->Draw();
   c36->cd();    
   TPad *pad2 = new TPad("pad2", "pad2", 0, 0.05, 1, 0.5);
  // pad2->SetTopMargin(0.4);
   pad2->SetBottomMargin(1);
   pad2->SetGridx(); // vertical grid
   pad2->Draw();
   pad2->cd();       // pad2 becomes the current pad
   c36->Update();
   
   //elettroni secondari
   
   TGraphErrors *prot = new TGraphErrors(datafileprot,"%lg %lg %lg %lg");
   prot->SetMarkerColor(kGreen);
   prot->SetLineColor(kGreen);
   prot->SetMarkerStyle(29);
   prot->SetMarkerSize(2);
   prot->GetHistogram()->SetTitle("");
    prot->Draw();
   //force drawing of canvas to generate the fit TPaveStats
  //mg->SetTitle("Energia- Quota di prima interazione dei #gamma e #font[12]{p}");
   prot->GetXaxis()->SetTitle("N_{e.m.}");
   prot->GetYaxis()->SetTitle("N_{#mu^{#pm}}");
   prot->GetYaxis()->SetTitleSize(26);
   prot->GetYaxis()->SetTitleFont(43);
   prot->GetYaxis()->SetTitleOffset(1.55);
   prot->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
   prot->GetYaxis()->SetLabelSize(15);
    prot->GetXaxis()->SetTitleSize(26);
   prot->GetXaxis()->SetTitleFont(43);
   prot->GetXaxis()->SetTitleOffset(1.55);
   prot->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
   prot->GetXaxis()->SetLabelSize(15);
  prot->GetYaxis()->SetTitleOffset(offy);
   c36->Modified();
   //force drawing of canvas to generate the fit TPaveStats
   c36->Update();
    TLegend* leg10 = new TLegend(0.7, 0.3, .55, .5);
    leg10->SetNColumns(1);
    leg10->AddEntry(gamma, " EAS indotti da #gamma", "lp");
    leg10->AddEntry(prot, "EAS indotti da #font[12]{p}", "lp");
    leg10->Draw("same");
   c36->Modified();
   gROOT->ProcessLine("gROOT->SetBatch()");
   c36->Print(dataout);delete c36;
   gROOT->ProcessLine("gROOT->SetBatch(kFALSE)");
   exit();
  
   }