#include <iostream>
#include <cmath>

#include "Rtypes.h"

#include "atlasstyle-00-03-01/AtlasUtils.h"
// #ifndef __CINT__
#include "atlasstyle-00-03-01/AtlasStyle.C"
// #endif

#include "TCanvas.h"
#include "TFile.h"
#include "TROOT.h"
#include "TH1F.h"
#include "TRandom.h"
#include "TGraphErrors.h"

using namespace std;

void threshold ()
{
#ifdef __CINT__
gROOT->LoadMacro("atlasstyle-00-03-01/AtlasUtils.C");
#endif

SetAtlasStyle();
Int_t icol1=5;
  Int_t icol2=5;

  // TFile * file10 = TFile::Open("Efficiency_EFMu4EFMu10_combcomb.root");
  TFile * file6 = TFile::Open("../LowPt_jpsi_note/Efficiency_EFMu4EFMu6_combcomb.root");
  TFile * file6Any = TFile::Open("../LowPt_jpsi_note/Efficiency_EFMu4EFMu6_combAny.root");
  TFile * file4 = TFile::Open("../LowPt_jpsi_note/Efficiency_EFMu4EFMu4_combcomb.root");
  TFile * file4Any = TFile::Open("../LowPt_jpsi_note/Efficiency_EFMu4EFMu4_combAny.root");

  
 
   TCanvas *c1 = new TCanvas("c1","",10,10,700,500);
 file4->cd();
  heff_bkgsub->SetMinimum(0.);
  heff_bkgsub->SetMaximum(1.1);
  heff_bkgsub->SetLineColor(1);
  heff_bkgsub->SetMarkerColor(1);
  heff_bkgsub->SetMarkerStyle(20);
  heff_bkgsub->Draw("AP"); 
  heff_bkgsub->GetXaxis()->SetTitle("p_{T} (GeV)");
  heff_bkgsub->GetYaxis()->SetTitle("efficiency w.r.t. to L1 and offline");

  file6->cd();
  heff->SetMinimum(0.);
  heff->SetMaximum(1.1);
  heff->SetLineColor(2);
  heff->SetMarkerColor(2);
  heff->SetMarkerStyle(21);
  heff->Draw("P");

  /*
 file10->cd();
  heff_bkgsub->SetMinimum(0.);
  heff_bkgsub->SetMaximum(1.1);
  heff_bkgsub->SetLineColor(4);
  heff_bkgsub->SetMarkerColor(4);
  heff_bkgsub->SetMarkerStyle(22);
  heff_bkgsub->Draw("P");*/

  ATLAS_LABEL(0.55,0.45);
  myText( 0.67,0.45,1,"Preliminary");

  myText(      0.6,0.4,1,"CB muons");
  myText(      0.6,0.35,1,"barrel");
  myMarkerText(0.6,0.3,1,20,"EFMu4",1.3);
  myMarkerText(0.6,0.25,2,21,"EFMu6",1.3);
  //  myMarkerText(0.6,0.2,4,22,"EFMu10",1.3);
  
//  myBoxText(   0.6,0.25,0.05,icol1,"MC: prompt J/#psi->#mu0#mu0");

    c1->Print("../LowPt_jpsi_note/figure/TurnOn_thresholds_data_combcomb.eps");

  
  /////ENDCAP

   TCanvas *c2 = new TCanvas("c2","",10,10,700,500);
   
 
file4->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(1);
  heff_ec->SetMarkerColor(1);
  heff_ec->SetMarkerStyle(20);
  heff_ec->Draw("AP");
   heff_ec->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_ec->GetYaxis()->SetTitle("EF_MU10 efficiency w.r.t. to L1 and offline");
   
  
  file6->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(2);
  heff_ec->SetMarkerColor(2);
  heff_ec->SetMarkerStyle(21);
  heff_ec->Draw("P");

  /*
 file10->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(4);
  heff_ec->SetMarkerColor(4);
  heff_ec->SetMarkerStyle(22);
  heff_ec->Draw("P");*/

  ATLAS_LABEL(0.55,0.45);
  myText( 0.67,0.45,1,"Preliminary");

  myText(      0.6,0.4,1,"CB muons");
  myText(      0.6,0.35,1,"endcaps");
  myMarkerText(0.6,0.3,1,20,"EFMu4",1.3);
  myMarkerText(0.6,0.25,2,21,"EFMu6",1.3);
  //  myMarkerText(0.6,0.2,4,22,"EFMu10",1.3);
  
  c2->Print("../LowPt_jpsi_note/figure/TurnOn_thresholds_data_endcap_combcomb.eps");


   /////// BARREL
   TCanvas *c3 = new TCanvas("c3","",10,10,700,500);
   file4->cd();
   heff_bkgsub->SetMinimum(0.);
   heff_bkgsub->SetMaximum(1.1);
   heff_bkgsub->SetLineColor(2);
   heff_bkgsub->SetMarkerColor(2);
   heff_bkgsub->SetMarkerStyle(21);
   heff_bkgsub->Draw("AP"); 
   heff_bkgsub->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_bkgsub->GetYaxis()->SetTitle("efficiency w.r.t. to L1 and offline");
     
   file6->cd();
   heff_bkgsub->SetMinimum(0.);
   heff_bkgsub->SetMaximum(1.1);
   heff_bkgsub->SetLineColor(4);
   heff_bkgsub->SetMarkerColor(4);
   heff_bkgsub->SetMarkerStyle(22);
   heff_bkgsub->Draw("P");
   
file4Any->cd();
   heff_bkgsub->SetMinimum(0.);
   heff_bkgsub->SetMaximum(1.1);
   heff_bkgsub->SetLineColor(2);
   heff_bkgsub->SetMarkerColor(2);
   heff_bkgsub->SetMarkerStyle(24);
   heff_bkgsub->Draw("P"); 
   heff_bkgsub->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_bkgsub->GetYaxis()->SetTitle("efficiency w.r.t. to L1 and offline");
     
   file6Any->cd();
   heff_bkgsub->SetMinimum(0.);
   heff_bkgsub->SetMaximum(1.1);
   heff_bkgsub->SetLineColor(4);
   heff_bkgsub->SetMarkerColor(4);
   heff_bkgsub->SetMarkerStyle(26);
   heff_bkgsub->Draw("P");

   ATLAS_LABEL(0.55,0.45);
   myText( 0.67,0.45,1,"Preliminary");
   
   // myText(      0.6,0.4,1,"CB muons");
   myText(      0.6,0.4,1,"barrel");
   //myMarkerText(0.6,0.3,1,20,"EFMu4",1.3);
   myMarkerText(0.6,0.35,2,21,"EFMu4/EFMu6 CB-CB",1.3);
   myMarkerText(0.6,0.3,4,22,"EFMu4/EFMu4 CB-Any",1.3);
   myMarkerText(0.6,0.25,2,21,"EFMu4/EFMu6 CB-CB",1.3);
   myMarkerText(0.6,0.2,4,22,"EFMu4/EFMu4 CB-Any",1.3);
  
   c3->Print("../LowPt_jpsi_note/figure/TurnOn_thresholdsEFMu4EFMu6-EFMu6EFMu6_data_combcomb.eps");

 /////ENDCAP

   TCanvas *c4 = new TCanvas("c4","",10,10,700,500);
   
   
   file4->cd();
   heff_bkgsub_ec->SetMinimum(0.);
   heff_bkgsub_ec->SetMaximum(1.1);
   heff_bkgsub_ec->SetLineColor(2);
   heff_bkgsub_ec->SetMarkerColor(2);
   heff_bkgsub_ec->SetMarkerStyle(21);
   heff_bkgsub_ec->Draw("AP");
   heff_bkgsub_ec->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_bkgsub_ec->GetYaxis()->SetTitle("EF_MU10 efficiency w.r.t. to L1 and offline");

   file6->cd();
   heff_bkgsub_ec->SetMinimum(0.);
   heff_bkgsub_ec->SetMaximum(1.1);
   heff_bkgsub_ec->SetLineColor(4);
   heff_bkgsub_ec->SetMarkerColor(4);
   heff_bkgsub_ec->SetMarkerStyle(22);
   heff_bkgsub_ec->Draw("P");
   
   file4Any->cd();
   heff_bkgsub_ec->SetMinimum(0.);
   heff_bkgsub_ec->SetMaximum(1.1);
   heff_bkgsub_ec->SetLineColor(2);
   heff_bkgsub_ec->SetMarkerColor(2);
   heff_bkgsub_ec->SetMarkerStyle(24);
   heff_bkgsub_ec->Draw("P");
   heff_bkgsub_ec->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_bkgsub_ec->GetYaxis()->SetTitle("EF_MU10 efficiency w.r.t. to L1 and offline");

   file6Any->cd();
   heff_bkgsub_ec->SetMinimum(0.);
   heff_bkgsub_ec->SetMaximum(1.1);
   heff_bkgsub_ec->SetLineColor(4);
   heff_bkgsub_ec->SetMarkerColor(4);
   heff_bkgsub_ec->SetMarkerStyle(26);
   heff_bkgsub_ec->Draw("P");
   
   ATLAS_LABEL(0.55,0.45);
   myText( 0.67,0.45,1,"Preliminary");
   
   //  myText(      0.6,0.4,1,"CB muons");
   myText(      0.6,0.4,1,"endcaps");
  myMarkerText(0.6,0.35,2,21,"EFMu4/EFMu6 CB-CB",1.3);
   myMarkerText(0.6,0.3,4,22,"EFMu4/EFMu4 CB-Any",1.3);
   myMarkerText(0.6,0.25,2,21,"EFMu4/EFMu6 CB-CB",1.3);
   myMarkerText(0.6,0.2,4,22,"EFMu4/EFMu4 CB-Any",1.3);
   
  c4->Print("../LowPt_jpsi_note/figure/TurnOn_thresholdsEFMu4EFMu6-EFMu6EFMu6_data_endcap_combcomb.eps");

  ///////////////////////////////

  //////CB ANY
  TCanvas *c5 = new TCanvas("c5","",10,10,700,500);
 file4Any->cd();
  heff_bkgsub->SetMinimum(0.);
  heff_bkgsub->SetMaximum(1.1);
  heff_bkgsub->SetLineColor(1);
  heff_bkgsub->SetMarkerColor(1);
  heff_bkgsub->SetMarkerStyle(24);
  heff_bkgsub->Draw("AP"); 
  heff_bkgsub->GetXaxis()->SetTitle("p_{T} (GeV)");
  heff_bkgsub->GetYaxis()->SetTitle("efficiency w.r.t. to L1 and offline");

  file6Any->cd();
  heff->SetMinimum(0.);
  heff->SetMaximum(1.1);
  heff->SetLineColor(2);
  heff->SetMarkerColor(2);
  heff->SetMarkerStyle(25);
  heff->Draw("P");

  /*file10->cd();
  heff_bkgsub->SetMinimum(0.);
  heff_bkgsub->SetMaximum(1.1);
  heff_bkgsub->SetLineColor(4);
  heff_bkgsub->SetMarkerColor(4);
  heff_bkgsub->SetMarkerStyle(22);
  heff_bkgsub->Draw("P");*/

  ATLAS_LABEL(0.55,0.45);
  myText( 0.67,0.45,1,"Preliminary");

  myText(      0.6,0.4,1,"CB+ST muons");
  myText(      0.6,0.35,1,"barrel");
  myMarkerText(0.6,0.3,1,24,"EFMu4",1.3);
  myMarkerText(0.6,0.25,2,25,"EFMu6",1.3);
  //  myMarkerText(0.6,0.2,4,26,"EFMu10",1.3);
  
//  myBoxText(   0.6,0.25,0.05,icol1,"MC: prompt J/#psi->#mu0#mu0");

    c5->Print("../LowPt_jpsi_note/figure/TurnOn_thresholds_data_combAny.eps");

  
  /////ENDCAP

   TCanvas *c6 = new TCanvas("c6","",10,10,700,500);
   
 
file4->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(1);
  heff_ec->SetMarkerColor(1);
  heff_ec->SetMarkerStyle(24);
  heff_ec->Draw("AP");
   heff_ec->GetXaxis()->SetTitle("p_{T} (GeV)");
   heff_ec->GetYaxis()->SetTitle("EF_MU10 efficiency w.r.t. to L1 and offline");
   
  
  file6->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(2);
  heff_ec->SetMarkerColor(2);
  heff_ec->SetMarkerStyle(25);
  heff_ec->Draw("P");
  /*
 file10->cd();
  heff_ec->SetMinimum(0.);
  heff_ec->SetMaximum(1.1);
  heff_ec->SetLineColor(4);
  heff_ec->SetMarkerColor(4);
  heff_ec->SetMarkerStyle(26);
  heff_ec->Draw("P");*/

  ATLAS_LABEL(0.55,0.45);
  myText( 0.67,0.45,1,"Preliminary");

  myText(      0.6,0.4,1,"CB+ST muons");
  myText(      0.6,0.35,1,"endcaps");
  myMarkerText(0.6,0.3,1,24,"EFMu4",1.3);
  myMarkerText(0.6,0.25,2,25,"EFMu6",1.3);
  // myMarkerText(0.6,0.2,4,26,"EFMu10",1.3);
  
  c6->Print("../LowPt_jpsi_note/figure/TurnOn_thresholds_data_endcap_combAny.eps");

}