#include <stdio.h>
#include <stdbool.h>
#include <stdlib.h>
#include <iostream>
#include <fstream>
#include <string>
#include <vector>

struct resultsStruc
{
    int bin;
    double epsilon;
    double branching_ratio;
};

struct dataStruc
{
    int b;
    int e;
    int g;
};


void efficiency()
{

ifstream filebin;  
 ifstream filedata; 
 ofstream results;

int bin;
int b,e,g;
double branching_ratio,epsilon;
int temp;

	//std::vector<int> bin;
    std::vector<resultsStruc> data_results;
	std::vector<dataStruc> data;
	filebin.open("bin.txt"); 
	filedata.open("EnergySpectrumResults.txt"); 	
		results.open("Results.txt",  ios::out);
	
	  while (filedata  >> b >> e >> g ){
	 	data.push_back(dataStruc{b,e,g});
	 }
	 

    while (filebin >> bin ){
        
        epsilon=data.at(bin).e;
        
        if(0<= bin<= data.size()-1)
            epsilon = epsilon - 0.5*data.at(bin-1).e - 0.5*data.at(bin+1).e;
        
        epsilon=epsilon/1.e6;
        
        branching_ratio=data.at(bin).g/1.e6;
        
        data_results.push_back({bin, epsilon, branching_ratio});
        
        
    }
    
    
    
    
    results << "**************************************** "<< endl;
   	results << "Energy" << "\t"<< "Efficiency" << "\t"<< "BR" << endl;
    for(int i=0;i<data_results.size();i++)
        results <<data_results.at(i).bin<< "\t "<< data_results.at(i).epsilon<< "\t "<< data_results.at(i).branching_ratio<< endl;
   	 		
	filebin.close(); 
	filedata.close(); 
	results.close(); 
	 
}