Changeset View
Changeset View
Standalone View
Standalone View
src/LauAbsFitModel.cc
| Show First 20 Lines • Show All 560 Lines • ▼ Show 20 Lines | |||||
| void LauAbsFitModel::fitExpt() | void LauAbsFitModel::fitExpt() | ||||
| { | { | ||||
| // Routine to perform the actual fit for the given experiment | // Routine to perform the actual fit for the given experiment | ||||
| // Update initial fit parameters if required (e.g. if using random numbers). | // Update initial fit parameters if required (e.g. if using random numbers). | ||||
| this->checkInitFitParams(); | this->checkInitFitParams(); | ||||
| // Initialise the fitter | // Initialise the fitter | ||||
| LauFitter::fitter()->useAsymmFitErrors( this->useAsymmFitErrors() ); | LauFitter::fitter().useAsymmFitErrors( this->useAsymmFitErrors() ); | ||||
| LauFitter::fitter()->twoStageFit( this->twoStageFit() ); | LauFitter::fitter().twoStageFit( this->twoStageFit() ); | ||||
| LauFitter::fitter()->initialise( this, fitVars_ ); | LauFitter::fitter().initialise( this, fitVars_ ); | ||||
| this->startNewFit( LauFitter::fitter()->nParameters(), LauFitter::fitter()->nFreeParameters() ); | this->startNewFit( LauFitter::fitter().nParameters(), LauFitter::fitter().nFreeParameters() ); | ||||
| // Now ready for minimisation step | // Now ready for minimisation step | ||||
| std::cout << "\nINFO in LauAbsFitModel::fitExpt : Start minimisation...\n"; | std::cout << "\nINFO in LauAbsFitModel::fitExpt : Start minimisation...\n"; | ||||
| LauAbsFitter::FitStatus fitResult = LauFitter::fitter()->minimise(); | LauAbsFitter::FitStatus fitResult = LauFitter::fitter().minimise(); | ||||
| // If we're doing a two stage fit we can now release (i.e. float) | // If we're doing a two stage fit we can now release (i.e. float) | ||||
| // the 2nd stage parameters and re-fit | // the 2nd stage parameters and re-fit | ||||
| if (this->twoStageFit()) { | if (this->twoStageFit()) { | ||||
| if ( fitResult.status != 3 ) { | if ( fitResult.status != 3 ) { | ||||
| std::cerr << "WARNING in LauAbsFitModel:fitExpt : Not running second stage fit since first stage failed." << std::endl; | std::cerr << "WARNING in LauAbsFitModel:fitExpt : Not running second stage fit since first stage failed." << std::endl; | ||||
| LauFitter::fitter()->releaseSecondStageParameters(); | LauFitter::fitter().releaseSecondStageParameters(); | ||||
| } else { | } else { | ||||
| LauFitter::fitter()->releaseSecondStageParameters(); | LauFitter::fitter().releaseSecondStageParameters(); | ||||
| this->startNewFit( LauFitter::fitter()->nParameters(), LauFitter::fitter()->nFreeParameters() ); | this->startNewFit( LauFitter::fitter().nParameters(), LauFitter::fitter().nFreeParameters() ); | ||||
| fitResult = LauFitter::fitter()->minimise(); | fitResult = LauFitter::fitter().minimise(); | ||||
| } | } | ||||
| } | } | ||||
| const TMatrixD& covMat = LauFitter::fitter()->covarianceMatrix(); | const TMatrixD& covMat = LauFitter::fitter().covarianceMatrix(); | ||||
| this->storeFitStatus( fitResult, covMat ); | this->storeFitStatus( fitResult, covMat ); | ||||
| // Store the final fit results and errors into protected internal vectors that | // Store the final fit results and errors into protected internal vectors that | ||||
| // all sub-classes can use within their own finalFitResults implementation | // all sub-classes can use within their own finalFitResults implementation | ||||
| // used below (e.g. putting them into an ntuple in a root file) | // used below (e.g. putting them into an ntuple in a root file) | ||||
| LauFitter::fitter()->updateParameters(); | LauFitter::fitter().updateParameters(); | ||||
| } | } | ||||
| void LauAbsFitModel::calculateSPlotData() | void LauAbsFitModel::calculateSPlotData() | ||||
| { | { | ||||
| if (sPlotNtuple_ != 0) { | if (sPlotNtuple_ != 0) { | ||||
| sPlotNtuple_->addFriendTree(inputFitData_->fileName(), inputFitData_->treeName()); | sPlotNtuple_->addFriendTree(inputFitData_->fileName(), inputFitData_->treeName()); | ||||
| sPlotNtuple_->writeOutGenResults(); | sPlotNtuple_->writeOutGenResults(); | ||||
| LauSPlot splot(sPlotNtuple_->fileName(), sPlotNtuple_->treeName(), this->firstExpt(), this->nExpt(), | LauSPlot splot(sPlotNtuple_->fileName(), sPlotNtuple_->treeName(), this->firstExpt(), this->nExpt(), | ||||
| ▲ Show 20 Lines • Show All 484 Lines • Show Last 20 Lines | |||||