Generate systematic response functions

Currently possible for uncorrelated interaction and flux systematic parameters.

For nue, we have several samples (intrinsic nue, oscillated sample, numu + dirt and cosmic). Need to generates splines independently from each sample with the appropriate subset of the mode list. Once the splines have been generated from the different samples, we can hadd them back to a single file for convenience.

Generate the splines

Uncorrelated flux splines

valor_sbn-make_flux_splines \
  --config analysis/sbn_numu_disapp_2020a1_1_universes_nd.xml \
  --systematic horncurrent \ // See the application for the list of possible parameters
  --variations varations_text_file.txt \ // These will exist in the data archive soon and correspond to the relevant data release
  --store-spectra \
  --save-modelist \
  --detector sbn_nd

The config must contain universe dsts and sterile params:

  <parameter_definitions>
    <filename desc="sterile_physics_parameter_values"> params/sterile_physics_parameter_values.xml </filename>
  </parameter_definitions>

Interaction splines

valor_sbn-make_interaction_splines \
  --config analysis/config.xml \ // This must be a configuration setup with the universe dsts 
  --systematic qema \ // See the application for the list of possible parameters
  --variations varations_text_file.txt \ // These will exist in the data archive soon and correspond to the relevant data release
  --store-spectra \
  --save-modelist \
  --detector sbn_nd

A simple bash script can be setup to run the full list of parameters successively in a tmux or screen session.

NB1: Check config/apps/application_name.arg for additional arguments, a number are set to defaults therefore are left off the command line here

Construct relevant covariance matrix for the parameters

Uncorrelated flux parameters

python $VALOR_SBN/macros/systematic_construction/make_flux_spline_cov.py \
--output output_file.root

Interaction parameters

python $VALOR_SBN/macros/systematic_construction/make_int_spline_cov.py \
--output output_file.root

Modify the python scripts to ensure the relevant parameters you have constructed splines for are included.

Additional modifications needed for use in the analysis

Add the relevant parameter and spline information to the following configurations, examples of the file format can generally be found in the relevant directory.

1. config/analysis/config.xml
2. config/mc/config.xml
3. config/params/config.xml
4. config/parameter_matrices/config.xml

Creating corresponding data products

Please see here for the creation of new responses data products for distribution.