#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 ("effEtaMajorDATA2010.root", "RECREATE"); 
  //TFile *_file1 = TFile::Open("testDATA2010noemi.root");
   my_h1= new TFile ("effEtaMajorDATA2011_trk2011nuovo.root", "RECREATE"); 
   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;
  
  cout<<"ciao\n";
  for(int i=1;i<=nbinsEta;i++)
    {
      if(EtaMajor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	
	  XStep[i-1]=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLow=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  Eta[i-1]=XStepLow+(XStep[i-1]/2.);
	Eff[i-1]=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/EtaMajor_Jpsi_EFMU4->GetBinContent(i);
	  cout<<"efficienza "<<Eff[i-1]<<endl;
	  ///metti la formula per il calcolo dell'errore sull'efficienza
	  ErrEff[i-1]=0;
	  if((i-1)!=9)
	    	  ErrEff[i-1]=sqrt(Eff[i-1]*(1-Eff[i-1])/WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
	  cout<<"Entrie"<<WEtaMajor_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(EtaMajor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	
	  XStepMC[i-1]=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLowMC=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  EtaMC[i-1]=XStepLowMC+(XStepMC[i-1]/2.);
	EffMC[i-1]=WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/EtaMajor_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])/WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
	  cout<<"Entrie"<<WEtaMajor_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 ("effEtaMajorMC.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);
  /////////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("effEtaMajorMC.root");
  TMultiGraph *mg = new TMultiGraph();
mg->SetTitle("Efficienza;pT(GeV);Efficienza");
  mg->Add(heffMC);
  heffMC->SetLineColor(2);
  mg->Draw("e3");
  mg->Draw("AP");
  mg->Draw("AP");
  mg->Draw("e3");
  TFile *_file1 = TFile::Open("effEtaMajorDATA2010.root");
  mg->Add(heff);
  mg->Draw("AP");
  mg->Draw("e3");
TFile *_file2 = TFile::Open("effEtaMajorDATA2011_trk2011nuovo.root");
heff->SetLineColor(4);
  mg->Add(heff);
  mg->Draw("AP");
  mg->Draw("e3");
myMarkerText(0.7, 0.40,2,20, "MC",1.0);
  myMarkerText(0.7, 0.40,4,20, "DATA 2011",1.0);
myMarkerText(0.7, 0.40,1,20, "DATA 2010",1.0);

}