102workspaceName =
"combined"
103modelConfigName =
"ModelConfig"
107confidenceLevel = 0.95
110additionalToysFac = 0.5
119 filename =
"results/example_combined_GaussExample_model.root"
124 print(f
"will run standard hist2workspace example")
127 print(
"\n\n---------------------")
128 print(
"Done creating example input")
129 print(
"---------------------\n\n")
143mc =
w.obj(modelConfigName)
148print(f
"Found data and ModelConfig:")
194 print(f
"Not sure what to do about this model")
197 print(f
"will use global observables for unconditional ensemble")
207 f
"\n95% interval on {firstPOI.GetName()} is : [{interval.LowerLimit(firstPOI)}, {interval.UpperLimit(firstPOI)}] "
256print(
"\nWill use these parameter points to generate pseudo data for bkg only")
271for imc
in range(nToyMC):
283 print(f
"Not sure what to do about this model")
299 if obsTSatObsUL < toyTSatObsUL:
300 CLb += (1.0) / nToyMC
301 if obsTSatObsUL <= toyTSatObsUL:
302 CLbinclusive += (1.0) / nToyMC
330c1.SaveAs(
"two-sided_upper_limit_output.pdf")
364print(f
"-2 sigma band {band2sigDown}")
365print(f
"-1 sigma band {band1sigDown} [Power Constraint)]")
366print(f
"median of band {bandMedian}")
367print(f
"+1 sigma band {band1sigUp}")
368print(f
"+2 sigma band {band2sigUp}")
371print(f
"\nobserved 95% upper-limit {interval.UpperLimit(firstPOI)}")
372print(f
"CLb strict [P(toy>obs|0)] for observed 95% upper-limit {CLb}")
373print(f
"CLb inclusive [P(toy>=obs|0)] for observed 95% upper-limit {CLbinclusive}")
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
void Print(GNN_Data &d, std::string txt="")