pisa.stages.xsec package

Submodules

pisa.stages.xsec.dis_sys module

Stage to apply pre-calculated DIS uncertainties Study done by Maria Liubarska & Juan Pablo Yanez, more information available here: https://drive.google.com/open?id=1SRBgIyX6kleYqDcvop6m0SInToAVhSX6 ToDo: tech note being written, link here as soon as available

class pisa.stages.xsec.dis_sys.dis_sys(extrapolation_type='constant', extrapolation_energy_threshold=<Quantity(100, 'gigaelectron_volt')>, **std_kwargs)

Bases: Stage

Stage to apply pre-calculated DIS systematics.

Parameters:

• data

• params –

  Must contain ::

  dis_csms : quantity (dimensionless)

  Expected container keys are ..

  "true_energy"
  "bjorken_y"
  "dis"
  "nubar"
  "weights"
  

• extrapolation_type (string) – choice of [‘constant’, ‘linear’, ‘higher’]

• extrapolation_energy_threshold (float) – Below what energy (in GeV) to extrapolate Defaults to 100. CSMS not considered reliable below 50-100 GeV

Notes

Requires the events have the following keys ::

true_energy

Neutrino energy in GeV

bjorken_y

Inelasticity

dis

1 if event is DIS, else 0

apply_function()

Implement in services (subclasses of Stage)

setup_function()

Implement in services (subclasses of Stage)

pisa.stages.xsec.dis_sys.init_test(**param_kwargs)

Instantiation example

pisa.stages.xsec.genie_sys module

Stage to apply pre-calculated Genie uncertainties

pisa.stages.xsec.genie_sys.apply_genie_sys(genie_params, linear_fits, quad_fits, out)

class pisa.stages.xsec.genie_sys.genie_sys(interactions='Genie_Ma_QE, Genie_Ma_RES', names='maccqe, maccres', **std_kwargs)

Bases: Stage

Stage to apply pre-calculated Genie systematics.

Parameters:

params –

Must contain

parameters specified in interactions (dimensionless)

Notes

Requires the events have the following keys for each included interaction

linear_fit_{name}
    Genie linear coefficient for interaction {name}
quad_fit_{name}
    Genie quadratic coefficient for interaction {name}

apply_function()

Implement in services (subclasses of Stage)

setup_function()

Check the range of the parameter in the analysis. Send a warning if these are beyond +- 2sigma

pisa.stages.xsec.genie_sys.init_test(**param_kwargs)

Instantiation example

pisa.stages.xsec.nutau_xsec module

A stage to apply nutau cross-section uncertainties as implemented in https://github.com/marialiubarska/nutau_xsec It interpolates between different nutau CC cross section models as compared in this paper: https://arxiv.org/pdf/1008.2984.pdf?fname=cm&font=TypeI

pisa.stages.xsec.nutau_xsec.init_test(**param_kwargs)

Instantiation example

class pisa.stages.xsec.nutau_xsec.nutau_xsec(xsec_file='cross_sections/interp_nutau_xsec_protocol2.pckl', **std_kwargs)

Bases: Stage

Nu_tau cross-section correction to interpolate between different nutau CC cross-section models. This requires the interpolated file produced by Maria Liubarska: https://github.com/marialiubarska/nutau_xsec

Parameters:

• xsec_file ((string)) – Path to pickled interpolated function. Default is included in PISA in pisa_examples/resources/cross_sections/interp_nutau_xsec_protocol2.pckl

• params (ParamSet or sequence with which to instantiate a ParamSet.) –

  Expected params ..

  nutau_xsec_scale : quantity (dimensionless)
      Scaling between different cross-section models. The range [-1, 1]
      covers all models tested in the paper.
  

  Expected container keys are ..

  "true_energy"
  "weights"
  

apply_function()

Implement in services (subclasses of Stage)

compute_function()

Implement in services (subclasses of Stage)

setup_function()

Implement in services (subclasses of Stage)

Module contents