#include "Experiment.h"
#include <cmath>
#include <vector>
#include <iterator>
#include <armadillo>
using namespace arma;
using namespace std; 

Experiment::Experiment(vector< int > ipar_) : ipar(ipar_) //: r(new TRandom3)
{   
    par_n = 10;
}

double Experiment::expectation_value(double *t, double *par){
    parameters(par);
    
    return 0.5;
}

void ModelA::parameters(double *par){
    for(int i=0; i<par_n; i++){
      par_internal[i] = par[ipar[i]];
    }
}

void ModelB::parameters(double *par){
    double norm_meas=0.0;
    for(int i=0; i<par_n; i++){
      par_internal[i] = par[ipar[i]];
    }
    for(int i=0; i<3; i++){
      norm_meas += pow(par_internal[i+3],2);
    }
    norm_meas = sqrt(norm_meas);
    for(int i=0; i<3; i++){
      par_internal[i+3] = par_internal[i+3]/norm_meas;
    }
}