#include <iostream>
#include <stdio.h>
#include<locale.h>
#include "TCanvas.h"
#include "TFitResult.h"
#include "TDirectory.h"
//#include "TAxis.h"
//#include "TH1D.h"
//#include "TAxis.h"
//#include "TStyle.h"  //for gStyle decleration
//#include <fstream>
//#include "TF1.h"
//#include "TFormula.h"
#include "TMath.h"
#include "TROOT.h"
#include <vector>
#include"TString.h"
#include <cmath>
#define PI 3.14159265
//#define SIZEJ  60    // note it must be:
//#define SIZEK  240   //   4*SIZEJ = SIZEK > 4*jmax is a good choice

using namespace std;

int roottest(){
    
    //delete gROOT->FindObject("h1");
    //delete gROOT->FindObject("c1");
    //initializing a vector
    std::vector<float>T_energyLIST{20., 19., 18., 17., 16.};
    for (unsigned int i = 0; i < T_energyLIST.size(); i++ ){
        cout << "T_energy values: " << T_energyLIST[i] << endl;
    }
    float z_particle, M_particle, rho_material, A_material, Z_material, I_material;
    float area= TMath::Power(10.,-28.); //m^2, 4.pi.r^2
    float mc2= 0.511; //MeV
    float NA= 6.022*TMath::Power(10.,23.); //atoms / mole
    float c2= 931.5 ; // MeV/u
    float N_material =  (rho_material * NA )/ A_material ;
    //Arrays
    std::vector<float> gamma;
    std::vector<float> Beta;
    std::vector<float> V_particle;
    std::vector<float> AverageEnergyLossPerUnitDistance_O;
    std::vector<float> AverageEnergyLossPerUnitDistance_I;
    std::vector<float> result;
    /*
    for (unsigned int i = 0; i < T_energyLIST.size(); i++ ){
     gamma.push_back( (T_energyLIST[i] / (M_particle * c2) ) + 1. ) ;
     Beta.push_back( sqrt(1.- (1./(gamma[i]* gamma[i]) ) )  );
     V_particle.push_back( Beta[i] * sqrt(c2)  );
     AverageEnergyLossPerUnitDistance_O.push_back( area * TMath::Power(z_particle, 2.)* mc2 * N_material * Z_material * (1./ TMath::Power(Beta[i], 2.)) );
     AverageEnergyLossPerUnitDistance_I.push_back( TMath::Log(2.* mc2* TMath::Power(Beta[i], 2.)* TMath::Power(gamma[i], 2.)* (1./ I_material) )- TMath::Power(Beta[i], 2.)   );
     result.push_back( AverageEnergyLossPerUnitDistance_O[i] * AverageEnergyLossPerUnitDistance_I[i]  );
    }
     */
    cout<<"done"<<endl;
    return 0;
    
    
}