#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <string>
#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include "stack.h"
#include "trigger.h"

using namespace std;
int main()
{
  TH1F * h1 = new TH1F("h1","h1",100,8000,11000);
  TFile * my_h1= new TFile ("TestL1.root", "RECREATE");
  fstream f_input;
  f_input.open("Dumpdata.txt",ios::in);
  fstream f_output;
  f_output.open("copia.txt",ios::out);
  
  if(f_input.fail())
    {
      cout<<"Il file di input non è stato aperto\n";
      exit(1);
    }
  if(f_output.fail())
    {
      cout<<"Il file di output non è stato aperto\n";
      exit(1);
    }
  
  Stack*  s= new Stack;
  Trigger *t = new Trigger;
  int j,i=0;
  int k=0;
  float r;
  float  a,b,c;
  for(i=0;f_input.eof()==0;)
    {     
      f_input >> a; 
      f_input >> b; 
      f_input >> c;
      s->push(a,b,c);
      i++;
    }
  f_output<< "Massinv\tPt1\tPt2\n";   
  float e=0;
  float f=0;
  cout<< "Inserisci i valori da usare nel trigger.\nPt1 > ";
  cin>>e;
  cout<<"Pt2 > ";
  cin>>f;
  for(j=0;j<i;j++)
    {
      t=s->pop();
      if(t->compare(t,e,f)!=0)
	{ 
	  f_output<<t->getMassinv()<<"\t"<<t->getPt1()<<"\t"<<t->getPt2()<<"\n";
	  k++;
	  h1->Fill(t->getMassinv());
	  
	}
    }
  
  r=(float)(k)/(float)(i)*100;
  cout<<"La percentuale di dati presi è "<<r<<"%\n";
  f_output<<"La percentuale di dati presi è "<<r<<"%\n";
  f_input.close();
  f_output.close();
  
  
  h1->Write();
  my_h1->Write();
  return 0;
}