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