TChain *InputDataEgamma = new TChain("mini");
TChain *InputDataMuons  = new TChain("mini");
TChain *InputMonteCarlo = new TChain("mini");

/*
InputDataEgamma->Add("/Users/clitemnestra/Downloads/DataEgammaA.root");
*/

InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaA.root");
/*
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaB.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaC.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaD.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaE.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaG.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaH.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaI.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaJ.root");
InputDataEgamma->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataEgammaL.root");
*/

/*
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsA.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsB.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsC.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsD.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsE.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsG.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsH.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsI.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsJ.root");
InputDataMuons->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/data/DataMuonsL.root");*/

/*
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_105985.WW.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_105986.ZZ.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_105987.WZ.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110090.stop_tchan_top.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110091.stop_tchan_antitop.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110119.stop_schan.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110140.stop_wtchan.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110899.ZPrime400.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110901.ZPrime500.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110902.ZPrime750.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110903.ZPrime1000.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110904.ZPrime1250.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110905.ZPrime1500.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110906.ZPrime1750.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110907.ZPrime2000.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110908.ZPrime2250.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110909.ZPrime2500.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_110910.ZPrime3000.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_117049.ttbar_had.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_117050.ttbar_lep.root");

InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_147770.Zee.root");

InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_147771.Zmumu.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_147772.Ztautau.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_158020.ZPrime_ee.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_158026.ZPrime_mm.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_160155.ggH125_ZZ4lep.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_160205.VBFH125_ZZ4lep.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_161005.ggH125_WW2lep.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_161055.VBFH125_WW2lep.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167740.WenuWithB.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167741.WenuJetsBVeto.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167742.WenuNoJetsBVeto.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167743.WmunuWithB.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167744.WmunuJetsBVeto.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167745.WmunuNoJetsBVeto.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167746.WtaunuWithB.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167747.WtaunuJetsBVeto.root");
InputMonteCarlo->Add("http://arturos.web.cern.ch/arturos/napoli/OutReach/mc/mc_167748.WtaunuNoJetsBVeto.root");
*/

Int_t InputEgamma_lep_n;

Float_t InputEgamma_lep_pt[5],
        InputEgamma_lep_eta[5],
        InputEgamma_lep_phi[5],
        InputEgamma_lep_type[5];

Int_t InputMuons_lep_n;

Float_t InputMuons_lep_pt[5],
        InputMuons_lep_eta[5],
        InputMuons_lep_phi[5],
        InputMuons_lep_type[5];

Int_t InputMC_lep_n;

Float_t InputMC_lep_pt[5],
        InputMC_lep_eta[5],
        InputMC_lep_phi[5],
        InputMC_lep_type[5];

InputDataEgamma->SetBranchAddress("lep_n",    &InputEgamma_lep_n);
InputDataEgamma->SetBranchAddress("lep_pt",   &InputEgamma_lep_pt);
InputDataEgamma->SetBranchAddress("lep_eta",  &InputEgamma_lep_eta);
InputDataEgamma->SetBranchAddress("lep_phi",  &InputEgamma_lep_phi);
InputDataEgamma->SetBranchAddress("lep_type", &InputEgamma_lep_type);

InputDataMuons->SetBranchAddress("lep_n",  &InputMuons_lep_n);
InputDataMuons->SetBranchAddress("lep_pt", &InputMuons_lep_pt);
InputDataMuons->SetBranchAddress("lep_eta",&InputMuons_lep_eta);
InputDataMuons->SetBranchAddress("lep_phi",&InputMuons_lep_phi);
InputDataMuons->SetBranchAddress("lep_type",&InputMuons_lep_type);

InputMonteCarlo->SetBranchAddress("lep_n",  &InputMC_lep_n);
InputMonteCarlo->SetBranchAddress("lep_pt", &InputMC_lep_pt);
InputMonteCarlo->SetBranchAddress("lep_eta",&InputMC_lep_eta);
InputMonteCarlo->SetBranchAddress("lep_phi",&InputMC_lep_phi);
InputMonteCarlo->SetBranchAddress("lep_type",&InputMC_lep_type);

TH1F *h_lep_pt   = new TH1F("h_lep_pt","h_lep_pt",100,0,250000);
TH1F *h_lep_pt_2 = new TH1F("h_lep_pt_2","h_lep_pt",100,0,250000);
TH1F *h_lep12_IM = new TH1F("h_lep12_IM","h_lep12_IM",100,0,800000);
TH1I *h_lep12_type = new TH1I("h_lep12_type","h_lep12_type",40,-20,20);

h_lep_pt->SetMarkerStyle(8);
h_lep_pt->SetMarkerSize(0.5);
h_lep_pt_2->SetLineColor(0);

h_lep_pt_2->SetMarkerStyle(8);
h_lep_pt_2->SetMarkerSize(0.5);
h_lep_pt_2->SetLineColor(2);

int nentries,nbytes,i;
nentries=(Int_t)InputDataEgamma->GetEntries();

// Variables for...
float M;

for (i = 0;i<nentries;i++) {
    nbytes = InputDataEgamma->GetEntry(i);
        
    //Uncut histogram
    h_lep_pt->Fill(InputEgamma_lep_pt[0]);
    h_lep12_type->Fill(InputEgamma_lep_type[0]);
    

    //Cut histogram
    if (InputEgamma_lep_n>1 && InputEgamma_lep_pt[1]>15000.)
    {
        h_lep_pt_2->Fill(InputEgamma_lep_pt[1]);
        
        // Calculation of the invarriant mass of the two leading leptons, using:
        // InvMass^2 = 2.pt_{lep1}.pt_{lep2}.(\cosh{(eta_{lep1}.eta_{lep2})} - \cos{(phi_{lep1}.phi_{lep2})})
        M = sqrt(2*InputEgamma_lep_pt[0]*InputEgamma_lep_pt[1]*
                 (  TMath::CosH(InputEgamma_lep_eta[0]*InputEgamma_lep_eta[1]) 
                  - TMath::Cos(InputEgamma_lep_phi[0]*InputEgamma_lep_phi[1])));
        
        
      //Fill ...
      h_lep12_IM->Fill(M);
      h_lep12_type->Fill(InputEgamma_lep_type[1]);
    }
  }

TCanvas *c1 = new TCanvas("c1","c1",50,50,800,600);
c1->SetFillColor(0);
//gStyle->SetOptStat(0);
gStyle->SetPalette(1);

TCanvas *c2 = new TCanvas("c2","c2",50,50,800,600);
c2->SetFillColor(0);

  TLegend *leg = new TLegend(0.8,0.8,0.98,0.98);
  leg->AddEntry(h_lep_pt,"Uncut","P");
  leg->AddEntry(h_lep_pt_2,"Cut","P");

  c1->cd();
  h_lep_pt->Draw();
  h_lep_pt_2->Draw("same");
  //leg->Draw();

c1->Draw();

  c2->cd();
  h_lep12_IM->Draw();

c2->Draw();