HistFactory¶
A histfactory model is built as follows
- A channel is modeled by a stack of samples (signal, backgrounds)
- A measurement is composed of one or more channels
In [1]:
using namespace RooStats;
using namespace HistFactory;
Create the measurement¶
In [2]:
std::string inputFileName = "./input/example.root";
HistFactory::Measurement meas("meas", "meas");
meas.SetOutputFilePrefix( "./results/example_" );
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
Define the name of the POI
In [3]:
meas.SetPOI("mu");
Declare constants
In [4]:
meas.AddConstantParam("alpha_syst1");
meas.AddConstantParam("Lumi");
Set luminosity and its uncertainty
In [5]:
meas.SetLumi(1.0);
meas.SetLumiRelErr(0.10);
…and some technical configuration…
In [6]:
meas.SetExportOnly(false);
meas.SetBinHigh(2);
Create a channel¶
In [7]:
HistFactory::Channel chan("channel1");
chan.SetData("data",inputFileName) ;
chan.SetStatErrorConfig(0.05, "Poisson");
Populate this channel with samples
Signal sample - scaled by a normalization factor \(\mu\) and
with a 5% systematic scale uncertainty named syst1
In [8]:
HistFactory::Sample sig( "signal", "signal", inputFileName ) ;
sig.AddOverallSys("syst1",0.95,1.05) ;
sig.AddNormFactor("mu",1,0,3 ) ;
chan.AddSample(sig) ;
Background sample 1 - with a 5% systematic scale uncertainty named syst2 and with modeling of MC statistical, uncertainties
In [9]:
HistFactory::Sample background1("background1","background1",inputFileName) ;
background1.ActivateStatError("background1_statUncert",inputFileName) ;
background1.AddOverallSys("syst2",0.95,1.05) ;
chan.AddSample(background1) ;
Background sample 2 - with a 5% systematic scale uncertainty named syst3 - and with modeling of MC statistical uncertainties
In [10]:
HistFactory::Sample background2( "background2", "background2", inputFileName );
background2.ActivateStatError();
background2.AddOverallSys("syst3",0.95,1.05);
chan.AddSample( background2) ;
Add channel to the measurementa
In [11]:
meas.AddChannel( chan );
Collect all the required templated histograms from the input file
In [12]:
meas.CollectHistograms();
Getting histogram. InputFile ./input/example.root HistoPath HistoName data
Opened input file: ./input/example.root: 0x7f507c3030f0
Collecting Nominal Histogram
Getting histogram. InputFile ./input/example.root HistoPath HistoName signal
Opened input file: ./input/example.root: 0x7f507c510350
Collecting Nominal Histogram
Getting histogram. InputFile ./input/example.root HistoPath HistoName background1
Opened input file: ./input/example.root: 0x7f507c419120
Getting histogram. InputFile ./input/example.root HistoPath HistoName background1_statUncert
Opened input file: ./input/example.root: 0x7f507c5111b0
Collecting Nominal Histogram
Getting histogram. InputFile ./input/example.root HistoPath HistoName background2
Opened input file: ./input/example.root: 0x7f507c53c950
Now build a workspace with a pdf and a modelconfig
In [13]:
RooWorkspace* w = MakeModelAndMeasurementFast(meas) ;
Making Model and Measurements (Fast) for measurement: meas
using lumi = 1 and lumiError = 0.1 including bins between 0 and 2
fixing the following parameters:
alpha_syst1
Lumi
Checking if output directory : ./results - exists
Creating the output file: ./results/example__meas.root
Creating the table file: ./results/example__results.table
Creating the HistoToWorkspaceFactoryFast factory
Setting preprocess functions
Starting to process channel: channel1
-------------------
Starting to process channel1 channel with 1 observables
lumi str = [1,0,10]
lumi Error str = nominalLumi[1,0,2],0.1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon
making normFactor: mu
signal_channel1 has no variation histograms
processing hist signal
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(signal_channel1nominalDHist): fit range of variable obs_x_channel1 expanded to nearest bin boundaries: [1,2] --> [1,2]
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing dataset signal_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooHistFunc::signal_channel1_nominal
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon
background1_channel1 has no variation histograms
processing hist background1
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(background1_channel1nominalDHist): fit range of variable obs_x_channel1 expanded to nearest bin boundaries: [1,2] --> [1,2]
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing dataset background1_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooHistFunc::background1_channel1_nominal
Sample: background1 to be included in Stat Error for channel channel1
Using external histogram for Stat Errors for Channel: channel1 Sample: background1
Error Histogram: background1_statUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooRealVar::gamma_stat_channel1_bin_0
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooRealVar::gamma_stat_channel1_bin_1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing ParamHistFunc::mc_stat_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::obs_x_channel1 for import of ParamHistFunc::mc_stat_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_0 for import of ParamHistFunc::mc_stat_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_1 for import of ParamHistFunc::mc_stat_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of ParamHistFunc::mc_stat_channel1 for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::obs_x_channel1 for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_0 for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_1 for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background1_channel1_overallSyst_x_Exp for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooHistFunc::background1_channel1_nominal for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::alpha_syst2 for import of RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon
background2_channel1 has no variation histograms
processing hist background2
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(background2_channel1nominalDHist): fit range of variable obs_x_channel1 expanded to nearest bin boundaries: [1,2] --> [1,2]
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing dataset background2_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooHistFunc::background2_channel1_nominal
Sample: background2 to be included in Stat Error for channel channel1
Making Statistical Uncertainty Hist for Channel: channel1 Sample: background2
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of ParamHistFunc::mc_stat_channel1 for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::obs_x_channel1 for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_0 for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_1 for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background2_channel1_overallSyst_x_Exp for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooHistFunc::background2_channel1_nominal for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::alpha_syst3 for import of RooProduct::background2_channel1_overallSyst_x_StatUncert
Making Total Uncertainty for bin 1 Error = 5 Val = 100 RelativeError = 0.05
Making Total Uncertainty for bin 2 Error = 10 Val = 100 RelativeError = 0.1
About to create Constraint Terms from: mc_stat_channel1 params: (gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1)
Using Poisson StatErrors in channel: channel1
Creating constraint for: gamma_stat_channel1_bin_0. Type of constraint: 1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooPoisson::gamma_stat_channel1_bin_0_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooRealVar::nom_gamma_stat_channel1_bin_0
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooProduct::gamma_stat_channel1_bin_0_poisMean
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_0 for import of RooPoisson::gamma_stat_channel1_bin_0_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooConstVar::gamma_stat_channel1_bin_0_tau
Creating constraint for: gamma_stat_channel1_bin_1. Type of constraint: 1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooPoisson::gamma_stat_channel1_bin_1_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooRealVar::nom_gamma_stat_channel1_bin_1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooProduct::gamma_stat_channel1_bin_1_poisMean
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_1 for import of RooPoisson::gamma_stat_channel1_bin_1_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooConstVar::gamma_stat_channel1_bin_1_tau
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooRealSumPdf::channel1_model
-----------------------------------------
import model into workspace
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooGaussian::lumiConstraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooConstVar::0.1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::Lumi for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nominalLumi for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooGaussian::alpha_syst1Constraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooConstVar::1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::alpha_syst1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nom_alpha_syst1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooGaussian::alpha_syst2Constraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::alpha_syst2 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nom_alpha_syst2 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooGaussian::alpha_syst3Constraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::alpha_syst3 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nom_alpha_syst3 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooPoisson::gamma_stat_channel1_bin_0_constraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nom_gamma_stat_channel1_bin_0 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::gamma_stat_channel1_bin_0_poisMean for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooConstVar::gamma_stat_channel1_bin_0_tau for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_0 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooPoisson::gamma_stat_channel1_bin_1_constraint for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::nom_gamma_stat_channel1_bin_1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::gamma_stat_channel1_bin_1_poisMean for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooConstVar::gamma_stat_channel1_bin_1_tau for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::gamma_stat_channel1_bin_1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealSumPdf::channel1_model for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::L_x_signal_channel1_overallSyst_x_Exp for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::signal_channel1_overallSyst_x_Exp for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooHistFunc::signal_channel1_nominal for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::obs_x_channel1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::signal_channel1_overallNorm_x_sigma_epsilon for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::mu for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::binWidth_obs_x_channel1_0 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::L_x_background1_channel1_overallSyst_x_StatUncert for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background1_channel1_overallSyst_x_StatUncert for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of ParamHistFunc::mc_stat_channel1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background1_channel1_overallSyst_x_Exp for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooHistFunc::background1_channel1_nominal for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::binWidth_obs_x_channel1_1 for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::L_x_background2_channel1_overallSyst_x_StatUncert for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background2_channel1_overallSyst_x_StatUncert for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooProduct::background2_channel1_overallSyst_x_Exp for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooHistFunc::background2_channel1_nominal for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon for import of RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) using existing copy of RooRealVar::binWidth_obs_x_channel1_2 for import of RooProdPdf::model_channel1
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
RooDataSet::AsimovData[obs_x_channel1,weight:binWeightAsimov] = 2 entries (230 weighted)
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing dataset AsimovData
[#1] INFO:ObjectHandling -- RooWorkSpace::import(channel1) changing name of dataset from AsimovData to asimovData
[#1] INFO:ObjectHandling -- RooWorkspace::import(channel1) importing dataset obsData
RooWorkspace(channel1) channel1 workspace contents
variables
---------
(Lumi,alpha_syst1,alpha_syst2,alpha_syst3,binWidth_obs_x_channel1_0,binWidth_obs_x_channel1_1,binWidth_obs_x_channel1_2,gamma_stat_channel1_bin_0,gamma_stat_channel1_bin_1,mu,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/230
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
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) ] = 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) ] = 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) ] = 10
RooProduct::signal_channel1_overallNorm_x_sigma_epsilon[ mu * 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)
RooDataSet::obsData(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)
named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)
ModelConfig_Observables:(obs_x_channel1)
coefList:(binWidth_obs_x_channel1_0,binWidth_obs_x_channel1_1,binWidth_obs_x_channel1_2)
constraintTerms:(lumiConstraint,alpha_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
globalObservables:(nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)
likelihoodTerms:(channel1_model)
obsAndWeight:(weightVar,obs_x_channel1)
observables:(obs_x_channel1)
observablesSet:(obs_x_channel1)
shapeList:(L_x_signal_channel1_overallSyst_x_Exp,L_x_background1_channel1_overallSyst_x_StatUncert,L_x_background2_channel1_overallSyst_x_StatUncert)
generic objects
---------------
RooStats::ModelConfig::ModelConfig
Setting Parameter(s) of Interest as: mu
=== Using the following for ModelConfig ===
Observables: RooArgSet:: = (obs_x_channel1)
Parameters of Interest: RooArgSet:: = (mu)
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: 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
Opening File to hold channel: ./results/example__channel1_meas_model.root
About to write channel measurement to file
Writing sample: signal
Writing sample: background1
Writing sample: background2
Saved all histograms
Saved Measurement
Successfully wrote channel to file
---------------
---------------- Doing channel1 Fit
---------------
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Minization -- Including the following contraint terms in minimization: (lumiConstraint,alpha_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model_channel1_obsData_with_constr) Summation contains a RooNLLVar, using its error level
[#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization
[#1] INFO:Minization -- The following expressions have been identified as constant and will be precalculated and cached: (signal_channel1_nominal,background1_channel1_nominal,background2_channel1_nominal)
[#1] INFO:Minization -- The following expressions will be evaluated in cache-and-track mode: (mc_stat_channel1)
**********
** 1 **SET PRINT 1
**********
**********
** 2 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 alpha_syst2 0.00000e+00 1.00000e+00 -5.00000e+00 5.00000e+00
2 alpha_syst3 0.00000e+00 1.00000e+00 -5.00000e+00 5.00000e+00
3 gamma_stat_channel1_bin_0 1.00000e+00 1.25000e-01 0.00000e+00 1.25000e+00
4 gamma_stat_channel1_bin_1 1.00000e+00 1.50000e-01 0.00000e+00 1.50000e+00
5 mu 1.00000e+00 3.00000e-01 0.00000e+00 3.00000e+00
**********
** 3 **SET ERR 0.5
**********
**********
** 4 **SET PRINT 1
**********
**********
** 5 **SET STR 1
**********
NOW USING STRATEGY 1: TRY TO BALANCE SPEED AGAINST RELIABILITY
**********
** 6 **MIGRAD 2500 1
**********
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03
FCN=-1044.81 FROM MIGRAD STATUS=INITIATE 16 CALLS 17 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 alpha_syst2 0.00000e+00 1.00000e+00 2.01358e-01 -4.16828e-01
2 alpha_syst3 0.00000e+00 1.00000e+00 2.01358e-01 -4.54721e-01
3 gamma_stat_channel1_bin_0 1.00000e+00 1.25000e-01 2.57889e-01 -8.33599e-01
4 gamma_stat_channel1_bin_1 1.00000e+00 1.50000e-01 2.14402e-01 -1.28585e+00
5 mu 1.00000e+00 3.00000e-01 2.14402e-01 -7.28473e-01
ERR DEF= 0.5
MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-1044.84 FROM MIGRAD STATUS=CONVERGED 83 CALLS 84 TOTAL
EDM=3.27927e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 alpha_syst2 -7.65554e-03 9.82184e-01 4.06870e-03 7.64595e-03
2 alpha_syst3 1.30151e-02 9.47542e-01 4.03080e-03 -2.57978e-02
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93234e-02 2.03189e-03 1.64223e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.00861e-02 2.30372e-03 1.15819e-02
5 mu 1.11502e+00 5.86100e-01 7.54848e-03 7.38809e-04
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
9.774e-01 4.466e-02 -1.168e-03 8.943e-03 -2.108e-01
4.466e-02 9.088e-01 2.228e-03 -1.792e-02 -7.712e-02
-1.168e-03 2.228e-03 2.441e-03 4.462e-04 -1.052e-02
8.943e-03 -1.792e-02 4.462e-04 6.441e-03 -1.544e-02
-2.108e-01 -7.712e-02 -1.052e-02 -1.544e-02 3.641e-01
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.38745 1.000 0.047 -0.024 0.113 -0.353
2 0.32155 0.047 1.000 0.047 -0.234 -0.134
3 0.38693 -0.024 0.047 1.000 0.113 -0.353
4 0.42409 0.113 -0.234 0.113 1.000 -0.319
5 0.58365 -0.353 -0.134 -0.353 -0.319 1.000
**********
** 7 **SET ERR 0.5
**********
**********
** 8 **SET PRINT 1
**********
**********
** 9 **HESSE 2500
**********
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-1044.84 FROM HESSE STATUS=OK 31 CALLS 115 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 alpha_syst2 -7.65554e-03 9.82481e-01 8.13739e-04 -1.53111e-03
2 alpha_syst3 1.30151e-02 9.47601e-01 8.06160e-04 2.60302e-03
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 4.06379e-04 6.42748e-01
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 4.60745e-04 3.45220e-01
5 mu 1.11502e+00 5.86553e-01 1.50970e-03 -2.59558e-01
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
9.780e-01 4.481e-02 -1.142e-03 8.985e-03 -2.115e-01
4.481e-02 9.089e-01 2.237e-03 -1.792e-02 -7.728e-02
-1.142e-03 2.237e-03 2.442e-03 4.485e-04 -1.055e-02
8.985e-03 -1.792e-02 4.485e-04 6.445e-03 -1.549e-02
-2.115e-01 -7.728e-02 -1.055e-02 -1.549e-02 3.647e-01
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.38812 1.000 0.048 -0.023 0.113 -0.354
2 0.32173 0.048 1.000 0.047 -0.234 -0.134
3 0.38771 -0.023 0.047 1.000 0.113 -0.354
4 0.42463 0.113 -0.234 0.113 1.000 -0.320
5 0.58459 -0.354 -0.134 -0.354 -0.320 1.000
**********
** 10 **MINOS 2500 1
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 73 CALLS 188 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65554e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30151e-02 9.47601e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01
ERR DEF= 0.5
**********
** 11 **MINOS 2500 2
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 55 CALLS 243 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65554e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30151e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01
ERR DEF= 0.5
**********
** 12 **MINOS 2500 3
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 74 CALLS 317 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65554e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30151e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01
ERR DEF= 0.5
**********
** 13 **MINOS 2500 4
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 74 CALLS 391 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65554e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30151e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 -7.85700e-02 8.20232e-02
5 mu 1.11502e+00 5.86553e-01
ERR DEF= 0.5
**********
** 14 **MINOS 2500 5
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 82 CALLS 473 TOTAL
EDM=3.28068e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65554e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30151e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 -7.85700e-02 8.20232e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
[#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization
printing results for mu at 1.11502 high -0.596829 low 0.61159
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Minization -- Including the following contraint terms in minimization: (lumiConstraint,alpha_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_model_channel1_obsData_with_constr_Profile[mu]) Creating instance of MINUIT
[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model_channel1_obsData_with_constr) Summation contains a RooNLLVar, using its error level
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_model_channel1_obsData_with_constr_Profile[mu]) determining minimum likelihood for current configurations w.r.t all observable
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_model_channel1_obsData_with_constr_Profile[mu]) minimum found at (mu=1.11488)
..........................................................................................................................................................................................................full list of observables:
RooArgList:: = (obs_x_channel1)
------------------
Entering combination
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nom_alpha_syst2
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nom_alpha_syst3
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nom_gamma_stat_channel1_bin_0
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nom_gamma_stat_channel1_bin_1
-----------------------------------------
create toy data for channel1
RooDataSet::AsimovData0[obs_x_channel1,channelCat,weight:binWeightAsimov] = 2 entries (230 weighted)
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset asimovDataFullModel
[#1] INFO:ObjectHandling -- RooWorkSpace::import(combined) changing name of dataset from asimovDataFullModel to asimovData
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::obs_x_channel1
Merging data for channel channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset channel1
[#1] INFO:ObjectHandling -- RooWorkSpace::import(combined) changing name of dataset from channel1 to obsData
RooWorkspace(combined) combined contents
variables
---------
(channelCat,nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1,obs_x_channel1,weightVar)
datasets
--------
RooDataSet::asimovData(obs_x_channel1,weightVar,channelCat)
RooDataSet::obsData(channelCat,obs_x_channel1)
named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_syst2,nom_alpha_syst3,nom_gamma_stat_channel1_bin_0,nom_gamma_stat_channel1_bin_1)
ModelConfig_Observables:(obs_x_channel1,weightVar,channelCat)
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)
----------------
Importing combined model
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset signal_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset background1_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing dataset background2_channel1nominalDHist
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooCategory::channelCat for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProdPdf::model_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooGaussian::lumiConstraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooConstVar::0.1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::Lumi
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nominalLumi
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooGaussian::alpha_syst1Constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooConstVar::1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::alpha_syst1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::nom_alpha_syst1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooGaussian::alpha_syst2Constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::alpha_syst2
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooRealVar::nom_alpha_syst2 for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooGaussian::alpha_syst3Constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::alpha_syst3
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooRealVar::nom_alpha_syst3 for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooPoisson::gamma_stat_channel1_bin_0_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooRealVar::nom_gamma_stat_channel1_bin_0 for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::gamma_stat_channel1_bin_0_poisMean
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooConstVar::gamma_stat_channel1_bin_0_tau
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::gamma_stat_channel1_bin_0
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooPoisson::gamma_stat_channel1_bin_1_constraint
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooRealVar::nom_gamma_stat_channel1_bin_1 for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::gamma_stat_channel1_bin_1_poisMean
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooConstVar::gamma_stat_channel1_bin_1_tau
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::gamma_stat_channel1_bin_1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealSumPdf::channel1_model
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::L_x_signal_channel1_overallSyst_x_Exp
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::signal_channel1_overallSyst_x_Exp
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooHistFunc::signal_channel1_nominal
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) using existing copy of RooRealVar::obs_x_channel1 for import of RooSimultaneous::simPdf
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::signal_channel1_overallNorm_x_sigma_epsilon
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::mu
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::binWidth_obs_x_channel1_0
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::L_x_background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::background1_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing ParamHistFunc::mc_stat_channel1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::background1_channel1_overallSyst_x_Exp
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooHistFunc::background1_channel1_nominal
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::binWidth_obs_x_channel1_1
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::L_x_background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::background2_channel1_overallSyst_x_StatUncert
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooProduct::background2_channel1_overallSyst_x_Exp
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooHistFunc::background2_channel1_nominal
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon
[#1] INFO:ObjectHandling -- RooWorkspace::import(combined) importing RooRealVar::binWidth_obs_x_channel1_2
setting alpha_syst1 constant
setting Lumi constant
Setting Parameter(s) of Interest as: mu
=== Using the following for ModelConfig ===
Observables: RooArgSet:: = (obs_x_channel1,weightVar,channelCat)
Parameters of Interest: RooArgSet:: = (mu)
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
Writing combined workspace to file: ./results/example__combined_meas_model.root
Writing combined measurement to file: ./results/example__combined_meas_model.root
Info in <TCanvas::Print>: eps file ./results/example__channel1_meas_profileLR.eps has been created
Writing sample: signal
Writing sample: background1
Writing sample: background2
Saved all histograms
Saved Measurement
---------------
---------------- Doing combined Fit
---------------
[#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_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_simPdf_obsData_with_constr) Summation contains a RooNLLVar, using its error level
[#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization
RooAbsTestStatistic::initSimMode: creating slave calculator #0 for state channel1 (2 dataset entries)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Fitting -- RooAbsTestStatistic::initSimMode: created 1 slave calculators.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Minization -- The following expressions have been identified as constant and will be precalculated and cached: (signal_channel1_nominal,background1_channel1_nominal,background2_channel1_nominal)
[#1] INFO:Minization -- The following expressions will be evaluated in cache-and-track mode: (mc_stat_channel1)
**********
** 1842 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 alpha_syst2 0.00000e+00 1.00000e+00 -5.00000e+00 5.00000e+00
2 alpha_syst3 0.00000e+00 1.00000e+00 -5.00000e+00 5.00000e+00
3 gamma_stat_channel1_bin_0 1.00000e+00 1.25000e-01 0.00000e+00 1.25000e+00
4 gamma_stat_channel1_bin_1 1.00000e+00 1.50000e-01 0.00000e+00 1.50000e+00
5 mu 1.00000e+00 3.00000e-01 0.00000e+00 3.00000e+00
**********
** 1843 **SET ERR 0.5
**********
**********
** 1844 **SET PRINT 1
**********
**********
** 1845 **SET STR 1
**********
NOW USING STRATEGY 1: TRY TO BALANCE SPEED AGAINST RELIABILITY
**********
** 1846 **MIGRAD 2500 1
**********
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03
FCN=-1044.81 FROM MIGRAD STATUS=INITIATE 16 CALLS 17 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 alpha_syst2 0.00000e+00 1.00000e+00 2.01358e-01 -4.16836e-01
2 alpha_syst3 0.00000e+00 1.00000e+00 2.01358e-01 -4.54731e-01
3 gamma_stat_channel1_bin_0 1.00000e+00 1.25000e-01 2.57889e-01 -8.33599e-01
4 gamma_stat_channel1_bin_1 1.00000e+00 1.50000e-01 2.14402e-01 -1.28585e+00
5 mu 1.00000e+00 3.00000e-01 2.14402e-01 -7.28473e-01
ERR DEF= 0.5
MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-1044.84 FROM MIGRAD STATUS=CONVERGED 83 CALLS 84 TOTAL
EDM=3.27939e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 alpha_syst2 -7.65551e-03 9.82184e-01 4.06870e-03 7.64609e-03
2 alpha_syst3 1.30150e-02 9.47542e-01 4.03080e-03 -2.57984e-02
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93234e-02 2.03189e-03 1.64238e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.00861e-02 2.30372e-03 1.15810e-02
5 mu 1.11502e+00 5.86100e-01 7.54848e-03 7.38706e-04
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
9.774e-01 4.466e-02 -1.168e-03 8.943e-03 -2.108e-01
4.466e-02 9.088e-01 2.228e-03 -1.792e-02 -7.712e-02
-1.168e-03 2.228e-03 2.441e-03 4.462e-04 -1.052e-02
8.943e-03 -1.792e-02 4.462e-04 6.441e-03 -1.544e-02
-2.108e-01 -7.712e-02 -1.052e-02 -1.544e-02 3.641e-01
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.38745 1.000 0.047 -0.024 0.113 -0.353
2 0.32155 0.047 1.000 0.047 -0.234 -0.134
3 0.38693 -0.024 0.047 1.000 0.113 -0.353
4 0.42409 0.113 -0.234 0.113 1.000 -0.319
5 0.58365 -0.353 -0.134 -0.353 -0.319 1.000
**********
** 1847 **SET ERR 0.5
**********
**********
** 1848 **SET PRINT 1
**********
**********
** 1849 **HESSE 2500
**********
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-1044.84 FROM HESSE STATUS=OK 31 CALLS 115 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 alpha_syst2 -7.65551e-03 9.82481e-01 8.13739e-04 -1.53110e-03
2 alpha_syst3 1.30150e-02 9.47601e-01 8.06160e-04 2.60300e-03
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 4.06379e-04 6.42748e-01
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 4.60745e-04 3.45220e-01
5 mu 1.11502e+00 5.86553e-01 1.50970e-03 -2.59558e-01
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
9.780e-01 4.481e-02 -1.142e-03 8.985e-03 -2.115e-01
4.481e-02 9.089e-01 2.237e-03 -1.792e-02 -7.728e-02
-1.142e-03 2.237e-03 2.442e-03 4.485e-04 -1.055e-02
8.985e-03 -1.792e-02 4.485e-04 6.445e-03 -1.549e-02
-2.115e-01 -7.728e-02 -1.055e-02 -1.549e-02 3.647e-01
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.38812 1.000 0.048 -0.023 0.113 -0.354
2 0.32173 0.048 1.000 0.047 -0.234 -0.134
3 0.38771 -0.023 0.047 1.000 0.113 -0.354
4 0.42463 0.113 -0.234 0.113 1.000 -0.320
5 0.58459 -0.354 -0.134 -0.354 -0.320 1.000
**********
** 1850 **MINOS 2500 1
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 73 CALLS 188 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65551e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30150e-02 9.47601e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
**********
** 1851 **MINOS 2500 2
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 55 CALLS 243 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65551e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30150e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
**********
** 1852 **MINOS 2500 3
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 74 CALLS 317 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65551e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30150e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
**********
** 1853 **MINOS 2500 4
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 74 CALLS 391 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65551e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30150e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 -7.85700e-02 8.20232e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
**********
** 1854 **MINOS 2500 5
**********
FCN=-1044.84 FROM MINOS STATUS=SUCCESSFUL 82 CALLS 473 TOTAL
EDM=3.2808e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 alpha_syst2 -7.65551e-03 9.82481e-01 -9.90051e-01 9.87631e-01
2 alpha_syst3 1.30150e-02 9.47601e-01 -9.47046e-01 9.60252e-01
3 gamma_stat_channel1_bin_0 9.99623e-01 4.93408e-02 -4.86933e-02 5.01526e-02
4 gamma_stat_channel1_bin_1 1.00380e+00 8.01088e-02 -7.85700e-02 8.20232e-02
5 mu 1.11502e+00 5.86553e-01 -5.96829e-01 6.11590e-01
ERR DEF= 0.5
[#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization
printing results for mu at 1.11502 high -0.596829 low 0.61159
[#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_syst1Constraint,alpha_syst2Constraint,alpha_syst3Constraint,gamma_stat_channel1_bin_0_constraint,gamma_stat_channel1_bin_1_constraint)
RooAbsTestStatistic::initSimMode: creating slave calculator #0 for state channel1 (2 dataset entries)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Fitting -- RooAbsTestStatistic::initSimMode: created 1 slave calculators.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_simPdf_obsData_with_constr_Profile[mu]) 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[mu]) 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:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Fitting -- RooAbsTestStatistic::initSimMode: created 1 slave calculators.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_0_constraint_Int[gamma_stat_channel1_bin_0]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_0)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gamma_stat_channel1_bin_1_constraint_Int[gamma_stat_channel1_bin_1]) using numeric integrator RooIntegrator1D to calculate Int(gamma_stat_channel1_bin_1)
[#1] INFO:Minization -- RooProfileLL::evaluate(nll_simPdf_obsData_with_constr_Profile[mu]) minimum found at (mu=1.11488)
..........................................................................................................................................................................................................
RooAbsTestStatistic::initSimMode: creating slave calculator #0 for state channel1 (2 dataset entries)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
[#1] INFO:Fitting -- RooAbsTestStatistic::initSimMode: created 1 slave calculators.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(channel1_model_Int[obs_x_channel1]) using numeric integrator RooBinIntegrator to calculate Int(obs_x_channel1)
Info in <TCanvas::Print>: eps file ./results/example__combined_meas_profileLR.eps has been created
Make some adjustments to histfactory workspace to conform to naming conventions of this set of example macros
In [14]:
w->SetName("w") ;
w->data("obsData")->SetName("observed_data") ;
Print resulting workspace and save to file
In [15]:
w->Print("t") ;
w->writeToFile("model.root") ;
RooWorkspace(w) combined contents
variables
---------
(Lumi,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,mu,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
-------
RooSimultaneous::simPdf[ indexCat=channelCat channel1=model_channel1 ] = 0.174888
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
RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.1 ] = 1/0.250663
RooGaussian::alpha_syst1Constraint[ x=alpha_syst1 mean=nom_alpha_syst1 sigma=1 ] = 1/2.50663
RooGaussian::alpha_syst2Constraint[ x=alpha_syst2 mean=nom_alpha_syst2 sigma=1 ] = 1/2.50663
RooGaussian::alpha_syst3Constraint[ x=alpha_syst3 mean=nom_alpha_syst3 sigma=1 ] = 1/2.50663
RooPoisson::gamma_stat_channel1_bin_0_constraint[ x=nom_gamma_stat_channel1_bin_0 mean=gamma_stat_channel1_bin_0_poisMean ] = 0.019943/0.00250021
RooProduct::gamma_stat_channel1_bin_0_poisMean[ gamma_stat_channel1_bin_0 * gamma_stat_channel1_bin_0_tau ] = 400
RooPoisson::gamma_stat_channel1_bin_1_constraint[ x=nom_gamma_stat_channel1_bin_1 mean=gamma_stat_channel1_bin_1_poisMean ] = 0.039861/0.00999989
RooProduct::gamma_stat_channel1_bin_1_poisMean[ gamma_stat_channel1_bin_1 * gamma_stat_channel1_bin_1_tau ] = 100
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
RooProduct::L_x_signal_channel1_overallSyst_x_Exp[ Lumi * signal_channel1_overallSyst_x_Exp ] = 10
RooProduct::signal_channel1_overallSyst_x_Exp[ signal_channel1_nominal * signal_channel1_overallNorm_x_sigma_epsilon ] = 10
RooHistFunc::signal_channel1_nominal[ depList=(obs_x_channel1) ] = 10
RooProduct::signal_channel1_overallNorm_x_sigma_epsilon[ mu * signal_channel1_epsilon ] = 1
RooStats::HistFactory::FlexibleInterpVar::signal_channel1_epsilon[ paramList=(alpha_syst1) ] = 1
RooProduct::L_x_background1_channel1_overallSyst_x_StatUncert[ Lumi * background1_channel1_overallSyst_x_StatUncert ] = 0
RooProduct::background1_channel1_overallSyst_x_StatUncert[ mc_stat_channel1 * background1_channel1_overallSyst_x_Exp ] = 0
ParamHistFunc::mc_stat_channel1[ ] = 1
RooProduct::background1_channel1_overallSyst_x_Exp[ background1_channel1_nominal * background1_channel1_epsilon ] = 0
RooHistFunc::background1_channel1_nominal[ depList=(obs_x_channel1) ] = 0
RooStats::HistFactory::FlexibleInterpVar::background1_channel1_epsilon[ paramList=(alpha_syst2) ] = 1
RooProduct::L_x_background2_channel1_overallSyst_x_StatUncert[ Lumi * background2_channel1_overallSyst_x_StatUncert ] = 100
RooProduct::background2_channel1_overallSyst_x_StatUncert[ mc_stat_channel1 * background2_channel1_overallSyst_x_Exp ] = 100
ParamHistFunc::mc_stat_channel1[ ] = 1
RooProduct::background2_channel1_overallSyst_x_Exp[ background2_channel1_nominal * background2_channel1_epsilon ] = 100
RooHistFunc::background2_channel1_nominal[ depList=(obs_x_channel1) ] = 100
RooStats::HistFactory::FlexibleInterpVar::background2_channel1_epsilon[ paramList=(alpha_syst3) ] = 1
datasets
--------
RooDataSet::asimovData(obs_x_channel1,weightVar,channelCat)
RooDataSet::observed_data(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,mu=1,binWidth_obs_x_channel1_0=2[C],binWidth_obs_x_channel1_1=2[C],binWidth_obs_x_channel1_2=2[C])
InitialValues = (alpha_syst2=0,alpha_syst3=0,gamma_stat_channel1_bin_0=1,gamma_stat_channel1_bin_1=1,mu=1)
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:(mu)
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