Skip to content
Snippets Groups Projects
Commit f22668c8 authored by Lorenzo Moneta's avatar Lorenzo Moneta
Browse files

remove running of systematics for 1/s priors

git-svn-id: http://root.cern.ch/svn/root/trunk@31858 27541ba8-7e3a-0410-8455-c3a389f83636
parent fa05d15a
No related branches found
No related tags found
No related merge requests found
...@@ -49,6 +49,7 @@ void rs701_BayesianCalculator(bool useBkg = true, double confLevel = 0.90) ...@@ -49,6 +49,7 @@ void rs701_BayesianCalculator(bool useBkg = true, double confLevel = 0.90)
RooArgSet * nuisPar = 0; RooArgSet * nuisPar = 0;
if (useBkg) nuisPar = &nuisanceParameters; if (useBkg) nuisPar = &nuisanceParameters;
//if (!useBkg) ((RooRealVar *)w->var("b"))->setVal(0);
double size = 1.-confLevel; double size = 1.-confLevel;
std::cout << "\nBayesian Result using a Flat prior " << std::endl; std::cout << "\nBayesian Result using a Flat prior " << std::endl;
...@@ -65,9 +66,9 @@ void rs701_BayesianCalculator(bool useBkg = true, double confLevel = 0.90) ...@@ -65,9 +66,9 @@ void rs701_BayesianCalculator(bool useBkg = true, double confLevel = 0.90)
c1->cd(1); c1->cd(1);
plot->Draw(); plot->Draw();
// don;t run with background events - will take too long
std::cout << "\nBayesian Result using a 1/s prior " << std::endl; std::cout << "\nBayesian Result using a 1/s prior (no background) " << std::endl;
BayesianCalculator bcalc2(data,*model,RooArgSet(*POI),*priorPOI2, nuisPar); BayesianCalculator bcalc2(data,*model,RooArgSet(*POI),*priorPOI2);
bcalc2.SetTestSize(size); bcalc2.SetTestSize(size);
SimpleInterval* interval2 = bcalc2.GetInterval(); SimpleInterval* interval2 = bcalc2.GetInterval();
cl = bcalc2.ConfidenceLevel(); cl = bcalc2.ConfidenceLevel();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment