#include"TopFort.h"
#include<fstream>
#include<sstream>
#include<iostream>

using namespace std;

TopFort::~TopFort()
{

}


TopFort::TopFort()
{

}

vector<double>
TopFort::execute( vector<double> p_ele, vector<double> b0ele,
                  vector<double> p_muo, vector<double> b0muo )
{

	float p_ele_local[4];
	float p_muo_local[4];

	float b0ele_local[4];
        float b0muo_local[4];

	float result[200];

	vector<double> answer;

	for( unsigned i = 0; i < 200; ++i )
	{
		result[i] = 0.;
	}

	

	if( ( p_ele.size() != 4 ) || ( b0ele.size() != 4 ) || 
            ( p_muo.size() != 4 ) || ( b0muo.size() != 4 ) )
	{
		cout << "ERROR!  Input 4 vectors not valid!" << endl;
		return answer;
	}


        for( unsigned i = 0; i < 4; ++i )
        {
		p_ele_local[i] = p_ele[i];
		p_muo_local[i] = p_muo[i];
		b0ele_local[i] = b0ele[i];
		b0muo_local[i] = b0muo[i];	
        }


	dilep_( p_ele_local, b0ele_local, p_muo_local, b0muo_local, result );


	for( unsigned i = 0; i < 200; ++i )
	{
		answer.push_back( result[i] ); 
	}

	return answer;

}