#include <iostream>
#include <math>
using namespace std;

void testplot()
{
 TFile *_file1 = TFile::Open("testMCnoemi.root");
 TFile *_file2 = TFile::Open("testmc.root");
 const int nbinsPt=7;
  _file1->cd();
  TH1F *eff = new TH1F(*PtMajor_Jpsi_EFMU4);
  double N[nbinsPt];
  double NMC[nbinsPt];
  double epsilon[nbinsPt];
  double error[nbinsPt];
 double epsilonMC[nbinsPt];
  double errorMC[nbinsPt];
  double num;
  int i;
  for(i=0;i<=(nbinsPt-1);i++)
    {
      Double_t DenbinContent =eff->GetBinContent(i+1);
      cout<<DenbinContent;
      N[i]=DenbinContent;
      cout<<"N["<<i<<"] = "<<N[i]<<"\n";
    }
 
  eff->Divide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4);
  for(i=0;i<=(nbinsPt-1);i++)
    {
      Double_t DenbinContent =eff->GetBinContent(i+1);
      cout<<DenbinContent;
      epsilon[i]=DenbinContent;
      cout<<"epsilon["<<i<<"] = "<<epsilon[i]<<"\n";
    }
  error[0]=0;
cout<<"error[0] = "<<error[0]<<"\n";
 for(i=1;i<=(nbinsPt-1);i++)
    {
      cout<<DenbinContent;
      num=sqrt(epsilon[i]*(1-epsilon[i]));
      error[i]=num/N[i];
      cout<<"error["<<i<<"] = "<<error[i]<<"\n";
    }



  eff->SetLineColor(4);
  eff->Draw();
  _file2->cd();
 TH1F *eff = new TH1F(*PtMajor_Jpsi_EFMU4);
 eff->SetLineColor(2);
for(i=0;i<=(nbinsPt-1);i++)
    {
      Double_t DenbinContent =eff->GetBinContent(i+1);
      cout<<DenbinContent;
      NMC[i]=DenbinContent;
      cout<<"NMC["<<i<<"] = "<<NMC[i]<<"\n";
    }
 eff->Divide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4);
 eff->Draw("same");
for(i=0;i<=(nbinsPt-1);i++)
    {
      Double_t DenbinContent =eff->GetBinContent(i+1);
      cout<<DenbinContent;
      epsilonMC[i]=DenbinContent;
      cout<<"epsilonMC["<<i<<"] = "<<epsilonMC[i]<<"\n";
    }
  errorMC[0]=0;
cout<<"errorMC[0] = "<<errorMC[0]<<"\n";
 for(i=1;i<=(nbinsPt-1);i++)
    {
      cout<<DenbinContent;
      num=sqrt(epsilonMC[i]*(1-epsilonMC[i]));
      errorMC[i]=num/NMC[i];
      cout<<"errorMC["<<i<<"] = "<<errorMC[i]<<"\n";
    }
}