diff --git a/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_nu.cxx b/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_nu.cxx
index cc4a31c..60b226b 100644
--- a/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_nu.cxx
+++ b/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_nu.cxx
@@ -1,119 +1,119 @@
 // Copyright 2016 L. Pickering, P Stowell, R. Terri, C. Wilkinson, C. Wret
 
 /*******************************************************************************
 *    This file is part of NUISANCE.
 *
 *    NUISANCE 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.
 *
 *    NUISANCE 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 NUISANCE.  If not, see <http://www.gnu.org/licenses/>.
 *******************************************************************************/
 
 #include "MINERvA_SignalDef.h"
 #include "MINERvA_CCDIS_XSec_1DEnu_nu.h"
 
 
 //********************************************************************
 MINERvA_CCDIS_XSec_1DEnu_nu::MINERvA_CCDIS_XSec_1DEnu_nu(std::string name, std::string inputfile, std::string type){
 //********************************************************************
 
   // Sample overview ---------------------------------------------------
   std::string descrip = "MINERvA_CCDIS_XSec_1DEnu_nu sample.";
 
   // Setup common settings
   fSettings = LoadSampleSettings(name, inputfile, type);
   fSettings.SetDescription(descrip);
   fSettings.SetXTitle("E_{#nu} (GeV)");
   fSettings.SetYTitle("d#sigma/dE_{#nu} (cm^{2}/GeV/nucleon)");
   fSettings.SetAllowedTypes("FIX/DIAG/MASK", "FIX/DIAG");
   fSettings.SetEnuRange(5.0, 50.0);
   fSettings.DefineAllowedTargets("Pb,Fe,C,H");
   fSettings.DefineAllowedSpecies("numu");
   fSettings.SetTitle("MINERvA_CCDIS_XSec_1DEnu_nu");
 
   target = "";
   if      (name.find("C12")  != std::string::npos) target =   "C12";
   else if (name.find("Fe56")  != std::string::npos) target =  "Fe56";
   else if (name.find("Pb208") != std::string::npos) target = "Pb208";
   if      (name.find("DEN")   != std::string::npos) target =   "CH";
   if (target == "") ERR(WRN) << "target " << target << " was not found!" << std::endl;
 
   FinaliseSampleSettings();
 
   // Scaling Setup ---------------------------------------------------
   // ScaleFactor automatically setup for DiffXSec/cm2/Nucleon
   fScaleFactor =  (GetEventHistogram()->Integral("width")*1E-38/(fNEvents+0.))/this->TotalIntegratedFlux();
 
   // Plot Setup -------------------------------------------------------
   double binsx[8] = {5.00, 10.00, 15.00, 20.00, 25.00, 30.00, 40.00, 50.00};
   CreateDataHistogram(7, binsx);
 
   // Final setup  ---------------------------------------------------
   FinaliseMeasurement();
 
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1DEnu_nu::FillEventVariables(FitEvent *event){
 //********************************************************************
 
   if (event->NumFSParticle(13) == 0)
     return;
 
   TLorentzVector Pnu  = event->GetNeutrinoIn()->fP;
   TLorentzVector Pmu  = event->GetHMFSParticle(13)->fP;
 
   // Need to calculate Q2 and W using the MINERvA method
   double Enu   = Pnu.E()/1000.;
   double Emu   = Pmu.E()/1000.;
   double theta = FitUtils::th(Pnu, Pmu);
 
   Q2    = 4*Enu*Emu*sin(theta/2.)*sin(theta/2.);
   W     = sqrt(PhysConst::mass_nucleon*PhysConst::mass_nucleon + 2*PhysConst::mass_nucleon*(Enu-Emu) - Q2);
   fXVar = Enu;
 
   return;
 }
 
 
 //********************************************************************
 bool MINERvA_CCDIS_XSec_1DEnu_nu::isSignal(FitEvent *event){
 //*******************************************************************
 
   if (!SignalDef::isCCINC(event, 14, this->EnuMin, this->EnuMax)) return false;
 
   // Restrict the phase space to theta < 17 degrees
   if (!SignalDef::IsRestrictedAngle(event, 14, 13, 17)) return false;
 
   if (Q2 < 1.0) return false;
   if ( W < 2.0) return false;
 
   return true;
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1DEnu_nu::ScaleEvents(){
 //********************************************************************
 
   // Get rid of this because it causes odd behaviour
-  //  Measurement1D::ScaleEvents();
+  Measurement1D::ScaleEvents();
 
-  this->fMCHist->Scale(this->fScaleFactor, "width");
+  // this->fMCHist->Scale(this->fScaleFactor, "width");
   
-  // Proper error scaling - ROOT Freaks out with xsec weights sometimes
-  for(int i=0; i<this->fMCStat->GetNbinsX();i++) {
+  // // Proper error scaling - ROOT Freaks out with xsec weights sometimes
+  // for(int i=0; i<this->fMCStat->GetNbinsX();i++) {
     
-    if (this->fMCStat->GetBinContent(i+1) != 0)
-      this->fMCHist->SetBinError(i+1, this->fMCHist->GetBinContent(i+1) * this->fMCStat->GetBinError(i+1) / this->fMCStat->GetBinContent(i+1) );
-    else this->fMCHist->SetBinError(i+1, this->fMCHist->Integral());
-  }
+  //   if (this->fMCStat->GetBinContent(i+1) != 0)
+  //     this->fMCHist->SetBinError(i+1, this->fMCHist->GetBinContent(i+1) * this->fMCStat->GetBinError(i+1) / this->fMCStat->GetBinContent(i+1) );
+  //   else this->fMCHist->SetBinError(i+1, this->fMCHist->Integral());
+  // }
 
 }
diff --git a/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_ratio.cxx b/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_ratio.cxx
index fc31aa8..8ea7398 100644
--- a/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_ratio.cxx
+++ b/src/MINERvA/MINERvA_CCDIS_XSec_1DEnu_ratio.cxx
@@ -1,122 +1,129 @@
 // Copyright 2016 L. Pickering, P Stowell, R. Terri, C. Wilkinson, C. Wret
 
 /*******************************************************************************
 *    This file is part of NUISANCE.
 *
 *    NUISANCE 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.
 *
 *    NUISANCE 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 NUISANCE.  If not, see <http://www.gnu.org/licenses/>.
 *******************************************************************************/
 
 #include "MINERvA_SignalDef.h"
 #include "MINERvA_CCDIS_XSec_1DEnu_ratio.h"
 
 
 //********************************************************************
 MINERvA_CCDIS_XSec_1DEnu_ratio::MINERvA_CCDIS_XSec_1DEnu_ratio(nuiskey samplekey) {
 //********************************************************************
 
   // Sample overview ---------------------------------------------------
   std::string descrip = "MINERvA_CCDIS_XSec_1DEnu_ratio sample.";
 
   // Setup common settings
   fSettings = LoadSampleSettings(samplekey);
   fSettings.SetDescription(descrip);
   fSettings.SetXTitle("E_{#nu} (GeV)");
   fSettings.SetYTitle(" d#sigma/dE_{#nu} (cm^{2}/GeV/nucleon)");
   fSettings.SetAllowedTypes("FIX/DIAG,FULL/MASK", "FIX/FULL");
   fSettings.SetEnuRange(5.0, 50.0);
 
   // CCQELike plot information
   fSettings.SetTitle("MINERvA_CCDIS_XSec_1DEnu_ratio");
   
   fIsRatio = true;
   nBins = 7;
 
   target  = "";
   if      (fSettings.Found("name", "C12"))   target = "C12";
   else if (fSettings.Found("name", "Fe56"))  target = "Fe56";
   else if (fSettings.Found("name", "Pb208")) target = "Pb208";
   else {
     ERR(FTL) << "target " << target << " was not found!" << std::endl;
     exit(-1);
   }
 
   fSettings.SetDataInput(  GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_Enu_data.csv" );
   fSettings.SetCovarInput( GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_Enu_stat.csv;" +
 			   GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_Enu_syst.csv");
   FinaliseSampleSettings();
 
   // Get parsed input files
   if (fSubInFiles.size() != 2) ERR(FTL) << "MINERvA CCDIS ratio requires input files in format: NUMERATOR;DENOMINATOR" << std::endl;
   std::string inFileNUM = fSubInFiles.at(0);
   std::string inFileDEN = fSubInFiles.at(1);
 
   // Scaling Setup ---------------------------------------------------
   // Ratio of sub classes so non needed
 
   // Plot Setup -------------------------------------------------------
   SetDataFromTextFile( fSettings.GetDataInput() );
   SetCovarFromMultipleTextFiles(fSettings.GetCovarInput());
+  
+  // Need to overlay the sqrt covariance diagonals onto the data histogram
+  StatUtils::SetDataErrorFromCov(fDataHist, fFullCovar);
+
+  // Need to scale the covariance by 1E-76... this cancels with the factor of 1E76 introduced in StatUtils::GetChi2FromCov
+  // Who says two wrongs don't make a right
+  ScaleCovar(1E76);
 
   // Setup Experiments  -------------------------------------------------------
   std::string type = samplekey.GetS("type");
   NUM  = new MINERvA_CCDIS_XSec_1DEnu_nu("MINERvA_CCDIS_XSec_1DEnu_"+target+"_CH_NUM", inFileNUM, type);
   DEN  = new MINERvA_CCDIS_XSec_1DEnu_nu("MINERvA_CCDIS_XSec_1DEnu_"+target+"_CH_DEN", inFileDEN, type);
   NUM  ->SetNoData();
   DEN  ->SetNoData();
 
   // Add to chain for processing
   this->fSubChain.clear();
   this->fSubChain.push_back(NUM);
   this->fSubChain.push_back(DEN);
 
   // Final setup  ---------------------------------------------------
   FinaliseMeasurement();
 
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1DEnu_ratio::MakePlots() {
 //********************************************************************
 
   UInt_t sample = 0;
   for (std::vector<MeasurementBase*>::const_iterator expIter = this->fSubChain.begin(); expIter != this->fSubChain.end(); expIter++) {
     MeasurementBase* exp = static_cast<MeasurementBase*>(*expIter);
 
     if      (sample == 0) this->NUM = static_cast<MINERvA_CCDIS_XSec_1DEnu_nu*>(exp);
     else if (sample == 1) this->DEN = static_cast<MINERvA_CCDIS_XSec_1DEnu_nu*>(exp);
     else break;
     sample++;
   }
 
   // Now make the ratio histogram
   TH1D* NUM_MC = (TH1D*)this->NUM->GetMCList().at(0)->Clone();
   TH1D* DEN_MC = (TH1D*)this->DEN->GetMCList().at(0)->Clone();
 
   for (int i = 0; i < nBins; ++i) {
     double binVal = 0;
     double binErr = 0;
 
     if (DEN_MC->GetBinContent(i + 1) && NUM_MC->GetBinContent(i + 1)) {
       binVal = NUM_MC->GetBinContent(i + 1) / DEN_MC->GetBinContent(i + 1);
       double fractErrNUM = NUM_MC->GetBinError(i + 1) / NUM_MC->GetBinContent(i + 1);
       double fractErrDEN = DEN_MC->GetBinError(i + 1) / DEN_MC->GetBinContent(i + 1);
       binErr = binVal * sqrt(fractErrNUM * fractErrNUM + fractErrDEN * fractErrDEN);
     }
 
     this->fMCHist->SetBinContent(i + 1, binVal);
     this->fMCHist->SetBinError(i + 1, binErr);
   }
 
   return;
 }
diff --git a/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_nu.cxx b/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_nu.cxx
index f1cbd25..b0ee35d 100644
--- a/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_nu.cxx
+++ b/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_nu.cxx
@@ -1,120 +1,120 @@
 // Copyright 2016 L. Pickering, P Stowell, R. Terri, C. Wilkinson, C. Wret
 
 /*******************************************************************************
 *    This file is part of NUISANCE.
 *
 *    NUISANCE 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.
 *
 *    NUISANCE 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 NUISANCE.  If not, see <http://www.gnu.org/licenses/>.
 *******************************************************************************/
 
 #include "MINERvA_SignalDef.h"
 #include "MINERvA_CCDIS_XSec_1Dx_nu.h"
 
 
 //********************************************************************
 MINERvA_CCDIS_XSec_1Dx_nu::MINERvA_CCDIS_XSec_1Dx_nu(std::string name, std::string inputfile, std::string type){
 //********************************************************************
 
   // Sample overview ---------------------------------------------------
   std::string descrip = "MINERvA_CCDIS_XSec_1Dx_nu sample.";
 
   // Setup common settings
   fSettings = LoadSampleSettings(name, inputfile, type);
   fSettings.SetDescription(descrip);
   fSettings.SetXTitle("Reconstructed Bjorken x");
   fSettings.SetYTitle("d#sigma/dx (cm^{2}/nucleon)");
   fSettings.SetAllowedTypes("FIX/DIAG/MASK", "FIX/DIAG");
   fSettings.SetEnuRange(5.0, 50.0);
   fSettings.DefineAllowedTargets("Pb,Fe,C,H");
   fSettings.DefineAllowedSpecies("numu");
   fSettings.SetTitle("MINERvA_CCDIS_XSec_1Dx_nu");
 
   target = "";
   if      (name.find("C12")  != std::string::npos) target =   "C12";
   else if (name.find("Fe56")  != std::string::npos) target =  "Fe56";
   else if (name.find("Pb208") != std::string::npos) target = "Pb208";
   if      (name.find("DEN")   != std::string::npos) target =   "CH";
   if (target == "") ERR(WRN) << "target " << target << " was not found!" << std::endl;
 
   FinaliseSampleSettings();
 
   // Scaling Setup ---------------------------------------------------
   // ScaleFactor automatically setup for DiffXSec/cm2/Nucleon
   fScaleFactor =  (GetEventHistogram()->Integral("width")*1E-38/(fNEvents+0.))/this->TotalIntegratedFlux();
 
   // Plot Setup -------------------------------------------------------
   double binsx[6] = {0.00, 0.10, 0.20, 0.30, 0.40, 0.75};
   CreateDataHistogram(5, binsx);
 
   FinaliseMeasurement();
 
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1Dx_nu::FillEventVariables(FitEvent *event){
 //********************************************************************
 
   if (event->NumFSParticle(13) == 0)
     return;
 
   TLorentzVector Pnu  = event->GetNeutrinoIn()->fP;
   TLorentzVector Pmu  = event->GetHMFSParticle(13)->fP;
 
   // Need to calculate Q2 and W using the MINERvA method
   double Enu   = Pnu.E()/1000.;
   double Emu   = Pmu.E()/1000.;
   double theta = FitUtils::th(Pnu, Pmu);
 
   Q2    = 4*Enu*Emu*sin(theta/2.)*sin(theta/2.);
   W     = sqrt(PhysConst::mass_nucleon*PhysConst::mass_nucleon + 2*PhysConst::mass_nucleon*(Enu-Emu) - Q2);
   fXVar = Q2/(2*PhysConst::mass_nucleon*(Enu-Emu));
 
   return;
 }
 
 
 
 //********************************************************************
 bool MINERvA_CCDIS_XSec_1Dx_nu::isSignal(FitEvent *event){
 //*******************************************************************
 
   // Only look at numu events
   if (!SignalDef::isCCINC(event, 14, EnuMin, EnuMax)) return false;
 
   // Restrict the phase space to theta < 17 degrees
   if (!SignalDef::IsRestrictedAngle(event, 14, 13, 17)) return false;
 
   if (Q2 < 1.0) return false;
   if ( W < 2.0) return false;
 
   return true;
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1Dx_nu::ScaleEvents(){
 //********************************************************************
 
-  //Measurement1D::ScaleEvents();
-  this->fDataHist = (TH1D*)this->GetMCList().at(0)->Clone();
+  Measurement1D::ScaleEvents();
+  // this->fDataHist = (TH1D*)this->GetMCList().at(0)->Clone();
   
-  this->fMCHist->Scale(this->fScaleFactor, "width");
+  // this->fMCHist->Scale(this->fScaleFactor, "width");
   
-  // Proper error scaling - ROOT Freaks out with xsec weights sometimes
-  for(int i=0; i<this->fMCStat->GetNbinsX();i++) {
+  // // Proper error scaling - ROOT Freaks out with xsec weights sometimes
+  // for(int i=0; i<this->fMCStat->GetNbinsX();i++) {
     
-    if (this->fMCStat->GetBinContent(i+1) != 0)
-      this->fMCHist->SetBinError(i+1, this->fMCHist->GetBinContent(i+1) * this->fMCStat->GetBinError(i+1) / this->fMCStat->GetBinContent(i+1) );
-    else this->fMCHist->SetBinError(i+1, this->fMCHist->Integral());
-  }
+  //   if (this->fMCStat->GetBinContent(i+1) != 0)
+  //     this->fMCHist->SetBinError(i+1, this->fMCHist->GetBinContent(i+1) * this->fMCStat->GetBinError(i+1) / this->fMCStat->GetBinContent(i+1) );
+  //   else this->fMCHist->SetBinError(i+1, this->fMCHist->Integral());
+  // }
 
 }
diff --git a/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_ratio.cxx b/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_ratio.cxx
index 74acc1a..0981a58 100644
--- a/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_ratio.cxx
+++ b/src/MINERvA/MINERvA_CCDIS_XSec_1Dx_ratio.cxx
@@ -1,121 +1,128 @@
 // Copyright 2016 L. Pickering, P Stowell, R. Terri, C. Wilkinson, C. Wret
 
 /*******************************************************************************
 *    This file is part of NUISANCE.
 *
 *    NUISANCE 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.
 *
 *    NUISANCE 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 NUISANCE.  If not, see <http://www.gnu.org/licenses/>.
 *******************************************************************************/
 
 #include "MINERvA_SignalDef.h"
 #include "MINERvA_CCDIS_XSec_1Dx_ratio.h"
 
 //********************************************************************
 MINERvA_CCDIS_XSec_1Dx_ratio::MINERvA_CCDIS_XSec_1Dx_ratio(nuiskey samplekey) {
 //********************************************************************
 
   // Sample overview ---------------------------------------------------
   std::string descrip = "MINERvA_CCDIS_XSec_1Dx_ratio sample.";
 
   // Setup common settings
   fSettings = LoadSampleSettings(samplekey);
   fSettings.SetDescription(descrip);
   fSettings.SetXTitle("Reconstructed Bjorken x");
   fSettings.SetYTitle("d#sigma/dx (cm^{2}/nucleon)");
   fSettings.SetAllowedTypes("FIX/DIAG,FULL/MASK", "FIX/FULL");
   fSettings.SetEnuRange(5.0, 50.0);
 
   // CCQELike plot information
   fSettings.SetTitle("MINERvA_CCDIS_XSec_1Dx_ratio");
   fIsRatio = true;
   nBins = 5;
 
   target  = "";
   if      (fSettings.Found("name", "C12"))   target =   "C12";
   else if (fSettings.Found("name", "Fe56"))  target =   "Fe56";
   else if (fSettings.Found("name", "Pb208")) target =   "Pb208";
   else {
     ERR(FTL) << "target " << target << " was not found!" << std::endl;
     exit(-1);
   }
 
   fSettings.SetDataInput(  GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_x_data.csv" );
   fSettings.SetCovarInput( GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_x_stat.csv;" +
                            GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCDIS/CCDIS_"+target+"_CH_ratio_x_syst.csv");
   FinaliseSampleSettings();
 
 
   // Get parsed input files
   if (fSubInFiles.size() != 2) ERR(FTL) << "MINERvA CCDIS ratio requires input files in format: NUMERATOR;DENOMINATOR" << std::endl;
   std::string inFileNUM = fSubInFiles.at(0);
   std::string inFileDEN = fSubInFiles.at(1);
 
   // Scaling Setup ---------------------------------------------------
   // Ratio of sub classes so non needed
 
   // Plot Setup -------------------------------------------------------
   SetDataFromTextFile( fSettings.GetDataInput() );
   SetCovarFromMultipleTextFiles(fSettings.GetCovarInput());
 
+  // Need to overlay the sqrt covariance diagonals (*1E-38) onto the data histogram
+  StatUtils::SetDataErrorFromCov(fDataHist, fFullCovar);
+
+  // Need to scale the covariance by 1E-76... this cancels with the factor of 1E76 introduced in StatUtils::GetChi2FromCov
+  // Who says two wrongs don't make a right
+  ScaleCovar(1E76);
+
   // Setup Experiments  -------------------------------------------------------
   std::string type = samplekey.GetS("type");
   NUM  = new MINERvA_CCDIS_XSec_1Dx_nu("MINERvA_CCDIS_XSec_1Dx_" + target + "_CH_NUM", inFileNUM, type);
   DEN  = new MINERvA_CCDIS_XSec_1Dx_nu("MINERvA_CCDIS_XSec_1Dx_" + target + "_CH_DEN", inFileDEN, type);
   NUM  ->SetNoData();
   DEN  ->SetNoData();
 
   // Add to chain for processing
   this->fSubChain.clear();
   this->fSubChain.push_back(NUM);
   this->fSubChain.push_back(DEN);
 
   // Final setup  ---------------------------------------------------
   FinaliseMeasurement();
 
 };
 
 //********************************************************************
 void MINERvA_CCDIS_XSec_1Dx_ratio::MakePlots() {
 //********************************************************************
 
   UInt_t sample = 0;
   for (std::vector<MeasurementBase*>::const_iterator expIter = this->fSubChain.begin(); expIter != this->fSubChain.end(); expIter++) {
     MeasurementBase* exp = static_cast<MeasurementBase*>(*expIter);
 
     if      (sample == 0) this->NUM = static_cast<MINERvA_CCDIS_XSec_1Dx_nu*>(exp);
     else if (sample == 1) this->DEN = static_cast<MINERvA_CCDIS_XSec_1Dx_nu*>(exp);
     else break;
     sample++;
   }
 
   // Now make the ratio histogram
   TH1D* NUM_MC = (TH1D*)this->NUM->GetMCList().at(0)->Clone();
   TH1D* DEN_MC = (TH1D*)this->DEN->GetMCList().at(0)->Clone();
 
   for (int i = 0; i < nBins; ++i) {
     double binVal = 0;
     double binErr = 0;
 
     if (DEN_MC->GetBinContent(i + 1) && NUM_MC->GetBinContent(i + 1)) {
       binVal = NUM_MC->GetBinContent(i + 1) / DEN_MC->GetBinContent(i + 1);
       double fractErrNUM = NUM_MC->GetBinError(i + 1) / NUM_MC->GetBinContent(i + 1);
       double fractErrDEN = DEN_MC->GetBinError(i + 1) / DEN_MC->GetBinContent(i + 1);
       binErr = binVal * sqrt(fractErrNUM * fractErrNUM + fractErrDEN * fractErrDEN);
     }
 
     this->fMCHist->SetBinContent(i + 1, binVal);
     this->fMCHist->SetBinError(i + 1, binErr);
   }
 
   return;
 }