NPM/numass-framework
4
0
Fork 1
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed Monte-carlo integration

Browse Source
This commit is contained in:
Alexander Nozik 2017-07-07 15:12:21 +03:00
parent 9395ab5467
commit c63514999c
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
numass-main/src/main/groovy/inr/numass/scripts/OldTest.groovy
View File

@ -89,7 +89,7 @@ res.print(out);
//spectrum.counter.print(onComplete); //spectrum.counter.print(onComplete);
// //
//// fm.setPriorProb(new GaussianPrior("X", 0.47, 0.47*0.03)); //// fm.setPriorProb(new Gaussian("X", 0.47, 0.47*0.03));
//// fm.setPriorProb(new MultivariateGaussianPrior(allPars.getSubSet("X","trap"))); //// fm.setPriorProb(new MultivariateGaussianPrior(allPars.getSubSet("X","trap")));
//res = fm.runStage(res, "MINUIT", "run", "E0", "N", "bkg", "mnu2"); //res = fm.runStage(res, "MINUIT", "run", "E0", "N", "bkg", "mnu2");
//// ////

2
numass-main/src/main/groovy/inr/numass/scripts/TritiumTest.groovy
View File

@ -106,7 +106,7 @@ beta.setSuppressWarnings(true);
BayesianConfidenceLimit bm = new BayesianConfidenceLimit(); BayesianConfidenceLimit bm = new BayesianConfidenceLimit();
// bm.setPriorProb(new OneSidedUniformPrior("trap", 0, true)); // bm.setPriorProb(new OneSidedUniformPrior("trap", 0, true));
// bm.setPriorProb(new GaussianPrior("trap", 1d, 0.002)); // bm.setPriorProb(new Gaussian("trap", 1d, 0.002));
// bm.printMarginalLikelihood(Out.onComplete,"U2", res); // bm.printMarginalLikelihood(Out.onComplete,"U2", res);
FitState conf = bm.getConfidenceInterval("U2", res, ["U2", "N", "trap"]); FitState conf = bm.getConfidenceInterval("U2", res, ["U2", "N", "trap"]);