{
  TFile * file = TFile::Open("../TPout_GHI_efmu4efmu4_CombComb_staco_new.root");
  TFile * file0 = TFile::Open("PlotsEffJpsi.root");
  // TFile * file3 = TFile::Open("../TPout_EF_EFMu6EFMu10_new_CombComb.root");
  TFile *outfile= new TFile("Efmu4TurnOn.root", "RECREATE");
   
  /* file->cd();
  TGraphAsymmErrors* heff = new TGraphAsymmErrors();
  heff->BayesDivide(h_pt_num_ba,h_pt_den_ba);
  TGraphAsymmErrors* heffE = new TGraphAsymmErrors();
  heffE->BayesDivide(h_pt_num_ec,h_pt_den_ec);
  */ 
  TH1F     * h_probes_pass_all  =(TH1F*)    file->Get("hptmsselall");
  TH1F     * h_probes_all  =(TH1F*)    file->Get("hptmsall");

  TH1F     * h_probes_pass  =(TH1F*)    file->Get("hptmsselb");
  TH1F     * h_probes  =(TH1F*)    file->Get("hptmsb");

  TH1F     * h_probes_pass_ec  =(TH1F*)    file->Get("hptmssele");
  TH1F     * h_probes_ec  =(TH1F*)    file->Get("hptmse");

  
 TGraphAsymmErrors* heff_all = new TGraphAsymmErrors();
  heff_all->BayesDivide(h_probes_pass_all,h_probes_all);
  TGraphAsymmErrors* heff = new TGraphAsymmErrors();
  heff->BayesDivide(h_probes_pass,h_probes);
  TGraphAsymmErrors* heffE = new TGraphAsymmErrors();
  heffE->BayesDivide(h_probes_pass_ec,h_probes_ec);
    
  /////////////////////
 TCanvas *c3 = new TCanvas("c3","",10,10,700,500);
  TH1F *hr = c3->DrawFrame(0.,0.,22.,1.1);
  
  hr->GetXaxis()->SetTitle("#font[42]{p_{T} (GeV)}");
  hr->GetYaxis()->SetTitle("efficiency");
  //c1->Divide(1,2);
  // c1->cd(1);
  file->cd();
  heff_all->SetLineColor(2);
  heff_all->SetMarkerColor(2);
  heff_all->SetFillColor(3);
  heff_all->SetMarkerStyle(24);
  heff_all->SetFillStyle(3001);
  //  heff_all->Draw("P");

 
 file0->cd();
Eff_ptNoTrig_All->SetLineColor(4);
Eff_ptNoTrig_All->SetMarkerColor(4);
Eff_ptNoTrig_All->SetMarkerStyle(22);
  Eff_ptNoTrig_All->Draw("P");

  TLegend *tleg = new TLegend(0.55,0.25,0.9,0.35);
  tleg->SetFillColor(0);
  // tleg->AddEntry(heff,"J/psi","pl");
  // tleg->AddEntry(heff3,"EFMu10/EFMu10 Z","pl");
  //  tleg->AddEntry(heff_all,"efmu4","pl");
  tleg->AddEntry(Eff_ptNoTrig_All,"Eff_ptNoTrig_All","pl");
  //tleg->AddEntry(heff,"BARREL - 2GeV<invmass<100GeV","pl");
  tleg->Draw();
  
  c3->Print("TurnOn_DimuPt_periodGHI_ALL.eps");


  //////////
  

TCanvas *c1 = new TCanvas("c1","",10,10,700,500);
  TH1F *hr = c1->DrawFrame(0.,0.,22.,1.1);
  
  hr->GetXaxis()->SetTitle("#font[42]{p_{T} (GeV)}");
  hr->GetYaxis()->SetTitle("efficiency");
  //c1->Divide(1,2);
  // c1->cd(1);
  file->cd();
  heff->SetLineColor(2);
  heff->SetMarkerColor(2);
  heff->SetFillColor(3);
  heff->SetMarkerStyle(24);
  heff->SetFillStyle(3001);
  //  heff->Draw("P");

 
 file0->cd();
Eff_ptNoTrig_Barrel->SetLineColor(4);
Eff_ptNoTrig_Barrel->SetMarkerColor(4);
Eff_ptNoTrig_Barrel->SetMarkerStyle(22);
  Eff_ptNoTrig_Barrel->Draw("P");

  TLegend *tleg = new TLegend(0.55,0.25,0.9,0.35);
  tleg->SetFillColor(0);
  // tleg->AddEntry(heff,"J/psi","pl");
  // tleg->AddEntry(heff3,"EFMu10/EFMu10 Z","pl");
  //  tleg->AddEntry(heff,"efmu4","pl");
  tleg->AddEntry(Eff_ptNoTrig_Barrel,"Eff_ptNoTrig_Barrel","pl");
  //tleg->AddEntry(heff,"BARREL - 2GeV<invmass<100GeV","pl");
  tleg->Draw();
  
  c1->Print("TurnOn_DimuPt_periodGHI_barrel.eps");

 
  /////ENDCAP

 TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
  

file->cd();
  heffE->SetLineColor(2);
  heffE->SetMarkerColor(2);
  heffE->SetFillColor(3);
  heffE->SetMarkerStyle(24);
  heffE->SetFillStyle(3001);
  // heffE->Draw("AP");

 file0->cd();

Eff_ptNoTrig_Endcap->SetLineColor(4);
Eff_ptNoTrig_Endcap->SetMarkerColor(4);
Eff_ptNoTrig_Endcap->SetMarkerStyle(22);
  Eff_ptNoTrig_Endcap->Draw("AP");



 TLegend *tleg = new TLegend(0.55,0.25,0.9,0.35);
  tleg->SetFillColor(0);
  //  tleg->AddEntry(heff,"J/psi","pl");
  // tleg->AddEntry(heff3,"EFMu10/EFMu10 Z","pl");
  // tleg->AddEntry(heffE,"efmu4","pl");
  tleg->AddEntry(Eff_ptNoTrig_Endcap,"Eff_ptNoTrig_Endcap","pl");
  tleg->Draw();
 
     c2->Print("TurnOn_DimuPt_periodGHI_endcap.eps");

     outfile->cd();
     heff->Write("heffB");
     heffE->Write("heffE");
     heff_all->Write("heffAll");
     
     outfile->Write();
     outfile->Close();  
}