ana::GenieMultiverseParameters Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/XSec/GenieMultiverseSyst.h"

Classes
struct	id
struct	knob_sampling_mode
struct	name

Public Member Functions
	GenieMultiverseParameters (unsigned int nuniverses, std::string config_pathname="knob_config.txt")
	This class generates the parameters of the multi-universes. No more, no less. More...
	GenieMultiverseParameters (unsigned int nuniverses, std::vector< const ISyst * > systs)
std::vector< std::map< const ISyst *, double > >	ShiftTables ()
	Return the tables of knob shift values. More...
std::vector< std::map< const ISyst *, double > >	NuTruthShiftTables ()
std::vector< SystShifts >	GetSystShifts ()

Private Types
typedef multi_index_container< knob_sampling_mode, indexed_by< ordered_unique< tag< id >, BOOST_MULTI_INDEX_MEMBER(knob_sampling_mode, int, id)>, ordered_unique< tag< name >, BOOST_MULTI_INDEX_MEMBER(knob_sampling_mode, std::string, name)> > >	knob_sampling_mode_set

Private Member Functions
void	BuildKnobConfigTable (std::string config_pathname)
	Function for building knob configuration table. More...

Private Attributes
int	fUniverseSeed
bool	fIsNuTruth = false
std::vector< std::map< const ISyst *, double > >	fShiftTables
	Container of shift tables for all universes. More...
std::vector< std::map< const ISyst *, double > >	fNuTruthShiftTables
	Container of neutrino truth shift tables for all universes. More...
std::vector< SystShifts >	fSystShifts
knob_sampling_mode_set	fKnobConfigTable
	knob configuration table More...

Detailed Description

Definition at line 237 of file GenieMultiverseSyst.h.

Member Typedef Documentation

typedef multi_index_container< knob_sampling_mode, indexed_by< ordered_unique< tag<id>, BOOST_MULTI_INDEX_MEMBER(knob_sampling_mode,int,id)>, ordered_unique< tag<name>,BOOST_MULTI_INDEX_MEMBER(knob_sampling_mode,std::string,name)> > > ana::GenieMultiverseParameters::knob_sampling_mode_set

private

Definition at line 281 of file GenieMultiverseSyst.h.

Constructor & Destructor Documentation

ana::GenieMultiverseParameters::GenieMultiverseParameters	(	unsigned int	nuniverses,
		std::string	config_pathname = "knob_config.txt"
	)

This class generates the parameters of the multi-universes. No more, no less.

Definition at line 710 of file GenieMultiverseSyst.cxx.

References BuildKnobConfigTable(), fKnobConfigTable, fShiftTables, fSystShifts, fUniverseSeed, ana::getAllXsecNuTruthSysts_2017(), ana::GetGenieKnobSyst(), nuniverses, central_limit::rand, and systs.

                               :
    fUniverseSeed(1001)
  {
    // Retrieve knob configuration from file.
    BuildKnobConfigTable(config_pathname);
    // Generate tables of knobs for all universes.
    // The first entry corresponds to the nominal universe.
    for(unsigned int uIdx = 0; uIdx < nuniverses; uIdx++)
    {
      map<const ISyst*, double> knob_weight_table;
      for(auto& idIt: fKnobConfigTable.get<0>())
      {
        rwgt::ReweightLabel_t knoblabel = (rwgt::ReweightLabel_t)idIt.id;
        double nsigma = 0;
        if(uIdx > 0){
          TRandom3 rand(fUniverseSeed++);
          if(idIt.mode == 0) continue;
          else if(idIt.mode == 1) nsigma = rand.Gaus();
          else if(idIt.mode == 2) nsigma = rand.Integer(2);
        }
        // For tags later than 2017/08/14, when the underlying framework
        // underwent an overhaul, this conditional check is needed.
        vector<const ISyst*> systs = getAllXsecNuTruthSysts_2017();
        if(find(systs.begin(), systs.end(), GetGenieKnobSyst(knoblabel)) != systs.end())
          knob_weight_table[GetGenieKnobSyst(knoblabel)] = nsigma;
      }
      fShiftTables.push_back(knob_weight_table);
      fSystShifts.push_back(SystShifts(knob_weight_table));
    }
  }

ana::GenieMultiverseParameters::GenieMultiverseParameters	(	unsigned int	nuniverses,
		std::vector< const ISyst * >	systs
	)

Definition at line 744 of file GenieMultiverseSyst.cxx.

References fIsNuTruth, fShiftTables, fSystShifts, fUniverseSeed, nuniverses, and central_limit::rand.

                                   :
    fUniverseSeed(1001)
  {
    // Generate tables of knobs for all universes.
    // The first entry corresponds to the nominal universe.
    for(unsigned int uIdx = 0; uIdx < nuniverses; uIdx++)
    {
      map<const ISyst*, double> knob_weight_table;
      for(auto& knobIt: sel_systs)
      {
        double nsigma = 0;
        if(uIdx > 0){
          TRandom3 rand(fUniverseSeed++);
          nsigma = rand.Gaus();
        }
          knob_weight_table[knobIt] = nsigma;
      }
      fShiftTables.push_back(knob_weight_table);
      fSystShifts.push_back(SystShifts(knob_weight_table));
    }
    fIsNuTruth = false;
  }

Member Function Documentation

void ana::GenieMultiverseParameters::BuildKnobConfigTable ( std::string  config_pathname )

private

Function for building knob configuration table.

Definition at line 769 of file GenieMultiverseSyst.cxx.

References om::cerr, allTimeWatchdog::endl, exit(), MakeMiniprodValidationCuts::f, fKnobConfigTable, ana::GenieMultiverseNormalizedSpectra::GenieMultiverseNormalizedSpectra(), MECModelEnuComparisons::i, and allTimeWatchdog::index.

Referenced by GenieMultiverseParameters().

  {
    // check file existence
    ifstream f(fpn.c_str());
    if(!f.good()){
      cerr << "ERROR: GENIE knob configuration file \033[1;31m" << fpn << "\033[0m is not found." << endl;
      exit(0);
    }
    
    string knobname;
    int sample_mode;
    int index = 0;
    while(f >> knobname >> sample_mode)
      fKnobConfigTable.insert(knob_sampling_mode(index++,knobname,sample_mode));
  }

std::vector<SystShifts> ana::GenieMultiverseParameters::GetSystShifts ( )

inline

Definition at line 248 of file GenieMultiverseSyst.h.

Referenced by containment_optimization(), demo_cut_optimization(), demo_trivial_xsec_analysis(), fiducial_opt(), fiducial_optimization(), getCrossSectionAnalysis_Spectra_systematics(), getCrossSectionAnalysisSpectra(), getFitTemplates(), getSpectra_ForFitting(), muonid_opt(), muonid_optimization(), pion_multiverse(), prongcvn_optimization(), specprod_numuccinc(), vertex_optimize(), and wrong_sign_uncertainty().

{return fSystShifts;};

std::vector<std::map<const ISyst*, double> > ana::GenieMultiverseParameters::NuTruthShiftTables ( )

inline

Definition at line 247 of file GenieMultiverseSyst.h.

{return fNuTruthShiftTables;};

std::vector<std::map<const ISyst*, double> > ana::GenieMultiverseParameters::ShiftTables ( )

inline

Return the tables of knob shift values.

Definition at line 246 of file GenieMultiverseSyst.h.

Referenced by uncertainty1png(), and uncertainty2png().

{return fShiftTables;};

Member Data Documentation

bool ana::GenieMultiverseParameters::fIsNuTruth = false

private

Definition at line 254 of file GenieMultiverseSyst.h.

Referenced by GenieMultiverseParameters().

knob_sampling_mode_set ana::GenieMultiverseParameters::fKnobConfigTable

private

knob configuration table

Definition at line 283 of file GenieMultiverseSyst.h.

Referenced by BuildKnobConfigTable(), ana::GenieMultiverseNormalizedSpectra::CreateUniverse(), and GenieMultiverseParameters().

std::vector<std::map<const ISyst*, double> > ana::GenieMultiverseParameters::fNuTruthShiftTables

private

Container of neutrino truth shift tables for all universes.

Definition at line 257 of file GenieMultiverseSyst.h.

std::vector<std::map<const ISyst*, double> > ana::GenieMultiverseParameters::fShiftTables

private

Container of shift tables for all universes.

Definition at line 256 of file GenieMultiverseSyst.h.

Referenced by GenieMultiverseParameters().

std::vector<SystShifts> ana::GenieMultiverseParameters::fSystShifts

private

Definition at line 258 of file GenieMultiverseSyst.h.

Referenced by GenieMultiverseParameters().

int ana::GenieMultiverseParameters::fUniverseSeed

private

Random seed sequence. Need to have each universe with the same physical parameters. Start from 1001, giving way to the PPFX multi-universe.

Definition at line 248 of file GenieMultiverseSyst.h.

Referenced by ana::GenieMultiverseNormalizedSpectra::CreateUniverse(), and GenieMultiverseParameters().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/XSec/GenieMultiverseSyst.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/XSec/GenieMultiverseSyst.cxx