2 import SupernovaUtilities
as utils
8 print "\nrun : --- Parse and a flux file into a 2-d histogram" 11 from optparse
import OptionParser
12 parser = OptionParser()
13 parser.add_option(
"-m",
"--model",
14 help=
"which model to run, options are garching or livermore",
20 parser.add_option(
"-f",
"--fermi_dirac",
21 help=
"use a fermi-dirac parameterisaion",
26 parser.add_option(
"-v",
"--verbose_mode",
27 help=
"turn on verbose mode",
32 (options, args) = parser.parse_args()
34 print "run : --- Options:" 35 print "run : model: ", options.model
36 print "run : fermi-dirac: ", options.fermi
37 print "run : verbose mode: ", options.verbose
38 assert (options.model
in [
"garching",
"livermore"])
or (
39 "SND" in options.model),
"Unconfigured mode requested: %s" % options.model
47 if options.model ==
"garching":
48 raw_data = utils.ReadGarchingFluxes()
54 histograms = utils.ParameterisedDataToHistograms(raw_data,
55 parameterisation=param,
56 verbose=options.verbose)
58 elif options.model ==
"livermore":
59 raw_data = utils.ReadLivermoreFluxes()
60 histograms = utils.ParameterisedDataToHistograms(raw_data,
61 verbose=options.verbose)
65 elif "SND" in options.model:
66 if "integ" in options.model:
70 histograms = utils.ParseIntegratedSNDFile(options.model)
72 histograms = utils.ParseSNDFile(options.model, verbose=options.verbose)
76 root_output_name =
"./data/%s_FD_fluxes.root" % options.model
78 root_output_name =
"./data/%s_fluxes.root" % options.model
79 print "run : Writing histograms to file: %s" % root_output_name
80 root_file = ROOT.TFile(root_output_name,
"RECREATE")
84 histograms[
"Time_Energy_Fluence_%i" % i].
Write()
85 histograms[
"lumi_%i" % i].
Write()
86 histograms[
"lumi_mev_%i" % i].
Write()
87 histograms[
"E_%i" % i].
Write()
88 histograms[
"alpha_%i" % i].
Write()
89 histograms[
"eta_%i" % i].
Write()
90 histograms[
"T_%i" % i].
Write()