Variables | |
| addPdf0 = ROOT.RooAddPdf("addPdf0", "multiPdf0 + extra expo", ROOT.RooArgList(multiPdf0, gaussExtra), frac0) | |
| addPdf1 = ROOT.RooAddPdf("addPdf1", "multiPdf1 + extra gauss", ROOT.RooArgList(multiPdf1, expoExtra), frac1) | |
| c = ROOT.TCanvas("c_rf619", "NLL vs Mean for Different Discrete Combinations", 1200, 400) | |
| c0 = ROOT.RooRealVar("c0", "Cheby coeff 0", -1.0, -1.0, 1.0) | |
| c1 = ROOT.RooRealVar("c1", "Cheby coeff 1", 0.4, 0.05, 0.5) | |
| catIndex = ROOT.RooCategory("catIndex", "Category") | |
| cheb = ROOT.RooChebychev("cheb", "Chebyshev PDF", x, chebCoeffs) | |
| chebCoeffs = ROOT.RooArgList(c0, c1) | |
| list | colors = [ROOT.kRed, ROOT.kBlue, ROOT.kGreen + 2, ROOT.kMagenta, ROOT.kOrange + 7] |
| combined_data = ROOT.RooDataSet("data", "combined", ROOT.RooArgSet(x, catIndex)) | |
| combined_vars = ROOT.RooArgSet(x, catIndex) | |
| list | combosToPlot = [[i, j] for i in range(pdfIndex0.numTypes()) for j in range(pdfIndex1.numTypes())] |
| data0 = addPdf0.generate(ROOT.RooArgSet(x), 800) | |
| data1 = addPdf1.generate(ROOT.RooArgSet(x), 1000) | |
| expo1 = ROOT.RooExponential("expo1", "Exponential 1", x, lambda1) | |
| expoExtra = ROOT.RooExponential("expoExtra", "extra exponential", x, lambdaExtra) | |
| frac0 = ROOT.RooRealVar("frac0", "fraction for cat0", 0.7, 0.0, 1.0) | |
| frac1 = ROOT.RooRealVar("frac1", "fraction for cat1", 0.5, 0.0, 1.0) | |
| frame0 = x.frame() | |
| frame1 = x.frame() | |
| g = ROOT.TGraph(nMeanPoints) | |
| gauss1 = ROOT.RooGaussian("gauss1", "Gaussian", x, mean, sigmaG) | |
| gaussExtra = ROOT.RooGaussian("gaussExtra", "extra Gaussian", x, mean, sigmaExtra) | |
| list | graphs = [] |
| lambda1 = ROOT.RooRealVar("lambda1", "slope1", -0.025, -0.1, -0.02) | |
| lambdaExtra = ROOT.RooRealVar("lambdaExtra", "extra slope", -0.05, -1.0, -0.01) | |
| landau1 = ROOT.RooLandau("landau1", "Landau", x, mean, sigmaL) | |
| leg0 = ROOT.TLegend(0.6, 0.7, 0.9, 0.9) | |
| leg1 = ROOT.TLegend(0.6, 0.7, 0.9, 0.9) | |
| list | markers = [20, 21, 22, 23, 33] |
| mean = ROOT.RooRealVar("mean", "shared mean", 25, 0, 50) | |
| int | meanMax = 33 |
| int | meanMin = 17 |
| int | meanVal = meanMin + i * (meanMax - meanMin) / (nMeanPoints - 1) |
| mg = ROOT.TMultiGraph() | |
| minim = ROOT.RooMinimizer(nll) | |
| multiPdf0 = ROOT.RooMultiPdf("multiPdf0", "multiPdf0", pdfIndex0, ROOT.RooArgList(expo1, cheb)) | |
| multiPdf1 = ROOT.RooMultiPdf("multiPdf1", "multiPdf1", pdfIndex1, ROOT.RooArgList(gauss1, landau1)) | |
| nll = simPdf.createNLL(combined_data) | |
| int | nMeanPoints = 40 |
| pdfIndex0 = ROOT.RooCategory("pdfIndex0", "pdf index 0") | |
| pdfIndex1 = ROOT.RooCategory("pdfIndex1", "pdf index 1") | |
| profileGraph = ROOT.TGraph(nMeanPoints) | |
| sigmaExtra = ROOT.RooRealVar("sigmaExtra", "extra Gaussian width", 3.0, 1.0, 6.0) | |
| sigmaG = ROOT.RooRealVar("sigmaG", "Gaussian width", 2.0, 0.0, 5.0) | |
| sigmaL = ROOT.RooRealVar("sigmaL", "Landau width", 3.0, 1.0, 8.0) | |
| simPdf = ROOT.RooSimultaneous("simPdf", "simultaneous model", catIndex) | |
| x = ROOT.RooRealVar("x", "Observable", 0, 50) | |
| rf619_discrete_profiling.addPdf0 = ROOT.RooAddPdf("addPdf0", "multiPdf0 + extra expo", ROOT.RooArgList(multiPdf0, gaussExtra), frac0) |
Definition at line 51 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.addPdf1 = ROOT.RooAddPdf("addPdf1", "multiPdf1 + extra gauss", ROOT.RooArgList(multiPdf1, expoExtra), frac1) |
Definition at line 54 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.c = ROOT.TCanvas("c_rf619", "NLL vs Mean for Different Discrete Combinations", 1200, 400) |
Definition at line 170 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.c0 = ROOT.RooRealVar("c0", "Cheby coeff 0", -1.0, -1.0, 1.0) |
Definition at line 22 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.c1 = ROOT.RooRealVar("c1", "Cheby coeff 1", 0.4, 0.05, 0.5) |
Definition at line 23 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.catIndex = ROOT.RooCategory("catIndex", "Category") |
Definition at line 55 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.cheb = ROOT.RooChebychev("cheb", "Chebyshev PDF", x, chebCoeffs) |
Definition at line 25 of file rf619_discrete_profiling.py.
Definition at line 24 of file rf619_discrete_profiling.py.
| list rf619_discrete_profiling.colors = [ROOT.kRed, ROOT.kBlue, ROOT.kGreen + 2, ROOT.kMagenta, ROOT.kOrange + 7] |
Definition at line 128 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.combined_data = ROOT.RooDataSet("data", "combined", ROOT.RooArgSet(x, catIndex)) |
Definition at line 100 of file rf619_discrete_profiling.py.
Definition at line 102 of file rf619_discrete_profiling.py.
| list rf619_discrete_profiling.combosToPlot = [[i, j] for i in range(pdfIndex0.numTypes()) for j in range(pdfIndex1.numTypes())] |
Definition at line 126 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.data0 = addPdf0.generate(ROOT.RooArgSet(x), 800) |
Definition at line 64 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.data1 = addPdf1.generate(ROOT.RooArgSet(x), 1000) |
Definition at line 65 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.expo1 = ROOT.RooExponential("expo1", "Exponential 1", x, lambda1) |
Definition at line 20 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.expoExtra = ROOT.RooExponential("expoExtra", "extra exponential", x, lambdaExtra) |
Definition at line 32 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.frac0 = ROOT.RooRealVar("frac0", "fraction for cat0", 0.7, 0.0, 1.0) |
Definition at line 50 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.frac1 = ROOT.RooRealVar("frac1", "fraction for cat1", 0.5, 0.0, 1.0) |
Definition at line 53 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.frame0 = x.frame() |
Definition at line 68 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.frame1 = x.frame() |
Definition at line 84 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.g = ROOT.TGraph(nMeanPoints) |
Definition at line 133 of file rf619_discrete_profiling.py.
Definition at line 39 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.gaussExtra = ROOT.RooGaussian("gaussExtra", "extra Gaussian", x, mean, sigmaExtra) |
Definition at line 47 of file rf619_discrete_profiling.py.
| list rf619_discrete_profiling.graphs = [] |
Definition at line 131 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.lambda1 = ROOT.RooRealVar("lambda1", "slope1", -0.025, -0.1, -0.02) |
Definition at line 19 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.lambdaExtra = ROOT.RooRealVar("lambdaExtra", "extra slope", -0.05, -1.0, -0.01) |
Definition at line 31 of file rf619_discrete_profiling.py.
Definition at line 40 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.leg0 = ROOT.TLegend(0.6, 0.7, 0.9, 0.9) |
Definition at line 78 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.leg1 = ROOT.TLegend(0.6, 0.7, 0.9, 0.9) |
Definition at line 94 of file rf619_discrete_profiling.py.
| list rf619_discrete_profiling.markers = [20, 21, 22, 23, 33] |
Definition at line 129 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.mean = ROOT.RooRealVar("mean", "shared mean", 25, 0, 50) |
Definition at line 35 of file rf619_discrete_profiling.py.
| int rf619_discrete_profiling.meanMax = 33 |
Definition at line 124 of file rf619_discrete_profiling.py.
| int rf619_discrete_profiling.meanMin = 17 |
Definition at line 123 of file rf619_discrete_profiling.py.
| int rf619_discrete_profiling.meanVal = meanMin + i * (meanMax - meanMin) / (nMeanPoints - 1) |
Definition at line 150 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.mg = ROOT.TMultiGraph() |
Definition at line 187 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.minim = ROOT.RooMinimizer(nll) |
Definition at line 116 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.multiPdf0 = ROOT.RooMultiPdf("multiPdf0", "multiPdf0", pdfIndex0, ROOT.RooArgList(expo1, cheb)) |
Definition at line 28 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.multiPdf1 = ROOT.RooMultiPdf("multiPdf1", "multiPdf1", pdfIndex1, ROOT.RooArgList(gauss1, landau1)) |
Definition at line 43 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.nll = simPdf.createNLL(combined_data) |
Definition at line 115 of file rf619_discrete_profiling.py.
| int rf619_discrete_profiling.nMeanPoints = 40 |
Definition at line 122 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.pdfIndex0 = ROOT.RooCategory("pdfIndex0", "pdf index 0") |
Definition at line 27 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.pdfIndex1 = ROOT.RooCategory("pdfIndex1", "pdf index 1") |
Definition at line 42 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.profileGraph = ROOT.TGraph(nMeanPoints) |
Definition at line 140 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.sigmaExtra = ROOT.RooRealVar("sigmaExtra", "extra Gaussian width", 3.0, 1.0, 6.0) |
Definition at line 46 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.sigmaG = ROOT.RooRealVar("sigmaG", "Gaussian width", 2.0, 0.0, 5.0) |
Definition at line 36 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.sigmaL = ROOT.RooRealVar("sigmaL", "Landau width", 3.0, 1.0, 8.0) |
Definition at line 37 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.simPdf = ROOT.RooSimultaneous("simPdf", "simultaneous model", catIndex) |
Definition at line 59 of file rf619_discrete_profiling.py.
| rf619_discrete_profiling.x = ROOT.RooRealVar("x", "Observable", 0, 50) |
Definition at line 16 of file rf619_discrete_profiling.py.