NuFCFitter Class Reference

Fits a FC experiment. More...

#include <NuFCFitter.h>

List of all members.

Public Member Functions
	NuFCFitter (NuFCExperiment *experiment, NuMMHelperCPT *helper, NuMatrixSpectrum *ndNQData, NuMatrixSpectrum *ndPQData, const NuMMParameters &trueparams)
virtual	~NuFCFitter ()
Double_t	Fit (const NuMMParameters &fitParameters)
void	AlwaysArchive (std::string filename)
	Call if you want to write out every experiment.
void	ArchiveTo (std::string filename)
void	Archive (std::string filename)
Private Member Functions
NuMMParameters	GridSearch (TH2D &likesurf, NuMMParameters mmGrid, NuMMParameters *Ref=0)
	Does a grid search on an arbitrary surface.
NuMMParameters	CoarseGridSearch (NuMMParameters start)
	Do a coarse grid search, to seed minuit.
NuMMParameters	FineGridSearch (NuMMParameters start)
	Does a finer grid search, if minuit fails.
NuMMParameters	MinuitFit (NuMMParameters start)
	Does a minuit fit.
Double_t	DeltaChi (const NuMMParameters &bestfit)
	Calculates the Delta chi2 value.
NuMMParameters	RecoverNegativeDeltaChi (const NuMMParameters &bestfit)
	Recovers a negative delta chi2 value.
void	ClassifyFitPoint (const NuMMParameters &mmFit)
	Examines the fit point and classifies it.
void	DoOutput (const NuMMParameters &mmFit, Double_t deltachi)
	Dos the summary output form the fitting.
void	ShiftHistogram (TH2D &surface, Double_t shift)
	Shifts a histogram. Doesn't really belong here, but not sure where else.
	ClassDef (NuFCFitter, 0)
Private Attributes
NuMMRunFC *	fRun
	Pointer to the run object.
NuMatrixSpectrum *	ffdNQData
NuMatrixSpectrum *	ffdPQData
NuSystFitter	fsystFitter
	Fitter object.
NuMMParameters *	fBestFit
	Copy of the best fit point.
Bool_t	fFineGridSearch
Bool_t	fNegativeDeltaChi
Bool_t	fHighFitPoint
Bool_t	fZeroFitPoint
Bool_t	fArchive
std::string	fArchiveFilename
const NuMMParameters	fTrueParams

---

Detailed Description

Fits a FC experiment.

Author:

Nick Devenish, Last checkin

Author

nickd

Version:

Revision

1.5

Date:

Date

2008/09/12 14:19:44

Created on: Sun Aug 17, 2008

Definition at line 32 of file NuFCFitter.h.

---

Constructor & Destructor Documentation

NuFCFitter::NuFCFitter (  NuFCExperiment *  experiment, NuMMHelperCPT *  helper, NuMatrixSpectrum *  ndNQData, NuMatrixSpectrum *  ndPQData, const NuMMParameters &  trueparams )

Creates the fitter object. We have to accept pointers here, as everything else we want to use only accepts them, and we can't easily change without possibly breaking large amounts of existing code. Parameters: experiment  The experiment (i.e. far data) that we are fitting. helper  The NuMMHelperCPT object that contains the helper data ndNQData  The near detector neutrino data ndPQData  The enar detector antineutrino data Definition at line 23 of file NuFCFitter.cxx. References ffdNQData, ffdPQData, fRun, fsystFitter, NuFCExperiment::GetNQSpectrum(), NuFCExperiment::GetPQSpectrum(), NuMMRunFC::PredictNus(), NuSystFitter::push_back(), NuSystFitter::QuietModeOn(), and NuMMRun::SetMaximumEnergy(). : fRun(0), ffdNQData(0), ffdPQData(0), fsystFitter(), fBestFit(0), fFineGridSearch(false), fNegativeDeltaChi(false), fHighFitPoint(false), fZeroFitPoint(false), fArchive(false), fArchiveFilename("archive.root"), fTrueParams(trueparams) { // Firstly, get copies of the far detector spectrum //NuMatrixSpectrum *fdNQ = &experiment->GetSpectrum(); ffdPQData = new NuMatrixSpectrum(experiment->GetPQSpectrum()); ffdNQData = new NuMatrixSpectrum(experiment->GetNQSpectrum()); // Create the run! fRun = new NuMMRunFC(helper, ndNQData, ndPQData, ffdNQData, ffdPQData); fRun->PredictNus(false); fRun->SetMaximumEnergy(49.999); // fRun->QuietModeOn(); // Now add this run to the fitter fsystFitter.push_back(fRun); fsystFitter.QuietModeOn(); }

NuFCFitter::~NuFCFitter (   )  [virtual]

Definition at line 60 of file NuFCFitter.cxx. References fBestFit, ffdNQData, ffdPQData, and fRun. { // Delete our run object if (fRun) delete fRun; fRun = 0; // Delete our copies of the far detector data if (ffdNQData) delete ffdNQData; if (ffdPQData) delete ffdPQData; ffdNQData = ffdPQData = 0; // Delete the copy of the best fit point if (fBestFit) delete fBestFit; fBestFit = 0; }

---

Member Function Documentation

void NuFCFitter::AlwaysArchive (  std::string  filename  )  [inline]

Call if you want to write out every experiment. Definition at line 57 of file NuFCFitter.h. Referenced by NuFCGridPoint::Run(). { fArchiveFilename = filename; fArchive = true; }

void NuFCFitter::Archive (  std::string  filename  )

Definition at line 468 of file NuFCFitter.cxx. References ffdPQData, NuMatrixSpectrum::GetName(), MSG, and NuMatrixSpectrum::Spectrum(). Referenced by Fit(). { MSG("NuFCFitter",Msg::kInfo) << "Disabled all archiving." << endl; return; if (!ffdPQData || !ffdNQData) { MSG("NUFCFitter",Msg::kError) << "No far detector spectrum: Cannot archive" << endl; return; } if (!fBestFit) { MSG("NuFCFitter", Msg::kError) << "Error: Trying to archive before fit" << endl; return; } MSG("NuFCFitter", Msg::kInfo) << " Archiving to file " << filename << "..." << endl; // The actual experiment spectrum // TH1D fdSpec(*(ffdPQData->Spectrum())); NuMatrixSpectrum fdSpec(*ffdPQData); // // fdSpec.RebinTo4GeV(); // // Grab the name out of this string basename = fdSpec.Spectrum()->GetName(); // // // No longer do the likelihood surface - it takes too long. Instead just write out the spectrum, // // and we can regenerate the surface later if we wish // // Now we want the likelihood surfaces. First do the linear portion string surfname = basename + "_surf"; TH2D likesurface(surfname.c_str(), surfname.c_str(), 50, 0, 1, 100, 0, 20e-3); // TH2D likesurface(surfname.c_str(), surfname.c_str(), 100, 0.9, 1, 100, 0.002, 0.003); // Double_t minlin = fsystFitter.FillLikelihoodSurfaceBar(likesurface, *fBestFit); // // // Now do the logarithmic plot // vector<Double_t> bins; // for (Int_t i = -204; i <= 0; i++) { // bins.push_back(pow(10, (Double_t)i/100.0)); // } // Int_t bincount = bins.size()-1; // string logname = surfname + "_log"; // TH2D logsurf(logname.c_str(), logname.c_str(), 50, 0, 1, bincount, &(bins.front())); // Double_t minlog = fsystFitter.FillLikelihoodSurfaceBar(logsurf, *fBestFit); // // // Now, shift both histograms by the minimum // Double_t surfmin = minlin < minlog ? minlin : minlog; // ShiftHistogram(likesurface, -surfmin); // ShiftHistogram(logsurf, -surfmin); // Now, write these out! TFile tf(filename.c_str(), "UPDATE"); fdSpec.Spectrum()->Write(); likesurface.Write(); // logsurf.Write(); tf.Close(); }

void NuFCFitter::ArchiveTo (  std::string  filename  )  [inline]

Definition at line 62 of file NuFCFitter.h. Referenced by NuFCGridPoint::Run(). { fArchiveFilename = filename; }

NuFCFitter::ClassDef (  NuFCFitter  , 0  )  [private]

void NuFCFitter::ClassifyFitPoint (  const NuMMParameters &  mmFit  )  [private]

Examines the fit point and classifies it. Definition at line 443 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), fHighFitPoint, fTrueParams, fZeroFitPoint, NuMMParameters::LowerDm2BarLimit(), NuMMParameters::LowerSn2BarLimit(), MSG, and NuMMParameters::Sn2Bar(). Referenced by Fit(). { // cout << " Checking against minimum dm: " << fTrueParams.LowerDm2BarLimit() << ", sn: " << fTrueParams.LowerSn2BarLimit() << endl; // Check to see if the fit is to 'zero' if (mmFit.Sn2Bar() <= fTrueParams.LowerSn2BarLimit() + 1e-4 || mmFit.Dm2Bar() <= fTrueParams.LowerDm2BarLimit() + 1e-6 ) { MSG("NuFCFitter", Msg::kInfo) << " Best fit is to no oscillation case\n"; fZeroFitPoint = true; } // Now check if the best fit is maxing out // This message is a little out of date, but will do for now if (mmFit.Dm2Bar() >= 30e-3 - 0.25e-4) { MSG("NuFCFitter", Msg::kInfo) << " Best fit is maxing out on dm2bar (" << mmFit.Dm2Bar() << ")" << endl; fHighFitPoint = true; } }

NuMMParameters NuFCFitter::CoarseGridSearch (  NuMMParameters  start  )  [private]

Do a coarse grid search, to seed minuit. Definition at line 198 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), NuFCConfig::FixedSin(), fTrueParams, NuFCConfig::GetConfig(), GridSearch(), MSG, and NuMMParameters::Sn2Bar(). Referenced by Fit(). { // The width and height of the grid Int_t gridWidth = 5; Double_t sinMin = 0.2; Double_t sinMax = 1.0; const Int_t gridHeight = 15; // Check if we are doing one-parameter fit, and if so constrain the grid search // to one-dimension NuFCConfig &cfg = NuFCConfig::GetConfig(); if (cfg.FixedSin()) { gridWidth = 1; sinMin = sinMax = fTrueParams.Sn2Bar(); } // Generate the surface TH2D likesurf("likesurf", "likesurf", gridWidth, sinMin, sinMax, gridHeight, 2e-3, 30e-3); mmGrid = GridSearch(likesurf, mmGrid); // Now do a log search vector<Double_t> bins; for (Int_t i = -15; i <= -1; i++) { bins.push_back(pow(10, (Double_t)i/10.0)); } Int_t bincount = bins.size()-1; TH2D logsurf("logsurf", "logsurf", gridWidth, sinMin, sinMax, bincount, &(bins.front())); mmGrid = GridSearch(logsurf, mmGrid, &mmGrid); // Output the minimum MSG("NuFCFitter", Msg::kInfo) << " Coarse grid search found sn2: " << mmGrid.Sn2Bar() << ", dm: " << mmGrid.Dm2Bar() << '\n'; return mmGrid; }

Double_t NuFCFitter::DeltaChi (  const NuMMParameters &  bestfit  )  [private]

Calculates the Delta chi2 value. Definition at line 358 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), fsystFitter, fTrueParams, NuMMParameters::LowerDm2BarLimit(), NuMMParameters::LowerSn2BarLimit(), MSG, NuMMParameters::Sn2Bar(), and NuMMParameters::VectorParameters(). Referenced by Fit(). { // Calculate the two likelihood points Double_t fitlike = fsystFitter(bestfit.VectorParameters()); Double_t truelikelihood = 0; // cout << "bfdm: " << bestfit.Dm2Bar() << ", limit = " << fTrueParams.LowerDm2BarLimit() << endl; // Don't allow the parameters to go below constrained - this can // give us an artificially negative delta chi2 if ( (fTrueParams.Dm2Bar() < fTrueParams.LowerDm2BarLimit()) || (fTrueParams.Sn2Bar() < fTrueParams.LowerSn2BarLimit()) ) { MSG("NuFCFitter", Msg::kDebug) << " True point is lower than fit is allowed." << endl; NuMMParameters pc(bestfit); pc.Dm2Bar(fTrueParams.LowerDm2BarLimit()); pc.Sn2Bar(fTrueParams.LowerSn2BarLimit()); truelikelihood = fsystFitter(pc.VectorParameters()); } else { truelikelihood = fsystFitter(fTrueParams.VectorParameters()); } // fmmPars.Dm2Bar(fmmPars.LowerDm2BarLimit()); // MAXMSG("NuFCGriPoint", 5, Msg::kInfo) << "Moving true grid point " // fmmPars.Sn2Bar(fmmPars.LowerSn2BarLimit()); // // Double_t truelikelihood = fsystFitter(fTrueParams.VectorParameters()); // Return the difference return truelikelihood - fitlike; }

void NuFCFitter::DoOutput (  const NuMMParameters &  mmFit, Double_t  deltachi )  [private]

Dos the summary output form the fitting. Definition at line 120 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), fFineGridSearch, fsystFitter, MSG, NuMMParameters::Sn2Bar(), and NuMMParameters::VectorParameters(). Referenced by Fit(). { // Output the fit results Double_t fitlike = fsystFitter(mmFit.VectorParameters()); MSG("NuFCFitter", Msg::kInfo) << " Minimum is dm2bar: " << mmFit.Dm2Bar() << " sn2bar: " << mmFit.Sn2Bar() << '\n' << " Fit Likelyhood: " << fitlike << '\n' << " Delta Chi2: " << deltachi << '\n'; // Did we want to ask the control script to log this? if (fFineGridSearch || fNegativeDeltaChi) MSG("NuFCFitter", Msg::kInfo) << " Please Log this experiment." << endl; }

NuMMParameters NuFCFitter::FineGridSearch (  NuMMParameters  start  )  [private]

Does a finer grid search, if minuit fails. Definition at line 237 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), NuFCConfig::FixedSin(), fTrueParams, NuFCConfig::GetConfig(), GridSearch(), NuMMParameters::LowerDm2BarLimit(), NuMMParameters::LowerSn2BarLimit(), MSG, and NuMMParameters::Sn2Bar(). Referenced by MinuitFit(). { // Calculate the location of the fine grid search Double_t gridL, gridR, gridT, gridB; gridL = gridR = gridT = gridB = -1.0; // The width and height of the grid Int_t gridWidth = 19; const Int_t gridHeight = 19; // Configuration NuFCConfig &cfg = NuFCConfig::GetConfig(); // change behavior if we want a one-parameter fit if (!cfg.FixedSin()) { // Go +- 0.2 in sin space gridL = mmStart.Sn2Bar() - 0.2; gridR = gridL + 0.4; // Now limit it to the physical realm if (gridL <= fTrueParams.LowerSn2BarLimit()) gridL = fTrueParams.LowerSn2BarLimit(); if (gridR > 1.0) gridR = 1.0; } else { // One-parameter fit gridL = gridR = fTrueParams.Sn2Bar(); gridWidth = 1; } // And go +-2e-3 in mass space gridB = mmStart.Dm2Bar() - 2e-3; gridT = gridB + 4e-3; // We only really need to worry about lower limit here if (gridB <= fTrueParams.LowerDm2BarLimit()) gridB = fTrueParams.LowerSn2BarLimit(); MSG("NuFCFitter", Msg::kInfo) << " Doing fine grid search on sin = [ " << gridL << ", " << gridR << " ]\n" << " dm = [ " << gridB << ", " << gridT << " ]\n"; // Ensure that we haven't done anything stupid like set low above high if (gridL > gridR || gridT <= gridB) { // cout << "Error: Negative range in fine grid search" << endl; throw runtime_error("Negative range in fine grid search"); } // Make the surface // Generate the surface TH2D likesurf("likesurf", "likesurf", gridWidth, gridL, gridR, gridHeight, gridB, gridT); // And now do a grid search on this NuMMParameters mmGrid = GridSearch(likesurf, mmStart); // Temporarily, archive any fine grid surface // TFile tf("archives/finegrid.root", "UPDATE"); // likesurf.Write(); // tf.Close(); // Output the minimum MSG("NuFCFitter", Msg::kInfo) << " Fine grid search found sn2: " << mmGrid.Sn2Bar() << ", dm: " << mmGrid.Dm2Bar() << '\n'; return mmGrid; }

Double_t NuFCFitter::Fit (  const NuMMParameters &  fitParameters  )

Fits the previously passed in FC experiment. Returns an std::pair with the delta chi2, and the best fit Parameters: fitParameters  a parameters object, set to the true value Definition at line 78 of file NuFCFitter.cxx. References Archive(), ClassifyFitPoint(), CoarseGridSearch(), DeltaChi(), DoOutput(), fArchiveFilename, fBestFit, fFineGridSearch, MinuitFit(), and RecoverNegativeDeltaChi(). Referenced by NuFCGridPoint::Run(). { // Start with a coarse grid search, to seed minuit const NuMMParameters mmGrid = CoarseGridSearch(fitParameters); // Now do a minuit fit NuMMParameters mmFit = MinuitFit(mmGrid); Double_t dchi = DeltaChi(mmFit); // If this gave a negative value, rescue it (allow a tiny gap) if (dchi < 0) { if (dchi < -1e-5) { mmFit = RecoverNegativeDeltaChi(mmFit); dchi = DeltaChi(mmFit); } else { // It is below zero, but within error. Set it to zero! dchi = 0; } } // Make a copy of the best fit point if (fBestFit) delete fBestFit; fBestFit = new NuMMParameters(mmFit); // Classify the best fit point i.e. is it zero, high etc // this is mainly for informational purposes (i.e. the run script) ClassifyFitPoint(mmFit); // Do we want to archive? if (fFineGridSearch || fArchive) { // || fNegativeDeltaChi Archive(fArchiveFilename); } // Do the output DoOutput(mmFit, dchi); return dchi; }

NuMMParameters NuFCFitter::GridSearch (  TH2D &  likesurf, NuMMParameters  mmGrid, NuMMParameters *  Ref = 0 )  [private]

Does a grid search on an arbitrary surface. Definition at line 306 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), NuSystFitter::FillLikelihoodSurfaceBar(), NuFCConfig::FixedSin(), fsystFitter, fTrueParams, NuFCConfig::GetConfig(), MSG, NuMMParameters::Sn2Bar(), and NuMMParameters::VectorParameters(). Referenced by CoarseGridSearch(), and FineGridSearch(). { Double_t likelihood = fsystFitter.FillLikelihoodSurfaceBar(likesurf, mmGrid); // Find the point on this surface that is minimum Int_t minbinX = 0, minbinY = 0, minbinZ = 0; likesurf.GetMinimumBin(minbinX, minbinY, minbinZ); Double_t minsin = likesurf.GetXaxis()->GetBinCenter(minbinX); Double_t minmass = likesurf.GetYaxis()->GetBinCenter(minbinY); // Access the configuration NuFCConfig &cfg = NuFCConfig::GetConfig(); // Check the zero point, to see if it is lower than anything else on the grid mmGrid.Dm2Bar(0); // Different 'zero' point if one parameter if (cfg.FixedSin()) { mmGrid.Sn2Bar(fTrueParams.Sn2Bar()); } else { mmGrid.Sn2Bar(0); } Double_t zeropoint = fsystFitter(mmGrid.VectorParameters()); if (zeropoint < likelihood) { // Then the no-oscillation case is the minimum minsin = 0.0; minmass = 0.0; if (cfg.FixedSin()) { minsin = fTrueParams.Sn2Bar(); } } // If we have been given a reference point, check to see if this is lower if (Ref) { Double_t refpoint = fsystFitter(Ref->VectorParameters()); if (refpoint < likelihood) { MSG("NuFCFitter", Msg::kDebug) << " Reference point is lower than grid search" << endl; minsin = Ref->Sn2Bar(); minmass = Ref->Dm2Bar(); } } // Set mmPars to the newly found point, then return it mmGrid.Sn2Bar(minsin); mmGrid.Dm2Bar(minmass); return mmGrid; }

NuMMParameters NuFCFitter::MinuitFit (  NuMMParameters  start  )  [private]

Does a minuit fit. Definition at line 138 of file NuFCFitter.cxx. References NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), fFineGridSearch, FineGridSearch(), NuFCConfig::FixedSin(), NuMMParameters::FixSn2Bar(), fsystFitter, fTrueParams, NuFCConfig::GetConfig(), NuMMParameters::LowerDm2BarLimit(), NuMMParameters::LowerSn2BarLimit(), NuSystFitter::Minimise(), NuMMParameters::MinuitPass(), MSG, NuMMParameters::Sn2(), NuMMParameters::Sn2Bar(), and NuMMParameters::VectorParameters(). Referenced by Fit(). { // If it is trying to start us at zero, reset to the search minimum // fTrueParams.LowerDm2BarLimit()) { if (mmStart.Sn2Bar() == 0) { MSG("NuFCFitter", Msg::kDebug) << " Starting minuit at no-oscillation; shifting to lower limits" << endl; mmStart.Sn2Bar(fTrueParams.LowerDm2BarLimit()); mmStart.Dm2Bar(fTrueParams.LowerSn2BarLimit()); } // Fix the sin, if specified NuFCConfig &cfg = NuFCConfig::GetConfig(); if (cfg.FixedSin()) { mmStart.FixSn2Bar(); mmStart.Sn2Bar(fTrueParams.Sn2Bar()); } // Now run minuit itself NuMMParameters mmFit = fsystFitter.Minimise(mmStart); // Did minuit succeed? if (mmFit.MinuitPass()) { // Tell them what minuit found MSG("NuFCFitter", Msg::kInfo) << " Minuit found: dm2: " << mmFit.Dm2() << " sn2: " << mmFit.Sn2() << " dm2bar: " << mmFit.Dm2Bar() << " sn2bar: " << mmFit.Sn2Bar() << '\n'; } else { // Minuit failed to fit. Use a more intensive grid search, // starting on the minuit starting point (that we get from // a grid search) MSG("NuFCFitter", Msg::kInfo) << " Minuit fit failed." << endl; MSG("NuFCFitter", Msg::kInfo) << " Searching fine grid around sn2: " << mmStart.Sn2Bar() << ", dm2: " << mmStart.Dm2Bar() << endl; NuMMParameters mmFineGrid = FineGridSearch(mmStart); fFineGridSearch = true; // Test the fine grid search against the result from minuit = it could be possible // that minuit found a 'better' point, even if it failed Double_t minuit_chi = fsystFitter(mmFit.VectorParameters()); Double_t finegrid_chi = fsystFitter(mmFineGrid.VectorParameters()); if (minuit_chi > finegrid_chi) { // The new, fine grid point is better mmFit = mmFineGrid; } else { // The minuit fail point is actually better. Use that. MSG("NuFCFitter", Msg::kInfo) << " Failed Minuit fit is better than fine grid, using." << endl; } } return mmFit; }

NuMMParameters NuFCFitter::RecoverNegativeDeltaChi (  const NuMMParameters &  bestfit  )  [private]

Recovers a negative delta chi2 value. Definition at line 394 of file NuFCFitter.cxx. References NuMMParameters::Dm2Bar(), NuFCConfig::FixedSin(), NuMMParameters::FixSn2Bar(), fNegativeDeltaChi, fsystFitter, fTrueParams, NuFCConfig::GetConfig(), NuSystFitter::Minimise(), NuMMParameters::MinuitPass(), MSG, NuMMParameters::Sn2Bar(), and NuMMParameters::VectorParameters(). Referenced by Fit(). { // Set the failure flag fNegativeDeltaChi = true; // Recalculate the likelihoods Double_t truelike = fsystFitter(fTrueParams.VectorParameters()); Double_t oldlike = fsystFitter(bestfit.VectorParameters());; // Make a copy of the true parameters, to ensure that we do a 1D fit if specified // (this is assuming we cannot guarantee that the true parameters has the parameter fixed) NuMMParameters tpcopy(fTrueParams); NuFCConfig &cfg = NuFCConfig::GetConfig(); if (cfg.FixedSin()) { tpcopy.FixSn2Bar(); } // NuMMParameters newFit = fsystFitter.Minimise(fTrueParams); NuMMParameters newFit = fsystFitter.Minimise(tpcopy); MSG("NuFCFitter", Msg::kInfo) << " Negative delta chi2 likelihood (" << truelike - oldlike << ") found. Refitting from true (" << truelike << ")" << endl; // If this passed, and is lower than our old likelihood, then use it if (newFit.MinuitPass()) { MSG("NuFCFitter", Msg::kInfo) << " Refit passed." << endl; Double_t newlike = fsystFitter(newFit.VectorParameters()); if (newlike > oldlike) { MSG("NuFCFitter", Msg::kInfo) << " New fit point ( dm: " << newFit.Dm2Bar() << ", sn: " << newFit.Sn2Bar() << " )\n" << " at " << newlike << " is higher than before ( " << oldlike << " )! strange!\n" << " Returning true fit point, as is lowest we have seen." << endl; newFit = fTrueParams; } } else { MSG("NuFCFitter", Msg::kError) << " Unrecoverable error: Refit on negative deltachi failed." << endl; throw logic_error("NuFCFitter Refit on negative deltachi failed."); } return newFit; }

void NuFCFitter::ShiftHistogram (  TH2D &  surface, Double_t  shift )  [private]

Shifts a histogram. Doesn't really belong here, but not sure where else. Definition at line 527 of file NuFCFitter.cxx. { for (Int_t x = 1; x <= surface.GetNbinsX(); x++) { for (Int_t y = 1; y <= surface.GetNbinsY(); y++) { Double_t contents = surface.GetBinContent(x, y); surface.SetBinContent(x, y, contents + shift); } } }

---

Member Data Documentation

Bool_t NuFCFitter::fArchive [private]

Definition at line 116 of file NuFCFitter.h.

std::string NuFCFitter::fArchiveFilename [private]

Definition at line 117 of file NuFCFitter.h. Referenced by Fit().

NuMMParameters* NuFCFitter::fBestFit [private]

Copy of the best fit point. Definition at line 107 of file NuFCFitter.h. Referenced by Fit(), and ~NuFCFitter().

NuMatrixSpectrum* NuFCFitter::ffdNQData [private]

Pointers to copies of the far detector spectrum - the NuMMRunFC objects only accept pointers, but don't alter or copy them. Definition at line 100 of file NuFCFitter.h. Referenced by NuFCFitter(), and ~NuFCFitter().

NuMatrixSpectrum* NuFCFitter::ffdPQData [private]

Definition at line 101 of file NuFCFitter.h. Referenced by Archive(), NuFCFitter(), and ~NuFCFitter().

Bool_t NuFCFitter::fFineGridSearch [private]

Definition at line 110 of file NuFCFitter.h. Referenced by DoOutput(), Fit(), and MinuitFit().

Bool_t NuFCFitter::fHighFitPoint [private]

Definition at line 113 of file NuFCFitter.h. Referenced by ClassifyFitPoint().

Bool_t NuFCFitter::fNegativeDeltaChi [private]

Definition at line 111 of file NuFCFitter.h. Referenced by RecoverNegativeDeltaChi().

NuMMRunFC* NuFCFitter::fRun [private]

Pointer to the run object. Definition at line 96 of file NuFCFitter.h. Referenced by NuFCFitter(), and ~NuFCFitter().

NuSystFitter NuFCFitter::fsystFitter [private]

Fitter object. Definition at line 104 of file NuFCFitter.h. Referenced by DeltaChi(), DoOutput(), GridSearch(), MinuitFit(), NuFCFitter(), and RecoverNegativeDeltaChi().

const NuMMParameters NuFCFitter::fTrueParams [private]

Definition at line 120 of file NuFCFitter.h. Referenced by ClassifyFitPoint(), CoarseGridSearch(), DeltaChi(), FineGridSearch(), GridSearch(), MinuitFit(), and RecoverNegativeDeltaChi().

Bool_t NuFCFitter::fZeroFitPoint [private]

Definition at line 114 of file NuFCFitter.h. Referenced by ClassifyFitPoint().

---

The documentation for this class was generated from the following files:

• NuFCFitter.h
• NuFCFitter.cxx

---