#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 ("effDeltaRDATA2011_trk2011nuovo.root", "RECREATE"); 
  //my_h1= new TFile ("effDeltaRDATA2010.root", "RECREATE"); 
  //TFile *_file1 = TFile::Open("testDATA2010noemi.root");  
  //TFile *_file1 = TFile::Open("test2011_MODtrkJpsinuovo.root");
  TFile *_file2 = TFile::Open("testMCnoemi.root");
 
  int nbinsDeltaR=20;
  
  float Eff[20];
  float ErrEff[20];
  float DeltaR[20];
  float ErrDeltaR[20];
  float EffMC[20];
  float ErrEffMC[20];
  float DeltaRMC[20];
  float ErrDeltaRMC[20];
  Float_t XStep[20];
  Float_t XStepLow;
  Float_t XStepMC[20];
  Float_t XStepLowMC;
 
  cout<<"ciao\n";
  for(int i=1;i<=nbinsDeltaR;i++)
    {  
	
	  XStep[i-1]=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLow=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  DeltaR[i-1]=XStepLow+(XStep[i-1]/2.);
	   ErrDeltaR[i-1]=(XStep[i-1]/2.);
	 
	   if(DeltaR_Jpsi_EFMU4->GetBinContent(i)!=0)
	     { 
	       cout<<i<<endl;
	       Eff[i-1]=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/DeltaR_Jpsi_EFMU4->GetBinContent(i);
	       cout<<"efficienza "<<Eff[i-1]<<endl;
	       ///metti la formula per il calcolo dell'errore sull'efficienza
	       if(Eff[i-1]==0) 
		 {ErrEff[i-1]=0;}
	       else
		 {ErrEff[i-1]=sqrt(Eff[i-1]*(1-Eff[i-1])/WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));}
	       cout<<"Entrie"<<WDeltaR_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)
	     }
	   else
	     {
	       cout<<i<<endl;
	       Eff[i-1]=0;
	       ErrEff[i-1]=0;
	       cout<<"efficienza "<<Eff[i-1]<<endl;

	       cout<<"Entrie"<<WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	       cout<<"Errore"<<ErrEff[i-1]<<endl;
	     }
	  
	 
	
	}
    
my_h1= new TFile ("effDeltaRMC.root", "RECREATE");
 ///definizione del nuovo histos con gli errori
  TGraphErrors *heff = new TGraphErrors(nbinsDeltaR,DeltaR,Eff,ErrDeltaR,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<=nbinsDeltaR;i++)
    {
      if(DeltaR_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	
	  XStepMC[i-1]=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLowMC=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  DeltaRMC[i-1]=XStepLowMC+(XStepMC[i-1]/2.);
	EffMC[i-1]=WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/DeltaR_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])/WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
	  cout<<"Entrie"<<WDeltaR_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)
	  ErrDeltaRMC[i-1]=(XStepMC[i-1]/2.);

	}
    }
my_h2= new TFile ("effDeltaRMC.root", "RECREATE");
//TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
 ///definizione del nuovo histos con gli errori
  TGraphErrors *heffMC = new TGraphErrors(nbinsDeltaR,DeltaRMC,EffMC,ErrDeltaRMC,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("effDeltaRMC.root");
  TMultiGraph *mg = new TMultiGraph();
 mg->SetTitle("Efficienza;DeltaR;Efficienza");
mg->Add(heff);
heff->SetLineColor(2);
  mg->Draw("e3");
  mg->Draw("AP");
  mg->Draw("AP");
  mg->Draw("e3");
    TFile *_file1 = TFile::Open("effDeltaRDATA2010.root");
  mg->Add(heff);
  mg->Draw("AP");
   mg->Draw("e3");
   //TFile *_file2 = TFile::Open("effDeltaRDATA2011_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);

}