#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> > 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->TomLikehisdata();
    
    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->SavePlots();
    
    //SA->LoadPlots();
    //SA->GetHistSets()->SetInnateColors();
    //SA->ScaleHistSets();
    
    SA->DrawScaledPlots("HIST BOX");
    
    //SA->DrawPlots("HIST");
    SA->DrawRickPlots();

    SA->NN();  //TD at present just saves the trees for TMVA.
    //SA->TMVAnalysis();
    //SA->JetNet();
    SA->likelihood(true);//TD simulation
    SA->likelihood(false);//TD Data
    //SA->TomLikelihood7708();
    SA->TomLikelihoodTop();
    /*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->TomLikelihoodemu();
    //SA->TomLikelihoodDetailed(10);//detailed is a 'per-dilepton' analysis (depricated due to empty histograms)
    //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->PrintAcceptances(AcceptanceTexfile);
    if(start)
    {
      higgsacceptances.push_back(SA->getDatasets());//TD this fills in the first column of the vector once, they being the row lables.
      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

    delete SA;
//  }
  //}
  std::cout<<" end likelihoodvector.size()"<<likelihoodvector.size()<<std::endl; 
  TexTable(AcceptanceTexfile, higgsacceptances,"Acceptances for all masses"); //TD write the table
  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