Showing multiple gaussian histograms with GROOT
Code: "groot_multigauss.java". Programming language: Java
DMelt Version 1.9. Last modified: 07/16/1971. License: Free
https://datamelt.org/code/cache/groot_multigauss_6115.java
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.
import java.awt.Dimension;
import java.awt.Toolkit;
import java.util.Random;
import javax.swing.JFrame;
import org.jlab.groot.data.H1F;
import org.jlab.groot.fitter.DataFitter;
import org.jlab.groot.graphics.EmbeddedCanvas;
import org.jlab.groot.math.F1D;
import org.jlab.groot.math.RandomFunc;
public class groot_multigauss {
public static void main(String[] args) {
Dimension screensize = Toolkit.getDefaultToolkit().getScreenSize();
JFrame frame = new JFrame("GROOT DEMO");
frame.setSize((int)(screensize.getHeight()*.75*1.618), (int) (screensize.getHeight()*.75));
EmbeddedCanvas c1 = new EmbeddedCanvas();
//c1.showFPS(true);
c1.divide(4, 4);
Random rand = new Random();
H1F[] h1 = new H1F[16];
for (int i = 0; i < h1.length-1; i++) {
h1[i] = new H1F("h" + i, "", 200, -5.0, 5.0);
h1[i].setTitleX("Randomly Generated Function");
h1[i].setTitleY("Counts");
h1[i].setOptStat(0);
F1D f1 = new F1D("f1","[amp]*gaus(x,[mean],[sigma])", -5.0, 5.0);
f1.setParameter(0, 120.0);
f1.setParameter(1, (-3.0 + rand.nextDouble() * 6));
f1.setParameter(2, .4 + (rand.nextDouble() * 1));
RandomFunc rndm = new RandomFunc(f1);
for (int j = 0; j < 34000; j++) {
h1[i].fill(rndm.random());
}
h1[i].setLineWidth(2);
h1[i].setLineColor(21);
h1[i].setFillColor(30 + (i % 4) + 2);
c1.cd(i);
c1.getPad(i).setTitle("OptionsPanel Demo");
String optStatString = "";
for(int j=0; j<3-i%4; j++){
optStatString +="1";
}
optStatString+="00";
f1.setOptStat(Integer.parseInt(optStatString));
c1.getPad(i).getAxisX().setRange(-6.0, 5.0);
c1.draw(h1[i]);
f1.setParameter(0, h1[i].getEntries()); //Due to initial parameter estimates not existing
DataFitter.fit(f1, h1[i], "Q"); //No options uses error for sigma
c1.draw(f1,"same");
f1.setLineColor(30 + (i % 4) + 2);
f1.setLineWidth(3);
f1.setLineStyle(i%4);
c1.setFont("HanziPen TC");
c1.setTitleSize(16);
c1.setAxisTitleSize(14);
c1.setAxisLabelSize(12);
}
frame.add(c1);
frame.setLocationRelativeTo(null);
frame.setVisible(true);
}
}