CombinedGenesisAnalysis.h

Go to the documentation of this file.

#ifndef _COMBINED_GENESIS_ANALYSIS_H_
#define _COMBINED_GENESIS_ANALYSIS_H_
#include "GammaProductionAna.h"
#include "GenesisNeutronSinglesAna.h"
#include "TALYSUtils.h"
#include <random>
#include <Eigen/Dense>
class CombinedGenesisAnalysis
{
public:
  CombinedGenesisAnalysis();
  void addGammaProdAna(GammaProductionAna a_gammaProd);
  void addNeutronSinglesAna(GenesisNeutronSinglesAna a_neutronProd);
 
  double talysModelFit(const double* a_params);
  void combinedMinimum(bool a_findMinimum=false);
  void plotCorrelationMatrix();
  std::map<std::string, double> generateResampledParams();
  map<double, Eigen::MatrixXf> generateUncertainCrossSections(
                                        // vector<double> a_gammaEnergies,
                                        int a_numResamples=100);
 
  GammaProductionAna m_gammaProd;
  GenesisNeutronSinglesAna m_neutronProd;
// private:
  ROOT::Math::Minimizer* m_minimum;
  Eigen::MatrixXf m_covarianceMatrix;
  Eigen::VectorXf m_minimumParams;
  int m_numParams;
  std::mt19937_64 m_randGen;
private:
  // bool m_
  void accumulateUncGammaXSec();
  void accumulateUncNeutronXSec();
  vector<Eigen::MatrixXf> m_uncElasticAngDists;
  // Eigen::MatrixXf
  vector<Eigen::VectorXf> m_uncElasticXsec;
  Eigen::MatrixXf m_uncElasticCov;
  map<double, vector<Eigen::VectorXf>> m_uncGammaProd;
  map<double, Eigen::MatrixXf> m_uncGammaProdCov;
  map<int, vector<Eigen::VectorXf>> m_uncLevelPop;
  map<int, Eigen::MatrixXf> m_uncLevelPopCov;
};
 
#endif