NuFCFitter.h

Go to the documentation of this file.

#ifndef NUFCFITTER_H
#define NUFCFITTER_H

// class NuMMParamters;
class NuMatrixSpectrum;
class NuFCExperiment;
class NuMMHelperCPT;
class NuMMRunFC;

// #include "NtupleUtils/NuMMRunFC.h"
#include "NtupleUtils/NuSystFitter.h"
#include "NtupleUtils/NuMMParameters.h"

class NuFCFitter
{
public:
  NuFCFitter(NuFCExperiment *experiment, NuMMHelperCPT *helper,
      NuMatrixSpectrum *ndNQData, NuMatrixSpectrum *ndPQData,
      const NuMMParameters& trueparams);
  NuFCFitter(NuMMHelperCPT *helper,
             NuMatrixSpectrum *fdNQData, NuMatrixSpectrum *fdPQData,
             NuMatrixSpectrum *ndNQData, NuMatrixSpectrum *ndPQData,
             const NuMMParameters& trueparams);
  NuFCFitter(NuMMRunFC *_run,
             const NuMMParameters& trueparams);
  virtual ~NuFCFitter();
  
  Double_t Fit(const NuMMParameters &fitParameters);

  const NuMMParameters* BestFit() const;

  void AlwaysArchive(std::string filename) {
    fArchiveFilename = filename;
    fArchive = true;
  }

  void ArchiveTo(std::string filename) {
    fArchiveFilename = filename;
  }
  
  void Archive(std::string filename);
  
  void SetSilent(bool _silent = true) {silent = _silent;}
  
  
private:
  NuMMParameters GridSearch(TH2D &likesurf, NuMMParameters mmGrid, NuMMParameters *Ref = 0);

  NuMMParameters CoarseGridSearch(NuMMParameters start);
  NuMMParameters FineGridSearch(NuMMParameters start);
  
  NuMMParameters MinuitFit(NuMMParameters start);
  
  Double_t DeltaChi(const NuMMParameters& bestfit);
    
  NuMMParameters RecoverNegativeDeltaChi( const NuMMParameters& bestfit);
  
  void ClassifyFitPoint(const NuMMParameters& mmFit);
  
  void DoOutput(const NuMMParameters &mmFit, Double_t deltachi );
  
  void ShiftHistogram(TH2D &surface, Double_t shift);
  
  NuMMRunFC *fRun;
  
  NuMatrixSpectrum *ffdNQData;
  NuMatrixSpectrum *ffdPQData;
  
  NuSystFitter fsystFitter;
  
  NuMMParameters *fBestFit;
  
  // Flags to check if various stages of fit failed
  Bool_t fFineGridSearch;
  Bool_t fNegativeDeltaChi;
  // Classifications
  Bool_t fHighFitPoint;
  Bool_t fZeroFitPoint;
  
  Bool_t fArchive;
  std::string fArchiveFilename;

  // Copy of the true parameters
  const NuMMParameters fTrueParams;
  
  bool silent;
  
  ClassDef(NuFCFitter,0);
};

#endif

---