#ifndef GUARD_TDANA
#define GUARD_TDANA
//#include "TDHistSet.hh"
//#include "Histtest.hh"
#include "TDSuperAna.hh"
#include <ostream>
//#include <ofstream>
#include "TDExcelPlot.hh"
#include "TDCanvMan.hh"
#include <string>

std::ostream& TexTable(std::ostream& os, std::vector<std::vector<std::string> > thetable, char * caption)
{
  int cols=thetable.size();
  int rows=(thetable[0]).size();
  std::string colstring="|l||";
  for(int c=0;c!=cols-1;++c)colstring+="c|";
  
  os<<" \\begin{table}[H] \\begin{sideways}"<<std::endl;
  os<<"\\tiny"<<std::endl;
  os<<"\\begin{tabular}{"<<colstring<<"}"<<std::endl;
  
  os<<" \\hline \\hline ";
  for(int row=0; row!=rows; ++row)
  {
    if(row==1)
    os<<" \\hline ";
    for(std::vector<std::vector<std::string> >::iterator col=thetable.begin(); col!=thetable.end(); ++col)
    {
      if(col!=thetable.begin())os<<" & ";
      os<<(*col)[row];
    }
    os<<" \\\\ "<<std::endl;
  }
  os<<" \\hline \\hline "<<std::endl;
  os<<"\\end{tabular} \\end{sideways}"<<std::endl;
  os<<"\\caption{"<<caption<<"}"<<std::endl;
  os<<"\\label{table:XXXYYY}"<<std::endl;
  os<<"\\end{table}"<<std::endl;
  os<<std::endl;


  return os;
}

std::ostream& setupTex(std::ostream& os)
{
  os<<"\\documentclass[12pt]{article}"<<std::endl;
  os<<"\\usepackage{epsfig}"<<std::endl;
  os<<"\\usepackage{verbatim}"<<std::endl;
  os<<"\\usepackage{rotating}"<<std::endl;
  os<<"\\def\\rightboxfour#1#2#3#4{\\bigskip\\rightline{\\vbox{\\hbox{#1}"<<std::endl;
  os<<"        \\hbox{#2}\\hbox{#3}\\hbox{#4}}}\\bigskip}"<<std::endl;
  os<<"        \\parskip12pt plus 1pt minus 1pt"<<std::endl;
  os<<"        \\topsep0pt plus 1pt"<<std::endl;
  os<<"        \\setlength{\\textheight}{22.0cm}"<<std::endl;
  os<<"        \\setlength{\\topmargin}{0.0cm}"<<std::endl;
  os<<"        \\setlength{\\textwidth}{15.0cm}"<<std::endl;
  os<<"        \\setlength{\\oddsidemargin}{-0.5cm}"<<std::endl;
  os<<"        \\setlength{\\topsep}{1mm}"<<std::endl;
  os<<"        \\renewcommand{\\topfraction}{1.0}"<<std::endl;
  os<<"        \\renewcommand{\\bottomfraction}{1.0}"<<std::endl;
  os<<"        \\renewcommand{\\textfraction}{0.2}"<<std::endl;
  os<<"        \\renewcommand{\\floatpagefraction}{0.5}"<<std::endl;
  os<<"        \\pretolerance=10000"<<std::endl;
  os<<""<<std::endl;
  os<<""<<std::endl;
  os<<""<<std::endl;
  os<<"        \\begin{document}"<<std::endl;
  return os;
}

void runSuperAna(int diltype, std::ofstream& AnaTextFile, std::ofstream& ExcelTextFile, std::ofstream& ExcelDataFile, std::ofstream& ExcelMassFile, std::ofstream& CutTableFile,  std::ofstream& CutTableTexFile, int mass)
{
  std::ofstream AcceptanceTexfile;
  //char buffer[500];
  //sprintf(buffer,"results/Accept_%d.tex", (int)higgsmass);
  AcceptanceTexfile.open("results/Accept.tex");
  setupTex(AcceptanceTexfile);
  std::ofstream LikelihoodTexfile;
  LikelihoodTexfile.open("results/Likelihood.tex");
  setupTex(LikelihoodTexfile);

  //std::ofstream AllAcceptanceTexfile;
  //  AllAcceptanceTexfile.open("results/AllAccept.tex");
  //    setupTex(AllAcceptanceTexfile);
  
  std::vector<std::vector<std::string> > higgsacceptances;  //TD vectors to store table information until all the mass points are done.
  std::vector<std::vector<std::string> > higgsacceptances2;
  std::vector<std::vector<std::string> > likelihoodvector;//TD vectors to store table information until all the mass points are done.
  bool start=true;
  
//  for(int higgsmass=110; higgsmass!=210; higgsmass+=10) //TD to loop over mass points
//  {
  //for(int weightype=0; weightype<3; ++weightype) //TD to loop over weights - check that this still works, and that they are properly labled.
  //{
    //
    //int higgsmass=160; //td normally 160
    int higgsmass=mass;
    int weightype=2;
    TDSuperAna *SA=new TDSuperAna(false,"/data/nglas08/b/stdenis/hww_1fb_results/", diltype, AnaTextFile, weightype, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile, CutTableTexFile,  higgsmass);
    //SA->mirclikelihood(true);  // to test whether getting the limit from loaded histograms works
    
    SA->RunSubAnas();
    SA->XSroutine2(); //not clear, I know; corresponds to 20-20 et cuts
    SA->XSroutine3(); //not clear, I know; corresponds to 20-10 et cuts
    SA->XSroutine4(); //not clear, I know; corresponds to 20-10 BUT NOT 20-20 et cuts ( abit useless due to lack of stats.)
                                                       ////SA->XSroutine(); //old xs method (uses histograms and was more difficult to use.)
    SA->GetHistSets()->SetInnateColors();
    SA->DrawUnscaledPlots("HIST BOX");
    SA->ScaleHistSets();
   
    SA->GetSRB("MVA_TMlpANN_fine_ee", "MVA_TMlpANN_fine_ee_srb");
    SA->GetSRB("MVA_TMlpANN_fine_emu", "MVA_TMlpANN_fine_emu_srb");
    SA->GetSRB("MVA_TMlpANN_fine_mumu", "MVA_TMlpANN_fine_mumu_srb");
    SA->GetSRDYB("MVA_TMlpANN_fine_ee", "MVA_TMlpANN_fine_ee_srdyb");
    SA->GetSRDYB("MVA_TMlpANN_fine_emu", "MVA_TMlpANN_fine_emu_srdyb");
    SA->GetSRDYB("MVA_TMlpANN_fine_mumu", "MVA_TMlpANN_fine_mumu_srdyb");
    //SA->GetUCSD("MVA2_TMlpANN_cut1","MVA2_TMlpANN_cut1_ucsd");
    
    SA->SavePlots();
    
    /*SA->LoadPlots();
      SA->GetSRB("MVA_TMlpANN_fine_ee", "MVA_TMlpANN_fine_ee_srb");
      SA->GetSRB("MVA_TMlpANN_fine_emu", "MVA_TMlpANN_fine_emu_srb");
      SA->GetSRB("MVA_TMlpANN_fine_mumu", "MVA_TMlpANN_fine_mumu_srb");
      SA->GetSRDYB("MVA_TMlpANN_fine_ee", "MVA_TMlpANN_fine_ee_srdyb");
      SA->GetSRDYB("MVA_TMlpANN_fine_emu", "MVA_TMlpANN_fine_emu_srdyb");
      SA->GetSRDYB("MVA_TMlpANN_fine_mumu", "MVA_TMlpANN_fine_mumu_srdyb");
    SA->GetHistSets()->SetInnateColors();*/
    //    SA->ScaleHistSets();
   
    SA->DrawScaledPlots("HIST BOX");
    
    //SA->DrawPlots("HIST");
//    SA->DrawRickPlots();
    
    //SA->TomLikehisdata();

    //SA->NN();             //TD at present just saves the trees for TMVA.
    //SA->TMVAnalysis();    //this does not work - it was an attempt to run the TMVA training in the program.
    //SA->JetNet();         //have not checked whether this works for a while
    //SA->likelihood(true);//TD simulation 7708 code
    //SA->likelihood(false);//TD Data 7708 code 
    //SA->TomLikelihood7708();                      //uses tom's method with 7708 systematics
    //SA->TomLikelihoodTop();                       //uses tom's method with his systematics
    /*SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_0");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_1");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_2");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_3");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_4");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_5");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_6");*/
    //SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_7");
    /*SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_8");
    SA->TomLikelihoodTop("MVA2_CFMlpANN_cut2_9");*/
    //SA->TomLikelihoodTop("dilDR");
    //SA->TomLikelihoodTop("dilDphi_TMVA");
//    SA->TomLikelihoodTop("MVA2_TMlpANN_cut1", false);
    //SA->TomLikelihoodTop2D("MVA2_TMlpANN_2D", true); //"true" indicates that the lengthy frequentist method is to be used too.
    //SA->TomLikelihoodTop2D("MVA2_TMlpANN_2D_binexp");
    //SA->TomLikelihoodTop2D("MVA2_TMlpANN_2D_small", true);
    
    SA->TomLikelihoodemu_justtheerrorplots();   //draws the error plots used in TomLikelihoodemu.
    SA->TomLikelihoodemu_NN_justtheerrorplots();
    //SA->TomLikelihoodemu();
//    SA->TomLikelihoodDilDphi(false);
//    SA->TMLP2_justtheerrorplots();
//    SA->TomLikelihoodNNnp(false);
//    SA->TomLikelihoodNNshape(false); //isthis still available?
    SA->TomLikelihood_NN_shape_emu();
//    SA->TomLikelihood_NN_shape_emu_wgamma();
   	     
    //SA->TomLikelihoodDetailed(10);//detailed is a 'per-dilepton' analysis (depricated due to empty histograms (low stats))
    //SA->TomLikelihoodDetailed(20);//the number is the binning of the histograms - this was a test at one point.
    SA->WWXS(); // TD estimates the WW XS. gives a ratio of the integrals of the histograms
    SA->GetMaxEtCutError();
    SA->GetNNDYCutError();
    SA->GetMetSigCutEff();
 
    SA->PrintAcceptances(AcceptanceTexfile);
    SA->PrintAcceptances2(AcceptanceTexfile);
    SA->PrintAcceptances3(AcceptanceTexfile);
    SA->PrintNNAcceptances3(AcceptanceTexfile);
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "MVA2_TMlpANN_cut1_samesign", "MVA2_TMlpANN_cut1_samesign_ee", "MVA2_TMlpANN_cut1_samesign_emu", "MVA2_TMlpANN_cut1_samesign_mumu");
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "WWmetmag", "WWmetmag_ee", "WWmetmag_emu","WWmetmag_mumu");
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "MVA2_TMlpANN_cut1_samesign_cmupcmup", "MVA2_TMlpANN_cut1_samesign_cmupcmup", "MVA2_TMlpANN_cut1_samesign_cmupcmup", "MVA2_TMlpANN_cut1_samesign_cmupcmup");
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "MVA2_TMlpANN_cut1_samesign_tcetce", "MVA2_TMlpANN_cut1_samesign_tcetce", "MVA2_TMlpANN_cut1_samesign_tcetce", "MVA2_TMlpANN_cut1_samesign_tcetce");
    
   /* SA->PrintAcceptances2Generic(AcceptanceTexfile, "HighMllControlNNDYCut_dildphi", "HighMllControlNNDYCut_dildphi_ee", "HighMllControlNNDYCut_dildphi_emu", "HighMllControlNNDYCut_dildphi_mumu");  //hists1
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "ZPeakNNInputVar_dildphi", "ZPeakNNInputVar_dildphi_ee", "ZPeakNNInputVar_dildphi_emu", "ZPeakNNInputVar_dildphi_mumu"); //hists1
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "ZPeakNNDYCut_dildphi", "ZPeakNNDYCut_dildphi_ee", "ZPeakNNDYCut_dildphi_emu", "ZPeakNNDYCut_dildphi_mumu");   //hists1
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "HighMllControlNNInputVar_dildphi", "HighMllControlNNInputVar_dildphi_ee", "HighMllControlNNInputVar_dildphi_emu", "HighMllControlNNInputVar_dildphi_mumu"); //hists1
    */
   /* SA->PrintAcceptances2Generic(AcceptanceTexfile, "NNInputVar_dildphi", "NNInputVar_dildphi_ee", "NNInputVar_dildphi_emu", "NNInputVar_dildphi_mumu");//hists2
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "NNDYCut_dildphi", "NNDYCut_dildphi_ee", "NNDYCut_dildphi_emu", "NNDYCut_dildphi_mumu");//hists2
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LepEtControlNNInputVar_dildphi", "LepEtControlNNInputVar_dildphi_ee", "LepEtControlNNInputVar_dildphi_emu", "LepEtControlNNInputVar_dildphi_mumu"); //hists2
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LepEtControlNNDYCut_dildphi", "LepEtControlNNDYCut_dildphi_ee", "LepEtControlNNDYCut_dildphi_emu", "LepEtControlNNDYCut_dildphi_mumu"); //hists2
    */
/*    SA->PrintAcceptances2Generic(AcceptanceTexfile, "WWNNlowCut_dildphi", "WWNNlowCut_dildphi_ee", "WWNNlowCut_dildphi_emu", "WWNNlowCut_dildphi_mumu"); //hists3
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "WWNNmidCut_dildphi", "WWNNmidCut_dildphi_ee", "WWNNmidCut_dildphi_emu", "WWNNmidCut_dildphi_mumu"); //hists3
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "WWNNhighCut_dildphi", "WWNNhighCut_dildphi_ee", "WWNNhighCut_dildphi_emu", "WWNNhighCut_dildphi_mumu"); //hists3
*/
/*    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LowMllControlNNInputVar_dildphi", "LowMllControlNNInputVar_dildphi_ee", "LowMllControlNNInputVar_dildphi_emu", "LowMllControlNNInputVar_dildphi_mumu");//hists4
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LowMllControlNNDYCut_dildphi", "LowMllControlNNDYCut_dildphi_ee", "LowMllControlNNDYCut_dildphi_emu", "LowMllControlNNDYCut_dildphi_mumu");//hists4
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "NNDYCutMid_dildphi", "NNDYCutMid_dildphi_ee", "NNDYCutMid_dildphi_emu", "NNDYCutMid_dildphi_mumu");//hists4
*/
/*    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LMetHMetSNNInputVar_dildphi", "LMetHMetSNNInputVar_dildphi_ee", "LMetHMetSNNInputVar_dildphi_emu", "LMetHMetSNNInputVar_dildphi_mumu");//hists5
    SA->PrintAcceptances2Generic(AcceptanceTexfile, "LMetSHMetNNInputVar_dildphi", "LMetSHMetNNInputVar_dildphi_ee", "LMetSHMetNNInputVar_dildphi_emu", "LMetSHMetNNInputVar_dildphi_mumu");//hists5
*/
    SA->doKS();
    
    if(start)
    {
      higgsacceptances.push_back(SA->getDatasets());//TD this fills in the first column of the vector once, they being the row lables.
      higgsacceptances2.push_back(SA->getDatasets());
      likelihoodvector.push_back(SA->GetLikelihoodTypes());  //TD this fills in the first column of the vector, they being the row lables.
      start=false;
    }
    likelihoodvector.push_back(SA->GetLikelihoodValues());//TD these vectors build the columns up of the tables they suppy
    std::cout<<"likelihoodvector.size()"<<likelihoodvector.size()<<std::endl;
    higgsacceptances.push_back(SA->getAcceptanceVector());//TD these vectors build the columns up of the tables they suppy
    higgsacceptances2.push_back(SA->getAcceptanceVector2());
    //SA->TomLikelihoodemu();
    delete SA;
//  }
  //}
  std::cout<<" end likelihoodvector.size()"<<likelihoodvector.size()<<std::endl; 
  TexTable(AcceptanceTexfile, higgsacceptances,"Acceptances for all masses"); //TD write the table
  TexTable(AcceptanceTexfile, higgsacceptances2,"Acceptances for all masses");
  AcceptanceTexfile<<"        \\end{document}"<<std::endl;  //TD finish off the file
  AcceptanceTexfile.close();
  TexTable(LikelihoodTexfile, likelihoodvector,"likelihoods");  // TDwrite the table
  LikelihoodTexfile<<"        \\end{document}"<<std::endl; //TD finish off the file
  LikelihoodTexfile.close();
  
  //AllAcceptanceTexfile<<"        \\end{document}"<<std::endl;
  //AllAcceptanceTexfile.close();
    
  ExcelTextFile<<std::endl;
}


int main(int argc, const char * argv[])
{
   int mass=std::atoi(argv[1]);
   std::cout<<mass<<std::endl;
   if(mass<110 && mass>200) mass=160;
//   return 0;
    std::ofstream AnaTextFile;
    std::ofstream ExcelTextFile;
    std::ofstream ExcelDataFile;
    std::ofstream ExcelMassFile;
    std::ofstream CutTableFile;
    std::ofstream CutTableTexFile;
    AnaTextFile.open("results/AnaTextFile.txt");
    ExcelTextFile.open("results/ExcelTextFile.txt");
    ExcelDataFile.open("results/ExcelDataFile.txt");
    ExcelMassFile.open("results/ExcelMassFile.txt");
    CutTableFile.open("results/CutTableFile.txt");
    CutTableTexFile.open("results/CutTableTexFile.tex");
    setupTex(CutTableTexFile);

    //for(int diltype=0; diltype<31; ++diltype) runSuperAna(diltype, AnaTextFile, ExcelTextFile, ExcelDataFile);
  
runSuperAna(-1, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile, CutTableTexFile, mass);
    //runSuperAna(1, AnaTextFile, ExcelTextFile, ExcelDataFile);
   
    /*runSuperAna(0, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(1, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    :runSuperAna(17, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(18, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(19, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(20, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(21, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(22, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(23, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(24, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(25, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(29, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(30, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);
    runSuperAna(31, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile);*/
//runSuperAna(32, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile, CutTableTexFile);
//runSuperAna(33, AnaTextFile, ExcelTextFile, ExcelDataFile, ExcelMassFile, CutTableFile, CutTableTexFile);
  
  AnaTextFile.close();
  ExcelTextFile.close();
  ExcelDataFile.close();
  ExcelMassFile.close();
  CutTableFile.close();


CutTableTexFile<<"        \\end{document}"<<std::endl;
  CutTableTexFile.close();

 /* TDExcelPlot excelplot;
  std::cout<<"here"<<std::endl;
  excelplot.ReadFile("results/EDFmuon_0_1st4_combos.txt");
  std::cout<<"here"<<std::endl;
  TDCanvMan canvman("results/EDFmuon_0_1st4_combos.ps", "Ana_Plots", 1, 2);
  std::cout<<"here"<<std::endl;
  excelplot.plot(&canvman);*/
  return 0;
}


#endif