121 for (
int i = 0; i <
data.numEntries(); i++) {
122 value +=
data.get(i)->getRealValue(fColumnName.c_str());
126 virtual const TString GetVarName()
const {
return fColumnName; }
141void HybridStandardForm(
int ntoys = 6000)
177 w->factory(
"Uniform::f(m[0,1])");
178 w->factory(
"ExtendPdf::px(f,sum::splusb(s[0,0,100],b[100,0.1,300]))");
179 w->factory(
"Poisson::py(y[100,0.1,500],prod::taub(tau[1.],b))");
180 w->factory(
"PROD::model(px,py)");
181 w->factory(
"Uniform::prior_b(b)");
195 double p_Bi = NumberCountingUtils::BinomialWithTauObsP(150, 100, 1);
196 double Z_Bi = NumberCountingUtils::BinomialWithTauObsZ(150, 100, 1);
197 cout <<
"-----------------------------------------" << endl;
198 cout <<
"Part 3" << endl;
199 std::cout <<
"Z_Bi p-value (analytic): " <<
p_Bi << std::endl;
200 std::cout <<
"Z_Bi significance (analytic): " <<
Z_Bi << std::endl;
224 w->defineSet(
"obs",
"m");
225 w->defineSet(
"poi",
"s");
230 std::unique_ptr<RooDataSet>
data{
w->pdf(
"px")->generate(*
w->set(
"obs"), 150)};
237 b_model.SetObservables(*
w->set(
"obs"));
238 b_model.SetParametersOfInterest(*
w->set(
"poi"));
239 w->var(
"s")->setVal(0.0);
240 b_model.SetSnapshot(*
w->set(
"poi"));
246 sb_model.SetParametersOfInterest(*
w->set(
"poi"));
247 w->var(
"s")->setVal(50.0);
281 w->factory(
"Gaussian::gauss_prior(b,y, expr::sqrty('sqrt(y)',y))");
285 w->factory(
"Lognormal::lognorm_prior(b,y, expr::kappa('1+1./sqrt(y)',y))");
309 hc1.ForcePriorNuisanceAlt(*
w->pdf(
"py"));
310 hc1.ForcePriorNuisanceNull(*
w->pdf(
"py"));
327 cout <<
"-----------------------------------------" << endl;
328 cout <<
"Part 4" << endl;
359 hc2.ForcePriorNuisanceAlt(*
w->pdf(
"py"));
360 hc2.ForcePriorNuisanceNull(*
w->pdf(
"py"));
376 cout <<
"-----------------------------------------" << endl;
377 cout <<
"Part 5" << endl;
#define ClassDef(name, id)
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void data
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void value
Abstract base class for binned and unbinned datasets.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
static RooMsgService & instance()
Return reference to singleton instance.
Same purpose as HybridCalculatorOriginal, but different implementation.
This class provides the plots for the result of a study performed with any of the HypoTestCalculatorG...
HypoTestResult is a base class for results from hypothesis tests.
ModelConfig is a simple class that holds configuration information specifying how a model should be u...
NumEventsTestStat is a simple implementation of the TestStatistic interface used for simple number co...
TestStatistic class that returns -log(L[null] / L[alt]) where L is the likelihood.
TestStatistic is an interface class to provide a facility for construction test statistics distributi...
ToyMCSampler is an implementation of the TestStatSampler interface.
Persistable container for RooFit projects.
void Divide(Int_t nx=1, Int_t ny=1, Float_t xmargin=0.01, Float_t ymargin=0.01, Int_t color=0) override
Automatic pad generation by division.
void Start(Bool_t reset=kTRUE)
Start the stopwatch.
void Stop()
Stop the stopwatch.
void Print(Option_t *option="") const override
Print the real and cpu time passed between the start and stop events.
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
MsgLevel
Verbosity level for RooMsgService::StreamConfig in RooMsgService.
Namespace for the RooStats classes.