Execution Time83.55s

Test: tutorial-roostats-TwoSidedFrequentistUpperLimitWithBands (Passed)
Build: master-x86_64-centos7-gcc62-opt-no-rt-cxxmodules (olsnba08.cern.ch) on 2019-11-14 01:02:24
Repository revision: 32b17abcda23e44b64218a42d0ca69cb30cda7e0

Test Timing: Passed
Processors1

Show Command Line
Display graphs:

Test output
Processing /data/sftnight/workspace/root-benchmark-no-rt-cxxmodules/BUILDTYPE/Release/COMPILER/gcc62/LABEL/performance-sandy-cc7/root/tutorials/roostats/TwoSidedFrequentistUpperLimitWithBands.C...

[NON-XML-CHAR-0x1B][1mRooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby[NON-XML-CHAR-0x1B][0m 
                Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University
                All rights reserved, please read http://roofit.sourceforge.net/license.txt

Found data and ModelConfig:

=== Using the following for ModelConfig ===
Observables:             RooArgSet:: = (obs_x_channel1,weightVar,channelCat)
Parameters of Interest:  RooArgSet:: = (SigXsecOverSM)
Nuisance Parameters:     RooArgSet:: = (alpha_syst2,alpha_syst3,gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1)
Global Observables:      RooArgSet:: = (nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)
PDF:                     RooSimultaneous::simPdf[ indexCat=channelCat channel1=model_channel1 ] = 0.174888

FeldmanCousins: ntoys per point = 499
FeldmanCousins: nEvents per toy will fluctuate about  expectation
will use global observables for unconditional ensemble
RooArgSet:: = (nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)

=== Using the following for ModelConfig ===
Observables:             RooArgSet:: = (obs_x_channel1,weightVar,channelCat)
Parameters of Interest:  RooArgSet:: = (SigXsecOverSM)
Nuisance Parameters:     RooArgSet:: = (alpha_syst2,alpha_syst3,gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1)
Global Observables:      RooArgSet:: = (nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)
PDF:                     RooSimultaneous::simPdf[ indexCat=channelCat channel1=model_channel1 ] = 0.174888

FeldmanCousins: Model has nuisance parameters, will do profile construction
FeldmanCousins: # points to test = 20
lookup index = 0
NeymanConstruction: Prog: 1/20 total MC = 499 this test stat = 1.54009
 SigXsecOverSM=0.075 alpha_syst2=0.656048 alpha_syst3=0.244567 gamma_stat_channel1_bin_0=1.03397 gamma_stat_channel1_bin_1=1.04971 [-1e+30, 1.4945]  in interval = 0

NeymanConstruction: Prog: 2/20 total MC = 499 this test stat = 1.12274
 SigXsecOverSM=0.225 alpha_syst2=0.549579 alpha_syst3=0.211046 gamma_stat_channel1_bin_0=1.02849 gamma_stat_channel1_bin_1=1.04283 [-1e+30, 1.70069]  in interval = 1

NeymanConstruction: Prog: 3/20 total MC = 499 this test stat = 0.77243
 SigXsecOverSM=0.375 alpha_syst2=0.44775 alpha_syst3=0.177828 gamma_stat_channel1_bin_0=1.02318 gamma_stat_channel1_bin_1=1.03602 [-1e+30, 1.41504]  in interval = 1

NeymanConstruction: Prog: 4/20 total MC = 499 this test stat = 0.488545
 SigXsecOverSM=0.525 alpha_syst2=0.349863 alpha_syst3=0.144895 gamma_stat_channel1_bin_0=1.01804 gamma_stat_channel1_bin_1=1.02929 [-1e+30, 1.80371]  in interval = 1

NeymanConstruction: Prog: 5/20 total MC = 499 this test stat = 0.270334
 SigXsecOverSM=0.675 alpha_syst2=0.255203 alpha_syst3=0.11224 gamma_stat_channel1_bin_0=1.01307 gamma_stat_channel1_bin_1=1.02264 [-1e+30, 1.75591]  in interval = 1

NeymanConstruction: Prog: 6/20 total MC = 499 this test stat = 0.11689
 SigXsecOverSM=0.825 alpha_syst2=0.16608 alpha_syst3=0.0821439 gamma_stat_channel1_bin_0=1.00838 gamma_stat_channel1_bin_1=1.01547 [-1e+30, 1.89781]  in interval = 1

NeymanConstruction: Prog: 7/20 total MC = 499 this test stat = 0.0271562
 SigXsecOverSM=0.975 alpha_syst2=0.0745784 alpha_syst3=0.0487673 gamma_stat_channel1_bin_0=1.00375 gamma_stat_channel1_bin_1=1.00931 [-1e+30, 2.02136]  in interval = 1

NeymanConstruction: Prog: 8/20 total MC = 499 this test stat = 0.000125096
 SigXsecOverSM=1.125 alpha_syst2=-0.0140996 alpha_syst3=0.0151511 gamma_stat_channel1_bin_0=0.999297 gamma_stat_channel1_bin_1=1.00317 [-1e+30, 2.36523]  in interval = 1

NeymanConstruction: Prog: 9/20 total MC = 499 this test stat = 0.0347209
 SigXsecOverSM=1.275 alpha_syst2=-0.101327 alpha_syst3=-0.0176472 gamma_stat_channel1_bin_0=0.994967 gamma_stat_channel1_bin_1=0.99714 [-1e+30, 1.79773]  in interval = 1

NeymanConstruction: Prog: 10/20 total MC = 499 this test stat = 0.12986
 SigXsecOverSM=1.425 alpha_syst2=-0.185801 alpha_syst3=-0.0506339 gamma_stat_channel1_bin_0=0.990814 gamma_stat_channel1_bin_1=0.991141 [-1e+30, 2.11797]  in interval = 1

NeymanConstruction: Prog: 11/20 total MC = 499 this test stat = 0.284381
 SigXsecOverSM=1.575 alpha_syst2=-0.266353 alpha_syst3=-0.0784096 gamma_stat_channel1_bin_0=0.987025 gamma_stat_channel1_bin_1=0.984445 [-1e+30, 2.41483]  in interval = 1

NeymanConstruction: Prog: 12/20 total MC = 499 this test stat = 0.497424
 SigXsecOverSM=1.725 alpha_syst2=-0.3458 alpha_syst3=-0.10947 gamma_stat_channel1_bin_0=0.983247 gamma_stat_channel1_bin_1=0.978371 [-1e+30, 1.95432]  in interval = 1

NeymanConstruction: Prog: 13/20 total MC = 499 this test stat = 0.767852
 SigXsecOverSM=1.875 alpha_syst2=-0.422672 alpha_syst3=-0.140541 gamma_stat_channel1_bin_0=0.979595 gamma_stat_channel1_bin_1=0.972417 [-1e+30, 2.03922]  in interval = 1

NeymanConstruction: Prog: 14/20 total MC = 499 this test stat = 1.09468
 SigXsecOverSM=2.025 alpha_syst2=-0.496602 alpha_syst3=-0.171004 gamma_stat_channel1_bin_0=0.976094 gamma_stat_channel1_bin_1=0.966484 [-1e+30, 1.90696]  in interval = 1

NeymanConstruction: Prog: 15/20 total MC = 499 this test stat = 1.47704
 SigXsecOverSM=2.175 alpha_syst2=-0.567566 alpha_syst3=-0.201186 gamma_stat_channel1_bin_0=0.972718 gamma_stat_channel1_bin_1=0.960626 [-1e+30, 1.41513]  in interval = 0

NeymanConstruction: Prog: 16/20 total MC = 499 this test stat = 1.91389
 SigXsecOverSM=2.325 alpha_syst2=-0.63553 alpha_syst3=-0.231072 gamma_stat_channel1_bin_0=0.96946 gamma_stat_channel1_bin_1=0.954845 [-1e+30, 1.98699]  in interval = 1

NeymanConstruction: Prog: 17/20 total MC = 499 this test stat = 2.40432
 SigXsecOverSM=2.475 alpha_syst2=-0.700748 alpha_syst3=-0.261131 gamma_stat_channel1_bin_0=0.966278 gamma_stat_channel1_bin_1=0.949239 [-1e+30, 1.50323]  in interval = 0

NeymanConstruction: Prog: 18/20 total MC = 499 this test stat = 2.94737
 SigXsecOverSM=2.625 alpha_syst2=-0.763083 alpha_syst3=-0.290538 gamma_stat_channel1_bin_0=0.963225 gamma_stat_channel1_bin_1=0.943654 [-1e+30, 1.5594]  in interval = 0

NeymanConstruction: Prog: 19/20 total MC = 499 this test stat = 3.54205
 SigXsecOverSM=2.775 alpha_syst2=-0.822865 alpha_syst3=-0.319576 gamma_stat_channel1_bin_0=0.960275 gamma_stat_channel1_bin_1=0.938148 [-1e+30, 1.29457]  in interval = 0

NeymanConstruction: Prog: 20/20 total MC = 499 this test stat = 4.18752
 SigXsecOverSM=2.925 alpha_syst2=-0.880381 alpha_syst3=-0.348221 gamma_stat_channel1_bin_0=0.957426 gamma_stat_channel1_bin_1=0.932721 [-1e+30, 1.21599]  in interval = 0

[#1] INFO:Eval -- 14 points in interval

95% interval on SigXsecOverSM is : [0.225, 2.325] 
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:Minization -- createNLL picked up cached consraints from workspace with 6 entries
[#1] INFO:Minization --  Including the following constraint terms in minimization: (lumiConstraint,alpha_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_simPdf_obsData_with_constr_Profile[SigXsecOverSM]) Creating instance of MINUIT
[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_simPdf_obsData_with_constr) Summation contains a RooNLLVar, using its error level
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_simPdf_obsData_with_constr_Profile[SigXsecOverSM]) determining minimum likelihood for current configurations w.r.t all observable
RooAbsTestStatistic::initSimMode: creating slave calculator #0 for state channel1 (2 dataset entries)
[#1] INFO:Fitting -- RooAbsTestStatistic::initSimMode: created 1 slave calculators.
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_simPdf_obsData_with_constr_Profile[SigXsecOverSM]) minimum found at (SigXsecOverSM=1.11597)
.
Will use these parameter points to generate pseudo data for bkg only
  1) 0x99c76e0 RooRealVar::               alpha_syst2 = 0.711172 +/- 0.914106  L(-5 - 5)  "alpha_syst2"
  2) 0x112b64f0 RooRealVar::               alpha_syst3 = 0.26146 +/- 0.9291  L(-5 - 5)  "alpha_syst3"
  3) 0x112d3410 RooRealVar:: gamma_stat_channel1_bin_0 = 1.03677 +/- 0.0462899  L(0 - 1.25)  "gamma_stat_channel1_bin_0"
  4) 0x11342a20 RooRealVar:: gamma_stat_channel1_bin_1 = 1.05319 +/- 0.0761205  L(0 - 1.5)  "gamma_stat_channel1_bin_1"
  5) 0x11386ed0 RooRealVar::             SigXsecOverSM = 0 +/- 0  L(0 - 3) B(20)  "SigXsecOverSM"
Info in <TCanvas::Print>: pdf file two-sided_upper_limit_output.pdf has been created
-2 sigma  band 0
-1 sigma  band 0.495 [Power Constraint)]
median of band 1.095
+1 sigma  band 1.545
+2 sigma  band 1.995

observed 95% upper-limit 2.325
CLb strict [P(toy>obs|0)] for observed 95% upper-limit 0.96
CLb inclusive [P(toy>=obs|0)] for observed 95% upper-limit 0.96