type=str,
action='store',
default="event_rate.json",
help='''file to store the output''')
parser.add_argument('-v',
'--verbose',
action='count',
default=None,
help='''set verbosity level''')
args = parser.parse_args()
#Set verbosity level
set_verbosity(args.verbose)
#Check binning
ebins, czbins = check_binning(args.osc_flux_maps)
logging.info("Defining aeff_service...")
if args.mc_mode:
logging.info(" Using effective area from EVENT DATA...")
aeff_service = AeffServiceMC(ebins,
czbins,
aeff_weight_file=args.weighted_aeff_file)
else:
logging.info(" Using effective area from PARAMETRIZATION...")
aeff_settings = from_json(find_resource(args.settings_file))
aeff_service = AeffServicePar(ebins, czbins, **aeff_settings)
event_rate_maps = get_event_rates(args.osc_flux_maps, aeff_service,
args.livetime, args.aeff_scale)
{"nue_cc": {'czbins':[...], 'ebins':[...], 'map':[...]},
"numu_cc": {...},
"nutau_cc": {...},
"nuall_nc": {...} }''')
parser.add_argument('--settings',metavar='SETTINGS',type=from_json,
default=from_json(find_resource('pid/V15_pid.json')),
help='''json file containing parameterizations of the particle ID for each event type.''')
parser.add_argument('-o', '--outfile', dest='outfile', metavar='FILE', type=str,
action='store',default="pid.json",
help='''file to store the output''')
parser.add_argument('-v', '--verbose', action='count', default=0,
help='''set verbosity level''')
args = parser.parse_args()
#Set verbosity level
set_verbosity(args.verbose)
#Check binning
ebins, czbins = check_binning(args.reco_event_maps)
#Initialize the PID service
pid_service = PIDServicePar(args.settings,ebins,czbins)
#Galculate event rates after PID
event_rate_pid = get_pid_maps(args.reco_event_maps,pid_service)
logging.info("Saving output to: %s"%args.outfile)
to_json(event_rate_pid,args.outfile)
"--outfile",
dest="outfile",
metavar="FILE",
type=str,
action="store",
default="event_rate.json",
help="""file to store the output""",
)
parser.add_argument("-v", "--verbose", action="count", default=None, help="""set verbosity level""")
args = parser.parse_args()
# Set verbosity level
set_verbosity(args.verbose)
# Check binning
ebins, czbins = check_binning(args.osc_flux_maps)
logging.info("Defining aeff_service...")
if args.mc_mode:
logging.info(" Using effective area from EVENT DATA...")
aeff_service = AeffServiceMC(ebins, czbins, aeff_weight_file=args.weighted_aeff_file)
else:
logging.info(" Using effective area from PARAMETRIZATION...")
aeff_settings = from_json(find_resource(args.settings_file))
aeff_service = AeffServicePar(ebins, czbins, **aeff_settings)
event_rate_maps = get_event_rates(args.osc_flux_maps, aeff_service, args.livetime, args.aeff_scale)
logging.info("Saving output to: %s" % args.outfile)
to_json(event_rate_maps, args.outfile)
help='JSON file containing parameterizations '
'of the particle ID \nfor each event type.')
parser.add_argument('--kernel_file', metavar='JSON', type=str, default=None,
help='JSON file containing pre-calculated PID kernels')
parser.add_argument('-o', '--outfile', dest='outfile', metavar='FILE', type=str,
action='store',default="pid.json",
help='''file to store the output''')
parser.add_argument('-v', '--verbose', action='count', default=None,
help='''set verbosity level''')
args = parser.parse_args()
#Set verbosity level
set_verbosity(args.verbose)
#Check binning
ebins, czbins = check_binning(args.reco_event_maps)
#Initialize the PID service
if args.mode=='param':
pid_service = PIDServiceParam(ebins, czbins,
pid_paramfile=args.param_file, **vars(args))
elif args.mode=='stored':
pid_service = PIDServiceKernelFile(ebins, czbins,
pid_kernelfile=args.kernel_file, **vars(args))
#Calculate event rates after PID
event_rate_pid = get_pid_maps(args.reco_event_maps, pid_service=pid_service)
logging.info("Saving output to: %s"%args.outfile)
to_json(event_rate_pid,args.outfile)
parser.add_argument('--e_reco_scale', type=float, default=1.0,
help='''Reconstructed energy scaling.''')
parser.add_argument('--cz_reco_scale', type=float, default=1.0,
help='''Reconstructed coszen scaling.''')
parser.add_argument('-o', '--outfile', dest='outfile', metavar='JSON',
type=str, action='store', default="reco.json",
help='''file to store the output''')
parser.add_argument('-v', '--verbose', action='count', default=None,
help='''set verbosity level''')
args = parser.parse_args()
# Set verbosity level
set_verbosity(args.verbose)
# Check binning
ebins, czbins = check_binning(args.event_rate_maps)
logging.info("Defining RecoService...")
if args.mode=='MC':
reco_service = RecoServiceMC(ebins, czbins,
reco_mc_wt_file=args.mc_file,
**vars(args))
elif args.mode=='param':
reco_service = RecoServiceParam(ebins, czbins,
reco_param_file=args.param_file,
**vars(args))
elif args.mode=='stored':
reco_service = RecoServiceKernelFile(ebins, czbins,
reco_kernel_file=args.kernel_file,
**vars(args))
elif args.mode=='vbwkde':
parser.add_argument('--e_reco_scale', type=float, default=1.0,
help='''Reconstructed energy scaling.''')
parser.add_argument('--cz_reco_scale', type=float, default=1.0,
help='''Reconstructed coszen scaling.''')
parser.add_argument('-o', '--outfile', dest='outfile', metavar='JSON',
type=str, action='store', default="reco.json",
help='''file to store the output''')
parser.add_argument('-v', '--verbose', action='count', default=None,
help='''set verbosity level''')
args = parser.parse_args()
# Set verbosity level
set_verbosity(args.verbose)
# Check binning
ebins, czbins = check_binning(args.event_rate_maps)
logging.info("Defining RecoService...")
if args.mode=='MC':
reco_service = RecoServiceMC(ebins, czbins,
reco_mc_wt_file=args.mc_file,
**vars(args))
elif args.mode=='param':
reco_service = RecoServiceParam(ebins, czbins,
reco_param_file=args.param_file,
**vars(args))
elif args.mode=='stored':
reco_service = RecoServiceKernelFile(ebins, czbins,
reco_kernel_file=args.kernel_file,
**vars(args))
elif args.mode=='vbwkde':
