#include <iostream>
using std::cout;
using std::endl;

#include "PropSurfList.hpp"
#include "trfbase/SurfTest.h"
#include "trfbase/ETrack.h"
#include "trfbase/PropDirectedTest.h"

using namespace trf;

int main () { 

  cout << endl;
  cout << endl;
  cout << endl;
  cout << "Performing component test of PropSurfList" << endl;

  ///////////////////// initial setup
  double bfield = 0.00;
  PropagatorPtr prop1( new PropTest() );
  PropagatorPtr prop2( new PropTest() );

//make a surface
  SurfacePtr mysrf( new SurfTest(123) );

//make an ETrack
  TrackVector myvec;
  myvec(0)=1.0;myvec(0)=2.0;myvec(0)=3.0;myvec(0)=4.0;myvec(0)=5.0;

  TrackError myerr;
  myerr(0,0)=0.1; myerr(0,0)=0.1; myerr(0,1)=0.01;
  myerr(0,2)=0.02;myerr(0,3)=0.03;myerr(0,4)=0.04;
  myerr(1,1)=0.2;myerr(1,2)=0.02;myerr(1,3)=0.03;
  myerr(1,4)=0.04;myerr(2,2)=0.3;myerr(2,3)=0.03;
  myerr(2,4)=0.04;myerr(3,3)=0.4;myerr(3,4)=0.04;myerr(4,4)=0.5;

  TrackSurfaceDirection mydir=TSD_FORWARD;

  ETrack mytrack(mysrf,myvec,myerr,mydir);
  cout << mytrack << endl;

  ////////////

  PropSurfList a1;
  int i;
  i = a1.getLength();
  cout << "length of PropSurfList is " << i << endl;

  a1.push_back(prop1,mysrf);
  cout << "length is now " << a1.getLength() << endl;

  a1.push_back(prop2,mysrf);

  cout << "length is now " << a1.getLength() << endl;


return 0; }