Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

NuEZFitter Class Reference

Simplifies the process of setting up and running a fit. More...

#include <NuEZFitter.h>

List of all members.

Public Member Functions

 NuEZFitter (std::string Helper, std::string NeardetData)
 Creates the fitter for a single run.
 NuEZFitter (std::string Helper, std::string NeardetData, std::string FardetData)
 Creates a fitter, using far detector data.
 NuEZFitter (const NuEZFitter &)
 ~NuEZFitter ()
void MakeFDData (Double_t POT, Double_t dm2bar, Double_t sn2bar, Double_t dm2=2.4e-3, Double_t sn2=1.0)
 Creates a FD data spectrum This function generates a fake FD spectrum to use as FD data by extrapolating the near detector data.
void UseFDData (std::string FardetData)
 Uses an existing Far detector data spectrum, by loading from file.
void UseFDData (const NuMatrixSpectrum &fdNQdata, const NuMatrixSpectrum &fdPQdata)
 Uses an exiting spectrum as far detector data.
void UseFDFakeData (Double_t POT, std::string FardetData, std::string FarDetTauData)
 Uses an existing far detector fake-data spectrum.
void SetPoT (Double_t POT)
 Set's the value of the far detector POT.
void ScaleToPot (Double_t POT)
 Scale the value of the far detector POT.
Double_t Chi2 (Double_t dm2bar, Double_t sn2bar, Double_t dm2=2.4e-3, Double_t sn2=1.0) const
 Returns the Delta-chi squared of a point This function accepts a set of oscillation parameters, makes a prediction at those values and then does a log-likelihood comparison to calculate the chi2.
Double_t Chi2 (const NuMMParameters &params) const
 Returns the Delta-chi squared of a point Accepts a NuMMParameters, and returns the chi2 at the point represented.
Bool_t Fit (const NuMMParameters &Initial)
 Does a fit to the data Does a simple Minuit fit to the data. Returns a bool indicating how successful the fit was. Get the fit results through FitResults().
NuMatrixSpectrum PredictionNQ (Double_t dm2bar, Double_t sn2bar, Double_t dm2=2.4e-3, Double_t sn2=1.0) const
 Generates an NQ prediction.
NuMatrixSpectrum PredictionNQ (const NuMMParameters &params) const
 Generates an NQ prediction.
NuMatrixSpectrum PredictionPQ (Double_t dm2bar, Double_t sn2bar, Double_t dm2=2.4e-3, Double_t sn2=1.0) const
 Generates an PQ prediction.
NuMatrixSpectrum PredictionPQ (const NuMMParameters &params) const
 Generates an PQ prediction.
NuSystFitter * Fitter ()
 Returns a pointer to the fitter.
NuMMRunFC * Run ()
 Returns a pointer to the run object.
const NuMatrixSpectrum & NDDataNQ () const
 Returns the NQ Near detector spectrum.
const NuMatrixSpectrum & NDDataPQ () const
 Returns the PQ Near detector spectrum.
const NuMatrixSpectrum & FDDataNQ () const
 Returns the NQ Far detector spectrum.
const NuMatrixSpectrum & FDDataPQ () const
 Returns the PQ Far detector spectrum.
const NuMMParameters & FitResult () const
 Returns the results of the fit.

Private Member Functions

void Init (void)
 Initialises the class variables.
void Init_Helper_Near (std::string Helper, std::string NeardetData)
 Initialises the near detector and helper objects.
void Rebuild (void)
 Rebuilds the internal objects.

Private Attributes

NuMatrixSpectrum * fNDDataPQ
 PQ ND Spectrum.
NuMatrixSpectrum * fNDDataNQ
 NQ ND Spectrum.
NuMatrixSpectrum * fFDDataPQ
 PQ FD Spectrum.
NuMatrixSpectrum * fFDDataNQ
 NQ FD Spectrum.
NuMMHelperCPT * fHelper
 The MM helper object.
NuMMRunFC * fRun
 The run object that does the extrapolation.
NuSystFitter * fFitter
 The fitter object.
NuMMParameters * fFitResult

Detailed Description

Simplifies the process of setting up and running a fit.

Author:
Nick Devenish, Last checkin
Author
ahimmel
Version:
Revision
1.3
Date:
Date
2009/10/15 16:41:38
Created on: Sun Dec 14, 2008

Definition at line 32 of file NuEZFitter.h.

Constructor & Destructor Documentation

NuEZFitter::NuEZFitter (  std::string  Helper,
std::string  NeardetData
) 
 

Creates the fitter for a single run.

Parameters:
NeardetData The near detector data
Helper The helper file

Definition at line 41 of file NuEZFitter.cxx.

References Init(), Init_Helper_Near(), and Rebuild().

00041                                                            :
00042  //  fNDDataPQ(0), fNDDataNQ(0), fFDDataPQ(0), fFDDataNQ(0),
00043  //  fHelper(0), fRun(0), fFitter(0), fFitResult(0)
00044 {
00045   Init();
00046   Init_Helper_Near(Helper, NeardetData);
00047   Rebuild();
00048 }

NuEZFitter::NuEZFitter (  std::string  Helper,
std::string  NeardetData,
std::string  FardetData
) 
 

Creates a fitter, using far detector data.

Definition at line 50 of file NuEZFitter.cxx.

References Init(), Init_Helper_Near(), and UseFDData().

00051 {
00052   Init();
00053   Init_Helper_Near(Helper, NeardetData);
00054   UseFDData(FardetData);
00055 }

NuEZFitter::NuEZFitter (  const NuEZFitter &   ) 
 

Definition at line 57 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, fFitResult, fFitter, fHelper, fNDDataNQ, fNDDataPQ, fRun, and Init().

00058 { 
00059   Init();
00060   
00061   if (o.fNDDataPQ) fNDDataPQ = new NuMatrixSpectrum(*o.fNDDataPQ);
00062   if (o.fNDDataNQ) fNDDataNQ = new NuMatrixSpectrum(*o.fNDDataNQ);
00063   if (o.fFDDataPQ) fFDDataPQ = new NuMatrixSpectrum(*o.fFDDataPQ);
00064   if (o.fFDDataNQ) fFDDataNQ = new NuMatrixSpectrum(*o.fFDDataNQ);
00065   
00066   if (o.fHelper) fHelper = new NuMMHelperCPT(*o.fHelper);
00067   if (o.fRun) fRun = new NuMMRunFC(*o.fRun);
00068   if (o.fFitter) fFitter = new NuSystFitter(*o.fFitter);
00069   if (o.fFitResult) fFitResult = new NuMMParameters(*o.fFitResult);
00070 }

NuEZFitter::~NuEZFitter (   ) 
 

Definition at line 72 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, fFitResult, fFitter, fHelper, fNDDataNQ, fNDDataPQ, and fRun.

00073 {
00074   // Free all of our memory
00075   if (fNDDataPQ) delete fNDDataPQ; fNDDataPQ = 0;
00076   if (fNDDataNQ) delete fNDDataNQ; fNDDataNQ = 0;
00077   
00078   if (fFDDataPQ) delete fFDDataPQ;  fFDDataPQ = 0;
00079   if (fFDDataNQ) delete fFDDataNQ;  fFDDataNQ = 0;
00080 
00081   if (fHelper) delete fHelper; fHelper = 0;
00082   if (fRun) delete fRun; fRun = 0;
00083   if (fFitter) delete fFitter; fFitter = 0;
00084 
00085   if (fFitResult) delete fFitResult; fFitResult = 0;
00086 // if () delete ;  = 0;
00087 
00088 }

Member Function Documentation

Double_t NuEZFitter::Chi2 (  const NuMMParameters &  params  )  const
 

Returns the Delta-chi squared of a point Accepts a NuMMParameters, and returns the chi2 at the point represented.

Definition at line 259 of file NuEZFitter.cxx.

References MSG, and NuMMParameters::VectorParameters().

00260 {
00261   if (fFitter) {
00262     return (*fFitter)(params.VectorParameters());
00263   } else {
00264     MSG("NuEZFitter",Msg::kError) << "No Fitter exists. Have you added FD data?" << endl;
00265     return 0;
00266   }
00267 }

Double_t NuEZFitter::Chi2 (  Double_t  dm2bar,
Double_t  sn2bar,
Double_t  dm2 = 2.4e-3,
Double_t  sn2 = 1.0
)  const
 

Returns the Delta-chi squared of a point This function accepts a set of oscillation parameters, makes a prediction at those values and then does a log-likelihood comparison to calculate the chi2.

Definition at line 247 of file NuEZFitter.cxx.

References NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), NuMMParameters::Sn2(), and NuMMParameters::Sn2Bar().

00248 {
00249   NuMMParameters mmPars;
00250   mmPars.Dm2(dm2);
00251   mmPars.Sn2(sn2);
00252   mmPars.Dm2Bar(dm2bar);
00253   mmPars.Sn2Bar(sn2bar);
00254   
00255   // return (*fFitter)(mmPars.VectorParameters());
00256   return Chi2(mmPars);
00257 }

const NuMatrixSpectrum& NuEZFitter::FDDataNQ (   )  const [inline]
 

Returns the NQ Far detector spectrum.

Definition at line 111 of file NuEZFitter.h.

00111 { return *fFDDataNQ; }

const NuMatrixSpectrum& NuEZFitter::FDDataPQ (   )  const [inline]
 

Returns the PQ Far detector spectrum.

Definition at line 113 of file NuEZFitter.h.

00113 { return *fFDDataPQ; }

Bool_t NuEZFitter::Fit (  const NuMMParameters &  Initial  ) 
 

Does a fit to the data Does a simple Minuit fit to the data. Returns a bool indicating how successful the fit was. Get the fit results through FitResults().

Parameters:
Initial The initial parameters of the fit

Definition at line 313 of file NuEZFitter.cxx.

References fFitResult, fFitter, NuSystFitter::Minimise(), NuMMParameters::MinuitPass(), and MSG.

00314 {
00315   // Can we do a fit?
00316   if (!fFitter) {
00317     MSG("NuEZFitter",Msg::kError) << "No Fitter exists. Have you added FD data?" << endl;
00318     fFitResult = new NuMMParameters(Initial);
00319     return false;
00320   }
00321   
00322   // Do the fit
00323   // Bool_t minuitresult = false;
00324   NuMMParameters mmFit = fFitter->Minimise(Initial);
00325   // Copy the fit result to our class member
00326   if (fFitResult) delete fFitResult; fFitResult = 0;
00327   fFitResult = new NuMMParameters(mmFit);
00328   
00329   return fFitResult->MinuitPass();
00330 }

const NuMMParameters& NuEZFitter::FitResult (   )  const [inline]
 

Returns the results of the fit.

Definition at line 116 of file NuEZFitter.h.

00116 { return *fFitResult; }

NuSystFitter* NuEZFitter::Fitter (   )  [inline]
 

Returns a pointer to the fitter.

Definition at line 101 of file NuEZFitter.h.

00101 { return fFitter; }

void NuEZFitter::Init (  void   )  [private]
 

Initialises the class variables.

Definition at line 19 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, fFitResult, fFitter, fHelper, fNDDataNQ, fNDDataPQ, fRun, and NuMMParameters::MinuitPass().

Referenced by NuEZFitter().

00020 {
00021   fNDDataPQ = fNDDataNQ = fFDDataPQ = fFDDataNQ = 0;
00022   fHelper = 0; fRun = 0; fFitter = 0; fFitResult = 0;
00023   
00024   // Make a dummy fit result
00025   fFitResult = new NuMMParameters();
00026   fFitResult->MinuitPass(false);
00027 }

void NuEZFitter::Init_Helper_Near (  std::string  Helper,
std::string  NeardetData
)  [private]
 

Initialises the near detector and helper objects.

Definition at line 29 of file NuEZFitter.cxx.

References fHelper, fNDDataNQ, fNDDataPQ, and MSG.

Referenced by NuEZFitter().

00030 {
00031   // Reading near detector data
00032   MSG("NuEZFitter",Msg::kInfo) << "Reading Near detector data " << NeardetData << endl;
00033   fNDDataPQ = new NuMatrixSpectrum(NeardetData, "RecoEnergyPQ_ND");
00034   fNDDataNQ = new NuMatrixSpectrum(NeardetData, "RecoEnergy_ND");
00035   
00036   // Grab the helper file
00037   MSG("NuEZFitter",Msg::kInfo) << "Extrapolating with helper " << Helper << endl;
00038   string xsec = "$SRT_PRIVATE_CONTEXT/NtupleUtils/data/xsec_minos_modbyrs4_v3_5_0_mk.root";
00039   fHelper = new NuMMHelperCPT(Helper.c_str(),xsec);
00040 }

void NuEZFitter::MakeFDData (  Double_t  POT,
Double_t  dm2bar,
Double_t  sn2bar,
Double_t  dm2 = 2.4e-3,
Double_t  sn2 = 1.0
) 
 

Creates a FD data spectrum This function generates a fake FD spectrum to use as FD data by extrapolating the near detector data.

Parameters:
POT The exposure at the far detector

Definition at line 124 of file NuEZFitter.cxx.

References NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), fFDDataNQ, fFDDataPQ, fHelper, fNDDataNQ, fNDDataPQ, NuMMRunNuBar::MakeFDPredNoPair(), NuMMRun::QuietModeOn(), Rebuild(), NuMMRun::SetMaximumEnergy(), NuMMParameters::Sn2(), and NuMMParameters::Sn2Bar().

00125 {
00126     // Make the extrapolator
00127     NuMatrixSpectrum blankpot(POT);
00128     NuMMRunFC mmFakeFar(fHelper, fNDDataNQ, fNDDataPQ, &blankpot, &blankpot);
00129     mmFakeFar.QuietModeOn();
00130     mmFakeFar.SetMaximumEnergy(49.999);
00131  
00132     NuMMParameters mmPars;
00133     mmPars.Dm2(dm2);
00134     mmPars.Sn2(sn2);
00135     mmPars.Dm2Bar(dm2bar);
00136     mmPars.Sn2Bar(sn2bar);
00137     
00138     // Clear any existing far detector spectra
00139     if (fFDDataPQ) delete fFDDataPQ; fFDDataPQ = 0;
00140     if (fFDDataNQ) delete fFDDataNQ; fFDDataNQ = 0;
00141 
00142     // Grab the extrapolations
00143     // NuMatrixSpectrum *nqms = 0;
00144     // NuMatrixSpectrum *pqms = 0;
00145     mmFakeFar.MakeFDPredNoPair(mmPars, &fFDDataNQ, &fFDDataPQ);
00146     
00147     // Rebuild the internal objects
00148     Rebuild();
00149 }

const NuMatrixSpectrum& NuEZFitter::NDDataNQ (   )  const [inline]
 

Returns the NQ Near detector spectrum.

Definition at line 106 of file NuEZFitter.h.

00106 { return *fNDDataNQ; }

const NuMatrixSpectrum& NuEZFitter::NDDataPQ (   )  const [inline]
 

Returns the PQ Near detector spectrum.

Definition at line 108 of file NuEZFitter.h.

00108 { return *fNDDataPQ; }

NuMatrixSpectrum NuEZFitter::PredictionNQ (  const NuMMParameters &  params  )  const
 

Generates an NQ prediction.

Definition at line 293 of file NuEZFitter.cxx.

References fRun, and NuMMRunNuBar::MakeFDPredNuMu().

00294 {
00295   // Generate a prediction and return it in a memory-safe way
00296   NuMatrixSpectrum *pred = fRun->MakeFDPredNuMu(params);
00297   NuMatrixSpectrum ret(*pred);
00298   if (pred) delete pred; pred = 0;
00299   
00300   return ret;
00301 }

NuMatrixSpectrum NuEZFitter::PredictionNQ (  Double_t  dm2bar,
Double_t  sn2bar,
Double_t  dm2 = 2.4e-3,
Double_t  sn2 = 1.0
)  const
 

Generates an NQ prediction.

Definition at line 269 of file NuEZFitter.cxx.

References NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), NuMMParameters::Sn2(), and NuMMParameters::Sn2Bar().

00271 {
00272   NuMMParameters mmPars;
00273   mmPars.Dm2(dm2);
00274   mmPars.Sn2(sn2);
00275   mmPars.Dm2Bar(dm2bar);
00276   mmPars.Sn2Bar(sn2bar);
00277   
00278   return PredictionNQ(mmPars);
00279 }

NuMatrixSpectrum NuEZFitter::PredictionPQ (  const NuMMParameters &  params  )  const
 

Generates an PQ prediction.

Definition at line 303 of file NuEZFitter.cxx.

References fRun, and NuMMRunNuBar::MakeFDPredNuBar().

00304 {
00305   // Generate a prediction and return it in a memory-safe way
00306   NuMatrixSpectrum *pred = fRun->MakeFDPredNuBar(params);
00307   NuMatrixSpectrum ret(*pred);
00308   if (pred) delete pred; pred = 0;
00309   
00310   return ret;
00311 }

NuMatrixSpectrum NuEZFitter::PredictionPQ (  Double_t  dm2bar,
Double_t  sn2bar,
Double_t  dm2 = 2.4e-3,
Double_t  sn2 = 1.0
)  const
 

Generates an PQ prediction.

Definition at line 281 of file NuEZFitter.cxx.

References NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), NuMMParameters::Sn2(), and NuMMParameters::Sn2Bar().

00283 {
00284   NuMMParameters mmPars;
00285   mmPars.Dm2(dm2);
00286   mmPars.Sn2(sn2);
00287   mmPars.Dm2Bar(dm2bar);
00288   mmPars.Sn2Bar(sn2bar);
00289   
00290   return PredictionPQ(mmPars);
00291 }

void NuEZFitter::Rebuild (  void   )  [private]
 

Rebuilds the internal objects.

Definition at line 95 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, fFitter, fHelper, fNDDataNQ, fNDDataPQ, fRun, MSG, NuSystFitter::push_back(), and NuMMRun::SetMaximumEnergy().

Referenced by MakeFDData(), NuEZFitter(), ScaleToPot(), SetPoT(), UseFDData(), and UseFDFakeData().

00096 {
00097   // Check the needed plots exist
00098   if (!(fHelper || fNDDataNQ || fNDDataPQ || fFDDataNQ || fFDDataPQ)) {
00099     // We don't have all the required plots - just return without doing anything
00100     MSG("NuEZFitter",Msg::kWarning) << "Rebuild() was called without all the plots in place" << endl;
00101     return;
00102   }
00103   
00104   // Delete the old objects if necessary
00105   if (fRun) delete fRun; fRun = 0;
00106   if (fFitter) delete fFitter; fFitter = 0;
00107   
00108   // Build the run object
00109   // Use a different initialisation if we have no far detector data
00110   if (!(fFDDataPQ && fFDDataNQ)) {
00111     fRun = new NuMMRunFC(fHelper, fNDDataNQ, fNDDataPQ, 0);
00112   } else {
00113     fRun = new NuMMRunFC(fHelper, fNDDataNQ, fNDDataPQ, fFDDataNQ, fFDDataPQ);
00114   }
00115 
00116   // Limit the maximum fitting energy
00117   fRun->SetMaximumEnergy(49.99999);
00118   
00119   // And the fitter
00120   fFitter = new NuSystFitter();
00121   fFitter->push_back(fRun);
00122 }

NuMMRunFC* NuEZFitter::Run (   )  [inline]
 

Returns a pointer to the run object.

Definition at line 103 of file NuEZFitter.h.

00103 { return fRun; }

void NuEZFitter::ScaleToPot (  Double_t  POT  ) 
 

Scale the value of the far detector POT.

Definition at line 232 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, MSG, Rebuild(), and NuMatrixSpectrum::ScaleToPot().

00233 {
00234     // If we don't have far detector spectra,
00235     if (!fFDDataPQ || !fFDDataNQ) {
00236         MSG("NuEZFitter",Msg::kWarning) << "No far detector spectra exists. Cannot rescale pot." << endl;
00237         return;
00238     } else {
00239         fFDDataPQ->ScaleToPot(POT);
00240         fFDDataNQ->ScaleToPot(POT);
00241     }
00242     
00243     // Far detector spectra have changed... rebuild
00244     Rebuild();
00245 }

void NuEZFitter::SetPoT (  Double_t  POT  ) 
 

Set's the value of the far detector POT.

Definition at line 212 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, MSG, Rebuild(), and NuMatrixSpectrum::ResetPoT().

00213 {
00214   // If we don't have far detector spectra,
00215   if (!fFDDataPQ || !fFDDataNQ) {
00216     MSG("NuEZFitter",Msg::kWarning) << "No far detector spectra exists. Creating blank spectra" << endl;
00217     // In case we somehow have one of the two
00218     if (fFDDataPQ) delete fFDDataPQ; fFDDataPQ = 0;
00219     if (fFDDataNQ) delete fFDDataPQ; fFDDataPQ = 0;
00220     
00221     fFDDataPQ = new NuMatrixSpectrum(POT);
00222     fFDDataNQ = new NuMatrixSpectrum(POT);
00223   } else {
00224     fFDDataPQ->ResetPoT(POT);
00225     fFDDataNQ->ResetPoT(POT);
00226   }
00227 
00228   // Far detector spectra have changed... rebuild
00229   Rebuild();
00230 }

void NuEZFitter::UseFDData (  const NuMatrixSpectrum &  fdNQdata,
const NuMatrixSpectrum &  fdPQdata
) 
 

Uses an exiting spectrum as far detector data.

Definition at line 165 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, and Rebuild().

00166 {
00167   // Clear any existing far detector spectra
00168   if (fFDDataPQ) delete fFDDataPQ; fFDDataPQ = 0;
00169   if (fFDDataNQ) delete fFDDataNQ; fFDDataNQ = 0;
00170   
00171   fFDDataNQ = new NuMatrixSpectrum(fdNQdata);
00172   fFDDataPQ = new NuMatrixSpectrum(fdPQdata);
00173   
00174   // Now rebuild the fitters
00175   Rebuild();
00176 }

void NuEZFitter::UseFDData (  std::string  FardetData  ) 
 

Uses an existing Far detector data spectrum, by loading from file.

Definition at line 151 of file NuEZFitter.cxx.

References fFDDataNQ, fFDDataPQ, MSG, and Rebuild().

Referenced by NuEZFitter().

00152 {
00153   // Clear any existing far detector spectra
00154   if (fFDDataPQ) delete fFDDataPQ; fFDDataPQ = 0;
00155   if (fFDDataNQ) delete fFDDataNQ; fFDDataNQ = 0;
00156 
00157   MSG("NuEZFitter",Msg::kInfo) << "Reading Far detector data " << FardetData << endl;
00158   fFDDataPQ = new NuMatrixSpectrum(FardetData, "RecoEnergyPQ_FD");
00159   fFDDataNQ = new NuMatrixSpectrum(FardetData, "RecoEnergy_FD");
00160   
00161   // Rebuild the internal objects
00162   Rebuild();
00163 }

void NuEZFitter::UseFDFakeData (  Double_t  POT,
std::string  FardetData,
std::string  FarDetTauData
) 
 

Uses an existing far detector fake-data spectrum.

Definition at line 178 of file NuEZFitter.cxx.

References NuMatrixSpectrum::Add(), fFDDataNQ, fFDDataPQ, MSG, NuMatrixSpectrum::PoT(), Rebuild(), and NuMatrixSpectrum::ScaleToPot().

00179 {
00180   // Clear any existing far detector spectra
00181   if (fFDDataPQ) delete fFDDataPQ; fFDDataPQ = 0;
00182   if (fFDDataNQ) delete fFDDataNQ; fFDDataNQ = 0;
00183   
00184   // Combine the fake far data and taudata then use as the spectrum
00185   MSG("NuEZFitter",Msg::kInfo) << "Reading Far detector fake-data " << FardetData << endl;
00186   MSG("NuEZFitter",Msg::kInfo) << "Reading Far detector fake-Tau data " << FarDetTauData << endl;
00187   
00188   NuMatrixSpectrum *fdDataPQ = new NuMatrixSpectrum(FardetData, "RecoEnergyPQ_ND");
00189   NuMatrixSpectrum *fdDataNQ = new NuMatrixSpectrum(FardetData, "RecoEnergy_ND");
00190   NuMatrixSpectrum *fdTauDataPQ = new NuMatrixSpectrum(FarDetTauData, "RecoEnergyPQ_ND");
00191   NuMatrixSpectrum *fdTauDataNQ = new NuMatrixSpectrum(FarDetTauData, "RecoEnergy_ND");
00192 
00193   if (fdDataPQ->PoT() < POT || fdTauDataPQ->PoT() < POT) {
00194     MSG("NuEZFitter",Msg::kWarning) << "Fake Data files have lower POT than requested - scaling data up" << endl;
00195   }
00196   
00197   fdDataPQ->ScaleToPot(POT);
00198   fdDataNQ->ScaleToPot(POT);
00199   fdTauDataPQ->ScaleToPot(POT);
00200   fdTauDataNQ->ScaleToPot(POT);
00201   
00202   fdDataPQ->Add(*fdTauDataPQ);
00203   fdDataNQ->Add(*fdTauDataNQ);
00204   
00205   fFDDataNQ = fdDataNQ;
00206   fFDDataPQ = fdDataPQ;
00207   
00208   // Rebuild the internal objects
00209   Rebuild();
00210 }

Member Data Documentation

NuMatrixSpectrum* NuEZFitter::fFDDataNQ [private]
 

NQ FD Spectrum.

Definition at line 133 of file NuEZFitter.h.

Referenced by Init(), MakeFDData(), NuEZFitter(), Rebuild(), ScaleToPot(), SetPoT(), UseFDData(), UseFDFakeData(), and ~NuEZFitter().

NuMatrixSpectrum* NuEZFitter::fFDDataPQ [private]
 

PQ FD Spectrum.

Definition at line 132 of file NuEZFitter.h.

Referenced by Init(), MakeFDData(), NuEZFitter(), Rebuild(), ScaleToPot(), SetPoT(), UseFDData(), UseFDFakeData(), and ~NuEZFitter().

NuMMParameters* NuEZFitter::fFitResult [private]
 

Definition at line 144 of file NuEZFitter.h.

Referenced by Fit(), Init(), NuEZFitter(), and ~NuEZFitter().

NuSystFitter* NuEZFitter::fFitter [private]
 

The fitter object.

Definition at line 141 of file NuEZFitter.h.

Referenced by Fit(), Init(), NuEZFitter(), Rebuild(), and ~NuEZFitter().

NuMMHelperCPT* NuEZFitter::fHelper [private]
 

The MM helper object.

Definition at line 136 of file NuEZFitter.h.

Referenced by Init(), Init_Helper_Near(), MakeFDData(), NuEZFitter(), Rebuild(), and ~NuEZFitter().

NuMatrixSpectrum* NuEZFitter::fNDDataNQ [private]
 

NQ ND Spectrum.

Definition at line 129 of file NuEZFitter.h.

Referenced by Init(), Init_Helper_Near(), MakeFDData(), NuEZFitter(), Rebuild(), and ~NuEZFitter().

NuMatrixSpectrum* NuEZFitter::fNDDataPQ [private]
 

PQ ND Spectrum.

Definition at line 128 of file NuEZFitter.h.

Referenced by Init(), Init_Helper_Near(), MakeFDData(), NuEZFitter(), Rebuild(), and ~NuEZFitter().

NuMMRunFC* NuEZFitter::fRun [private]
 

The run object that does the extrapolation.

Definition at line 139 of file NuEZFitter.h.

Referenced by Init(), NuEZFitter(), PredictionNQ(), PredictionPQ(), Rebuild(), and ~NuEZFitter().

The documentation for this class was generated from the following files:
Generated on Mon Feb 15 11:09:50 2010 for loon by  doxygen 1.3.9.1