Using P1D data holder and calculate weighted average
Code: "p1d_average.py". Programming language: Python
DMelt Version 1. Last modified: 12/08/2015. License: Pro
https://datamelt.org/code/cache/p1d_average_6147.py
To run this script using the DMelt IDE,
copy the above URL link to the menu [File]→[Read script from URL] of the DMelt IDE.
from java.awt import Color,Font
from java.util import Random
from jhplot import HPlot,P1D
from math import sqrt;
c1 = HPlot("Canvas",600,400,2, 1)
c1.visible(1)
c1.setAutoRange()
# keep 2 sets of the data
data=[];
c1.setGTitle("Average of experimental data", Color.blue) #put title
c1.viewHisto(0) # make sure starts from 0
# fill first data set
p1= P1D("data1")
p1.add(10,100,5,5)
p1.add(20,50,5,5)
p1.add(30,80,5,5)
c1.draw(p1)
# fill second data set
p2= P1D("data2")
p2.add(10,40,5,5)
p2.add(20,40,5,5)
p2.add(30,40,5,5)
data.append(p2)
c1.draw(p2)
# fill 3rd data set
p3= P1D("data3")
p3.add(10,50,5,5)
p3.add(20,45,5,5)
p3.add(30,20,5,5)
data.append(p3)
c1.draw(p3)
# add p2,p3 to p1 using weighted average
p4=p1.addAndAverage(data)
p4.setTitle("Weighted Average")
p4.setColor(Color.blue)
# show points with statistical errors in y
c1.cd(2,1)
c1.setAutoRange()
c1.draw(p4)
You see the box below because you did not login.