def apply_nue_numu_ratio(flux_maps, nue_numu_ratio):
"""
Applies the nue_numu_ratio systematic to the flux maps
and returns the scaled maps. The actual calculation is
done by apply_ratio_scale.
"""
# keep both nu and nubar flux constant
scaled_nue_flux, scaled_numu_flux = apply_ratio_scale(
orig_maps=flux_maps,
key1='nue',
key2='numu',
ratio_scale=nue_numu_ratio,
is_flux_scale=True)
scaled_nue_bar_flux, scaled_numu_bar_flux = apply_ratio_scale(
orig_maps=flux_maps,
key1='nue_bar',
key2='numu_bar',
ratio_scale=nue_numu_ratio,
is_flux_scale=True)
flux_maps['nue']['map'] = scaled_nue_flux
flux_maps['nue_bar']['map'] = scaled_nue_bar_flux
flux_maps['numu']['map'] = scaled_numu_flux
flux_maps['numu_bar']['map'] = scaled_numu_bar_flux
return flux_maps
def apply_nue_numu_ratio(flux_maps, nue_numu_ratio):
"""
Applies the nue_numu_ratio systematic to the flux maps
and returns the scaled maps. The actual calculation is
done by apply_ratio_scale.
"""
# keep both nu and nubar flux constant
scaled_nue_flux, scaled_numu_flux = apply_ratio_scale(
orig_maps = flux_maps,
key1 = 'nue', key2 = 'numu',
ratio_scale = nue_numu_ratio,
is_flux_scale = True
)
scaled_nue_bar_flux, scaled_numu_bar_flux = apply_ratio_scale(
orig_maps = flux_maps,
key1 = 'nue_bar', key2 = 'numu_bar',
ratio_scale = nue_numu_ratio,
is_flux_scale = True
)
flux_maps['nue']['map'] = scaled_nue_flux
flux_maps['nue_bar']['map'] = scaled_nue_bar_flux
flux_maps['numu']['map'] = scaled_numu_flux
flux_maps['numu_bar']['map'] = scaled_numu_bar_flux
return flux_maps
def apply_nu_nubar_ratio(event_rate_maps, nu_nubar_ratio):
'''
Applies the nu_nubar_ratio systematic to the event rate
maps and returns the scaled maps. The actual calculation is
done by apply_ratio_scale.
'''
flavours = event_rate_maps.keys()
if 'params' in flavours: flavours.remove('params')
for flavour in flavours:
# process nu and nubar in one go
if not 'bar' in flavour:
# do this for each interaction channel (cc and nc)
for int_type in event_rate_maps[flavour].keys():
scaled_nu_rates, scaled_nubar_rates = apply_ratio_scale(
orig_maps=event_rate_maps,
key1=flavour,
key2=flavour + '_bar',
ratio_scale=nu_nubar_ratio,
is_flux_scale=False,
int_type=int_type)
event_rate_maps[flavour][int_type]['map'] = scaled_nu_rates
event_rate_maps[flavour +
'_bar'][int_type]['map'] = scaled_nubar_rates
return event_rate_maps
def apply_nu_nubar_ratio(event_rate_maps, nu_nubar_ratio):
'''
Applies the nu_nubar_ratio systematic to the event rate
maps and returns the scaled maps. The actual calculation is
done by apply_ratio_scale.
'''
flavours = event_rate_maps.keys()
if 'params' in flavours: flavours.remove('params')
for flavour in flavours:
# process nu and nubar in one go
if not 'bar' in flavour:
# do this for each interaction channel (cc and nc)
for int_type in event_rate_maps[flavour].keys():
scaled_nu_rates, scaled_nubar_rates = apply_ratio_scale(
orig_maps = event_rate_maps,
key1 = flavour, key2 = flavour+'_bar',
ratio_scale = nu_nubar_ratio,
is_flux_scale = False,
int_type = int_type
)
event_rate_maps[flavour][int_type]['map'] = scaled_nu_rates
event_rate_maps[flavour+'_bar'][int_type]['map'] = scaled_nubar_rates
return event_rate_maps
