#include <iostream>
#include <math>
#include "AtlasUtils.h"
// #ifndef __CINT__
#include "AtlasStyle.C"
// #endif

using namespace std;

void testplot()
{
#ifdef __CINT__
  gROOT->LoadMacro("AtlasUtils.C");
#endif
  SetAtlasStyle();
  
  //my_h1= new TFile ("effEtaMinorDATA2011_trk2011nuovo.root", "RECREATE"); 
  //my_h1= new TFile ("effEtaMinorDATA2010.root", "RECREATE"); 
  //TFile *_file1 = TFile::Open("testDATA2010noemi.root");  
  //TFile *_file1 = TFile::Open("test2011_MODtrkJpsinuovo.root");
  TFile *_file2 = TFile::Open("testMCnoemi.root");
  int nbinsEta=18;
  
  float Eff[18];
  float ErrEff[18];
  float Eta[18];
  float ErrEta[18];
  float EffMC[18];
  float ErrEffMC[18];
  float EtaMC[18];
  float ErrEtaMC[18];
  Float_t XStep[18];
  Float_t XStepLow;
  Float_t XStepMC[18];
  Float_t XStepLowMC;
  /* _file1->cd();
  cout<<"ciao\n";
  for(int i=1;i<=nbinsEta;i++)
    {
      if(EtaMinor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	
	  XStep[i-1]=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLow=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  Eta[i-1]=XStepLow+(XStep[i-1]/2.);
	Eff[i-1]=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/EtaMinor_Jpsi_EFMU4->GetBinContent(i);
	  cout<<"efficienza "<<Eff[i-1]<<endl;
	  ///metti la formula per il calcolo dell'errore sull'efficienza
	  ErrEff[i-1]=sqrt(Eff[i-1]*(1-Eff[i-1])/WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
	  cout<<"Entrie"<<WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	  cout<<"Errore"<<ErrEff[i-1]<<endl;
	  //metti la formula per il calcolo dell'errore sul pT (larghezza del bin/2)
	  ErrEta[i-1]=(XStep[i-1]/2.);

	}
    }

 ///definizione del nuovo histos con gli errori
  TGraphErrors *heff = new TGraphErrors(nbinsEta,Eta,Eff,ErrEta,ErrEff);
TCanvas *c1 = new TCanvas("c1","",10,10,700,500);
  /////////PLOT
 heff->SetMarkerColor(1);
 //heff->Draw("AP");

 my_h1->cd();
 heff->Write("heff");

 my_h1->Write(); 
 my_h1->Close();
  

  */
for(int i=1;i<=nbinsEta;i++)
    {
      if(EtaMinor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	
	  XStepMC[i-1]=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLowMC=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  EtaMC[i-1]=XStepLowMC+(XStepMC[i-1]/2.);
	EffMC[i-1]=WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/EtaMinor_Jpsi_EFMU4->GetBinContent(i);
	  cout<<"efficienza "<<EffMC[i-1]<<endl;
	  ///metti la formula per il calcolo dell'errore sull'efficienza
	  ErrEffMC[i-1]=sqrt(EffMC[i-1]*(1-EffMC[i-1])/WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
	  cout<<"Entrie"<<WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	  cout<<"Errore"<<ErrEffMC[i-1]<<endl;
	  //metti la formula per il calcolo dell'errore sul pT (larghezza del bin/2)
	  ErrEtaMC[i-1]=(XStepMC[i-1]/2.);

	}
    }
my_h2= new TFile ("effEtaMinorMC.root", "RECREATE");
TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
 ///definizione del nuovo histos con gli errori
  TGraphErrors *heffMC = new TGraphErrors(nbinsEta,EtaMC,EffMC,ErrEtaMC,ErrEffMC);
TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
  /////////PLOT
 heffMC->SetMarkerColor(2);
 //se invece di "AP" metti "e3" ti disegna una banda colorata larga quanto gli errori 
 my_h2->cd();
 heffMC->Write("heffMC");

 my_h2->Write(); 
 my_h2->Close();
  
}


void show(){

#ifdef __CINT__
  gROOT->LoadMacro("AtlasUtils.C");
#endif
  SetAtlasStyle();
TCanvas *c1 = new TCanvas("c1","",10,10,700,500);
  TFile *_file0 = TFile::Open("effEtaMinorMC.root");
  TMultiGraph *mg = new TMultiGraph();
  mg->Add(heffMC);
  heffMC->SetLineColor(2);
  mg->Draw("e3");
  mg->Draw("AP");
  mg->Draw("AP");
  mg->Draw("e3");
  TFile *_file1 = TFile::Open("effEtaMinorDATA2010.root");
  mg->Add(heff);
  mg->Draw("AP");
  mg->Draw("e3");
TFile *_file2 = TFile::Open("effEtaMinorDATA2011_trk2011nuovo.root");
heff->SetLineColor(4);
  mg->Add(heff);
  mg->Draw("AP");
  mg->Draw("e3");

}