{
 
  TFile *F1= new TFile("../test.root","read");
  TFile *outfile= new TFile("test_divide.root", "RECREATE");
  //histo useful for efficiency
  //numerators
  TH1D *PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                    ; 
  TH1D *PtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                  ; 
  TH1D *EtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("EtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                   ; 
  TH1D *EtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("EtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                   ; 
  TH1D *DeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("DeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                    ; 
  
  //Barrel
  //numerators
  TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b")                    ; 
  TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b")                  ; 
 
   
  
  //denominators
  TH1D *PtMajor_Jpsi_EFMU4b=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4b")                                          ; 
  TH1D *PtMinor_Jpsi_EFMU4b=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4b")                                        ; 
 
  
 //EndCap
  //numerators
  TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e")                    ; 
  TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e")                  ;  
 
  //denominators
  TH1D *PtMajor_Jpsi_EFMU4e=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4e")                                          ; 
  TH1D *PtMinor_Jpsi_EFMU4e=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4e")                                        ; 
  
 
  //numerators
  TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                    ; 
  TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                  ; 
  TH1D *WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                   ; 
  TH1D *WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                   ; 
  TH1D *WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4")                    ; 

  TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis")                    ; 
  TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis")                  ; 
  //denominators
  TH1D *PtMajor_Jpsi_EFMU4=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4")                                          ; 
  TH1D *PtMinor_Jpsi_EFMU4=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4")                                        ; 
  TH1D *EtaMajor_Jpsi_EFMU4=(TH1D*)F1->Get("EtaMajor_Jpsi_EFMU4")                                         ; 
  TH1D *EtaMinor_Jpsi_EFMU4=(TH1D*)F1->Get("EtaMinor_Jpsi_EFMU4")                                         ; 
  TH1D *DeltaR_Jpsi_EFMU4=(TH1D*)F1->Get("DeltaR_Jpsi_EFMU4")                                          ; 
  TH1D *PtMajor_Jpsi_EFMU4_mis=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4_mis")                                          ; 
  TH1D *PtMinor_Jpsi_EFMU4_mis=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4_mis")                                        ; 

 
  TGraphAsymmErrors *Eff_ptMajor_All_noW=new TGraphAsymmErrors();
  Eff_ptMajor_All_noW->BayesDivide(PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4);
  TGraphAsymmErrors *Eff_ptMajor_All=new TGraphAsymmErrors();
  Eff_ptMajor_All->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4);
  TGraphAsymmErrors *Eff_ptMajor_Barrel=new TGraphAsymmErrors();  
  Eff_ptMajor_Barrel->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b,PtMajor_Jpsi_EFMU4b);
  TGraphAsymmErrors *Eff_ptMajor_Endcap=new TGraphAsymmErrors();  
  Eff_ptMajor_Endcap->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e,PtMajor_Jpsi_EFMU4e);
  //TGraphAsymmErrors *Eff_ptMajor_mis=new TGraphAsymmErrors();
  //Eff_ptMajor_mis->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis,PtMajor_Jpsi_EFMU4_mis);  
  TGraphAsymmErrors *Eff_ptMinor_All=new TGraphAsymmErrors();
  Eff_ptMinor_All->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMinor_Jpsi_EFMU4);
  TGraphAsymmErrors *Eff_ptMinor_Barrel=new TGraphAsymmErrors();  
  Eff_ptMinor_Barrel->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b,PtMinor_Jpsi_EFMU4b);
  TGraphAsymmErrors *Eff_ptMinor_Endcap=new TGraphAsymmErrors();  
  Eff_ptMinor_Endcap->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e,PtMinor_Jpsi_EFMU4e);
  //TGraphAsymmErrors *Eff_ptMinor_mis=new TGraphAsymmErrors();
  //Eff_ptMinor_mis->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis,PtMinor_Jpsi_EFMU4_mis);
  TGraphAsymmErrors *Eff_etaMajor_All=new TGraphAsymmErrors();
  Eff_etaMajor_All->BayesDivide(WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,EtaMajor_Jpsi_EFMU4);
  TGraphAsymmErrors *Eff_etaMinor_All=new TGraphAsymmErrors();
  Eff_etaMinor_All->BayesDivide(WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,EtaMinor_Jpsi_EFMU4);
  TGraphAsymmErrors *Eff_DeltaR_All=new TGraphAsymmErrors();
  Eff_DeltaR_All->BayesDivide(WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,DeltaR_Jpsi_EFMU4);
 

  Eff_ptMajor_All->Write("Eff_ptMajor_All");
  //Eff_ptMajor_mis->Write("Eff_ptMajor_mis");
  Eff_ptMajor_Barrel->Write("Eff_ptMajor_Barrel");  
  Eff_ptMajor_Endcap->Write("Eff_ptMajor_Endcap");
  Eff_ptMajor_All_noW->Write("Eff_ptMajor_All_noW");
  
  Eff_ptMinor_All->Write("Eff_ptMinor_All");
  //Eff_ptMinor_mis->Write("Eff_ptMinor_mis");
  Eff_ptMinor_Barrel->Write("Eff_ptMinor_Barrel");  
  Eff_ptMinor_Endcap->Write("Eff_ptMinor_Endcap");  
  Eff_etaMajor_All->Write("Eff_etaMajor_All");
  Eff_etaMinor_All->Write("Eff_etaMinor_All");
  Eff_DeltaR_All->Write("Eff_DeltaR_All");
  
  
  outfile->Close();
}