#include "JppPdfTest.hh"
#include "JppShowerPdf.hh"
#include "JppShowerFit.hh"

using namespace JPP;

int main(){

  JppShowerPdf* my_pdf = new JppShowerPdf( "/sps/km3net/users/jseneca/jpp/data/J14p.dat",
					   "/sps/km3net/users/jseneca/jpp/data/J13p.dat",
					   true,
					   1.0 );

  //JppShowerFit<JppShowerNpe>* my_fit_npe = new JppShowerFit<JppShowerNpe>( "name", my_npe );
  //my_fit_npe->_init();
  //JppShowerFit<JppShowerPdf> my_fit_pdf = JppShowerFit<JppShowerPdf>( "name", my_pdf );
  //my_fit_pdf._init();

  cout << "Created all objects" << endl;

  int N = 0;

  cout << "Call npe from pdf class:" << endl;
  cout << my_pdf->npe<1>( N, 1.0, 50.0, 0.1, 0.1, 0.1, 0.1) << endl;

  cout << "Call npe from base class:" << endl;
  cout << my_pdf->npe_tot<1>( N, 0.1, 0.1, 0.1, 0.1, 0.1 ) << endl;

  cout << "npe from base worked!" << endl;



  cout << "Call npe from pdf class:" << endl;
  cout << my_pdf->dPfirst_dt<1>( N, 6.0, 50.0, 0.1, 0.1, 0.1, 0.1) << endl;

  cout << "Press enter to test JppShowerNpe vs. JppShowerPdf" << endl;
  cin.ignore();
  JppShowerNpe* my_npe = new JppShowerNpe( *my_pdf );

  //compare NPE object to npe from PDF object
  cout << "Npe object:        " << my_pdf->npe_tot<0>( 0, 1.0, 50.0, 0.1, 0.1, 0.1 ) << endl;
  cout << "Npe in pdf object: " << my_npe->npe_tot<0>( 0, 1.0, 50.0, 0.1, 0.1, 0.1 ) << endl;

  my_npe->test_npef( *my_pdf );

//my_pdf.test_params( trk );
//my_npe.test_params( trk );
}