lxplus009:/tmp/ganis/roofitcrash $ root -b ------------------------------------------------------------ | Welcome to ROOT 6.04/14 http://root.cern.ch | | (c) 1995-2014, The ROOT Team | | Built for linuxx8664gcc | | From tag v6-04-14, 3 February 2016 | | Try '.help', '.demo', '.license', '.credits', '.quit'/'.q' | ------------------------------------------------------------ root [0] .L tutorials/roostats/StandardHypoTestInvDemo.C+ Info in : creating shared library /tmp/ganis/roofitcrash/./tutorials/roostats/StandardHypoTestInvDemo_C.so RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University All rights reserved, please read http://roofit.sourceforge.net/license.txt root [1] StandardHypoTestInvDemo() 0x2d20cc0 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,binWidth_obs_x_channel1_0,binWidth_obs_x_channel1_1,binWidth_obs_x_channel1_2,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,weightVar) 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[ binWidth_obs_x_channel1_0 * L_x_signal_channel1_overallSyst_x_Exp + binWidth_obs_x_channel1_1 * L_x_background1_channel1_overallSyst_x_StatUncert + binWidth_obs_x_channel1_2 * L_x_background2_channel1_overallSyst_x_StatUncert ] = 220 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.174888 RooSimultaneous::simPdf[ indexCat=channelCat channel1=model_channel1 ] = 0.174888 functions -------- RooProduct::L_x_background1_channel1_overallSyst_x_StatUncert[ Lumi * background1_channel1_overallSyst_x_StatUncert ] = 0 RooProduct::L_x_background2_channel1_overallSyst_x_StatUncert[ Lumi * background2_channel1_overallSyst_x_StatUncert ] = 100 RooProduct::L_x_signal_channel1_overallSyst_x_Exp[ Lumi * signal_channel1_overallSyst_x_Exp ] = 10 RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon[ paramList=(alpha_syst2) ] = 1 RooHistFunc::background1_channel1_nominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 0 RooProduct::background1_channel1_overallSyst_x_Exp[ background1_channel1_nominal * background1_channel1_epsilon ] = 0 RooProduct::background1_channel1_overallSyst_x_StatUncert[ mc_stat_channel1 * background1_channel1_overallSyst_x_Exp ] = 0 RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon[ paramList=(alpha_syst3) ] = 1 RooHistFunc::background2_channel1_nominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 100 RooProduct::background2_channel1_overallSyst_x_Exp[ background2_channel1_nominal * background2_channel1_epsilon ] = 100 RooProduct::background2_channel1_overallSyst_x_StatUncert[ mc_stat_channel1 * background2_channel1_overallSyst_x_Exp ] = 100 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 RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon[ paramList=(alpha_syst1) ] = 1 RooHistFunc::signal_channel1_nominal[ depList=(obs_x_channel1) depList=(obs_x_channel1) ] = 10 RooProduct::signal_channel1_overallNorm_x_sigma_epsilon[ SigXsecOverSM * signal_channel1_epsilon ] = 1 RooProduct::signal_channel1_overallSyst_x_Exp[ signal_channel1_nominal * signal_channel1_overallNorm_x_sigma_epsilon ] = 10 datasets -------- RooDataSet::asimovData(obs_x_channel1,weightVar,channelCat) RooDataSet::obsData(channelCat,obs_x_channel1) embedded datasets (in pdfs and functions) ----------------------------------------- RooDataHist::signal_channel1nominalDHist(obs_x_channel1) RooDataHist::background1_channel1nominalDHist(obs_x_channel1) RooDataHist::background2_channel1nominalDHist(obs_x_channel1) parameter snapshots ------------------- NominalParamValues = (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],weightVar=0,obs_x_channel1=1.75,Lumi=1[C],nominalLumi=1[C],alpha_syst1=0[C],nom_alpha_syst1=0[C],alpha_syst2=0,alpha_syst3=0,gamma_stat_channel1_bin_0=1,gamma_stat_channel1_bin_1=1,SigXsecOverSM=1,binWidth_obs_x_channel1_0=2[C],binWidth_obs_x_channel1_1=2[C],binWidth_obs_x_channel1_2=2[C]) named sets ---------- ModelConfig_GlobalObservables:(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,weightVar,channelCat) ModelConfig_POI:(SigXsecOverSM) globalObservables:(nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1) observables:(obs_x_channel1,weightVar,channelCat) generic objects --------------- RooStats::ModelConfig::ModelConfig Using data set obsData Info in : Model ModelConfig has no snapshot - make one using model poi StandardHypoTestInvDemo - Switch off all systematics by setting them constant to their initial values Info in : The background model does not exist Info in : Copy it from ModelConfig ModelConfig and set POI to zero StandardHypoTestInvDemo : POI initial value: SigXsecOverSM = 1 Info in : Using Minuit as minimizer for computing the test statistic Info in : Using LEP test statistic - an initial fit is not done and the TS will use the nuisances at the model value Warning in : generate binned is activated but the number of ovservable is 3. Too much memory could be needed for allocating all the bins [#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, use xmax = 3 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.158989 Snapshot: 1) 0x4424d10 RooRealVar:: SigXsecOverSM = 0 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.158989 Snapshot: 1) 0x4436f80 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session +++ Starting PROOF-Lite with 8 workers +++ Opening connections to workers: OK (8 workers) Setting up worker servers: OK (8 workers) PROOF set to parallel mode (8 workers) [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 1 Info in : nwrks: 8 [TProof::Progress] Total 8 events. | (3 |====================| 100.00 % [6.1 evts/s, time left: 0.0 s] Query processing time: 1.3 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 2 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (3 |====================| 100.00 % [7.9 evts/s, time left: 0.0 s] Query processing time: 1.0 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 0 CLs = 1 +/- 0 CLb = 1 +/- 0 CLsplusb = 1 +/- 0 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.168529 Snapshot: 1) 0x443e4d0 RooRealVar:: SigXsecOverSM = 0.6 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.168529 Snapshot: 1) 0x443e4d0 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 3 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [4.3 evts/s, time left: 0.0 s] Query processing time: 1.9 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 4 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (4 |====================| 100.00 % [11.1 evts/s, time left: 0.0 s] Query processing time: 0.7 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 0.6 CLs = 0.875422 +/- 0.0117442 CLb = 0.988095 +/- 0.00483108 CLsplusb = 0.865 +/- 0.0108062 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.178068 Snapshot: 1) 0x4f62300 RooRealVar:: SigXsecOverSM = 1.2 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.178068 Snapshot: 1) 0x4f62300 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 5 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [2.8 evts/s, time left: 0.0 s] Query processing time: 2.9 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 6 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [5.6 evts/s, time left: 0.0 s] Query processing time: 1.4 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 1.2 CLs = 0.45441 +/- 0.0160727 CLb = 0.988095 +/- 0.00483108 CLsplusb = 0.449 +/- 0.0157289 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.187607 Snapshot: 1) 0x54e7950 RooRealVar:: SigXsecOverSM = 1.8 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.187607 Snapshot: 1) 0x54e7950 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 7 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [3.8 evts/s, time left: 0.0 s] Query processing time: 2.1 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 8 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [7.1 evts/s, time left: 0.0 s] Query processing time: 1.1 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 1.8 CLs = 0.105253 +/- 0.00978303 CLb = 0.988095 +/- 0.00483108 CLsplusb = 0.104 +/- 0.00965319 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.197147 Snapshot: 1) 0x5bb7ef0 RooRealVar:: SigXsecOverSM = 2.4 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.197147 Snapshot: 1) 0x5bb7ef0 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 9 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (2 |====================| 100.00 % [4.9 evts/s, time left: 0.0 s] Query processing time: 1.6 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 10 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (3 |====================| 100.00 % [8.2 evts/s, time left: 0.0 s] Query processing time: 1.0 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 2.4 CLs = 0.00708434 +/- 0.00266846 CLb = 0.988095 +/- 0.00483108 CLsplusb = 0.007 +/- 0.00263647 [#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(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,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.206686 Snapshot: 1) 0x5b87900 RooRealVar:: SigXsecOverSM = 3 L(0 - 3) "SigXsecOverSM" === Using the following for ModelConfig_with_poi_0 === 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.206686 Snapshot: 1) 0x5b87900 RooRealVar:: SigXsecOverSM = 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 -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 11 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (3 |====================| 100.00 % [3.0 evts/s, time left: 0.0 s] Query processing time: 2.6 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 1000 [#1] INFO:InputArguments -- Profiling conditional MLEs for Alt. [#1] INFO:InputArguments -- Using a ToyMCSampler. Now configuring for Alt. [#0] PROGRESS:Generation -- RooStudyManager::runProof() opening PROOF session [#0] PROGRESS:Generation -- RooStudyManager::runProof() sending work package to PROOF servers [#0] PROGRESS:Generation -- RooStudyManager::runProof() starting PROOF processing of 8 experiments Info in : starting query: 12 Info in : nwrks: 8 [TProof::Progress] Total 8 events. / (4 |====================| 100.00 % [6.3 evts/s, time left: 0.0 s] Query processing time: 1.3 s Lite-0: all output objects have been merged [#0] PROGRESS:Generation -- RooStudyManager::runProof() aggregating results data [#0] PROGRESS:Generation -- RooStudyManager::runProof() cleaning up input list [#0] PROGRESS:Generation -- Merged data from nworkers # 8- merged data size is 504 [#0] PROGRESS:Eval -- P values for SigXsecOverSM = 3 CLs = 0 +/- 0 CLb = 0.988095 +/- 0.00483108 CLsplusb = 0 +/- 0 Time to perform limit scan Real time 0:00:28, CP time 1.360 The computed upper limit is: 2.1377 +/- 0.0529227 Expected upper limits, using the B (alternate) model : expected limit (median) 1.13336 expected limit (-1 sig) 0.898732 expected limit (+1 sig) 1.49524 expected limit (-2 sig) 0.59299 expected limit (+2 sig) 1.90383 Info in : detailed output will be written in output result file Info in : HypoTestInverterResult has been written in the file Freq_CLs_grid_ts0_ Warning in : Replacing existing TH1: CLs_observed (Potential memory leak). [#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. Warning in : Replacing existing TH1: CLb_observed (Potential memory leak). Warning in : Replacing existing TH1: CLs+b_observed (Potential memory leak). root [2]