gROOT->Reset();

#include "TMatrix.h"
#include "TSQLResult.h"
#include "TSystem.h"
#include "TROOT.h"
#include "TChain.h"
#include "TTree.h"

#include <iostream>
#include <fstream>
#include <vector>
#include <string>
#include <stdio.h>

#include <stdlib.h>

using namespace std;
bool RunSusy3LepCutFlow8TeV(){
    
    
    string filename = "";
    bool applyDPhicut = false, applyDRcut=false;
    string met_cut ="",mt_cut="";
    //string signalRegion = "";
    filename = string(getenv ("FILENAME_AKXNZ"));
    //signalRegion = string(getenv ("SIGNAL_REGION_AKXNZ"));
    if (string(getenv ("DPhi_AKXNZ"))=="TRUE") applyDPhicut = true;
    else if (string(getenv ("DR_AKXNZ"))=="TRUE") applyDRcut = true;

    if (string(getenv ("MET_AKXNZ"))!="") met_cut = string(getenv ("MET_AKXNZ"));
    if (string(getenv ("MT_AKXNZ"))!="") mt_cut = string(getenv ("MT_AKXNZ"));

    
    cout<<"filename: "<<filename<<endl;

    gROOT->ProcessLine (".x $ROOTCOREDIR/scripts/load_packages.C");
    
    
    gSystem->CompileMacro("Susy3LepCutflow8TeV.C","gk");
    
    cout<<"Susy3LepCutflow8TeV.C compiled"<<endl;
    TChain* chain = new TChain("chain","chain");
    
    chain->Add("root://eosatlas//eos/atlas/user/a/akourkou/NewPtCuts/user.akourkou.164323.Herwigpp_wA_noslep_lepW_50.merge.NTUP_SUSY.e1461_s1499_s1504_r3658_r3549_p1512_NewPtMll.root/susy3lep");

    Susy3LepCutflow8TeV t(chain, 1, 0, 1, 0, "output.root",filename,applyDPhicut,applyDRcut,met_cut,mt_cut);
    
    t.Loop();
    return 0;
}