No OneTemporary
Actions

Size

67 KB

Subscribers

None

View Options

	diff --git a/src/FCN/SampleList.cxx b/src/FCN/SampleList.cxx
	index a5ebb36..3db2359 100644
	--- a/src/FCN/SampleList.cxx
	+++ b/src/FCN/SampleList.cxx
	@@ -1,1006 +1,1016 @@
	#include "SampleList.h"

	#ifndef __NO_ANL__
	#include "ANL_CCQE_Evt_1DQ2_nu.h"
	#include "ANL_CCQE_XSec_1DEnu_nu.h"

	// ANL CC1ppip
	#include "ANL_CC1ppip_Evt_1DQ2_nu.h"
	#include "ANL_CC1ppip_Evt_1DcosmuStar_nu.h"
	#include "ANL_CC1ppip_Evt_1DcosmuStar_nu.h"
	#include "ANL_CC1ppip_Evt_1DcosthAdler_nu.h"
	#include "ANL_CC1ppip_Evt_1Dphi_nu.h"
	#include "ANL_CC1ppip_Evt_1Dppi_nu.h"
	#include "ANL_CC1ppip_Evt_1Dthpr_nu.h"
	#include "ANL_CC1ppip_XSec_1DEnu_nu.h"
	#include "ANL_CC1ppip_XSec_1DQ2_nu.h"
	// ANL CC1npip
	#include "ANL_CC1npip_Evt_1DQ2_nu.h"
	#include "ANL_CC1npip_Evt_1DcosmuStar_nu.h"
	#include "ANL_CC1npip_Evt_1Dppi_nu.h"
	#include "ANL_CC1npip_XSec_1DEnu_nu.h"
	// ANL CC1pi0
	#include "ANL_CC1pi0_Evt_1DQ2_nu.h"
	#include "ANL_CC1pi0_Evt_1DcosmuStar_nu.h"
	#include "ANL_CC1pi0_XSec_1DEnu_nu.h"
	// ANL NC1npip (mm, exotic!)
	#include "ANL_NC1npip_Evt_1Dppi_nu.h"
	// ANL NC1ppim (mm, exotic!)
	#include "ANL_NC1ppim_Evt_1DcosmuStar_nu.h"
	#include "ANL_NC1ppim_XSec_1DEnu_nu.h"
	// ANL CC2pi 1pim1pip (mm, even more exotic!)
	#include "ANL_CC2pi_1pim1pip_Evt_1Dpmu_nu.h"
	#include "ANL_CC2pi_1pim1pip_Evt_1Dppim_nu.h"
	#include "ANL_CC2pi_1pim1pip_Evt_1Dppip_nu.h"
	#include "ANL_CC2pi_1pim1pip_Evt_1Dpprot_nu.h"
	#include "ANL_CC2pi_1pim1pip_XSec_1DEnu_nu.h"
	// ANL CC2pi 1pip1pip (mm, even more exotic!)
	#include "ANL_CC2pi_1pip1pip_Evt_1Dpmu_nu.h"
	#include "ANL_CC2pi_1pip1pip_Evt_1Dpneut_nu.h"
	#include "ANL_CC2pi_1pip1pip_Evt_1DppipHigh_nu.h"
	#include "ANL_CC2pi_1pip1pip_Evt_1DppipLow_nu.h"
	#include "ANL_CC2pi_1pip1pip_XSec_1DEnu_nu.h"
	// ANL CC2pi 1pip1pi0 (mm, even more exotic!)
	#include "ANL_CC2pi_1pip1pi0_Evt_1Dpmu_nu.h"
	#include "ANL_CC2pi_1pip1pi0_Evt_1Dppi0_nu.h"
	#include "ANL_CC2pi_1pip1pi0_Evt_1Dppip_nu.h"
	#include "ANL_CC2pi_1pip1pi0_Evt_1Dpprot_nu.h"
	#include "ANL_CC2pi_1pip1pi0_XSec_1DEnu_nu.h"
	#endif

	#ifndef __NO_ArgoNeuT__
	// ArgoNeuT CC-inclusive
	#include "ArgoNeuT_CCInc_XSec_1Dpmu_antinu.h"
	#include "ArgoNeuT_CCInc_XSec_1Dpmu_nu.h"
	#include "ArgoNeuT_CCInc_XSec_1Dthetamu_antinu.h"
	#include "ArgoNeuT_CCInc_XSec_1Dthetamu_nu.h"
	#endif

	#ifndef __NO_BNL__
	// BNL CCQE
	#include "BNL_CCQE_Evt_1DQ2_nu.h"
	#include "BNL_CCQE_XSec_1DEnu_nu.h"
	// BNL CC1ppip
	#include "BNL_CC1ppip_Evt_1DQ2_nu.h"
	#include "BNL_CC1ppip_Evt_1DQ2_nu.h"
	#include "BNL_CC1ppip_Evt_1DcosthAdler_nu.h"
	#include "BNL_CC1ppip_Evt_1Dphi_nu.h"
	#include "BNL_CC1ppip_XSec_1DEnu_nu.h"
	// BNL CC1npip
	#include "BNL_CC1npip_Evt_1DQ2_nu.h"
	#include "BNL_CC1npip_XSec_1DEnu_nu.h"
	// BNL CC1pi0
	#include "BNL_CC1pi0_Evt_1DQ2_nu.h"
	#include "BNL_CC1pi0_XSec_1DEnu_nu.h"
	#endif

	#ifndef __NO_FNAL__
	// FNAL CCQE
	#include "FNAL_CCQE_Evt_1DQ2_nu.h"
	// FNAL CC1ppip
	#include "FNAL_CC1ppip_Evt_1DQ2_nu.h"
	#include "FNAL_CC1ppip_XSec_1DEnu_nu.h"
	#include "FNAL_CC1ppip_XSec_1DQ2_nu.h"
	// FNAL CC1ppim
	#include "FNAL_CC1ppim_XSec_1DEnu_antinu.h"
	#endif

	#ifndef __NO_BEBC__
	// BEBC CCQE
	#include "BEBC_CCQE_XSec_1DQ2_nu.h"
	// BEBC CC1ppip
	#include "BEBC_CC1ppip_XSec_1DEnu_nu.h"
	#include "BEBC_CC1ppip_XSec_1DQ2_nu.h"
	// BEBC CC1npip
	#include "BEBC_CC1npip_XSec_1DEnu_nu.h"
	#include "BEBC_CC1npip_XSec_1DQ2_nu.h"
	// BEBC CC1pi0
	#include "BEBC_CC1pi0_XSec_1DEnu_nu.h"
	#include "BEBC_CC1pi0_XSec_1DQ2_nu.h"
	// BEBC CC1npim
	#include "BEBC_CC1npim_XSec_1DEnu_antinu.h"
	#include "BEBC_CC1npim_XSec_1DQ2_antinu.h"
	// BEBC CC1ppim
	#include "BEBC_CC1ppim_XSec_1DEnu_antinu.h"
	#include "BEBC_CC1ppim_XSec_1DQ2_antinu.h"
	#endif

	#ifndef __NO_GGM__
	// GGM CC1ppip
	#include "GGM_CC1ppip_Evt_1DQ2_nu.h"
	#include "GGM_CC1ppip_XSec_1DEnu_nu.h"
	#endif

	#ifndef __NO_MiniBooNE__
	// MiniBooNE CCQE
	#include "MiniBooNE_CCQE_XSec_1DQ2_antinu.h"
	#include "MiniBooNE_CCQE_XSec_1DQ2_nu.h"
	#include "MiniBooNE_CCQE_XSec_2DTcos_antinu.h"
	#include "MiniBooNE_CCQE_XSec_2DTcos_antinu.h"
	#include "MiniBooNE_CCQE_XSec_2DTcos_nu.h"

	// MiniBooNE CC1pi+ 1D
	#include "MiniBooNE_CC1pip_XSec_1DEnu_nu.h"
	#include "MiniBooNE_CC1pip_XSec_1DQ2_nu.h"
	#include "MiniBooNE_CC1pip_XSec_1DTpi_nu.h"
	#include "MiniBooNE_CC1pip_XSec_1DTu_nu.h"
	// MiniBooNE CC1pi+ 2D
	#include "MiniBooNE_CC1pip_XSec_2DQ2Enu_nu.h"
	#include "MiniBooNE_CC1pip_XSec_2DTpiCospi_nu.h"
	#include "MiniBooNE_CC1pip_XSec_2DTpiEnu_nu.h"
	#include "MiniBooNE_CC1pip_XSec_2DTuCosmu_nu.h"
	#include "MiniBooNE_CC1pip_XSec_2DTuEnu_nu.h"

	// MiniBooNE CC1pi0
	#include "MiniBooNE_CC1pi0_XSec_1DEnu_nu.h"
	#include "MiniBooNE_CC1pi0_XSec_1DQ2_nu.h"
	#include "MiniBooNE_CC1pi0_XSec_1DTu_nu.h"
	#include "MiniBooNE_CC1pi0_XSec_1Dcosmu_nu.h"
	#include "MiniBooNE_CC1pi0_XSec_1Dcospi0_nu.h"
	#include "MiniBooNE_CC1pi0_XSec_1Dppi0_nu.h"
	#include "MiniBooNE_NC1pi0_XSec_1Dcospi0_antinu.h"
	#include "MiniBooNE_NC1pi0_XSec_1Dcospi0_nu.h"
	#include "MiniBooNE_NC1pi0_XSec_1Dppi0_antinu.h"
	#include "MiniBooNE_NC1pi0_XSec_1Dppi0_nu.h"

	// MiniBooNE NC1pi0
	//#include "MiniBooNE_NCpi0_XSec_1Dppi0_nu.h"

	// MiniBooNE NCEL
	#include "MiniBooNE_NCEL_XSec_Treco_nu.h"
	#endif

	#ifndef __NO_MINERvA__
	// MINERvA CCQE
	#include "MINERvA_CCQE_XSec_1DQ2_antinu.h"
	#include "MINERvA_CCQE_XSec_1DQ2_joint.h"
	#include "MINERvA_CCQE_XSec_1DQ2_nu.h"

	// MINERvA CC0pi
	#include "MINERvA_CC0pi_XSec_1DEe_nue.h"
	#include "MINERvA_CC0pi_XSec_1DQ2_nu_proton.h"
	#include "MINERvA_CC0pi_XSec_1DQ2_nue.h"
	#include "MINERvA_CC0pi_XSec_1DThetae_nue.h"

	// MINERvA CC1pi+
	#include "MINERvA_CC1pip_XSec_1DTpi_20deg_nu.h"
	#include "MINERvA_CC1pip_XSec_1DTpi_nu.h"
	#include "MINERvA_CC1pip_XSec_1Dth_20deg_nu.h"
	#include "MINERvA_CC1pip_XSec_1Dth_nu.h"

	// MINERvA CCNpi+
	#include "MINERvA_CCNpip_XSec_1DEnu_nu.h"
	#include "MINERvA_CCNpip_XSec_1DQ2_nu.h"
	#include "MINERvA_CCNpip_XSec_1DTpi_nu.h"
	#include "MINERvA_CCNpip_XSec_1Dpmu_nu.h"
	#include "MINERvA_CCNpip_XSec_1Dth_nu.h"
	#include "MINERvA_CCNpip_XSec_1Dthmu_nu.h"

	// MINERvA CC1pi0
	#include "MINERvA_CC1pi0_XSec_1DEnu_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1DQ2_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1DTpi0_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1Dpmu_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1Dppi0_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1Dth_antinu.h"
	#include "MINERvA_CC1pi0_XSec_1Dthmu_antinu.h"

	// MINERvA CCINC
	#include "MINERvA_CCinc_XSec_1DEnu_ratio.h"
	#include "MINERvA_CCinc_XSec_1Dx_ratio.h"
	#include "MINERvA_CCinc_XSec_2DEavq3_nu.h"

	// MINERvA CCDIS
	#include "MINERvA_CCDIS_XSec_1DEnu_ratio.h"
	#include "MINERvA_CCDIS_XSec_1Dx_ratio.h"

	// MINERvA CCCOH pion
	#include "MINERvA_CCCOHPI_XSec_1DEnu_antinu.h"
	#include "MINERvA_CCCOHPI_XSec_1DEnu_nu.h"
	#include "MINERvA_CCCOHPI_XSec_1DEpi_antinu.h"
	+#include "MINERvA_CCCOHPI_XSec_1DQ2_antinu.h"
	+
	#include "MINERvA_CCCOHPI_XSec_1DEpi_nu.h"
	#include "MINERvA_CCCOHPI_XSec_1Dth_antinu.h"
	#include "MINERvA_CCCOHPI_XSec_1Dth_nu.h"
	+#include "MINERvA_CCCOHPI_XSec_1DQ2_nu.h"

	#include "MINERvA_CC0pi_XSec_1DQ2_TgtRatio_nu.h"
	#include "MINERvA_CC0pi_XSec_1DQ2_Tgt_nu.h"
	#include "MINERvA_CC1pip_XSec_1D_2017Update.h"

	+
	+
	#endif

	#ifndef __NO_T2K__
	// T2K CC0pi
	#include "T2K_CC0pi_XSec_2DPcos_nu.h"

	// T2K CC1pi+ on CH
	#include "T2K_CC1pip_CH_XSec_1DQ2_nu.h"
	#include "T2K_CC1pip_CH_XSec_1DWrec_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dpmu_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dppi_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dq3_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dthmupi_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dthpi_nu.h"
	#include "T2K_CC1pip_CH_XSec_1Dthq3pi_nu.h"

	// T2K CC1pi+ on H2O
	#include "T2K_CC1pip_H2O_XSec_1DEnuDelta_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1DEnuMB_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1Dcosmu_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1Dcosmupi_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1Dcospi_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1Dpmu_nu.h"
	#include "T2K_CC1pip_H2O_XSec_1Dppi_nu.h"

	// T2K STV CC0pi
	#include "T2K_CC0pinp_STV_XSec_1Ddpt_nu.h"

	#include "T2K_CC0pi_XSec_2DPcos_nu_nonuniform.h"

	#endif

	#ifndef __NO_SciBooNE__

	// SciBooNE COH studies
	#include "SciBooNE_CCCOH_1TRK_1DQ2_nu.h"
	#include "SciBooNE_CCCOH_MuPiNoVA_1DQ2_nu.h"
	#include "SciBooNE_CCCOH_MuPiNoVA_1Dthetapi_nu.h"
	#include "SciBooNE_CCCOH_MuPiNoVA_1Dthetapr_nu.h"
	#include "SciBooNE_CCCOH_MuPiVA_1DQ2_nu.h"
	#include "SciBooNE_CCCOH_MuPr_1DQ2_nu.h"
	#include "SciBooNE_CCCOH_STOPFINAL_1DQ2_nu.h"
	#include "SciBooNE_CCCOH_STOP_NTrks_nu.h"
	#endif

	#ifndef __NO_K2K__
	// K2K NC1pi0
	#include "K2K_NC1pi0_Evt_1Dppi0_nu.h"
	#endif

	// MC Studies
	#include "ExpMultDist_CCQE_XSec_1DVar_FakeStudy.h"
	#include "ExpMultDist_CCQE_XSec_2DVar_FakeStudy.h"
	#include "MCStudy_CCQEHistograms.h"

	#include "GenericFlux_Tester.h"
	#include "GenericFlux_Vectors.h"

	#include "ElectronFlux_FlatTree.h"
	#include "ElectronScattering_DurhamData.h"
	#include "MCStudy_KaonPreSelection.h"
	#include "MCStudy_MuonValidation.h"

	#include "OfficialNIWGPlots.h"
	#include "T2K2017_FakeData.h"

	#include "Simple_Osc.h"

	#include "FitWeight.h"

	#include "NuisConfig.h"
	#include "NuisKey.h"



	//! Functions to make it easier for samples to be created and handled.
	namespace SampleUtils {

	//! Create a given sample given its name, file, type, fakdata(fkdt) file and the
	//! current rw engine and push it back into the list fChain.
	MeasurementBase* CreateSample(std::string name, std::string file,
	std::string type, std::string fkdt,
	FitWeight* rw) {
	nuiskey samplekey = Config::CreateKey("sample");
	samplekey.AddS("name", name);
	samplekey.AddS("input", file);
	samplekey.AddS("type", type);

	return CreateSample(samplekey);
	}

	MeasurementBase* CreateSample(nuiskey samplekey) {
	FitWeight* rw = FitBase::GetRW();
	std::string name = samplekey.GetS("name");
	std::string file = samplekey.GetS("input");
	std::string type = samplekey.GetS("type");
	std::string fkdt = "";

	/*
	ANL CCQE Samples
	*/

	#ifndef __NO_ANL__
	if (!name.compare("ANL_CCQE_XSec_1DEnu_nu") \|\|
	!name.compare("ANL_CCQE_XSec_1DEnu_nu_PRD26") \|\|
	!name.compare("ANL_CCQE_XSec_1DEnu_nu_PRL31") \|\|
	!name.compare("ANL_CCQE_XSec_1DEnu_nu_PRD16")) {
	return (new ANL_CCQE_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CCQE_Evt_1DQ2_nu") \|\|
	!name.compare("ANL_CCQE_Evt_1DQ2_nu_PRL31") \|\|
	!name.compare("ANL_CCQE_Evt_1DQ2_nu_PRD26") \|\|
	!name.compare("ANL_CCQE_Evt_1DQ2_nu_PRD16")) {
	return (new ANL_CCQE_Evt_1DQ2_nu(samplekey));
	/*
	ANL CC1ppip samples
	*/
	} else if (!name.compare("ANL_CC1ppip_XSec_1DEnu_nu") \|\|
	!name.compare("ANL_CC1ppip_XSec_1DEnu_nu_W14Cut") \|\|
	!name.compare("ANL_CC1ppip_XSec_1DEnu_nu_Uncorr") \|\|
	!name.compare("ANL_CC1ppip_XSec_1DEnu_nu_W14Cut_Uncorr") \|\|
	!name.compare("ANL_CC1ppip_XSec_1DEnu_nu_W16Cut_Uncorr")) {
	return (new ANL_CC1ppip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_XSec_1DQ2_nu")) {
	return (new ANL_CC1ppip_XSec_1DQ2_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1DQ2_nu") \|\|
	!name.compare("ANL_CC1ppip_Evt_1DQ2_nu_W14Cut")) {
	return (new ANL_CC1ppip_Evt_1DQ2_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1Dppi_nu")) {
	return (new ANL_CC1ppip_Evt_1Dppi_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1Dthpr_nu")) {
	return (new ANL_CC1ppip_Evt_1Dthpr_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1DcosmuStar_nu")) {
	return (new ANL_CC1ppip_Evt_1DcosmuStar_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1DcosthAdler_nu")) {
	return (new ANL_CC1ppip_Evt_1DcosthAdler_nu(samplekey));
	} else if (!name.compare("ANL_CC1ppip_Evt_1Dphi_nu")) {
	return (new ANL_CC1ppip_Evt_1Dphi_nu(samplekey));
	/*
	ANL CC1npip sample
	*/
	} else if (!name.compare("ANL_CC1npip_XSec_1DEnu_nu") \|\|
	!name.compare("ANL_CC1npip_XSec_1DEnu_nu_W14Cut") \|\|
	!name.compare("ANL_CC1npip_XSec_1DEnu_nu_Uncorr") \|\|
	!name.compare("ANL_CC1npip_XSec_1DEnu_nu_W14Cut_Uncorr") \|\|
	!name.compare("ANL_CC1npip_XSec_1DEnu_nu_W16Cut_Uncorr")) {
	return (new ANL_CC1npip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC1npip_Evt_1DQ2_nu") \|\|
	!name.compare("ANL_CC1npip_Evt_1DQ2_nu_W14Cut")) {
	return (new ANL_CC1npip_Evt_1DQ2_nu(samplekey));
	} else if (!name.compare("ANL_CC1npip_Evt_1Dppi_nu")) {
	return (new ANL_CC1npip_Evt_1Dppi_nu(samplekey));
	} else if (!name.compare("ANL_CC1npip_Evt_1DcosmuStar_nu")) {
	return (new ANL_CC1npip_Evt_1DcosmuStar_nu(samplekey));
	/*
	ANL CC1pi0 sample
	*/
	} else if (!name.compare("ANL_CC1pi0_XSec_1DEnu_nu") \|\|
	!name.compare("ANL_CC1pi0_XSec_1DEnu_nu_W14Cut") \|\|
	!name.compare("ANL_CC1pi0_XSec_1DEnu_nu_Uncorr") \|\|
	!name.compare("ANL_CC1pi0_XSec_1DEnu_nu_W14Cut_Uncorr") \|\|
	!name.compare("ANL_CC1pi0_XSec_1DEnu_nu_W16Cut_Uncorr")) {
	return (new ANL_CC1pi0_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC1pi0_Evt_1DQ2_nu") \|\|
	!name.compare("ANL_CC1pi0_Evt_1DQ2_nu_W14Cut")) {
	return (new ANL_CC1pi0_Evt_1DQ2_nu(samplekey));
	} else if (!name.compare("ANL_CC1pi0_Evt_1DcosmuStar_nu")) {
	return (new ANL_CC1pi0_Evt_1DcosmuStar_nu(samplekey));
	/*
	ANL NC1npip sample
	*/
	} else if (!name.compare("ANL_NC1npip_Evt_1Dppi_nu")) {
	return (new ANL_NC1npip_Evt_1Dppi_nu(samplekey));
	/*
	ANL NC1ppim sample
	*/
	} else if (!name.compare("ANL_NC1ppim_XSec_1DEnu_nu")) {
	return (new ANL_NC1ppim_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_NC1ppim_Evt_1DcosmuStar_nu")) {
	return (new ANL_NC1ppim_Evt_1DcosmuStar_nu(samplekey));
	/*
	ANL CC2pi sample
	*/
	} else if (!name.compare("ANL_CC2pi_1pim1pip_XSec_1DEnu_nu")) {
	return (new ANL_CC2pi_1pim1pip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pim1pip_Evt_1Dpmu_nu")) {
	return (new ANL_CC2pi_1pim1pip_Evt_1Dpmu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pim1pip_Evt_1Dppip_nu")) {
	return (new ANL_CC2pi_1pim1pip_Evt_1Dppip_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pim1pip_Evt_1Dppim_nu")) {
	return (new ANL_CC2pi_1pim1pip_Evt_1Dppim_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pim1pip_Evt_1Dpprot_nu")) {
	return (new ANL_CC2pi_1pim1pip_Evt_1Dpprot_nu(samplekey));

	} else if (!name.compare("ANL_CC2pi_1pip1pip_XSec_1DEnu_nu")) {
	return (new ANL_CC2pi_1pip1pip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pip_Evt_1Dpmu_nu")) {
	return (new ANL_CC2pi_1pip1pip_Evt_1Dpmu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pip_Evt_1Dpneut_nu")) {
	return (new ANL_CC2pi_1pip1pip_Evt_1Dpneut_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pip_Evt_1DppipHigh_nu")) {
	return (new ANL_CC2pi_1pip1pip_Evt_1DppipHigh_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pip_Evt_1DppipLow_nu")) {
	return (new ANL_CC2pi_1pip1pip_Evt_1DppipLow_nu(samplekey));

	} else if (!name.compare("ANL_CC2pi_1pip1pi0_XSec_1DEnu_nu")) {
	return (new ANL_CC2pi_1pip1pi0_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pi0_Evt_1Dpmu_nu")) {
	return (new ANL_CC2pi_1pip1pi0_Evt_1Dpmu_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pi0_Evt_1Dppip_nu")) {
	return (new ANL_CC2pi_1pip1pi0_Evt_1Dppip_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pi0_Evt_1Dppi0_nu")) {
	return (new ANL_CC2pi_1pip1pi0_Evt_1Dppi0_nu(samplekey));
	} else if (!name.compare("ANL_CC2pi_1pip1pi0_Evt_1Dpprot_nu")) {
	return (new ANL_CC2pi_1pip1pi0_Evt_1Dpprot_nu(samplekey));

	/*
	ArgoNeut Samples
	*/
	} else
	#endif
	#ifndef __NO_ArgoNeuT__
	if (!name.compare("ArgoNeuT_CCInc_XSec_1Dpmu_antinu")) {
	return (new ArgoNeuT_CCInc_XSec_1Dpmu_antinu(samplekey));
	} else if (!name.compare("ArgoNeuT_CCInc_XSec_1Dpmu_nu")) {
	return (new ArgoNeuT_CCInc_XSec_1Dpmu_nu(samplekey));
	} else if (!name.compare("ArgoNeuT_CCInc_XSec_1Dthetamu_antinu")) {
	return (new ArgoNeuT_CCInc_XSec_1Dthetamu_antinu(samplekey));
	} else if (!name.compare("ArgoNeuT_CCInc_XSec_1Dthetamu_nu")) {
	return (new ArgoNeuT_CCInc_XSec_1Dthetamu_nu(samplekey));

	/*
	BNL Samples
	*/
	} else
	#endif
	#ifndef __NO_BNL__
	if (!name.compare("BNL_CCQE_XSec_1DEnu_nu")) {
	return (new BNL_CCQE_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BNL_CCQE_Evt_1DQ2_nu")) {
	return (new BNL_CCQE_Evt_1DQ2_nu(samplekey));

	/*
	BNL CC1ppip samples
	*/
	} else if (!name.compare("BNL_CC1ppip_XSec_1DEnu_nu") \|\|
	!name.compare("BNL_CC1ppip_XSec_1DEnu_nu_Uncorr") \|\|
	!name.compare("BNL_CC1ppip_XSec_1DEnu_nu_W14Cut") \|\|
	!name.compare("BNL_CC1ppip_XSec_1DEnu_nu_W14Cut_Uncorr")) {
	return (new BNL_CC1ppip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BNL_CC1ppip_Evt_1DQ2_nu") \|\|
	!name.compare("BNL_CC1ppip_Evt_1DQ2_nu_W14Cut")) {
	return (new BNL_CC1ppip_Evt_1DQ2_nu(samplekey));
	} else if (!name.compare("BNL_CC1ppip_Evt_1DcosthAdler_nu")) {
	return (new BNL_CC1ppip_Evt_1DcosthAdler_nu(samplekey));
	} else if (!name.compare("BNL_CC1ppip_Evt_1Dphi_nu")) {
	return (new BNL_CC1ppip_Evt_1Dphi_nu(samplekey));

	/*
	BNL CC1npip samples
	*/
	} else if (!name.compare("BNL_CC1npip_XSec_1DEnu_nu") \|\|
	!name.compare("BNL_CC1npip_XSec_1DEnu_nu_Uncorr")) {
	return (new BNL_CC1npip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BNL_CC1npip_Evt_1DQ2_nu")) {
	return (new BNL_CC1npip_Evt_1DQ2_nu(samplekey));
	/*
	BNL CC1pi0 samples
	*/
	} else if (!name.compare("BNL_CC1pi0_XSec_1DEnu_nu")) {
	return (new BNL_CC1pi0_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BNL_CC1pi0_Evt_1DQ2_nu")) {
	return (new BNL_CC1pi0_Evt_1DQ2_nu(samplekey));

	/*
	FNAL Samples
	*/
	} else
	#endif
	#ifndef __NO_FNAL__
	if (!name.compare("FNAL_CCQE_Evt_1DQ2_nu")) {
	return (new FNAL_CCQE_Evt_1DQ2_nu(samplekey));
	/*
	FNAL CC1ppip
	*/
	} else if (!name.compare("FNAL_CC1ppip_XSec_1DEnu_nu")) {
	return (new FNAL_CC1ppip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("FNAL_CC1ppip_XSec_1DQ2_nu")) {
	return (new FNAL_CC1ppip_XSec_1DQ2_nu(samplekey));
	} else if (!name.compare("FNAL_CC1ppip_Evt_1DQ2_nu")) {
	return (new FNAL_CC1ppip_Evt_1DQ2_nu(samplekey));
	/*
	FNAL CC1ppim
	*/
	} else if (!name.compare("FNAL_CC1ppim_XSec_1DEnu_antinu")) {
	return (new FNAL_CC1ppim_XSec_1DEnu_antinu(samplekey));

	/*
	BEBC Samples
	*/
	} else
	#endif
	#ifndef __NO_BEBC__
	if (!name.compare("BEBC_CCQE_XSec_1DQ2_nu")) {
	return (new BEBC_CCQE_XSec_1DQ2_nu(samplekey));
	/*
	BEBC CC1ppip samples
	*/
	} else if (!name.compare("BEBC_CC1ppip_XSec_1DEnu_nu")) {
	return (new BEBC_CC1ppip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BEBC_CC1ppip_XSec_1DQ2_nu")) {
	return (new BEBC_CC1ppip_XSec_1DQ2_nu(samplekey));
	/*
	BEBC CC1npip samples
	*/
	} else if (!name.compare("BEBC_CC1npip_XSec_1DEnu_nu")) {
	return (new BEBC_CC1npip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BEBC_CC1npip_XSec_1DQ2_nu")) {
	return (new BEBC_CC1npip_XSec_1DQ2_nu(samplekey));
	/*
	BEBC CC1pi0 samples
	*/
	} else if (!name.compare("BEBC_CC1pi0_XSec_1DEnu_nu")) {
	return (new BEBC_CC1pi0_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("BEBC_CC1pi0_XSec_1DQ2_nu")) {
	return (new BEBC_CC1pi0_XSec_1DQ2_nu(samplekey));
	/*
	BEBC CC1npim samples
	*/
	} else if (!name.compare("BEBC_CC1npim_XSec_1DEnu_antinu")) {
	return (new BEBC_CC1npim_XSec_1DEnu_antinu(samplekey));
	} else if (!name.compare("BEBC_CC1npim_XSec_1DQ2_antinu")) {
	return (new BEBC_CC1npim_XSec_1DQ2_antinu(samplekey));
	/*
	BEBC CC1ppim samples
	*/
	} else if (!name.compare("BEBC_CC1ppim_XSec_1DEnu_antinu")) {
	return (new BEBC_CC1ppim_XSec_1DEnu_antinu(samplekey));
	} else if (!name.compare("BEBC_CC1ppim_XSec_1DQ2_antinu")) {
	return (new BEBC_CC1ppim_XSec_1DQ2_antinu(samplekey));

	/*
	GGM CC1ppip samples
	*/
	} else
	#endif
	#ifndef __NO_GGM__
	if (!name.compare("GGM_CC1ppip_XSec_1DEnu_nu")) {
	return (new GGM_CC1ppip_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("GGM_CC1ppip_Evt_1DQ2_nu")) {
	return (new GGM_CC1ppip_Evt_1DQ2_nu(samplekey));

	/*
	MiniBooNE Samples
	*/
	/*
	CCQE
	*/
	} else
	#endif
	#ifndef __NO_MiniBooNE__
	if (!name.compare("MiniBooNE_CCQE_XSec_1DQ2_nu") \|\|
	!name.compare("MiniBooNE_CCQELike_XSec_1DQ2_nu")) {
	return (new MiniBooNE_CCQE_XSec_1DQ2_nu(samplekey));
	} else if (!name.compare("MiniBooNE_CCQE_XSec_1DQ2_antinu") \|\|
	!name.compare("MiniBooNE_CCQELike_XSec_1DQ2_antinu") \|\|
	!name.compare("MiniBooNE_CCQE_CTarg_XSec_1DQ2_antinu")) {
	return (new MiniBooNE_CCQE_XSec_1DQ2_antinu(samplekey));

	} else if (!name.compare("MiniBooNE_CCQE_XSec_2DTcos_nu") \|\|
	!name.compare("MiniBooNE_CCQELike_XSec_2DTcos_nu")) {
	return (new MiniBooNE_CCQE_XSec_2DTcos_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CCQE_XSec_2DTcos_antinu") \|\|
	!name.compare("MiniBooNE_CCQELike_XSec_2DTcos_antinu")) {
	return (new MiniBooNE_CCQE_XSec_2DTcos_antinu(samplekey));

	/*
	MiniBooNE CC1pi+
	*/
	// 1D
	} else if (!name.compare("MiniBooNE_CC1pip_XSec_1DEnu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_1DEnu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_1DQ2_nu")) {
	return (new MiniBooNE_CC1pip_XSec_1DQ2_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_1DTpi_nu")) {
	return (new MiniBooNE_CC1pip_XSec_1DTpi_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_1DTu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_1DTu_nu(samplekey));

	// 2D
	} else if (!name.compare("MiniBooNE_CC1pip_XSec_2DQ2Enu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_2DQ2Enu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_2DTpiCospi_nu")) {
	return (new MiniBooNE_CC1pip_XSec_2DTpiCospi_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_2DTpiEnu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_2DTpiEnu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_2DTuCosmu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_2DTuCosmu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pip_XSec_2DTuEnu_nu")) {
	return (new MiniBooNE_CC1pip_XSec_2DTuEnu_nu(samplekey));

	/*
	MiniBooNE CC1pi0
	*/
	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1DEnu_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1DEnu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1DQ2_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1DQ2_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1DTu_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1DTu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1Dcosmu_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1Dcosmu_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1Dcospi0_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1Dcospi0_nu(samplekey));

	} else if (!name.compare("MiniBooNE_CC1pi0_XSec_1Dppi0_nu")) {
	return (new MiniBooNE_CC1pi0_XSec_1Dppi0_nu(samplekey));

	} else if (!name.compare("MiniBooNE_NC1pi0_XSec_1Dcospi0_antinu") \|\|
	!name.compare("MiniBooNE_NC1pi0_XSec_1Dcospi0_rhc")) {
	return (new MiniBooNE_NC1pi0_XSec_1Dcospi0_antinu(samplekey));

	} else if (!name.compare("MiniBooNE_NC1pi0_XSec_1Dcospi0_nu") \|\|
	!name.compare("MiniBooNE_NC1pi0_XSec_1Dcospi0_fhc")) {
	return (new MiniBooNE_NC1pi0_XSec_1Dcospi0_nu(samplekey));

	} else if (!name.compare("MiniBooNE_NC1pi0_XSec_1Dppi0_antinu") \|\|
	!name.compare("MiniBooNE_NC1pi0_XSec_1Dppi0_rhc")) {
	return (new MiniBooNE_NC1pi0_XSec_1Dppi0_antinu(samplekey));

	} else if (!name.compare("MiniBooNE_NC1pi0_XSec_1Dppi0_nu") \|\|
	!name.compare("MiniBooNE_NC1pi0_XSec_1Dppi0_fhc")) {
	return (new MiniBooNE_NC1pi0_XSec_1Dppi0_nu(samplekey));

	/*
	MiniBooNE NCEL
	*/
	} else if (!name.compare("MiniBooNE_NCEL_XSec_Treco_nu")) {
	return (new MiniBooNE_NCEL_XSec_Treco_nu(samplekey));

	/*
	MINERvA Samples
	*/
	} else
	#endif
	#ifndef __NO_MINERvA__
	if (!name.compare("MINERvA_CCQE_XSec_1DQ2_nu") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_nu_20deg") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_nu_oldflux") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_nu_20deg_oldflux")) {
	return (new MINERvA_CCQE_XSec_1DQ2_nu(samplekey));

	} else if (!name.compare("MINERvA_CCQE_XSec_1DQ2_antinu") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_antinu_20deg") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_antinu_oldflux") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_antinu_20deg_oldflux")) {
	return (new MINERvA_CCQE_XSec_1DQ2_antinu(samplekey));

	} else if (!name.compare("MINERvA_CCQE_XSec_1DQ2_joint_oldflux") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_joint_20deg_oldflux") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_joint") \|\|
	!name.compare("MINERvA_CCQE_XSec_1DQ2_joint_20deg")) {
	return (new MINERvA_CCQE_XSec_1DQ2_joint(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DEe_nue")) {
	return (new MINERvA_CC0pi_XSec_1DEe_nue(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DQ2_nue")) {
	return (new MINERvA_CC0pi_XSec_1DQ2_nue(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DThetae_nue")) {
	return (new MINERvA_CC0pi_XSec_1DThetae_nue(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DQ2_nu_proton")) {
	return (new MINERvA_CC0pi_XSec_1DQ2_nu_proton(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtC_nu") \|\|
	!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtCH_nu") \|\|
	!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtFe_nu") \|\|
	!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtPb_nu")) {
	return (new MINERvA_CC0pi_XSec_1DQ2_Tgt_nu(samplekey));

	} else if (!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtRatioC_nu") \|\|
	!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtRatioFe_nu") \|\|
	!name.compare("MINERvA_CC0pi_XSec_1DQ2_TgtRatioPb_nu")) {
	return (new MINERvA_CC0pi_XSec_1DQ2_TgtRatio_nu(samplekey));

	/*
	CC1pi+
	*/
	// DONE
	} else if (!name.compare("MINERvA_CC1pip_XSec_1DTpi_nu") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1DTpi_nu_20deg") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1DTpi_nu_fluxcorr") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1DTpi_nu_20deg_fluxcorr")) {
	return (new MINERvA_CC1pip_XSec_1DTpi_nu(samplekey));

	// DONE
	} else if (!name.compare("MINERvA_CC1pip_XSec_1Dth_nu") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dth_nu_20deg") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dth_nu_fluxcorr") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dth_nu_20deg_fluxcorr")) {
	return (new MINERvA_CC1pip_XSec_1Dth_nu(samplekey));

	} else if (!name.compare("MINERvA_CC1pip_XSec_1DTpi_nu_2017") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dth_nu_2017") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dpmu_nu_2017") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1Dthmu_nu_2017") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1DQ2_nu_2017") \|\|
	!name.compare("MINERvA_CC1pip_XSec_1DEnu_nu_2017")) {
	return (new MINERvA_CC1pip_XSec_1D_2017Update(samplekey));
	/*
	CCNpi+
	*/
	} else if (!name.compare("MINERvA_CCNpip_XSec_1Dth_nu") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1Dth_nu_2015") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1Dth_nu_2016") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1Dth_nu_2015_20deg") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1Dth_nu_2015_fluxcorr") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1Dth_nu_2015_20deg_fluxcorr")) {
	return (new MINERvA_CCNpip_XSec_1Dth_nu(samplekey));

	} else if (!name.compare("MINERvA_CCNpip_XSec_1DTpi_nu") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1DTpi_nu_2015") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1DTpi_nu_2016") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1DTpi_nu_2015_20deg") \|\|
	!name.compare("MINERvA_CCNpip_XSec_1DTpi_nu_2015_fluxcorr") \|\|
	!name.compare(
	"MINERvA_CCNpip_XSec_1DTpi_nu_2015_20deg_fluxcorr")) {
	return (new MINERvA_CCNpip_XSec_1DTpi_nu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CCNpip_XSec_1Dthmu_nu")) {
	return (new MINERvA_CCNpip_XSec_1Dthmu_nu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CCNpip_XSec_1Dpmu_nu")) {
	return (new MINERvA_CCNpip_XSec_1Dpmu_nu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CCNpip_XSec_1DQ2_nu")) {
	return (new MINERvA_CCNpip_XSec_1DQ2_nu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CCNpip_XSec_1DEnu_nu")) {
	return (new MINERvA_CCNpip_XSec_1DEnu_nu(samplekey));

	/*
	CC1pi0
	*/
	// Done
	} else if (!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu_2015") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu_2016") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu_fluxcorr") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu_2015_fluxcorr") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dth_antinu_2016_fluxcorr")) {
	return (new MINERvA_CC1pi0_XSec_1Dth_antinu(samplekey));

	} else if (!name.compare("MINERvA_CC1pi0_XSec_1Dppi0_antinu") \|\|
	!name.compare("MINERvA_CC1pi0_XSec_1Dppi0_antinu_fluxcorr")) {
	return (new MINERvA_CC1pi0_XSec_1Dppi0_antinu(samplekey));

	} else if (!name.compare("MINERvA_CC1pi0_XSec_1DTpi0_antinu")) {
	return (new MINERvA_CC1pi0_XSec_1DTpi0_antinu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CC1pi0_XSec_1DQ2_antinu")) {
	return (new MINERvA_CC1pi0_XSec_1DQ2_antinu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CC1pi0_XSec_1Dthmu_antinu")) {
	return (new MINERvA_CC1pi0_XSec_1Dthmu_antinu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CC1pi0_XSec_1Dpmu_antinu")) {
	return (new MINERvA_CC1pi0_XSec_1Dpmu_antinu(samplekey));

	// Done
	} else if (!name.compare("MINERvA_CC1pi0_XSec_1DEnu_antinu")) {
	return (new MINERvA_CC1pi0_XSec_1DEnu_antinu(samplekey));

	/*
	CCINC
	*/
	} else if (!name.compare("MINERvA_CCinc_XSec_2DEavq3_nu")) {
	return (new MINERvA_CCinc_XSec_2DEavq3_nu(samplekey));

	} else if (!name.compare("MINERvA_CCinc_XSec_1Dx_ratio_C12_CH") \|\|
	!name.compare("MINERvA_CCinc_XSec_1Dx_ratio_Fe56_CH") \|\|
	!name.compare("MINERvA_CCinc_XSec_1Dx_ratio_Pb208_CH")) {
	return (new MINERvA_CCinc_XSec_1Dx_ratio(samplekey));

	} else if (!name.compare("MINERvA_CCinc_XSec_1DEnu_ratio_C12_CH") \|\|
	!name.compare("MINERvA_CCinc_XSec_1DEnu_ratio_Fe56_CH") \|\|
	!name.compare("MINERvA_CCinc_XSec_1DEnu_ratio_Pb208_CH")) {
	return (new MINERvA_CCinc_XSec_1DEnu_ratio(samplekey));
	/*
	CCDIS
	*/
	} else if (!name.compare("MINERvA_CCDIS_XSec_1Dx_ratio_C12_CH") \|\|
	!name.compare("MINERvA_CCDIS_XSec_1Dx_ratio_Fe56_CH") \|\|
	!name.compare("MINERvA_CCDIS_XSec_1Dx_ratio_Pb208_CH")) {
	return (new MINERvA_CCDIS_XSec_1Dx_ratio(samplekey));

	} else if (!name.compare("MINERvA_CCDIS_XSec_1DEnu_ratio_C12_CH") \|\|
	!name.compare("MINERvA_CCDIS_XSec_1DEnu_ratio_Fe56_CH") \|\|
	!name.compare("MINERvA_CCDIS_XSec_1DEnu_ratio_Pb208_CH")) {
	return (new MINERvA_CCDIS_XSec_1DEnu_ratio(samplekey));

	/*
	CC-COH
	*/
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1DEnu_nu")) {
	return (new MINERvA_CCCOHPI_XSec_1DEnu_nu(samplekey));
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1DEpi_nu")) {
	return (new MINERvA_CCCOHPI_XSec_1DEpi_nu(samplekey));
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1Dth_nu")) {
	return (new MINERvA_CCCOHPI_XSec_1Dth_nu(samplekey));
	+ } else if (!name.compare("MINERvA_CCCOHPI_XSec_1DQ2_nu")) {
	+ return (new MINERvA_CCCOHPI_XSec_1DQ2_nu(samplekey));
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1DEnu_antinu")) {
	return (new MINERvA_CCCOHPI_XSec_1DEnu_antinu(samplekey));
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1DEpi_antinu")) {
	return (new MINERvA_CCCOHPI_XSec_1DEpi_antinu(samplekey));
	} else if (!name.compare("MINERvA_CCCOHPI_XSec_1Dth_antinu")) {
	return (new MINERvA_CCCOHPI_XSec_1Dth_antinu(samplekey));
	+ } else if (!name.compare("MINERvA_CCCOHPI_XSec_1DQ2_antinu")) {
	+ return (new MINERvA_CCCOHPI_XSec_1DQ2_antinu(samplekey));
	+
	/*
	T2K Samples
	*/

	} else
	#endif
	#ifndef __NO_T2K__
	if (!name.compare("T2K_CC0pi_XSec_2DPcos_nu") \|\|
	!name.compare("T2K_CC0pi_XSec_2DPcos_nu_I") \|\|
	!name.compare("T2K_CC0pi_XSec_2DPcos_nu_II")) {
	return (new T2K_CC0pi_XSec_2DPcos_nu(samplekey));

	} else if (!name.compare("T2K_CC0pi_XSec_2DPcos_nu_nonuniform")) {
	return (new T2K_CC0pi_XSec_2DPcos_nu_nonuniform(samplekey));

	/*
	T2K CC1pi+ CH samples
	*/
	// Comment these out for now because we don't have the proper data
	/*

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dpmu_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dpmu_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dppi_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dppi_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1DQ2_nu")) {
	return (new T2K_CC1pip_CH_XSec_1DQ2_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dq3_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dq3_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dthmupi_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dthmupi_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dthpi_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dthpi_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1Dthq3pi_nu")) {
	return (new T2K_CC1pip_CH_XSec_1Dthq3pi_nu(file, rw, type, fkdt));

	} else if (!name.compare("T2K_CC1pip_CH_XSec_1DWrec_nu")) {
	return (new T2K_CC1pip_CH_XSec_1DWrec_nu(file, rw, type, fkdt));
	*/

	/*
	T2K CC1pi+ H2O samples
	*/
	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1DEnuDelta_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1DEnuDelta_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1DEnuMB_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1DEnuMB_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1Dcosmu_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1Dcosmu_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1Dcosmupi_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1Dcosmupi_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1Dcospi_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1Dcospi_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1Dpmu_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1Dpmu_nu(samplekey));

	} else if (!name.compare("T2K_CC1pip_H2O_XSec_1Dppi_nu")) {
	return (new T2K_CC1pip_H2O_XSec_1Dppi_nu(samplekey));

	/*
	T2K CC0pi + np CH samples
	*/
	} else if (!name.compare("T2K_CC0pinp_STV_XSec_1Ddpt_nu")) {
	return (new T2K_CC0pinp_STV_XSec_1Ddpt_nu(samplekey));

	// SciBooNE COH studies
	} else
	#endif
	#ifndef __NO_SciBooNE__
	if (!name.compare("SciBooNE_CCCOH_STOP_NTrks_nu")) {
	return (new SciBooNE_CCCOH_STOP_NTrks_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_1TRK_1DQ2_nu")) {
	return (new SciBooNE_CCCOH_1TRK_1DQ2_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_MuPr_1DQ2_nu")) {
	return (new SciBooNE_CCCOH_MuPr_1DQ2_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_MuPiVA_1DQ2_nu")) {
	return (new SciBooNE_CCCOH_MuPiVA_1DQ2_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_MuPiNoVA_1DQ2_nu")) {
	return (new SciBooNE_CCCOH_MuPiNoVA_1DQ2_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_MuPiNoVA_1Dthetapr_nu")) {
	return (new SciBooNE_CCCOH_MuPiNoVA_1Dthetapr_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_MuPiNoVA_1Dthetapi_nu")) {
	return (new SciBooNE_CCCOH_MuPiNoVA_1Dthetapi_nu(samplekey));
	} else if (!name.compare("SciBooNE_CCCOH_STOPFINAL_1DQ2_nu")) {
	return (new SciBooNE_CCCOH_STOPFINAL_1DQ2_nu(samplekey));

	/*
	K2K Samples
	*/
	/*
	NC1pi0
	*/
	} else
	#endif
	#ifndef __NO_K2K__
	if (!name.compare("K2K_NC1pi0_Evt_1Dppi0_nu")) {
	return (new K2K_NC1pi0_Evt_1Dppi0_nu(samplekey));

	/*
	Fake Studies
	*/

	} else
	#endif
	if (name.find("ExpMultDist_CCQE_XSec_1D") != std::string::npos &&
	name.find("_FakeStudy") != std::string::npos) {
	return (
	new ExpMultDist_CCQE_XSec_1DVar_FakeStudy(name, file, rw, type, fkdt));

	} else if (name.find("ExpMultDist_CCQE_XSec_2D") != std::string::npos &&
	name.find("_FakeStudy") != std::string::npos) {
	return (
	new ExpMultDist_CCQE_XSec_2DVar_FakeStudy(name, file, rw, type, fkdt));

	} else if (name.find("GenericFlux_") != std::string::npos) {
	return (new GenericFlux_Tester(name, file, rw, type, fkdt));

	} else if (name.find("GenericVectors_") != std::string::npos) {
	return (new GenericFlux_Vectors(name, file, rw, type, fkdt));

	} else if (!name.compare("T2K2017_FakeData")) {
	return (new T2K2017_FakeData(samplekey));

	} else if (!name.compare("MCStudy_CCQE")) {
	return (new MCStudy_CCQEHistograms(name, file, rw, type, fkdt));

	} else if (!name.compare("ElectronFlux_FlatTree")) {
	return (new ElectronFlux_FlatTree(name, file, rw, type, fkdt));

	} else if (name.find("ElectronData_") != std::string::npos) {
	return new ElectronScattering_DurhamData(samplekey);

	} else if (name.find("MuonValidation_") != std::string::npos) {
	return (new MCStudy_MuonValidation(name, file, rw, type, fkdt));

	} else if (!name.compare("NIWGOfficialPlots")) {
	return (new OfficialNIWGPlots(samplekey));

	}else if (!name.compare("Simple_Osc")) {
	return (new Simple_Osc(samplekey));

	} else {
	ERR(FTL) << "Error: No such sample: " << name << std::endl;
	exit(-1);
	return NULL;
	}

	// Return NULL if no sample loaded.
	return NULL;
	}
	}
	diff --git a/src/MINERvA/CMakeLists.txt b/src/MINERvA/CMakeLists.txt
	index 91889ab..d41a492 100644
	--- a/src/MINERvA/CMakeLists.txt
	+++ b/src/MINERvA/CMakeLists.txt
	@@ -1,153 +1,156 @@
	# 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/>.
	################################################################################
	set(IMPLFILES
	MINERvA_CCQE_XSec_1DQ2_antinu.cxx
	MINERvA_CCQE_XSec_1DQ2_joint.cxx
	MINERvA_CCQE_XSec_1DQ2_nu.cxx

	MINERvA_CC0pi_XSec_1DEe_nue.cxx
	MINERvA_CC0pi_XSec_1DQ2_nue.cxx
	MINERvA_CC0pi_XSec_1DQ2_nu_proton.cxx
	MINERvA_CC0pi_XSec_1DThetae_nue.cxx

	MINERvA_CC1pi0_XSec_1DEnu_antinu.cxx
	MINERvA_CC1pi0_XSec_1DQ2_antinu.cxx
	MINERvA_CC1pi0_XSec_1Dpmu_antinu.cxx
	MINERvA_CC1pi0_XSec_1Dppi0_antinu.cxx
	MINERvA_CC1pi0_XSec_1DTpi0_antinu.cxx
	MINERvA_CC1pi0_XSec_1Dth_antinu.cxx
	MINERvA_CC1pi0_XSec_1Dthmu_antinu.cxx

	MINERvA_CC1pip_XSec_1DTpi_20deg_nu.cxx
	MINERvA_CC1pip_XSec_1DTpi_nu.cxx
	MINERvA_CC1pip_XSec_1Dth_20deg_nu.cxx
	MINERvA_CC1pip_XSec_1Dth_nu.cxx

	MINERvA_CCNpip_XSec_1DEnu_nu.cxx
	MINERvA_CCNpip_XSec_1DQ2_nu.cxx
	MINERvA_CCNpip_XSec_1DTpi_nu.cxx
	MINERvA_CCNpip_XSec_1Dpmu_nu.cxx
	MINERvA_CCNpip_XSec_1Dth_nu.cxx
	MINERvA_CCNpip_XSec_1Dthmu_nu.cxx

	MINERvA_CCinc_XSec_2DEavq3_nu.cxx
	MINERvA_CCinc_XSec_1Dx_ratio.cxx
	MINERvA_CCinc_XSec_1DEnu_ratio.cxx
	MINERvA_CCinc_XSec_1Dx_nu.cxx
	MINERvA_CCinc_XSec_1DEnu_nu.cxx

	MINERvA_CCDIS_XSec_1Dx_ratio.cxx
	MINERvA_CCDIS_XSec_1DEnu_ratio.cxx
	MINERvA_CCDIS_XSec_1Dx_nu.cxx
	MINERvA_CCDIS_XSec_1DEnu_nu.cxx

	MINERvA_CC0pi_XSec_1DQ2_Tgt_nu.cxx
	MINERvA_CC0pi_XSec_1DQ2_TgtRatio_nu.cxx
	MINERvA_CC1pip_XSec_1D_2017Update.cxx

	MINERvA_CCCOHPI_XSec_1DEnu_nu.cxx
	MINERvA_CCCOHPI_XSec_1DEpi_nu.cxx
	MINERvA_CCCOHPI_XSec_1Dth_nu.cxx
	+MINERvA_CCCOHPI_XSec_1DQ2_nu.cxx
	+
	MINERvA_CCCOHPI_XSec_1DEnu_antinu.cxx
	MINERvA_CCCOHPI_XSec_1DEpi_antinu.cxx
	MINERvA_CCCOHPI_XSec_1Dth_antinu.cxx
	+MINERvA_CCCOHPI_XSec_1DQ2_antinu.cxx

	MINERvAUtils.cxx
	MINERvA_SignalDef.cxx
	)

	set(HEADERFILES
	MINERvA_CCQE_XSec_1DQ2_antinu.h
	MINERvA_CCQE_XSec_1DQ2_joint.h
	MINERvA_CCQE_XSec_1DQ2_nu.h

	MINERvA_CC0pi_XSec_1DEe_nue.h
	MINERvA_CC0pi_XSec_1DQ2_nue.h
	MINERvA_CC0pi_XSec_1DQ2_nu_proton.h
	MINERvA_CC0pi_XSec_1DThetae_nue.h

	MINERvA_CC1pi0_XSec_1DEnu_antinu.h
	MINERvA_CC1pi0_XSec_1DQ2_antinu.h
	MINERvA_CC1pi0_XSec_1Dpmu_antinu.h
	MINERvA_CC1pi0_XSec_1Dppi0_antinu.h
	MINERvA_CC1pi0_XSec_1DTpi0_antinu.h
	MINERvA_CC1pi0_XSec_1Dth_antinu.h
	MINERvA_CC1pi0_XSec_1Dthmu_antinu.h

	MINERvA_CC1pip_XSec_1DTpi_20deg_nu.h
	MINERvA_CC1pip_XSec_1DTpi_nu.h
	MINERvA_CC1pip_XSec_1Dth_20deg_nu.h
	MINERvA_CC1pip_XSec_1Dth_nu.h

	MINERvA_CCNpip_XSec_1DEnu_nu.h
	MINERvA_CCNpip_XSec_1DQ2_nu.h
	MINERvA_CCNpip_XSec_1DTpi_nu.h
	MINERvA_CCNpip_XSec_1Dpmu_nu.h
	MINERvA_CCNpip_XSec_1Dth_nu.h
	MINERvA_CCNpip_XSec_1Dthmu_nu.h

	MINERvA_CCinc_XSec_2DEavq3_nu.h
	MINERvA_CCinc_XSec_1Dx_ratio.h
	MINERvA_CCinc_XSec_1DEnu_ratio.h
	MINERvA_CCinc_XSec_1Dx_nu.h
	MINERvA_CCinc_XSec_1DEnu_nu.h

	MINERvA_CCDIS_XSec_1Dx_ratio.h
	MINERvA_CCDIS_XSec_1DEnu_ratio.h
	MINERvA_CCDIS_XSec_1Dx_nu.h
	MINERvA_CCDIS_XSec_1DEnu_nu.h

	MINERvA_CC0pi_XSec_1DQ2_Tgt_nu.h
	MINERvA_CC0pi_XSec_1DQ2_TgtRatio_nu.h
	MINERvA_CC1pip_XSec_1D_2017Update.h

	MINERvA_CCCOHPI_XSec_1DEnu_nu.h
	MINERvA_CCCOHPI_XSec_1DEpi_nu.h
	MINERvA_CCCOHPI_XSec_1Dth_nu.h
	MINERvA_CCCOHPI_XSec_1DEnu_antinu.h
	MINERvA_CCCOHPI_XSec_1DEpi_antinu.h
	MINERvA_CCCOHPI_XSec_1Dth_antinu.h

	MINERvAUtils.h
	MINERvA_SignalDef.h
	MINERvAVariableBoxes.h
	)

	set(LIBNAME expMINERvA)

	if(CMAKE_BUILD_TYPE MATCHES DEBUG)
	add_library(${LIBNAME} STATIC ${IMPLFILES})
	else(CMAKE_BUILD_TYPE MATCHES RELEASE)
	add_library(${LIBNAME} SHARED ${IMPLFILES})
	endif()

	include_directories(${MINIMUM_INCLUDE_DIRECTORIES})

	set_target_properties(${LIBNAME} PROPERTIES VERSION
	"${NUISANCE_VERSION_MAJOR}.${NUISANCE_VERSION_MINOR}.${NUISANCE_VERSION_REVISION}")
	#set_target_properties(${LIBNAME} PROPERTIES LINK_FLAGS ${ROOT_LD_FLAGS})
	if(DEFINED PROJECTWIDE_EXTRA_DEPENDENCIES)
	add_dependencies(${LIBNAME} ${PROJECTWIDE_EXTRA_DEPENDENCIES})
	endif()

	install(TARGETS ${LIBNAME} DESTINATION lib)
	#Can uncomment this to install the headers... but is it really neccessary?
	#install(FILES ${HEADERFILES} DESTINATION include)

	set(MODULETargets ${MODULETargets} ${LIBNAME} PARENT_SCOPE)
	diff --git a/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.cxx b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.cxx
	new file mode 100644
	index 0000000..ae6ac5f
	--- /dev/null
	+++ b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.cxx
	@@ -0,0 +1,79 @@
	+// 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_CCCOHPI_XSec_1DQ2_antinu.h"
	+
	+
	+//********************************************************************
	+MINERvA_CCCOHPI_XSec_1DQ2_antinu::MINERvA_CCCOHPI_XSec_1DQ2_antinu(nuiskey samplekey) {
	+//********************************************************************
	+
	+ // Sample overview ---------------------------------------------------
	+ std::string descrip = "MINERvA_CCCOHPI_XSec_1DQ2_antinu sample. \n" \
	+ "Target: CH \n" \
	+ "Flux: MINERvA Reverse Horn Current numu \n" \
	+ "Signal: Any event with 1 mu+, 1pi-, and no other FS particles \n";
	+
	+ // Setup common settings
	+ fSettings = LoadSampleSettings(samplekey);
	+ fSettings.SetDescription(descrip);
	+ fSettings.SetXTitle("E_{#pi} (GeV)");
	+ fSettings.SetYTitle("d#sigma/dE_{#pi} (cm^{2}/GeV/C^{12})");
	+ fSettings.SetAllowedTypes("FIX,FREE,SHAPE/DIAG,FULL/NORM/MASK", "FIX/FULL");
	+ fSettings.SetEnuRange(1.5, 20.0);
	+ fSettings.DefineAllowedTargets("C,H");
	+ fSettings.DefineAllowedSpecies("numub");
	+ fSettings.SetTitle("MINERvA_CCCOHPI_XSec_1DQ2_antinu");
	+
	+ fSettings.SetDataInput( GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCcoh/Q2_nubar_data.csv");
	+ fSettings.SetCovarInput(GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCcoh/Q2_nubar_cov.csv");
	+
	+ FinaliseSampleSettings();
	+
	+ // Scaling Setup ---------------------------------------------------
	+ // ScaleFactor automatically setup for DiffXSec/cm2/Nucleon
	+ fScaleFactor = GetEventHistogram()->Integral("width") * double(1E-38) / double(fNEvents)/TotalIntegratedFlux("width")*13;
	+
	+ // Plot Setup -------------------------------------------------------
	+ SetDataFromTextFile( fSettings.GetDataInput() );
	+ SetCovarFromMultipleTextFiles(fSettings.GetCovarInput());
	+
	+ // Final setup ---------------------------------------------------
	+ FinaliseMeasurement();
	+
	+};
	+
	+
	+void MINERvA_CCCOHPI_XSec_1DQ2_antinu::FillEventVariables(FitEvent *event) {
	+
	+ if (event->NumFSParticle(-211) == 0) return;
	+
	+ TLorentzVector Ppi = event->GetHMFSParticle(-211)->fP;
	+ fXVar = 0.0;
	+
	+ return;
	+};
	+
	+//********************************************************************
	+bool MINERvA_CCCOHPI_XSec_1DQ2_antinu::isSignal(FitEvent *event) {
	+//********************************************************************
	+ // Signal: numu + 12C --> mu- + pi+ + 12C'
	+ return SignalDef::isCCCOH(event, -14, -211, EnuMin, EnuMax);
	+}
	+
	diff --git a/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.h b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.h
	new file mode 100644
	index 0000000..dde6ce6
	--- /dev/null
	+++ b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_antinu.h
	@@ -0,0 +1,36 @@
	+// 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/>.
	+*******************************************************************************/
	+
	+#ifndef MINERVA_CCCOHPI_XSEC_1DQ2_ANTINU_H_SEEN
	+#define MINERVA_CCCOHPI_XSEC_1DQ2_ANTINU_H_SEEN
	+
	+#include "Measurement1D.h"
	+
	+class MINERvA_CCCOHPI_XSec_1DQ2_antinu : public Measurement1D {
	+public:
	+ MINERvA_CCCOHPI_XSec_1DQ2_antinu(nuiskey samplekey);
	+ virtual ~MINERvA_CCCOHPI_XSec_1DQ2_antinu() {};
	+
	+ void FillEventVariables(FitEvent *event);
	+ bool isSignal(FitEvent *event);
	+
	+ private:
	+};
	+
	+#endif
	diff --git a/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.cxx b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.cxx
	new file mode 100644
	index 0000000..0cf7db8
	--- /dev/null
	+++ b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.cxx
	@@ -0,0 +1,79 @@
	+// 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_CCCOHPI_XSec_1DQ2_nu.h"
	+
	+//********************************************************************
	+MINERvA_CCCOHPI_XSec_1DQ2_nu::MINERvA_CCCOHPI_XSec_1DQ2_nu(nuiskey samplekey) {
	+//********************************************************************
	+
	+ // Sample overview ---------------------------------------------------
	+ std::string descrip = "MINERvA_CCCOHPI_XSec_1DQ2_nu sample. \n" \
	+ "Target: CH \n" \
	+ "Flux: MINERvA Forward Horn Current numu \n" \
	+ "Signal: Any event with 1 mu-, 1pi+, and no other FS particles \n";
	+
	+ // Setup common settings
	+ fSettings = LoadSampleSettings(samplekey);
	+ fSettings.SetDescription(descrip);
	+ fSettings.SetXTitle("E_{#pi} (GeV)");
	+ fSettings.SetYTitle("d#sigma/dE_{#pi} (cm^{2}/GeV/C^{12})");
	+ fSettings.SetAllowedTypes("FIX,FREE,SHAPE/DIAG,FULL/NORM/MASK", "FIX/FULL");
	+ fSettings.SetEnuRange(1.5, 20.0);
	+ fSettings.DefineAllowedTargets("C,H");
	+ fSettings.DefineAllowedSpecies("numu");
	+ fSettings.SetTitle("MINERvA_CCCOHPI_XSec_1DQ2_nu");
	+
	+ fSettings.SetDataInput( GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCcoh/Q2_nu_data.csv");
	+ fSettings.SetCovarInput(GeneralUtils::GetTopLevelDir() + "/data/MINERvA/CCcoh/Q2_nu_cov.csv");
	+
	+ FinaliseSampleSettings();
	+
	+ // Scaling Setup ---------------------------------------------------
	+ // ScaleFactor automatically setup for DiffXSec/cm2/Nucleon
	+ fScaleFactor = GetEventHistogram()->Integral("width") * double(1E-38) / double(fNEvents)/TotalIntegratedFlux("width")*13;
	+
	+ // Plot Setup -------------------------------------------------------
	+ SetDataFromTextFile( fSettings.GetDataInput() );
	+ SetCovarFromMultipleTextFiles(fSettings.GetCovarInput());
	+
	+ // Final setup ---------------------------------------------------
	+ FinaliseMeasurement();
	+};
	+
	+
	+void MINERvA_CCCOHPI_XSec_1DQ2_nu::FillEventVariables(FitEvent *event) {
	+
	+ if (event->NumFSParticle(211) == 0) return;
	+
	+ TLorentzVector Ppi = event->GetHMFSParticle(211)->fP;
	+
	+ double Q2 = 0.0;
	+ fXVar = Q2;
	+
	+ return;
	+};
	+
	+//********************************************************************
	+bool MINERvA_CCCOHPI_XSec_1DQ2_nu::isSignal(FitEvent *event) {
	+//********************************************************************
	+ // Signal: numu + 12C --> mu- + pi+ + 12C'
	+ return SignalDef::isCCCOH(event, 14, 211, EnuMin, EnuMax);
	+}
	+
	diff --git a/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.h b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.h
	new file mode 100644
	index 0000000..d22b181
	--- /dev/null
	+++ b/src/MINERvA/MINERvA_CCCOHPI_XSec_1DQ2_nu.h
	@@ -0,0 +1,36 @@
	+// 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/>.
	+*******************************************************************************/
	+
	+#ifndef MINERVA_CCCOHPI_XSEC_1DQ2_NU_H_SEEN
	+#define MINERVA_CCCOHPI_XSEC_1DQ2_NU_H_SEEN
	+
	+#include "Measurement1D.h"
	+
	+class MINERvA_CCCOHPI_XSec_1DQ2_nu : public Measurement1D {
	+public:
	+ MINERvA_CCCOHPI_XSec_1DQ2_nu(nuiskey samplekey);
	+ virtual ~MINERvA_CCCOHPI_XSec_1DQ2_nu() {};
	+
	+ void FillEventVariables(FitEvent *event);
	+ bool isSignal(FitEvent *event);
	+
	+ private:
	+};
	+
	+#endif
	diff --git a/src/Reweight/MINERvAWeightCalcs.cxx b/src/Reweight/MINERvAWeightCalcs.cxx
	new file mode 100644
	index 0000000..5fefde9
	--- /dev/null
	+++ b/src/Reweight/MINERvAWeightCalcs.cxx
	@@ -0,0 +1,345 @@
	+#ifdef __MINERVA_RW_ENABLED__
	+#ifdef __GENIE_ENABLED__
	+
	+#include "MINERvAWeightCalcs.h"
	+#include "BaseFitEvt.h"
	+
	+namespace nuisance {
	+namespace reweight {
	+
	+MINERvAReWeight_MEC::MINERvAReWeight_MEC() {
	+ fTwk_NormCCMEC = 0.0;
	+ fDef_NormCCMEC = 1.0;
	+ fCur_NormCCMEC = fDef_NormCCMEC;
	+}
	+
	+virtual ~MINERvAReWeight_MEC::MINERvAReWeight_MEC(){};
	+
	+double MINERvAReWeight_MEC::CalcWeight(BaseFitEvt* evt) {
	+ // Check GENIE
	+ if (evt->fType != kGENIE) return 1.0;
	+
	+ // Extract the GENIE Record
	+ GHepRecord* ghep = static_cast<GHepRecord*>(evt->genie_event->event);
	+ const Interaction* interaction = ghep->Summary();
	+ const InitialState& init_state = interaction->InitState();
	+ const ProcessInfo& proc_info = interaction->ProcInfo();
	+ const Target& tgt = init_state.Tgt();
	+
	+ // If the event is not MEC this Calc doesn't handle it
	+ if (!proc_info.IsMEC()) return 1.0;
	+
	+ // WEIGHT CALCULATIONS -------------
	+ double w = 1.0;
	+
	+ // CCMEC Dial
	+ if (!proc_info.IsCC()) w *= fCur_NormCCMEC;
	+
	+ // Return Combined Weight
	+ return w;
	+}
	+
	+void MINERvAReWeight_MEC::SetDialValue(std::string name, double val) {
	+ SetDialValue(Reweight::ConvDial(name, kCUSTOM), val);
	+}
	+
	+void MINERvAReWeight_MEC::SetDialValue(int rwenum, double val) {
	+ // Check Handled
	+ int curenum = rwenum % 1000;
	+ if (!IsHandled(curenum)) return;
	+
	+ // Set Values
	+ if (curenum == kMINERvARW_NormCCMEC) fTwk_NormCCMEC = val;
	+
	+ // Define Tweaked
	+ fTweaked = ((fTwk_NormCCMEC != 0.0));
	+}
	+
	+bool MINERvAReWeight_MEC::IsHandled(int rwenum) {
	+ int curenum = rwenum % 1000;
	+
	+ switch (curenum) {
	+ case kMINERvARW_NormCCMEC:
	+ return true;
	+ default:
	+ return false;
	+ }
	+}
	+} // namespace reweight
	+} // namespace nuisance
	+#endif
	+#endif
	+
	+
	+
	+
	+
	+
	+// RikRPA::RikRPA() {
	+
	+// // - Syst : kMINERvA_RikRPA_ApplyRPA
	+// // - Type : Binary
	+// // - Limits : 0.0 (false) -> 1.0 (true)
	+// // - Default : 0.0
	+// fApplyDial_RPACorrection = false;
	+
	+// // - Syst : kMINERvA_RikRPA_LowQ2
	+// // - Type : Absolute
	+// // - Limits : 1.0 -> 1.0
	+// // - Default : 0.0
	+// // - Frac Error : 100%
	+// fDefDial_RPALowQ2 = 0.0;
	+// fCurDial_RPALowQ2 = fDefDial_RPALowQ2;
	+// fErrDial_RPALowQ2 = 0.0;
	+
	+// // - Syst : kMINERvA_RikRPA_HighQ2
	+// // - Type : Absolute
	+// // - Limits : 1.0 -> 1.0
	+// // - Default : 0.0
	+// // - Frac Error : 100%
	+// fDefDial_RPAHighQ2 = 0.0;
	+// fCurDial_RPAHighQ2 = fDefDial_RPAHighQ2;
	+// fErrDial_RPAHighQ2 = 1.0;
	+
	+// fEventWeights = new double[5];
	+
	+// for (size_t i = 0; i < kMaxCalculators; i++) {
	+// fRPACalculators[i] = NULL;
	+// }
	+
	+// fTweaked = false;
	+// }
	+
	+// RikRPA::~RikRPA() {
	+
	+// delete fEventWeights;
	+
	+// for (size_t i = 0; i < kMaxCalculators; i++) {
	+// if (fRPACalculators[i]) delete fRPACalculators[i];
	+// fRPACalculators[i] = NULL;
	+// }
	+
	+// }
	+
	+// double RikRPA::CalcWeight(BaseFitEvt* evt) {
	+// // LOG(FIT) << "Calculating RikRPA" << std::endl;
	+// // Return 1.0 if not tweaked
	+// if (!fTweaked) return 1.0;
	+
	+// double w = 1.0;
	+
	+// // Extract the GENIE Record
	+// GHepRecord* ghep = static_cast<GHepRecord*>(evt->genie_event->event);
	+// const Interaction* interaction = ghep->Summary();
	+// const InitialState& init_state = interaction->InitState();
	+// const ProcessInfo& proc_info = interaction->ProcInfo();
	+// // const Kinematics & kine = interaction->Kine();
	+// // const XclsTag & xcls = interaction->ExclTag();
	+// const Target& tgt = init_state.Tgt();
	+
	+// // If not QE return 1.0
	+// // LOG(FIT) << "RikRPA : Event QE = " << proc_info.IsQuasiElastic() <<
	+// std::endl;
	+// if (!tgt.IsNucleus()) return 1.0;
	+// if (!proc_info.IsQuasiElastic()) return 1.0;
	+
	+// // Extract Beam and Target PDG
	+// GHepParticle* neutrino = ghep->Probe();
	+// int bpdg = neutrino->Pdg();
	+
	+// GHepParticle* target = ghep->Particle(1);
	+// assert(target);
	+// int tpdg = target->Pdg();
	+
	+// // Find the enum we need
	+// int calcenum = GetRPACalcEnum(bpdg, tpdg);
	+// if (calcenum == -1) return 1.0;
	+
	+// // Check we have the RPA Calc setup for this enum
	+// // if not, set it up at that point
	+// if (!fRPACalculators[calcenum]) SetupRPACalculator(calcenum);
	+// weightRPA* rpacalc = fRPACalculators[calcenum];
	+// if (!rpacalc) {
	+// THROW("Failed to grab the RPA Calculator : " << calcenum);
	+// }
	+
	+// // Extract Q0-Q3
	+// GHepParticle* fsl = ghep->FinalStatePrimaryLepton();
	+// const TLorentzVector& k1 = *(neutrino->P4());
	+// const TLorentzVector& k2 = *(fsl->P4());
	+// double q0 = fabs((k1 - k2).E());
	+// double q3 = fabs((k1 - k2).Vect().Mag());
	+
	+// // Now use q0-q3 and RPA Calculator to fill fWeights
	+// // LOG(FIT) << "Getting Weights = " << q0 << " " << q3 << std::endl;
	+// rpacalc->getWeight(q0, q3, fEventWeights);
	+
	+// // Apply Interpolation (for the time being simple linear)
	+
	+// // Syst Application : kMINERvA_RikRPA_ApplyRPA
	+// if (fApplyDial_RPACorrection) {
	+// w *= fEventWeights[0]; // CV
	+// }
	+
	+// LOG(FIT) << " fCurDial_RPALowQ2 = " << fCurDial_RPALowQ2
	+// << " fCurDial_RPAHighQ2 = " << fCurDial_RPAHighQ2 << " Weights "
	+// << fEventWeights[0] << " " << fEventWeights[1] << " "
	+// << fEventWeights[2] << " " << fEventWeights[3] << " "
	+// << fEventWeights[4] << std::endl;
	+
	+// // Syst Application : kMINERvA_RikRPA_LowQ2
	+// if (fabs(fCurDial_RPALowQ2) > 0.0) {
	+// double interpw = fEventWeights[0];
	+
	+// if (fCurDial_RPALowQ2 > 0.0) {
	+// interpw = fEventWeights[0] -
	+// (fEventWeights[0] - fEventWeights[1]) *
	+// fCurDial_RPALowQ2; // WLow+ } else if
	+// (fCurDial_RPALowQ2 < 0.0) {
	+// interpw = fEventWeights[0] - (fEventWeights[2] - fEventWeights[0]) *
	+// fCurDial_RPALowQ2; // WLow-
	+// }
	+// w *= interpw / fEventWeights[0]; // Div by CV again
	+// }
	+
	+// // Syst Application : kMINERvA_RikRPA_HighQ2
	+// if (fabs(fCurDial_RPAHighQ2) > 0.0) {
	+// double interpw = fEventWeights[0];
	+// if (fCurDial_RPAHighQ2 > 0.0) {
	+// interpw = fEventWeights[0] -
	+// (fEventWeights[0] - fEventWeights[3]) *
	+// fCurDial_RPAHighQ2; // WHigh+ } else
	+// if
	+// (fCurDial_RPAHighQ2 < 0.0) {
	+// interpw = fEventWeights[0] - (fEventWeights[4] - fEventWeights[0])
	+// *
	+// fCurDial_RPAHighQ2; // WHigh-
	+// }
	+// w *= interpw / fEventWeights[0]; // Div by CV again
	+// }
	+
	+// // LOG(FIT) << "RPA Weight = " << w << std::endl;
	+
	+// return w;
	+// }
	+
	+// void RikRPA::SetDialValue(std::string name, double val) {
	+// SetDialValue(Reweight::ConvDial(name, kCUSTOM), val);
	+// }
	+
	+// void RikRPA::SetDialValue(int rwenum, double val) {
	+// int curenum = rwenum % 1000;
	+
	+// // Check Handled
	+// if (!IsHandled(curenum)) return;
	+// if (curenum == kMINERvA_RikRPA_ApplyRPA)
	+// fApplyDial_RPACorrection = (val > 0.5);
	+// if (curenum == kMINERvA_RikRPA_LowQ2) fCurDial_RPALowQ2 = val;
	+// if (curenum == kMINERvA_RikRPA_HighQ2) fCurDial_RPAHighQ2 = val;
	+
	+// // Assign flag to say stuff has changed
	+// fTweaked = (fApplyDial_RPACorrection \|\|
	+// fabs(fCurDial_RPAHighQ2 - fDefDial_RPAHighQ2) > 0.0 \|\|
	+// fabs(fCurDial_RPALowQ2 - fDefDial_RPALowQ2) > 0.0);
	+// }
	+
	+// bool RikRPA::IsHandled(int rwenum) {
	+// int curenum = rwenum % 1000;
	+// switch (curenum) {
	+// case kMINERvA_RikRPA_ApplyRPA:
	+// return true;
	+// case kMINERvA_RikRPA_LowQ2:
	+// return true;
	+// case kMINERvA_RikRPA_HighQ2:
	+// return true;
	+// default:
	+// return false;
	+// }
	+// }
	+
	+// void RikRPA::SetupRPACalculator(int calcenum) {
	+// std::string rwdir = FitPar::GetDataBase() +
	+// "/reweight/MINERvA/RikRPA/"; std::string fidir = ""; switch (calcenum)
	+// {
	+// case kNuMuC12:
	+// fidir = "outNievesRPAratio-nu12C-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuO16:
	+// fidir = "outNievesRPAratio-nu16O-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuAr40:
	+// fidir = "outNievesRPAratio-nu40Ar-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuCa40:
	+// fidir = "outNievesRPAratio-nu40Ca-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuFe56:
	+// fidir = "outNievesRPAratio-nu56Fe-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuBarC12:
	+// fidir = "outNievesRPAratio-anu12C-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuBarO16:
	+// fidir = "outNievesRPAratio-anu16O-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuBarAr40:
	+// fidir = "outNievesRPAratio-anu40Ar-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuBarCa40:
	+// fidir = "outNievesRPAratio-anu40Ca-20GeV-20170202.root";
	+// break;
	+
	+// case kNuMuBarFe56:
	+// fidir = "outNievesRPAratio-anu56Fe-20GeV-20170202.root";
	+// break;
	+// }
	+
	+// LOG(FIT) << "Loading RPA CALC : " << fidir << std::endl;
	+// TDirectory* olddir = gDirectory;
	+
	+// fRPACalculators[calcenum] = new weightRPA(rwdir + "/" + fidir);
	+// olddir->cd();
	+// return;
	+// }
	+
	+// int RikRPA::GetRPACalcEnum(int bpdg, int tpdg) {
	+// if (bpdg == 14 && tpdg == 1000060120)
	+// return kNuMuC12;
	+// else if (bpdg == 14 && tpdg == 1000080160)
	+// return kNuMuO16;
	+// else if (bpdg == 14 && tpdg == 1000180400)
	+// return kNuMuAr40;
	+// else if (bpdg == 14 && tpdg == 1000200400)
	+// return kNuMuCa40;
	+// else if (bpdg == 14 && tpdg == 1000280560)
	+// return kNuMuFe56;
	+// else if (bpdg == -14 && tpdg == 1000060120)
	+// return kNuMuBarC12;
	+// else if (bpdg == -14 && tpdg == 1000080160)
	+// return kNuMuBarO16;
	+// else if (bpdg == -14 && tpdg == 1000180400)
	+// return kNuMuBarAr40;
	+// else if (bpdg == -14 && tpdg == 1000200400)
	+// return kNuMuBarCa40;
	+// else if (bpdg == -14 && tpdg == 1000280560)
	+// return kNuMuBarFe56;
	+// else {
	+// ERROR(WRN, "Unknown beam and target combination for RPA Calcs! "
	+// << bpdg << " " << tpdg);
	+// }
	+
	+// return -1;
	+// }
	+// } // namespace reweight
	+// } // namespace nuisance
	+
	+// #endif // __GENIE_ENABLED__
	+// #endif //__MINERVA_RW_ENABLED__
	diff --git a/src/Reweight/MINERvAWeightCalcs.h b/src/Reweight/MINERvAWeightCalcs.h
	new file mode 100644
	index 0000000..ec059dc
	--- /dev/null
	+++ b/src/Reweight/MINERvAWeightCalcs.h
	@@ -0,0 +1,96 @@
	+#ifndef MINERVA_WEIGHT_CALCS
	+#define MINERVA_WEIGHT_CALCS
	+
	+#include <string>
	+
	+// #ifdef __MINERVA_RW_ENABLED__
	+// #ifdef __GENIE_ENABLED__
	+#include "Conventions/Units.h"
	+#include "EVGCore/EventRecord.h"
	+#include "FitEvent.h"
	+#include "FitParameters.h"
	+#include "GHEP/GHepParticle.h"
	+#include "GHEP/GHepRecord.h"
	+#include "GHEP/GHepUtils.h"
	+#include "GeneralUtils.h"
	+#include "NUISANCESyst.h"
	+#include "NUISANCEWeightCalcs.h"
	+#include "Ntuple/NtpMCEventRecord.h"
	+#include "PDG/PDGUtils.h"
	+#include "WeightUtils.h"
	+// #include "weightRPA.h"
	+using namespace genie;
	+class BaseFitEvt;
	+
	+
	+namespace nuisance {
	+namespace reweight {
	+
	+// MEC Dials
	+class MINERvAReWeight_MEC : public NUISANCEWeightCalc {
	+ public:
	+ MINERvAReWeight_MEC();
	+ virtual ~MINERvAReWeight_MEC();
	+
	+ double CalcWeight(BaseFitEvt* evt);
	+ void SetDialValue(std::string name, double val);
	+ void SetDialValue(int rwenum, double val);
	+ bool IsHandled(int rwenum);
	+
	+ double fTwk_NormCCMEC;
	+ double fCur_NormCCMEC;
	+ double fDef_NormCCMEC;
	+}
	+
	+// /// RPA Weight Calculator that applies RPA systematics
	+// /// to GENIE events. GENIE EVENTS ONLY!
	+// class RikRPA : public NUISANCEWeightCalc {
	+// public:
	+// RikRPA();
	+// ~RikRPA();
	+
	+// double CalcWeight(BaseFitEvt* evt);
	+// void SetDialValue(std::string name, double val);
	+// void SetDialValue(int rwenum, double val);
	+// bool IsHandled(int rwenum);
	+
	+// void SetupRPACalculator(int calcenum);
	+// int GetRPACalcEnum(int bpdg, int tpdg);
	+
	+// bool fApplyDial_RPACorrection;
	+
	+// double fTwkDial_RPALowQ2;
	+// double fDefDial_RPALowQ2;
	+// double fCurDial_RPALowQ2;
	+// double fErrDial_RPALowQ2;
	+
	+// double fTwkDial_RPAHighQ2;
	+// double fDefDial_RPAHighQ2;
	+// double fCurDial_RPAHighQ2;
	+// double fErrDial_RPAHighQ2;
	+
	+// double* fEventWeights;
	+// bool fTweaked;
	+
	+// const static int kMaxCalculators = 10;
	+// enum rpacalcenums {
	+// kNuMuC12,
	+// kNuMuO16,
	+// kNuMuAr40,
	+// kNuMuCa40,
	+// kNuMuFe56,
	+// kNuMuBarC12,
	+// kNuMuBarO16,
	+// kNuMuBarAr40,
	+// kNuMuBarCa40,
	+// kNuMuBarFe56
	+// };
	+// weightRPA* fRPACalculators[kMaxCalculators];
	+// };
	+}; // namespace reweight
	+}; // namespace nuisance
	+
	+// #endif // __GENIE_ENABLED__
	+// #endif //__MINERVA_RW_ENABLED__
	+
	+#endif

File Metadata

Mime Type: text/x-diff
Expires: Sat, Dec 21, 5:07 PM (14 h, 36 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4018384
Default Alt Text: (67 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions