#include <sstream>
#include <string>
#include <iostream.h>
#include <fstream.h>
#include <stdlib.h>  
#include <time.h>
#include <math.h>
//ROOT-includes
#include <TROOT.h>
#include <TApplication.h>
#include <TH1.h>
#include <TFile.h> 
#include <TTree.h>
#include <TBranch.h>
#include <TCanvas.h>
#include <TLegend.h> 
#include <TLeaf.h> 
#include <TClonesArray.h>
#include <TGraph.h>
#include <TAxis.h>
//CLHEP-includes
#include "TNT2pulseshape.h"

TROOT hej("hej","hej"); 

int main(int argc, char** argv)

{

  Double_t fit_param_middle;
  Double_t fit_param_middle_err;
  Double_t fit_param_left;
  Double_t fit_param_left_err;
  Double_t fit_param_right;
  Double_t fit_param_right_err;
  
  Double_t amplitude_difference;
    
  Int_t event_counter = 0.0;
  
  char infile[80];
  
  int argcROOT = 1; 
  char* argvROOT = argv[0]; 
  TApplication theApp("App", &argcROOT, &argvROOT);

  for(int argc_up=1;argc_up<argc;argc_up++)
  {
  
    strcpy(infile,argv[argc_up]);
    TNT2pulseshape ps(infile,"TNT2");
    
    cout << "Processing file " << infile << "\n";
    
    while(ps.getNextEvent())
    {
      
      TF1 *pol0_fit_middle = new TF1("pol0_fit_middle","pol0",2.5,5);
      ps.getPS(3)->Fit("pol0_fit_middle","QRM");
      
      TF1 *pol0_fit_left = new TF1("pol0_fit_left","pol0",0,2.5);
      ps.getPS(3)->Fit("pol0_fit_left","QRM");
          
      TF1 *pol0_fit_right = new TF1("pol0_fit_right","pol0",5,20);
      ps.getPS(3)->Fit("pol0_fit_right","QRM");
      
    }
  
  }

  return 0;
}