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 the following for ModelConfigB_only === 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.190787 Snapshot: 1) 0x79a3280 RooRealVar:: SigXsecOverSM = 0 +/- 0 L(0 - 3) "SigXsecOverSM" === 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.190787 Snapshot: 1) 0x79a3890 RooRealVar:: SigXsecOverSM = 1 +/- 0 L(0 - 3) "SigXsecOverSM" [#0] PROGRESS:Generation -- Test Statistic on data: 1.77388 [#1] INFO:InputArguments -- Profiling conditional MLEs for Null. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Null. [#0] PROGRESS:Generation -- generated toys: 500 / 5000 [#0] PROGRESS:Generation -- generated toys: 1000 / 5000 [#0] PROGRESS:Generation -- generated toys: 1500 / 5000 [#0] PROGRESS:Generation -- generated toys: 2000 / 5000 [#0] PROGRESS:Generation -- generated toys: 2500 / 5000 [#0] PROGRESS:Generation -- generated toys: 3000 / 5000 [#0] PROGRESS:Generation -- generated toys: 3500 / 5000 [#0] PROGRESS:Generation -- generated toys: 4000 / 5000 [#0] PROGRESS:Generation -- generated toys: 4500 / 5000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- generated toys: 500 / 1250 [#0] PROGRESS:Generation -- generated toys: 1000 / 1250 Results HypoTestCalculator_result: - Null p-value = 0.0306 +/- 0.00243572 - Significance = 1.87205 +/- 0.0352146 sigma - Number of Alt toys: 1250 - Number of Null toys: 5000 - Test statistic evaluated on data: 1.77388 - CL_b: 0.0306 +/- 0.00243572 - CL_s+b: 0.4312 +/- 0.0140076 - CL_s: 14.0915 +/- 1.21148 Expected p -value and significance at -2 sigma = 0.7916 significance -0.811985 sigma Expected p -value and significance at -1 sigma = 0.2472 significance 0.683327 sigma Expected p -value and significance at 0 sigma = 0.0434 significance 1.71252 sigma Expected p -value and significance at 1 sigma = 0.0034 significance 2.70648 sigma Expected p -value and significance at 2 sigma = 0.0004 significance 3.35279 sigma