#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 ("effPtMajorDATA2011_trk2011nuovo.root", "RECREATE"); 
  TFile *_file1 = TFile::Open("test2011_MODtrkJpsinuovo.root");
  //my_h1= new TFile ("effPtMajorDATA2010.root", "RECREATE"); 
  //TFile *_file1 = TFile::Open("testDATA2010noemi.root");
   
  //TFile *_file2 = TFile::Open("testMCnoemi.root");
  int nbinsPt=7;
  
  float Eff[7];
  float ErrEff[7];
  float Pt[7];
  float ErrPt[7];
  float EffMC[7];
  float ErrEffMC[7];
  float PtMC[7];
  float ErrPtMC[7];
  Float_t XStep[7];
  Float_t XStepLow;
  Float_t XStepMC[7];
  Float_t XStepLowMC;
  
  cout<<"ciao\n";
  
  for(int i=1;i<=nbinsPt;i++)
    {
      
	
	  XStep[i-1]=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLow=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  Pt[i-1]=XStepLow+(XStep[i-1]/2.);
 //metti la formula per il calcolo dell'errore sul pT (larghezza del bin/2)
	  ErrPt[i-1]=(XStep[i-1]/2.);
///metti la formula per il calcolo dell'errore sull'efficienza
if(PtMajor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	  cout<<endl;
	Eff[i-1]=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/PtMajor_Jpsi_EFMU4->GetBinContent(i);
	  cout<<"efficienza "<<Eff[i-1]<<endl;
	   if(Eff[i-1]==0) 
		 {ErrEff[i-1]=0;}
	       else
		 {
		   ErrEff[i-1]=sqrt(Eff[i-1]*(1-Eff[i-1])/WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
		   if(Eff[i-1]<0.01)
		     ErrEff[i-1]=0;
		 }
	  cout<<"Entrie"<<WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	  cout<<"Errore"<<ErrEff[i-1]<<endl;
	}
 else
	     {
	       cout<<i<<endl;
	       Eff[i-1]=0;
	       ErrEff[i-1]=0;
	       cout<<"efficienza "<<Eff[i-1]<<endl;
	 cout<<"Entrie"<<WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	  cout<<"Errore"<<ErrEff[i-1]<<endl;

	}
    }

 ///definizione del nuovo histos con gli errori
  TGraphErrors *heff = new TGraphErrors(nbinsPt,Pt,Eff,ErrPt,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<=nbinsPt;i++)
  {
	  XStepMC[i-1]=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinWidth(i);
	  XStepLowMC=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetXaxis()->GetBinLowEdge(i);

	  PtMC[i-1]=XStepLowMC+(XStepMC[i-1]/2.);
 ErrPtMC[i-1]=(XStepMC[i-1]/2.);
 if(PtMajor_Jpsi_EFMU4->GetBinContent(i)!=0)
	{
	  cout<<endl;
	EffMC[i-1]=WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)/PtMajor_Jpsi_EFMU4->GetBinContent(i);
	  cout<<"efficienza "<<EffMC[i-1]<<endl;
	  ///metti la formula per il calcolo dell'errore sull'efficienza
	  if(Eff[i-1]==0) 
		 {ErrEff[i-1]=0;}
	       else
		 {
		 ErrEffMC[i-1]=sqrt(EffMC[i-1]*(1-EffMC[i-1])/WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i));
		 
		 }
	  cout<<"Entrie"<<WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	  cout<<"Errore"<<ErrEffMC[i-1]<<endl;
	 

	}
 else
   {cout<<i<<endl;
	       Eff[i-1]=0;
	       ErrEff[i-1]=0;
	       cout<<"efficienza "<<Eff[i-1]<<endl;
	       cout<<"Entrie"<<WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4->GetBinContent(i)<<endl;
	       cout<<"Errore"<<ErrEff[i-1]<<endl;
	     }
    }
my_h2= new TFile ("effPtMajorMC.root", "RECREATE");
//TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
 ///definizione del nuovo histos con gli errori
  TGraphErrors *heffMC = new TGraphErrors(nbinsPt,PtMC,EffMC,ErrPtMC,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("effPtMajorMC.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("effPtMajorDATA2010.root");

mg->Add(heff);
mg->Draw("AP");
  mg->Draw("e3");
  //  TFile *_file2 = TFile::Open("effPtMajorDATA2011_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);

}