0x6a88000 results/example_combined_GaussExample_model.root Running HypoTestInverter on the workspace combined RooWorkspace(combined) combined contents variables --------- (Lumi,SigXsecOverSM,alpha_syst1,alpha_syst2,alpha_syst3,channelCat,gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1,nom_alpha_syst1,nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1,nominalLumi,obs_x_channel1) p.d.f.s ------- RooGaussian::alpha_syst1Constraint[ x=alpha_syst1 mean=nom_alpha_syst1 sigma=1 ] = 1 RooGaussian::alpha_syst2Constraint[ x=alpha_syst2 mean=nom_alpha_syst2 sigma=1 ] = 1 RooGaussian::alpha_syst3Constraint[ x=alpha_syst3 mean=nom_alpha_syst3 sigma=1 ] = 1 RooRealSumPdf::channel1_model[ signal_channel1_scaleFactors * signal_channel1_shapes + background1_channel1_scaleFactors * background1_channel1_shapes + background2_channel1_scaleFactors * background2_channel1_shapes ] = 240 RooPoisson::gamma_stat_channel1_bin_0_constraint[ x=nom_gamma_stat_channel1_bin_0 mean=gamma_stat_channel1_bin_0_poisMean ] = 0.019943 RooPoisson::gamma_stat_channel1_bin_1_constraint[ x=nom_gamma_stat_channel1_bin_1 mean=gamma_stat_channel1_bin_1_poisMean ] = 0.039861 RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.1 ] = 1 RooProdPdf::model_channel1[ lumiConstraint * alpha_syst1Constraint * alpha_syst2Constraint * alpha_syst3Constraint * gamma_stat_channel1_bin_0_constraint * gamma_stat_channel1_bin_1_constraint * channel1_model(obs_x_channel1) ] = 0.190787 RooSimultaneous::simPdf[ indexCat=channelCat channel1=model_channel1 ] = 0.190787 functions -------- RooHistFunc::background1_channel1_Hist_alphanominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 100 RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon[ paramList=(alpha_syst2) ] = 1 RooProduct::background1_channel1_scaleFactors[ background1_channel1_epsilon * Lumi ] = 1 RooProduct::background1_channel1_shapes[ background1_channel1_Hist_alphanominal * mc_stat_channel1 * channel1_model_binWidth ] = 200 RooHistFunc::background2_channel1_Hist_alphanominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 0 RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon[ paramList=(alpha_syst3) ] = 1 RooProduct::background2_channel1_scaleFactors[ background2_channel1_epsilon * Lumi ] = 1 RooProduct::background2_channel1_shapes[ background2_channel1_Hist_alphanominal * mc_stat_channel1 * channel1_model_binWidth ] = 0 RooBinWidthFunction::channel1_model_binWidth[ HistFuncForBinWidth=signal_channel1_Hist_alphanominal HistFuncForBinWidth=signal_channel1_Hist_alphanominal ] = 2 RooProduct::gamma_stat_channel1_bin_0_poisMean[ gamma_stat_channel1_bin_0 * gamma_stat_channel1_bin_0_tau ] = 400 RooProduct::gamma_stat_channel1_bin_1_poisMean[ gamma_stat_channel1_bin_1 * gamma_stat_channel1_bin_1_tau ] = 100 ParamHistFunc::mc_stat_channel1[ ] = 1 RooHistFunc::signal_channel1_Hist_alphanominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 20 RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon[ paramList=(alpha_syst1) ] = 1 RooProduct::signal_channel1_scaleFactors[ signal_channel1_epsilon * SigXsecOverSM * Lumi ] = 1 RooProduct::signal_channel1_shapes[ signal_channel1_Hist_alphanominal * channel1_model_binWidth ] = 40 datasets -------- RooDataSet::obsData(obs_x_channel1,channelCat) RooDataSet::asimovData(obs_x_channel1,channelCat) embedded datasets (in pdfs and functions) ----------------------------------------- RooDataHist::signal_channel1_Hist_alphanominalDHist(obs_x_channel1) RooDataHist::background1_channel1_Hist_alphanominalDHist(obs_x_channel1) RooDataHist::background2_channel1_Hist_alphanominalDHist(obs_x_channel1) parameter snapshots ------------------- NominalParamValues = (nominalLumi=1[C],nom_alpha_syst1=0[C],nom_alpha_syst2=0[C],nom_alpha_syst3=0[C],nom_gamma_stat_channel1_bin_0=400[C],nom_gamma_stat_channel1_bin_1=100[C],obs_x_channel1=1.25,Lumi=1 +/- 0.1[C],alpha_syst1=0 +/- 1[C],alpha_syst2=0 +/- 1,alpha_syst3=0 +/- 1,gamma_stat_channel1_bin_0=1 +/- 0.05,gamma_stat_channel1_bin_1=1 +/- 0.1,SigXsecOverSM=1 +/- 0) named sets ---------- ModelConfig_GlobalObservables:(nominalLumi,nom_alpha_syst1,nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1) ModelConfig_NuisParams:(alpha_syst2,alpha_syst3,gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1) ModelConfig_Observables:(obs_x_channel1,channelCat) ModelConfig_POI:(SigXsecOverSM) globalObservables:(nominalLumi,nom_alpha_syst1,nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1) observables:(obs_x_channel1,channelCat) generic objects --------------- RooStats::ModelConfig::ModelConfig Using data set obsData StandardHypoTestInvDemo : POI initial value: SigXsecOverSM = 1 [#1] INFO:InputArguments -- HypoTestInverter ---- Input models: using as S+B (null) model : ModelConfig using as B (alternate) model : ModelConfig_with_poi_0 Doing a fixed scan in interval : 0 , 5 [#1] INFO:Eval -- HypoTestInverter::GetInterval - run a fixed scan [#0] WARNING:InputArguments -- HypoTestInverter::RunFixedScan - xMax > upper bound, using xmax = 3 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 0 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.158989 Snapshot: 1) 0x79a1ec0 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.158989 Snapshot: 1) 0x6aa4950 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: 0 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 0 CLs = 1 +/- 0 CLb = 1 +/- 0 CLsplusb = 1 +/- 0 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 0.6 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.178068 Snapshot: 1) 0x7ab93a0 RooRealVar:: SigXsecOverSM = 0.6 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.178068 Snapshot: 1) 0x79a2ae0 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: -2.35222 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 0.6 CLs = 0.826039 +/- 0.0187037 CLb = 0.914 +/- 0.0125383 CLsplusb = 0.755 +/- 0.0136006 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 1.2 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.197147 Snapshot: 1) 0x7b1a490 RooRealVar:: SigXsecOverSM = 1.2 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.197147 Snapshot: 1) 0x7b0b400 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: -2.9364 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 1.2 CLs = 0.495652 +/- 0.018325 CLb = 0.92 +/- 0.0121326 CLsplusb = 0.456 +/- 0.01575 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 1.8 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.216225 Snapshot: 1) 0x7ab4050 RooRealVar:: SigXsecOverSM = 1.8 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.216225 Snapshot: 1) 0x62fe000 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: -2.05075 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 1.8 CLs = 0.209052 +/- 0.0137241 CLb = 0.928 +/- 0.0115599 CLsplusb = 0.194 +/- 0.0125046 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 2.4 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.235304 Snapshot: 1) 0x7b9da70 RooRealVar:: SigXsecOverSM = 2.4 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.235304 Snapshot: 1) 0x7b26150 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: 0.0783908 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 2.4 CLs = 0.0503282 +/- 0.00728062 CLb = 0.914 +/- 0.0125383 CLsplusb = 0.046 +/- 0.0066245 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnapshot(combined) replacing previous snapshot with name ModelConfig__snapshot [#0] PROGRESS:Eval -- Running for SigXsecOverSM = 3 === Using the following for ModelConfig === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.254383 Snapshot: 1) 0x7bb24d0 RooRealVar:: SigXsecOverSM = 3 +/- 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === Observables: RooArgSet:: = (obs_x_channel1,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:: = (nominalLumi,nom_alpha_syst1,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.254383 Snapshot: 1) 0x7ade650 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: 3.27476 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 3 CLs = 0.00534188 +/- 0.0023838 CLb = 0.936 +/- 0.0109457 CLsplusb = 0.005 +/- 0.00223047 Time to perform limit scan Real time 0:00:03, CP time 3.990 The computed upper limit is: 2.40438 +/- 0.0566658 Expected upper limits, using the B (alternate) model : expected limit (median) 1.61927 expected limit (-1 sig) 1.09318 expected limit (+1 sig) 2.24198 expected limit (-2 sig) 0.787047 expected limit (+2 sig) 2.86153 [#0] WARNING:Plotting -- Could not determine xmin and xmax of sampling distribution that was given to plot. [#0] WARNING:Plotting -- Could not determine xmin and xmax of sampling distribution that was given to plot.