MaCh3DUNEFactory.h File Reference

#include "manager/manager.h"
#include "samplePDF/samplePDFFDBase.h"

Include dependency graph for MaCh3DUNEFactory.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
void	MakeMaCh3DuneInstance (manager *fitMan, std::vector< samplePDFFDBase * > &sample_vec, covarianceXsec *&xsec, covarianceOsc *&osc)
	Factory function that generates MaCh3 DUNE instance including configured samples.
samplePDFFDBase *	GetMaCh3DuneInstance (std::string SampleType, std::string SampleConfig, covarianceXsec *&xsec)
	Gets MaCh3 DUNE samplePDF instance.

Function Documentation

GetMaCh3DuneInstance()

samplePDFFDBase * GetMaCh3DuneInstance	(	std::string	SampleType,
		std::string	SampleConfig,
		covarianceXsec *&	xsec )

Gets MaCh3 DUNE samplePDF instance.

Parameters

SampleType	Sample type (BeamFD, BeamND, Atm)
SampleConfig	Configuration file
xsec	Cross-section covariance matrix

Returns

samplePDF instance

MakeMaCh3DuneInstance()

void MakeMaCh3DuneInstance	(	manager *	fitMan,
		std::vector< samplePDFFDBase * > &	sample_vec,
		covarianceXsec *&	xsec,
		covarianceOsc *&	osc )

Factory function that generates MaCh3 DUNE instance including configured samples.

Parameters

fitMan	Configuration Manager
sample_vec	Vector of samplePDF objects
xsec	Cross-section covariance matrix
osc	Oscillation covariance matrix

Definition at line 27 of file MaCh3DUNEFactory.cpp.

                                                                                                                                  {
 
  // there's a check inside the manager class that does this; left here for demonstrative purposes
  if (FitManager == nullptr) {
    MACH3LOG_ERROR("Didn't find a good config in input configuration");
    throw MaCh3Exception(__FILE__, __LINE__);
  }
  
  //Check that you have specified some DUNE samples
  if(!FitManager->raw()["General"]["DUNESamples"]){
    MACH3LOG_ERROR("You didn't specify any DUNESample configs to create samples from. Please add General:DUNESamples to your config");
    throw MaCh3Exception(__FILE__, __LINE__);
  }
  
  // Get inputted systematic parameters covariance matrices
  std::vector<std::string> xsecCovMatrixFile;
  if (CheckNodeExists(FitManager->raw(), "General", "Systematics", "XsecCovFile") ){
    xsecCovMatrixFile = FitManager->raw()["General"]["Systematics"]["XsecCovFile"].as<std::vector<std::string>>();
  } else {
    MACH3LOG_ERROR("Require General:Systematics:XsecCovFile node in {}, please add this to the file!", FitManager->GetFileName());
    throw MaCh3Exception(__FILE__, __LINE__);
  }
  
  // Setup the covariance matrices
  if(xsec == nullptr){
    xsec = new covarianceXsec(xsecCovMatrixFile, "xsec_cov");
  }
  else{
    MACH3LOG_INFO("covariance Xsec has already been created so I am not re-initialising the object"); 
  }
 
  //Setup the cross-section parameters
  //This should get the prior values.
  std::vector<double> XsecParVals = xsec->getNominalArray();
 
  xsec->setParameters(XsecParVals);  
  xsec->setStepScale(FitManager->raw()["General"]["Systematics"]["XsecStepScale"].as<double>());
  
  MACH3LOG_INFO("cov xsec setup");
  MACH3LOG_INFO("------------------------------");
  
  std::vector<std::string> OscMatrixFile = FitManager->raw()["General"]["Systematics"]["OscCovFile"].as<std::vector<std::string>>();
  std::string  OscMatrixName = FitManager->raw()["General"]["Systematics"]["OscCovName"].as<std::string>(); 
  std::vector<double> oscpars = FitManager->raw()["General"]["OscillationParameters"].as<std::vector<double>>();
 
  std::string OscPars = "";
 
  for (unsigned int i=0;i<oscpars.size();i++) {
    OscPars+=std::to_string(oscpars[i]);
    OscPars+=", ";
  }
  MACH3LOG_INFO("Oscillation Parameters being used: {} ", OscPars);
  
  osc = new covarianceOsc(OscMatrixFile,OscMatrixName.c_str());
  osc->setName("osc_cov");
 
  std::vector<std::string> OscFixParams   = GetFromManager<std::vector<std::string>>(FitManager->raw()["General"]["Systematics"]["OscFix"], {""});
  // Fixed xsec parameters loop
  if (OscFixParams.size() == 1 && OscFixParams.at(0) == "All") {
    for (int j = 0; j < osc->getNpars(); j++) {
      osc->toggleFixParameter(j);
    }
  } else {
    for (unsigned int j = 0; j < OscFixParams.size(); j++) {
      osc->toggleFixParameter(OscFixParams.at(j));
    }
  }
  
  MACH3LOG_INFO("Osc cov setup");
  MACH3LOG_INFO("------------------------------");
  
  // ==========================================================
  //read flat prior, fixed paramas from the config file
  std::vector<std::string> XsecFixParams   = GetFromManager<std::vector<std::string>>(FitManager->raw()["General"]["Systematics"]["XsecFix"], {""});
  
  // Fixed xsec parameters loop
  if (XsecFixParams.size() == 1 && XsecFixParams.at(0) == "All") {
    for (int j = 0; j < xsec->getNpars(); j++) {
      xsec->toggleFixParameter(j);
    }
  } else {
    for (unsigned int j = 0; j < XsecFixParams.size(); j++) {
      xsec->toggleFixParameter(XsecFixParams.at(j));
    }
  }
  MACH3LOG_INFO("xsec parameters loop done");
  
  // Fill the parameter values with their nominal values
  // should _ALWAYS_ be done before overriding with fix or flat
  xsec->setParameters();
  osc->setParameters(oscpars); 
  
  //####################################################################################
  //Create samplePDFDUNE Objs
  MACH3LOG_INFO("-------------------------------------------------------------------");
  MACH3LOG_INFO("Loading DUNE samples..");
  std::vector<std::string> DUNESampleConfigs = FitManager->raw()["General"]["DUNESamples"].as<std::vector<std::string>>();
  
  for(unsigned int Sample_i = 0 ; Sample_i < DUNESampleConfigs.size() ; Sample_i++){
 
    manager* tempSampleManager = new manager(DUNESampleConfigs[Sample_i].c_str());
    std::string SampleType = tempSampleManager->raw()["SampleType"].as<std::string>();
    
    DUNEPdfs.push_back(GetMaCh3DuneInstance(SampleType, DUNESampleConfigs[Sample_i], xsec, osc));
 
    // Pure for debugging, lets us set which weights we don't want via the manager
#if DEBUG_DUNE_WEIGHTS==1
    DUNEPdfs.back()->setWeightSwitchOffVector(FitManager->getWeightSwitchOffVector());
    DUNEPdfs.back()->setXsecWeightSwitchOffVector(FitManager->getXsecWeightSwitchOffVector());
#endif
  }
  
  return;
}