#include <math.h> #include <stdio.h> #include <boost/timer/timer.hpp> #include "Minuit2/MnUserParameters.h" #include "PwaUtils/PwaFcnBase.hh" #include "PwaUtils/AbsLh.hh" #include "PwaUtils/NetworkServer.hh" #include "ErrLogger/ErrLogger.hh" using namespace ROOT::Minuit2; boost::timer::cpu_timer theTimer; PwaFcnBase::PwaFcnBase(boost::shared_ptr<AbsLh> absLh, boost::shared_ptr<FitParamsBase> fitParamsBase, std::string suffix) : _absLhPtr(absLh) , _fitParamsBasePtr(fitParamsBase) , _fcnCounter(0) , _currentResFileName("currentResult"+suffix+".dat") , _serverMode(false) { if (0==_absLhPtr) { Alert << "AbsLh* _absLhPtr pointer is 0 !!!!" << endmsg; exit(1); } _absLhPtr->getDefaultParams(_defaultFitValParms, _defaultFitErrParms); } PwaFcnBase::~PwaFcnBase() { } double PwaFcnBase::operator()(const std::vector<double>& par) const { double result=0; fitParams theFitParmValTmp=_defaultFitValParms; _fitParamsBasePtr->getFitParamVal(par, theFitParmValTmp); if(_serverMode){ LHData theLHData; _networkServerPtr->BroadcastParams(par); if(!_networkServerPtr->WaitForLH(theLHData.logLH_data, theLHData.weightSum, theLHData.LH_mc)) result = 0; else result = _absLhPtr->mergeLogLhData(theLHData); } else{ result = _absLhPtr->calcLogLh(theFitParmValTmp); } _fcnCounter++; if(_fcnCounter%20 == 0){ theTimer.stop(); boost::timer::cpu_times elapsed(theTimer.elapsed()); if(elapsed.wall > 0){ Info << "Wall time: " << elapsed.wall / 1E9 << "s User: " << elapsed.user/1E9 << "s System: " << elapsed.system/1E9 << "s\n" << endmsg; } theTimer.start(); } if ( _fcnCounter%100 == 0) { _fitParamsBasePtr->printParams(theFitParmValTmp); } if ( _fcnCounter%200 == 0) { std::ofstream theStream (_currentResFileName.c_str()); _fitParamsBasePtr->dumpParams(theStream, theFitParmValTmp, (fitParams&)_defaultFitErrParms); } return result; } double PwaFcnBase::Up() const { return .5; } void PwaFcnBase::SetServerMode(boost::shared_ptr<NetworkServer> networkServerPtr){ _serverMode = true; _networkServerPtr = networkServerPtr; }