LauFormulaPar.cc

Go to the documentation of this file.

// Copyright University of Warwick 2006 - 2013.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

// Authors:
// Thomas Latham
// John Back
// Paul Harrison

#include <iostream>
#include <iomanip>
#include <vector>
#include <cstdlib>

using std::cout;
using std::cerr;
using std::endl;
using std::vector;

#include "TRandom.h"
#include "TMessage.h"
#include "TSystem.h"

#include "LauFormulaPar.hh"
#include "LauParameter.hh"
#include "LauRandom.hh"

ClassImp(LauFormulaPar)


LauFormulaPar::LauFormulaPar(const TString& forName, const TString& formula, const std::vector<LauParameter*>& params) :
        name_(forName),
        formula_(forName,formula),
        paramVec_(params),
        dummy_(0),
        paramArray_(0)

{
        // Dummy array for TFormula
        dummy_ = new Double_t[1];

        // Array of input parameters
        Int_t pars = paramVec_.size();

        // Check length of vector matches number of parameter in the formula
        if (formula_.GetNpar()!=pars){
                cerr<<"ERROR in LauFormulaPar::evaluate : Number of parameters in the formula is : "<<formula_.GetNpar()<< " and the number of LauParameters is : "<<pars<<endl;
                gSystem->Exit(EXIT_FAILURE);
        }

        paramArray_ = new Double_t[pars];

        if (formula_.GetNdim()!=1){
                cerr<<"ERROR in LauFormulaPar::evaluate : Given formula of dimension: "<<formula_.GetNdim()<<" and not 1"<<endl;
                gSystem->Exit(EXIT_FAILURE);
        }

}

LauFormulaPar::~LauFormulaPar()
{
        delete[] dummy_;
        delete[] paramArray_;
}

Double_t LauFormulaPar::value() const
{

        //Assign vector values to array
        Int_t pars = paramVec_.size();
        
        for(Int_t i=0; i<pars; ++i){
                paramArray_[i] = paramVec_[i]->value();
        }

        return formula_.EvalPar(dummy_,paramArray_);
}