  *******************************************
  *                                         *
  *        W E L C O M E  to  R O O T       *
  *                                         *
  *   Version   5.32/03        9 May 2012   *
  *                                         *
  *  You are welcome to visit our Web site  *
  *          http://root.cern.ch            *
  *                                         *
  *******************************************

ROOT 5.32/03 (tags/v5-32-03@44194, May 09 2012, 13:00:14 on linuxx8664gcc)

CINT/ROOT C/C++ Interpreter version 5.18.00, July 2, 2010
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.

Processing runSig_Factory.C+("/phys/groups/tev/scratch3/users/rompotis/CombinationResults/resultsV13updated/mumu_mA170tb35_comb_combined_AllSyst_model.root","mA170tb35")...

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

  1) 0x7f3b8f0 RooRealVar::                       Lumi = 1  L(0 - 10)  "Lumi"
  2) 0x8046f30 RooRealVar::                alpha_B_EFF = 0  L(-5 - 5)  "alpha_B_EFF"
  3) 0x84442c0 RooRealVar::          alpha_BkgFit_bjet = 0  L(-5 - 5)  "alpha_BkgFit_bjet"
  4) 0x83453a0 RooRealVar::        alpha_BkgFit_nobjet = 0  L(-5 - 5)  "alpha_BkgFit_nobjet"
  5) 0x804eaf0 RooRealVar::                alpha_C_EFF = 0  L(-5 - 5)  "alpha_C_EFF"
  6) 0x8059900 RooRealVar::                  alpha_JES = 0  L(-5 - 5)  "alpha_JES"
  7) 0x805aef0 RooRealVar::                alpha_L_EFF = 0  L(-5 - 5)  "alpha_L_EFF"
  8) 0x8070400 RooRealVar::                alpha_M_EFF = 0  L(-5 - 5)  "alpha_M_EFF"
  9) 0x84d4350 RooRealVar::   alpha_SignalFit_bbA_bjet = 0  L(-5 - 5)  "alpha_SignalFit_bbA_bjet"
 10) 0x81ba390 RooRealVar:: alpha_SignalFit_bbA_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_bbA_nobjet"
 11) 0x80ca7c0 RooRealVar::   alpha_SignalFit_bbH_bjet = 0  L(-5 - 5)  "alpha_SignalFit_bbH_bjet"
 12) 0x8233840 RooRealVar:: alpha_SignalFit_bbH_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_bbH_nobjet"
 13) 0x8040dd0 RooRealVar::   alpha_SignalFit_bbh_bjet = 0  L(-5 - 5)  "alpha_SignalFit_bbh_bjet"
 14) 0x81c7ec0 RooRealVar:: alpha_SignalFit_bbh_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_bbh_nobjet"
 15) 0x8077320 RooRealVar::   alpha_SignalFit_ggA_bjet = 0  L(-5 - 5)  "alpha_SignalFit_ggA_bjet"
 16) 0x81970c0 RooRealVar:: alpha_SignalFit_ggA_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_ggA_nobjet"
 17) 0x814e840 RooRealVar::   alpha_SignalFit_ggH_bjet = 0  L(-5 - 5)  "alpha_SignalFit_ggH_bjet"
 18) 0x82e7a00 RooRealVar:: alpha_SignalFit_ggH_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_ggH_nobjet"
 19) 0x8140000 RooRealVar::   alpha_SignalFit_ggh_bjet = 0  L(-5 - 5)  "alpha_SignalFit_ggh_bjet"
 20) 0x8293f30 RooRealVar:: alpha_SignalFit_ggh_nobjet = 0  L(-5 - 5)  "alpha_SignalFit_ggh_nobjet"
 21) 0x81c0eb0 RooRealVar::               alpha_XS_bbA = 0  L(-5 - 5)  "alpha_XS_bbA"
 22) 0x8147a10 RooRealVar::               alpha_XS_ggA = 0  L(-5 - 5)  "alpha_XS_ggA"
Setting max function calls
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:Minization --  Including the following contraint terms in minimization: (lumiConstraint,alpha_B_EFFConstraint,alpha_C_EFFConstraint,alpha_JESConstraint,alpha_L_EFFConstraint,alpha_M_EFFConstraint,alpha_SignalFit_bbA_bjetConstraint,alpha_XS_bbAConstraint,alpha_SignalFit_bbh_bjetConstraint,alpha_SignalFit_bbH_bjetConstraint,alpha_SignalFit_ggA_bjetConstraint,alpha_XS_ggAConstraint,alpha_SignalFit_ggh_bjetConstraint,alpha_SignalFit_ggH_bjetConstraint,alpha_BkgFit_bjetConstraint,alpha_SignalFit_bbA_nobjetConstraint,alpha_SignalFit_bbh_nobjetConstraint,alpha_SignalFit_bbH_nobjetConstraint,alpha_SignalFit_ggA_nobjetConstraint,alpha_SignalFit_ggh_nobjetConstraint,alpha_SignalFit_ggH_nobjetConstraint,alpha_BkgFit_nobjetConstraint)
DEBUG: making asimovData_0
Creating asimov data at mu = 0, profiling at mu = 0
[#0] ERROR:ObjectHandling -- RooWorkspace::loadSnapshot(combined) no snapshot with name nominalGlobs is available
nominalGlobs doesn't exist. Saving snapshot.
[#0] ERROR:ObjectHandling -- RooWorkspace::loadSnapshot(combined) no snapshot with name nominalNuis is available
nominalNuis doesn't exist. Saving snapshot.
makeAsimovData: run minimize()
makeAsimovData: run minimize() and got status 0
Saving conditional snapshots
Glob snapshot name = conditionalGlobs_0
Nuis snapshot name = conditionalNuis_0
[#1] INFO:NumericIntegration -- RooRealIntegral::init(MM_nobjet_model_Int[obs_x_MM_nobjet]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_MM_nobjet)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(MM_bjet_model_Int[obs_x_MM_bjet]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_MM_bjet)
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovData_0
DEBUG: making asimovData_0 DONE
DEBUG:-STA0 printing out conditionalNuis_0 set of parameters
  1) 0x7f3b8f0 RooRealVar::                       Lumi = 0.999961 +/- 0.0389982  L(0 - 10)  "Lumi"
  2) 0x8046f30 RooRealVar::                alpha_B_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_B_EFF"
  3) 0x84442c0 RooRealVar::          alpha_BkgFit_bjet = 0.386581 +/- 0.581017  L(-5 - 5)  "alpha_BkgFit_bjet"
  4) 0x83453a0 RooRealVar::        alpha_BkgFit_nobjet = 0.0447474 +/- 0.305682  L(-5 - 5)  "alpha_BkgFit_nobjet"
  5) 0x804eaf0 RooRealVar::                alpha_C_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_C_EFF"
  6) 0x8059900 RooRealVar::                  alpha_JES = 0 +/- 0.99374  L(-5 - 5)  "alpha_JES"
  7) 0x805aef0 RooRealVar::                alpha_L_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_L_EFF"
  8) 0x8070400 RooRealVar::                alpha_M_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_M_EFF"
  9) 0x84d4350 RooRealVar::   alpha_SignalFit_bbA_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbA_bjet"
 10) 0x81ba390 RooRealVar:: alpha_SignalFit_bbA_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbA_nobjet"
 11) 0x80ca7c0 RooRealVar::   alpha_SignalFit_bbH_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbH_bjet"
 12) 0x8233840 RooRealVar:: alpha_SignalFit_bbH_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbH_nobjet"
 13) 0x8040dd0 RooRealVar::   alpha_SignalFit_bbh_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbh_bjet"
 14) 0x81c7ec0 RooRealVar:: alpha_SignalFit_bbh_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbh_nobjet"
 15) 0x8077320 RooRealVar::   alpha_SignalFit_ggA_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggA_bjet"
 16) 0x81970c0 RooRealVar:: alpha_SignalFit_ggA_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggA_nobjet"
 17) 0x814e840 RooRealVar::   alpha_SignalFit_ggH_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggH_bjet"
 18) 0x82e7a00 RooRealVar:: alpha_SignalFit_ggH_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggH_nobjet"
 19) 0x8140000 RooRealVar::   alpha_SignalFit_ggh_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggh_bjet"
 20) 0x8293f30 RooRealVar:: alpha_SignalFit_ggh_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggh_nobjet"
 21) 0x81c0eb0 RooRealVar::               alpha_XS_bbA = 0 +/- 0.99374  L(-5 - 5)  "alpha_XS_bbA"
 22) 0x8147a10 RooRealVar::               alpha_XS_ggA = 0 +/- 0.99374  L(-5 - 5)  "alpha_XS_ggA"
DEBUG:-END0 printing out conditionalNuis_0 set of parameters
Asimov data doesn't exist! Please, allow me to build one for you...
DEBUG: making asimovData1
Creating asimov data at mu = 1, profiling at mu = 1
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpGlobs
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpNuis
makeAsimovData: run minimize()
makeAsimovData: run minimize() and got status 0
Saving conditional snapshots
Glob snapshot name = conditionalGlobs_1
Nuis snapshot name = conditionalNuis_1
[#0] WARNING:InputArguments -- RooWorkspace::defineSet(combined) WARNING redefining previously defined named set obsAndWeight
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovData_1
DEBUG: making asimovData1 DONE
DEBUG: making asimovData_a
Creating asimov data at mu = 0.1, profiling at mu = 0.1
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpGlobs
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpNuis
makeAsimovData: run minimize()
makeAsimovData: run minimize() and got status 0
Saving conditional snapshots
Glob snapshot name = conditionalGlobs_0.1
Nuis snapshot name = conditionalNuis_0.1
[#0] WARNING:InputArguments -- RooWorkspace::defineSet(combined) WARNING redefining previously defined named set obsAndWeight
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovData_0.1
DEBUG: making asimovData_a DONE
DEBUG: making asimovData_b
Creating asimov data at mu = 0.5, profiling at mu = 0.5
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpGlobs
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpNuis
makeAsimovData: run minimize()
makeAsimovData: run minimize() and got status 0
Saving conditional snapshots
Glob snapshot name = conditionalGlobs_0.5
Nuis snapshot name = conditionalNuis_0.5
[#0] WARNING:InputArguments -- RooWorkspace::defineSet(combined) WARNING redefining previously defined named set obsAndWeight
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovData_0.5
DEBUG: making asimovData_b DONE
DEBUG: making asimovData_c
Creating asimov data at mu = 0.8, profiling at mu = 0.8
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpGlobs
[#1] INFO:ObjectHandling -- RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name tmpNuis
makeAsimovData: run minimize()
makeAsimovData: run minimize() and got status 0
Saving conditional snapshots
Glob snapshot name = conditionalGlobs_0.8
Nuis snapshot name = conditionalNuis_0.8
[#0] WARNING:InputArguments -- RooWorkspace::defineSet(combined) WARNING redefining previously defined named set obsAndWeight
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovData_0.8
DEBUG: making asimovData_c DONE
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:Minization --  Including the following contraint terms in minimization: (lumiConstraint,alpha_B_EFFConstraint,alpha_C_EFFConstraint,alpha_JESConstraint,alpha_L_EFFConstraint,alpha_M_EFFConstraint,alpha_SignalFit_bbA_bjetConstraint,alpha_XS_bbAConstraint,alpha_SignalFit_bbh_bjetConstraint,alpha_SignalFit_bbH_bjetConstraint,alpha_SignalFit_ggA_bjetConstraint,alpha_XS_ggAConstraint,alpha_SignalFit_ggh_bjetConstraint,alpha_SignalFit_ggH_bjetConstraint,alpha_BkgFit_bjetConstraint,alpha_SignalFit_bbA_nobjetConstraint,alpha_SignalFit_bbh_nobjetConstraint,alpha_SignalFit_bbH_nobjetConstraint,alpha_SignalFit_ggA_nobjetConstraint,alpha_SignalFit_ggh_nobjetConstraint,alpha_SignalFit_ggH_nobjetConstraint,alpha_BkgFit_nobjetConstraint)
  1) 0x85256e0 RooRealVar::                    nominalLumi = 0.98829 C  L(0 - 1.39)  "nominalLumi"
  2) 0x84dbd70 RooRealVar::                nom_alpha_B_EFF = 0.0497461 C  L(-10 - 10)  "nom_alpha_B_EFF"
  3) 0x8051c50 RooRealVar::                nom_alpha_C_EFF = 0.015235 C  L(-10 - 10)  "nom_alpha_C_EFF"
  4) 0x84c86d0 RooRealVar::                  nom_alpha_JES = -0.127584 C  L(-10 - 10)  "nom_alpha_JES"
  5) 0x805e0b0 RooRealVar::                nom_alpha_L_EFF = 0.0222461 C  L(-10 - 10)  "nom_alpha_L_EFF"
  6) 0x80735c0 RooRealVar::                nom_alpha_M_EFF = -0.0731132 C  L(-10 - 10)  "nom_alpha_M_EFF"
  7) 0x81bd400 RooRealVar:: nom_alpha_SignalFit_bbA_nobjet = -0.242671 C  L(-10 - 10)  "nom_alpha_SignalFit_bbA_nobjet"
  8) 0x81c4070 RooRealVar::               nom_alpha_XS_bbA = -0.704408 C  L(-10 - 10)  "nom_alpha_XS_bbA"
  9) 0x81cb0e0 RooRealVar:: nom_alpha_SignalFit_bbh_nobjet = -0.000369848 C  L(-10 - 10)  "nom_alpha_SignalFit_bbh_nobjet"
 10) 0x82368f0 RooRealVar:: nom_alpha_SignalFit_bbH_nobjet = -0.244759 C  L(-10 - 10)  "nom_alpha_SignalFit_bbH_nobjet"
 11) 0x819a2e0 RooRealVar:: nom_alpha_SignalFit_ggA_nobjet = -0.104354 C  L(-10 - 10)  "nom_alpha_SignalFit_ggA_nobjet"
 12) 0x814abd0 RooRealVar::               nom_alpha_XS_ggA = -0.373992 C  L(-10 - 10)  "nom_alpha_XS_ggA"
 13) 0x8297010 RooRealVar:: nom_alpha_SignalFit_ggh_nobjet = -0.00565278 C  L(-10 - 10)  "nom_alpha_SignalFit_ggh_nobjet"
 14) 0x82eaae0 RooRealVar:: nom_alpha_SignalFit_ggH_nobjet = -0.105088 C  L(-10 - 10)  "nom_alpha_SignalFit_ggH_nobjet"
 15) 0x8348480 RooRealVar::        nom_alpha_BkgFit_nobjet = -0.245095 C  L(-10 - 10)  "nom_alpha_BkgFit_nobjet"
 16) 0x84d7570 RooRealVar::   nom_alpha_SignalFit_bbA_bjet = 0.0274456 C  L(-10 - 10)  "nom_alpha_SignalFit_bbA_bjet"
 17) 0x84dfa10 RooRealVar::   nom_alpha_SignalFit_bbh_bjet = -5.82153e-06 C  L(-10 - 10)  "nom_alpha_SignalFit_bbh_bjet"
 18) 0x80cd9e0 RooRealVar::   nom_alpha_SignalFit_bbH_bjet = 0.0265662 C  L(-10 - 10)  "nom_alpha_SignalFit_bbH_bjet"
 19) 0x807a540 RooRealVar::   nom_alpha_SignalFit_ggA_bjet = 0.00166464 C  L(-10 - 10)  "nom_alpha_SignalFit_ggA_bjet"
 20) 0x8143220 RooRealVar::   nom_alpha_SignalFit_ggh_bjet = 0.00118857 C  L(-10 - 10)  "nom_alpha_SignalFit_ggh_bjet"
 21) 0x8151a60 RooRealVar::   nom_alpha_SignalFit_ggH_bjet = 0.00163554 C  L(-10 - 10)  "nom_alpha_SignalFit_ggH_bjet"
 22) 0x84474e0 RooRealVar::          nom_alpha_BkgFit_bjet = 0.176819 C  L(-10 - 10)  "nom_alpha_BkgFit_bjet"
Trying to minimize asimov lnL(mu=0,theta^^)
Using conditionalNuis_0 Success with status: 0
Trying to minimize asimov lnL(mu=1,theta^^)
Using conditionalNuis_1 Success with status: 0
DEBUG (asimov) lnL(mu=0,theta^^) =  -124738.1268 lnL(mu=1,theta^^) =  -124742.5154 q0 =   8.777165
Trying to minimize observed lnL(mu=0,theta^^)
DEBUG:-STA printing out conditionalNuis_0 set of parameters
  1) 0x7f3b8f0 RooRealVar::                       Lumi = 0.999961 +/- 0.0389982  L(0 - 10)  "Lumi"
  2) 0x8046f30 RooRealVar::                alpha_B_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_B_EFF"
  3) 0x84442c0 RooRealVar::          alpha_BkgFit_bjet = 0.386581 +/- 0.581017  L(-5 - 5)  "alpha_BkgFit_bjet"
  4) 0x83453a0 RooRealVar::        alpha_BkgFit_nobjet = 0.0447474 +/- 0.305682  L(-5 - 5)  "alpha_BkgFit_nobjet"
  5) 0x804eaf0 RooRealVar::                alpha_C_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_C_EFF"
  6) 0x8059900 RooRealVar::                  alpha_JES = 0 +/- 0.99374  L(-5 - 5)  "alpha_JES"
  7) 0x805aef0 RooRealVar::                alpha_L_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_L_EFF"
  8) 0x8070400 RooRealVar::                alpha_M_EFF = 0 +/- 0.99374  L(-5 - 5)  "alpha_M_EFF"
  9) 0x84d4350 RooRealVar::   alpha_SignalFit_bbA_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbA_bjet"
 10) 0x81ba390 RooRealVar:: alpha_SignalFit_bbA_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbA_nobjet"
 11) 0x80ca7c0 RooRealVar::   alpha_SignalFit_bbH_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbH_bjet"
 12) 0x8233840 RooRealVar:: alpha_SignalFit_bbH_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbH_nobjet"
 13) 0x8040dd0 RooRealVar::   alpha_SignalFit_bbh_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbh_bjet"
 14) 0x81c7ec0 RooRealVar:: alpha_SignalFit_bbh_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_bbh_nobjet"
 15) 0x8077320 RooRealVar::   alpha_SignalFit_ggA_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggA_bjet"
 16) 0x81970c0 RooRealVar:: alpha_SignalFit_ggA_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggA_nobjet"
 17) 0x814e840 RooRealVar::   alpha_SignalFit_ggH_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggH_bjet"
 18) 0x82e7a00 RooRealVar:: alpha_SignalFit_ggH_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggH_nobjet"
 19) 0x8140000 RooRealVar::   alpha_SignalFit_ggh_bjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggh_bjet"
 20) 0x8293f30 RooRealVar:: alpha_SignalFit_ggh_nobjet = 0 +/- 0.99374  L(-5 - 5)  "alpha_SignalFit_ggh_nobjet"
 21) 0x81c0eb0 RooRealVar::               alpha_XS_bbA = 0 +/- 0.99374  L(-5 - 5)  "alpha_XS_bbA"
 22) 0x8147a10 RooRealVar::               alpha_XS_ggA = 0 +/- 0.99374  L(-5 - 5)  "alpha_XS_ggA"
DEBUG:-END printing out conditionalNuis_0 set of parameters
Using conditionalNuis_1 Success with status: 0
Trying to minimize observed lnL(mu^,theta^)
Using conditionalNuis_0 Success with status: 0
DEBUG (data)   lnL(mu=0,theta^^) =  -124631.0671 lnL(mu^,theta^)   =  -124631.1205 q0 =   0.106877 mu^ =   0.096431 Z0 =   0.326920 Z0d =   0.326920
DEBUG: Status-data   0 0
DEBUG: Status-asimov 0 0
mA170tb35 obs: 0.32692 obs-p-value: 0.371864
mA170tb35 obs-raw:0.32692 obs-raw-p-value: 0.371864
mA170tb35 obs-def:0.32692 obs-def-p-value: 0.371864
Observed significance: 0.32692
mA170tb35 Median test stat val: 8.77716
mA170tb35 med:   2.96263 med-p-value: 0.00152512
Real time 0:00:05, CP time 5.420
