//************************************************************************// // // // Copyright 2013 Bertram Kopf (bertram@ep1.rub.de) // // Julian Pychy (julian@ep1.rub.de) // // - Ruhr-Universität Bochum // // // // This file is part of Pawian. // // // // Pawian is free software: you can redistribute it and/or modify // // it under the terms of the GNU General Public License as published by // // the Free Software Foundation, either version 3 of the License, or // // (at your option) any later version. // // // // Pawian is distributed in the hope that it will be useful, // // but WITHOUT ANY WARRANTY; without even the implied warranty of // // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // // GNU General Public License for more details. // // // // You should have received a copy of the GNU General Public License // // along with Pawian. If not, see <http://www.gnu.org/licenses/>. // // // //************************************************************************// #include <math.h> #include <stdio.h> #include <boost/timer/timer.hpp> #include "MinFunctions/AbsFcn.hh" #include "Minuit2/MnUserParameters.h" #include "PwaUtils/GlobalEnv.hh" #include "PwaUtils/AbsLh.hh" #include "PwaUtils/NetworkServer.hh" #include "ErrLogger/ErrLogger.hh" using namespace ROOT::Minuit2; boost::timer::cpu_timer theTimer1; AbsFcn::AbsFcn() : _fcnCounter(0) , _currentResFileName("currentResult"+GlobalEnv::instance()->outputFileNameSuffix()+".dat") { } AbsFcn::~AbsFcn() { } double AbsFcn::Up() const { return .5; } void AbsFcn::printTimer() const{ theTimer1.stop(); boost::timer::cpu_times elapsed(theTimer1.elapsed()); if(elapsed.wall > 0){ Info << "Wall time: " << elapsed.wall / 1E9 << "s User: " << elapsed.user/1E9 << "s System: " << elapsed.system/1E9 << "s\n" << endmsg; } theTimer1.start(); } void AbsFcn::printFitParams(std::shared_ptr<AbsPawianParameters> par) const{ if ( _fcnCounter%1000 == 0) { par->print(std::cout, true); } } void AbsFcn::dumpFitParams(std::shared_ptr<AbsPawianParameters> par) const{ if ( _fcnCounter%200 == 0) { std::ofstream theStream (_currentResFileName.c_str()); par->print(theStream); } }