from ROOT import *
from array import array
from math import sqrt
from get_p import get_p
f = TFile("Navon.root", "RECREATE")
Ts = [50.]
mom = array("d", [get_p(T) for T in Ts])
print mom
moms = {
"C_xsec_abscx_piplus": mom,
"C_xsec_abscx_piminus": mom,
"Fe_xsec_abscx_piplus": mom,
"Fe_xsec_abscx_piminus": mom
}
xsecs = {
"C_xsec_abscx_piplus": array("d", [106.]),
"C_xsec_abscx_piminus": array("d", [238.]),
"Fe_xsec_abscx_piplus": array("d", [478.]),
"Fe_xsec_abscx_piminus": array("d", [1232.])
}
errs = {
"C_xsec_abscx_piplus": array("d", [21.]),
"C_xsec_abscx_piminus": array("d", [24.]),
"Fe_xsec_abscx_piplus": array("d", [48.]),
"Fe_xsec_abscx_piminus": array("d", [92.])
from ROOT import *
from array import array
from math import sqrt
from get_p import get_p
f = TFile("Gelderloos.root", "RECREATE")
moms = array("d", [get_p(410.), get_p(464.), get_p(492.)])
exs = array("d", [0.] * 3)
reacs = {
"C": array("d", [255., 244., 224.]),
"Al": array("d", [456., 451., 420.]),
"Cu": array("d", [827., 825., 847.]),
"Pb": array("d", [1922., 1888., 2005.])
}
'''errs = {
"C": array("d", [11., 10., 9.]),
"Al": array("d", [23., 27., 24.]),
"Cu": array("d", [35., 36., 34.]),
"Pb": array("d", [90., 86., 140.])
}'''
#Alternate errors. See Gelderloos
errs = {
"C": array("d", [15., 14., 12.]),
"Al": array("d", [26., 44., 29.]),
"Cu": array("d", [37., 42., 35.]),
"Pb": array("d", [95., 90., 236.])
}
f.cd()
