Page MenuHomeHEPForge
Files F7879025

No OneTemporary
Actions

View Options

diff --git a/data/priors/philcallums-bctune/bctuneprior_WithNormRES.root b/data/priors/philcallums-bctune/bctuneprior_WithNormRES.root
new file mode 100644
index 0000000..5f1f4e8
Binary files /dev/null and b/data/priors/philcallums-bctune/bctuneprior_WithNormRES.root differ
diff --git a/data/priors/philcallums-bctune/makeprior_normres.py b/data/priors/philcallums-bctune/makeprior_normres.py
new file mode 100644
index 0000000..bc33751
--- /dev/null
+++ b/data/priors/philcallums-bctune/makeprior_normres.py
@@ -0,0 +1,71 @@
+from ROOT import *
+
+def ConvertToGENIEDial(name, val, err):
+ if name == "MaRES":
+ convval = ((val / 1.12) - 1.0) / 0.20
+ converr = ((err / 1.12) / 0.20)
+ return "MaCCRES", convval, converr
+
+ elif name == "DISNorm":
+ convval = ((val - 100.0) / 100.0)/0.5
+ converr = (err / 100.0)/0.5
+ return "NonRESBGvnCC1pi,NonRESBGvpCC1pi,NonRESBGvbarnCC1pi,NonRESBGvbarpCC1pi", convval, converr
+
+ elif name == "RESNorm":
+ convval = ((val - 100.0)/100.0)
+ converr = err / 100.0
+ return "MINERvARW_NormCCRES", convval, converr
+
+def MakePriorROOTFile(name, fit, cor):
+
+ out = TFile("bctuneprior_" + name + ".root","RECREATE")
+ out.cd()
+
+ # Make Dial Plot
+ dials = TH1D("dials", "dials;Dials;Dial Value",
+ len(fit),0.0,float(len(fit)))
+ i = 0
+ for result in fit:
+ name, val, err = ConvertToGENIEDial( result[0], result[1], result[2] )
+ print name, val, err
+ dials.GetXaxis().SetBinLabel(i+1, name)
+ dials.SetBinContent(i+1, val)
+ dials.SetBinError(i+1, err)
+ i += 1
+ dials.Write()
+
+ # Make Cover Plot
+ covars = TH2D("covariance","covariance;Dials;Dials;Covar",
+ len(fit),0.0,float(len(fit)),
+ len(fit),0.0,float(len(fit)))
+ for i in range(len(fit)):
+ for j in range(len(fit)):
+
+ namei, vali, erri = ConvertToGENIEDial( fit[i][0], fit[i][1], fit[i][2] )
+ namej, valj, errj = ConvertToGENIEDial( fit[j][0], fit[j][1], fit[j][2] )
+
+ covars.GetXaxis().SetBinLabel(i+1,namei)
+ covars.GetYaxis().SetBinLabel(j+1,namej)
+ covars.SetBinContent(i+1,j+1, erri*errj*cor[i][j])
+ covars.Write()
+
+ # Close
+ out.Close()
+
+
+
+
+
+# Fit results [val,err]
+# Correlation should be in same order
+name_res = "WithNormRES"
+fit_res = [ ["MaRES", 0.94,0.05],
+ ["DISNorm", 46, 4],
+ ["RESNorm", 115,7] ]
+cor_res = [ [1.0, -0.3, -0.9],
+ [-0.3, 1.0, 0.2],
+ [-0.9, 0.2, 1.0] ]
+MakePriorROOTFile(name_res, fit_res, cor_res)
+
+
+

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:18 PM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3801379
Default Alt Text
(2 KB)

Event Timeline