//==================================================================
//    Automatic script for doing analysis with ana.cc
//==================================================================
{
  //
  // 07/10/03 new instructions to load jet corrections libraries
  //
  //  gSystem->SetIncludePath(" -I./ -I./include -I$CDFSOFT2_DIR/include -D__STANDARD_CPLUSPLUS -DDEFECT_NO_EXCEPTIONS");

  // "k" = keep the generated shared library
  // "k" simply loads the shared library if the files weren't
  // modified from the last time
  // "f" = force compilation of the shared library every time
  // use "f" when you modify the headers that have names different
  // from the main .cc file

  #include<vector> 
 
  gSystem->Load("libTopFortran.so");
  gSystem->CompileMacro("TopFort.cc","kO");

  vector<double> a,b,c,d;

  vector<double> myResult;


  a.push_back( 23.0877 );
  a.push_back( 5.06046 );
  a.push_back( -7.71082 );
  a.push_back( 24.8617 );

  b.push_back( 32.5073 );
  b.push_back( 70.7031 );
  b.push_back( 16.1516 );
  b.push_back( 80.2933 );

  c.push_back( 42.2484 );
  c.push_back( 4.41263 );
  c.push_back( 14.1839 );
  c.push_back( 44.7837 );

  d.push_back( 17.6892 );
  d.push_back( 18.9048 );
  d.push_back( 60.5532 );
  d.push_back( 65.9271 );

  TopFort test;

  myResult = test.execute( a, b, c, d );

  cout<<"myResult.size() = "<<myResult.size()<<endl;

  for( unsigned i = 0; i < myResult.size(); ++i )
  {
        cout << "MyResult[" << i <<"] = " << myResult[i] << endl;
  }


  //TopFort test;

  //test.go();

}