notebook.community

Edit and run 





In [8]:



    


import ROOT
ROOT.RooMsgService.instance().setGlobalKillBelow(5)



    











In [9]:



    


from Builder import get_workspace



    











In [3]:



    


workspace = get_workspace(nchannels = 3, events = 1000, nbins = 1, nnps = 3)



    


















    






this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2


-------------------
Starting to process Region0 channel with 1 observables
lumi str = [1,0,10]
lumi Error str = nominalLumi[1,0,1.2],0.02
Gaussian::alpha_background_0_shape_0Constraint(alpha_background_0_shape_0,nom_alpha_background_0_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_1Constraint(alpha_background_0_shape_1,nom_alpha_background_0_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_2Constraint(alpha_background_0_shape_2,nom_alpha_background_0_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_0Constraint(alpha_background_1_shape_0,nom_alpha_background_1_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_1Constraint(alpha_background_1_shape_1,nom_alpha_background_1_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_2Constraint(alpha_background_1_shape_2,nom_alpha_background_1_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_0Constraint(alpha_background_2_shape_0,nom_alpha_background_2_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_1Constraint(alpha_background_2_shape_1,nom_alpha_background_2_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_2Constraint(alpha_background_2_shape_2,nom_alpha_background_2_shape_2[0.,-10,10],1.)
making normFactor: SignalStrength
Gaussian::alpha_signal_shape_0Constraint(alpha_signal_shape_0,nom_alpha_signal_shape_0[0.,-10,10],1.)
Gaussian::alpha_signal_shape_1Constraint(alpha_signal_shape_1,nom_alpha_signal_shape_1[0.,-10,10],1.)
Gaussian::alpha_signal_shape_2Constraint(alpha_signal_shape_2,nom_alpha_signal_shape_2[0.,-10,10],1.)
-----------------------------------------
import model into workspace
RooDataSet::AsimovData[obs_x_Region0,weight:binWeightAsimov] = 1 entries (990 weighted)

RooWorkspace(Region0) Region0 workspace contents

variables
---------
(Lumi,SignalStrength,alpha_background_0_shape_0,alpha_background_0_shape_1,alpha_background_0_shape_2,alpha_background_1_shape_0,alpha_background_1_shape_1,alpha_background_1_shape_2,alpha_background_2_shape_0,alpha_background_2_shape_1,alpha_background_2_shape_2,alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2,alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2,alpha_signal_shape_0,alpha_signal_shape_1,alpha_signal_shape_2,binWidth_obs_x_Region0_0,binWidth_obs_x_Region0_1,binWidth_obs_x_Region0_2,binWidth_obs_x_Region0_3,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2,nominalLumi,obs_x_Region0,weightVar)

p.d.f.s
-------
RooRealSumPdf::Region0_model[ binWidth_obs_x_Region0_0 * L_x_model0_Region0_overallSyst_x_HistSyst + binWidth_obs_x_Region0_1 * L_x_model1_Region0_overallSyst_x_HistSyst + binWidth_obs_x_Region0_2 * L_x_model2_Region0_overallSyst_x_HistSyst + binWidth_obs_x_Region0_3 * L_x_signal_Region0_overallSyst_x_HistSyst ] = 990/990
RooGaussian::alpha_background_0_shape_0Constraint[ x=alpha_background_0_shape_0 mean=nom_alpha_background_0_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_1Constraint[ x=alpha_background_0_shape_1 mean=nom_alpha_background_0_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_2Constraint[ x=alpha_background_0_shape_2 mean=nom_alpha_background_0_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_0Constraint[ x=alpha_background_1_shape_0 mean=nom_alpha_background_1_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_1Constraint[ x=alpha_background_1_shape_1 mean=nom_alpha_background_1_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_2Constraint[ x=alpha_background_1_shape_2 mean=nom_alpha_background_1_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_0Constraint[ x=alpha_background_2_shape_0 mean=nom_alpha_background_2_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_1Constraint[ x=alpha_background_2_shape_1 mean=nom_alpha_background_2_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_2Constraint[ x=alpha_background_2_shape_2 mean=nom_alpha_background_2_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_0Constraint[ x=alpha_background_norm_uncertainty_0 mean=nom_alpha_background_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_1Constraint[ x=alpha_background_norm_uncertainty_1 mean=nom_alpha_background_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_2Constraint[ x=alpha_background_norm_uncertainty_2 mean=nom_alpha_background_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_0Constraint[ x=alpha_signal_norm_uncertainty_0 mean=nom_alpha_signal_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_1Constraint[ x=alpha_signal_norm_uncertainty_1 mean=nom_alpha_signal_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_2Constraint[ x=alpha_signal_norm_uncertainty_2 mean=nom_alpha_signal_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_0Constraint[ x=alpha_signal_shape_0 mean=nom_alpha_signal_shape_0 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_1Constraint[ x=alpha_signal_shape_1 mean=nom_alpha_signal_shape_1 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_2Constraint[ x=alpha_signal_shape_2 mean=nom_alpha_signal_shape_2 sigma=1 ] = 1
RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.02 ] = 1
RooProdPdf::model_Region0[ lumiConstraint * alpha_background_norm_uncertainty_0Constraint * alpha_background_norm_uncertainty_1Constraint * alpha_background_norm_uncertainty_2Constraint * alpha_background_0_shape_0Constraint * alpha_background_0_shape_1Constraint * alpha_background_0_shape_2Constraint * alpha_background_1_shape_0Constraint * alpha_background_1_shape_1Constraint * alpha_background_1_shape_2Constraint * alpha_background_2_shape_0Constraint * alpha_background_2_shape_1Constraint * alpha_background_2_shape_2Constraint * alpha_signal_norm_uncertainty_0Constraint * alpha_signal_norm_uncertainty_1Constraint * alpha_signal_norm_uncertainty_2Constraint * alpha_signal_shape_0Constraint * alpha_signal_shape_1Constraint * alpha_signal_shape_2Constraint * Region0_model(obs_x_Region0) ] = 990

functions
--------
RooProduct::L_x_model0_Region0_overallSyst_x_HistSyst[ 1 * model0_Region0_overallSyst_x_HistSyst ] = 300
RooProduct::L_x_model1_Region0_overallSyst_x_HistSyst[ 1 * model1_Region0_overallSyst_x_HistSyst ] = 300
RooProduct::L_x_model2_Region0_overallSyst_x_HistSyst[ 1 * model2_Region0_overallSyst_x_HistSyst ] = 300
RooProduct::L_x_signal_Region0_overallSyst_x_HistSyst[ 1 * signal_Region0_overallSyst_x_HistSyst ] = 90
PiecewiseInterpolation::model0_Region0_Hist_alpha[ ] = 300
RooHistFunc::model0_Region0_Hist_alpha_0high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model0_Region0_Hist_alpha_0low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model0_Region0_Hist_alpha_1high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model0_Region0_Hist_alpha_1low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model0_Region0_Hist_alpha_2high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model0_Region0_Hist_alpha_2low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model0_Region0_Hist_alphanominal[ depList=(obs_x_Region0) ] = 300
RooStats::HistFactory::FlexibleInterpVar::model0_Region0_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model0_Region0_overallSyst_x_HistSyst[ model0_Region0_Hist_alpha * model0_Region0_epsilon ] = 300
PiecewiseInterpolation::model1_Region0_Hist_alpha[ ] = 300
RooHistFunc::model1_Region0_Hist_alpha_0high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model1_Region0_Hist_alpha_0low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model1_Region0_Hist_alpha_1high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model1_Region0_Hist_alpha_1low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model1_Region0_Hist_alpha_2high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model1_Region0_Hist_alpha_2low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model1_Region0_Hist_alphanominal[ depList=(obs_x_Region0) ] = 300
RooStats::HistFactory::FlexibleInterpVar::model1_Region0_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model1_Region0_overallSyst_x_HistSyst[ model1_Region0_Hist_alpha * model1_Region0_epsilon ] = 300
PiecewiseInterpolation::model2_Region0_Hist_alpha[ ] = 300
RooHistFunc::model2_Region0_Hist_alpha_0high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model2_Region0_Hist_alpha_0low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model2_Region0_Hist_alpha_1high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model2_Region0_Hist_alpha_1low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model2_Region0_Hist_alpha_2high[ depList=(obs_x_Region0) ] = 330
RooHistFunc::model2_Region0_Hist_alpha_2low[ depList=(obs_x_Region0) ] = 270
RooHistFunc::model2_Region0_Hist_alphanominal[ depList=(obs_x_Region0) ] = 300
RooStats::HistFactory::FlexibleInterpVar::model2_Region0_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model2_Region0_overallSyst_x_HistSyst[ model2_Region0_Hist_alpha * model2_Region0_epsilon ] = 300
PiecewiseInterpolation::signal_Region0_Hist_alpha[ ] = 90
RooHistFunc::signal_Region0_Hist_alpha_0high[ depList=(obs_x_Region0) ] = 99
RooHistFunc::signal_Region0_Hist_alpha_0low[ depList=(obs_x_Region0) ] = 81
RooHistFunc::signal_Region0_Hist_alpha_1high[ depList=(obs_x_Region0) ] = 99
RooHistFunc::signal_Region0_Hist_alpha_1low[ depList=(obs_x_Region0) ] = 81
RooHistFunc::signal_Region0_Hist_alpha_2high[ depList=(obs_x_Region0) ] = 99
RooHistFunc::signal_Region0_Hist_alpha_2low[ depList=(obs_x_Region0) ] = 81
RooHistFunc::signal_Region0_Hist_alphanominal[ depList=(obs_x_Region0) ] = 90
RooStats::HistFactory::FlexibleInterpVar::signal_Region0_epsilon[ paramList=(alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2) ] = 1
RooProduct::signal_Region0_overallNorm_x_sigma_epsilon[ SignalStrength * signal_Region0_epsilon ] = 1
RooProduct::signal_Region0_overallSyst_x_HistSyst[ signal_Region0_Hist_alpha * signal_Region0_overallNorm_x_sigma_epsilon ] = 90

datasets
--------
RooDataSet::asimovData(obs_x_Region0)
RooDataSet::obsData(obs_x_Region0)

embedded datasets (in pdfs and functions)
-----------------------------------------
RooDataHist::model0_Region0_Hist_alphanominalDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_0lowDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_1lowDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_2lowDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_0highDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_1highDHist(obs_x_Region0)
RooDataHist::model0_Region0_Hist_alpha_2highDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alphanominalDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_0lowDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_1lowDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_2lowDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_0highDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_1highDHist(obs_x_Region0)
RooDataHist::model1_Region0_Hist_alpha_2highDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alphanominalDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_0lowDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_1lowDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_2lowDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_0highDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_1highDHist(obs_x_Region0)
RooDataHist::model2_Region0_Hist_alpha_2highDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alphanominalDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_0lowDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_1lowDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_2lowDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_0highDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_1highDHist(obs_x_Region0)
RooDataHist::signal_Region0_Hist_alpha_2highDHist(obs_x_Region0)

named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
ModelConfig_Observables:(obs_x_Region0)
coefList:(binWidth_obs_x_Region0_0,binWidth_obs_x_Region0_1,binWidth_obs_x_Region0_2,binWidth_obs_x_Region0_3)
constraintTerms:(lumiConstraint,alpha_background_norm_uncertainty_0Constraint,alpha_background_norm_uncertainty_1Constraint,alpha_background_norm_uncertainty_2Constraint,alpha_background_0_shape_0Constraint,alpha_background_0_shape_1Constraint,alpha_background_0_shape_2Constraint,alpha_background_1_shape_0Constraint,alpha_background_1_shape_1Constraint,alpha_background_1_shape_2Constraint,alpha_background_2_shape_0Constraint,alpha_background_2_shape_1Constraint,alpha_background_2_shape_2Constraint,alpha_signal_norm_uncertainty_0Constraint,alpha_signal_norm_uncertainty_1Constraint,alpha_signal_norm_uncertainty_2Constraint,alpha_signal_shape_0Constraint,alpha_signal_shape_1Constraint,alpha_signal_shape_2Constraint)
globalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
likelihoodTerms:(Region0_model)
obsAndWeight:(weightVar,obs_x_Region0)
observables:(obs_x_Region0)
observablesSet:(obs_x_Region0)
shapeList:(L_x_model0_Region0_overallSyst_x_HistSyst,L_x_model1_Region0_overallSyst_x_HistSyst,L_x_model2_Region0_overallSyst_x_HistSyst,L_x_signal_Region0_overallSyst_x_HistSyst)

generic objects
---------------
RooStats::ModelConfig::ModelConfig

Setting Parameter(s) of Interest as: SignalStrength 


-------------------
Starting to process Region1 channel with 1 observables
lumi str = [1,0,10]
lumi Error str = nominalLumi[1,0,1.2],0.02
Gaussian::alpha_background_0_shape_0Constraint(alpha_background_0_shape_0,nom_alpha_background_0_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_1Constraint(alpha_background_0_shape_1,nom_alpha_background_0_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_2Constraint(alpha_background_0_shape_2,nom_alpha_background_0_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_0Constraint(alpha_background_1_shape_0,nom_alpha_background_1_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_1Constraint(alpha_background_1_shape_1,nom_alpha_background_1_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_2Constraint(alpha_background_1_shape_2,nom_alpha_background_1_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_0Constraint(alpha_background_2_shape_0,nom_alpha_background_2_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_1Constraint(alpha_background_2_shape_1,nom_alpha_background_2_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_2Constraint(alpha_background_2_shape_2,nom_alpha_background_2_shape_2[0.,-10,10],1.)
making normFactor: SignalStrength
Gaussian::alpha_signal_shape_0Constraint(alpha_signal_shape_0,nom_alpha_signal_shape_0[0.,-10,10],1.)
Gaussian::alpha_signal_shape_1Constraint(alpha_signal_shape_1,nom_alpha_signal_shape_1[0.,-10,10],1.)
Gaussian::alpha_signal_shape_2Constraint(alpha_signal_shape_2,nom_alpha_signal_shape_2[0.,-10,10],1.)
-----------------------------------------
import model into workspace
RooDataSet::AsimovData[obs_x_Region1,weight:binWeightAsimov] = 1 entries (1089 weighted)

RooWorkspace(Region1) Region1 workspace contents

variables
---------
(Lumi,SignalStrength,alpha_background_0_shape_0,alpha_background_0_shape_1,alpha_background_0_shape_2,alpha_background_1_shape_0,alpha_background_1_shape_1,alpha_background_1_shape_2,alpha_background_2_shape_0,alpha_background_2_shape_1,alpha_background_2_shape_2,alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2,alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2,alpha_signal_shape_0,alpha_signal_shape_1,alpha_signal_shape_2,binWidth_obs_x_Region1_0,binWidth_obs_x_Region1_1,binWidth_obs_x_Region1_2,binWidth_obs_x_Region1_3,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2,nominalLumi,obs_x_Region1,weightVar)

p.d.f.s
-------
RooRealSumPdf::Region1_model[ binWidth_obs_x_Region1_0 * L_x_model0_Region1_overallSyst_x_HistSyst + binWidth_obs_x_Region1_1 * L_x_model1_Region1_overallSyst_x_HistSyst + binWidth_obs_x_Region1_2 * L_x_model2_Region1_overallSyst_x_HistSyst + binWidth_obs_x_Region1_3 * L_x_signal_Region1_overallSyst_x_HistSyst ] = 1089/1089
RooGaussian::alpha_background_0_shape_0Constraint[ x=alpha_background_0_shape_0 mean=nom_alpha_background_0_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_1Constraint[ x=alpha_background_0_shape_1 mean=nom_alpha_background_0_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_2Constraint[ x=alpha_background_0_shape_2 mean=nom_alpha_background_0_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_0Constraint[ x=alpha_background_1_shape_0 mean=nom_alpha_background_1_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_1Constraint[ x=alpha_background_1_shape_1 mean=nom_alpha_background_1_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_2Constraint[ x=alpha_background_1_shape_2 mean=nom_alpha_background_1_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_0Constraint[ x=alpha_background_2_shape_0 mean=nom_alpha_background_2_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_1Constraint[ x=alpha_background_2_shape_1 mean=nom_alpha_background_2_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_2Constraint[ x=alpha_background_2_shape_2 mean=nom_alpha_background_2_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_0Constraint[ x=alpha_background_norm_uncertainty_0 mean=nom_alpha_background_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_1Constraint[ x=alpha_background_norm_uncertainty_1 mean=nom_alpha_background_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_2Constraint[ x=alpha_background_norm_uncertainty_2 mean=nom_alpha_background_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_0Constraint[ x=alpha_signal_norm_uncertainty_0 mean=nom_alpha_signal_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_1Constraint[ x=alpha_signal_norm_uncertainty_1 mean=nom_alpha_signal_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_2Constraint[ x=alpha_signal_norm_uncertainty_2 mean=nom_alpha_signal_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_0Constraint[ x=alpha_signal_shape_0 mean=nom_alpha_signal_shape_0 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_1Constraint[ x=alpha_signal_shape_1 mean=nom_alpha_signal_shape_1 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_2Constraint[ x=alpha_signal_shape_2 mean=nom_alpha_signal_shape_2 sigma=1 ] = 1
RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.02 ] = 1
RooProdPdf::model_Region1[ lumiConstraint * alpha_background_norm_uncertainty_0Constraint * alpha_background_norm_uncertainty_1Constraint * alpha_background_norm_uncertainty_2Constraint * alpha_background_0_shape_0Constraint * alpha_background_0_shape_1Constraint * alpha_background_0_shape_2Constraint * alpha_background_1_shape_0Constraint * alpha_background_1_shape_1Constraint * alpha_background_1_shape_2Constraint * alpha_background_2_shape_0Constraint * alpha_background_2_shape_1Constraint * alpha_background_2_shape_2Constraint * alpha_signal_norm_uncertainty_0Constraint * alpha_signal_norm_uncertainty_1Constraint * alpha_signal_norm_uncertainty_2Constraint * alpha_signal_shape_0Constraint * alpha_signal_shape_1Constraint * alpha_signal_shape_2Constraint * Region1_model(obs_x_Region1) ] = 1089

functions
--------
RooProduct::L_x_model0_Region1_overallSyst_x_HistSyst[ 1 * model0_Region1_overallSyst_x_HistSyst ] = 90
RooProduct::L_x_model1_Region1_overallSyst_x_HistSyst[ 1 * model1_Region1_overallSyst_x_HistSyst ] = 900
RooProduct::L_x_model2_Region1_overallSyst_x_HistSyst[ 1 * model2_Region1_overallSyst_x_HistSyst ] = 90
RooProduct::L_x_signal_Region1_overallSyst_x_HistSyst[ 1 * signal_Region1_overallSyst_x_HistSyst ] = 9
PiecewiseInterpolation::model0_Region1_Hist_alpha[ ] = 90
RooHistFunc::model0_Region1_Hist_alpha_0high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model0_Region1_Hist_alpha_0low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model0_Region1_Hist_alpha_1high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model0_Region1_Hist_alpha_1low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model0_Region1_Hist_alpha_2high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model0_Region1_Hist_alpha_2low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model0_Region1_Hist_alphanominal[ depList=(obs_x_Region1) ] = 90
RooStats::HistFactory::FlexibleInterpVar::model0_Region1_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model0_Region1_overallSyst_x_HistSyst[ model0_Region1_Hist_alpha * model0_Region1_epsilon ] = 90
PiecewiseInterpolation::model1_Region1_Hist_alpha[ ] = 900
RooHistFunc::model1_Region1_Hist_alpha_0high[ depList=(obs_x_Region1) ] = 990
RooHistFunc::model1_Region1_Hist_alpha_0low[ depList=(obs_x_Region1) ] = 810
RooHistFunc::model1_Region1_Hist_alpha_1high[ depList=(obs_x_Region1) ] = 990
RooHistFunc::model1_Region1_Hist_alpha_1low[ depList=(obs_x_Region1) ] = 810
RooHistFunc::model1_Region1_Hist_alpha_2high[ depList=(obs_x_Region1) ] = 990
RooHistFunc::model1_Region1_Hist_alpha_2low[ depList=(obs_x_Region1) ] = 810
RooHistFunc::model1_Region1_Hist_alphanominal[ depList=(obs_x_Region1) ] = 900
RooStats::HistFactory::FlexibleInterpVar::model1_Region1_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model1_Region1_overallSyst_x_HistSyst[ model1_Region1_Hist_alpha * model1_Region1_epsilon ] = 900
PiecewiseInterpolation::model2_Region1_Hist_alpha[ ] = 90
RooHistFunc::model2_Region1_Hist_alpha_0high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model2_Region1_Hist_alpha_0low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model2_Region1_Hist_alpha_1high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model2_Region1_Hist_alpha_1low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model2_Region1_Hist_alpha_2high[ depList=(obs_x_Region1) ] = 99
RooHistFunc::model2_Region1_Hist_alpha_2low[ depList=(obs_x_Region1) ] = 81
RooHistFunc::model2_Region1_Hist_alphanominal[ depList=(obs_x_Region1) ] = 90
RooStats::HistFactory::FlexibleInterpVar::model2_Region1_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model2_Region1_overallSyst_x_HistSyst[ model2_Region1_Hist_alpha * model2_Region1_epsilon ] = 90
PiecewiseInterpolation::signal_Region1_Hist_alpha[ ] = 9
RooHistFunc::signal_Region1_Hist_alpha_0high[ depList=(obs_x_Region1) ] = 9.9
RooHistFunc::signal_Region1_Hist_alpha_0low[ depList=(obs_x_Region1) ] = 8.1
RooHistFunc::signal_Region1_Hist_alpha_1high[ depList=(obs_x_Region1) ] = 9.9
RooHistFunc::signal_Region1_Hist_alpha_1low[ depList=(obs_x_Region1) ] = 8.1
RooHistFunc::signal_Region1_Hist_alpha_2high[ depList=(obs_x_Region1) ] = 9.9
RooHistFunc::signal_Region1_Hist_alpha_2low[ depList=(obs_x_Region1) ] = 8.1
RooHistFunc::signal_Region1_Hist_alphanominal[ depList=(obs_x_Region1) ] = 9
RooStats::HistFactory::FlexibleInterpVar::signal_Region1_epsilon[ paramList=(alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2) ] = 1
RooProduct::signal_Region1_overallNorm_x_sigma_epsilon[ SignalStrength * signal_Region1_epsilon ] = 1
RooProduct::signal_Region1_overallSyst_x_HistSyst[ signal_Region1_Hist_alpha * signal_Region1_overallNorm_x_sigma_epsilon ] = 9

datasets
--------
RooDataSet::asimovData(obs_x_Region1)
RooDataSet::obsData(obs_x_Region1)

embedded datasets (in pdfs and functions)
-----------------------------------------
RooDataHist::model0_Region1_Hist_alphanominalDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_0lowDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_1lowDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_2lowDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_0highDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_1highDHist(obs_x_Region1)
RooDataHist::model0_Region1_Hist_alpha_2highDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alphanominalDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_0lowDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_1lowDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_2lowDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_0highDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_1highDHist(obs_x_Region1)
RooDataHist::model1_Region1_Hist_alpha_2highDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alphanominalDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_0lowDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_1lowDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_2lowDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_0highDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_1highDHist(obs_x_Region1)
RooDataHist::model2_Region1_Hist_alpha_2highDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alphanominalDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_0lowDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_1lowDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_2lowDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_0highDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_1highDHist(obs_x_Region1)
RooDataHist::signal_Region1_Hist_alpha_2highDHist(obs_x_Region1)

named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
ModelConfig_Observables:(obs_x_Region1)
coefList:(binWidth_obs_x_Region1_0,binWidth_obs_x_Region1_1,binWidth_obs_x_Region1_2,binWidth_obs_x_Region1_3)
constraintTerms:(lumiConstraint,alpha_background_norm_uncertainty_0Constraint,alpha_background_norm_uncertainty_1Constraint,alpha_background_norm_uncertainty_2Constraint,alpha_background_0_shape_0Constraint,alpha_background_0_shape_1Constraint,alpha_background_0_shape_2Constraint,alpha_background_1_shape_0Constraint,alpha_background_1_shape_1Constraint,alpha_background_1_shape_2Constraint,alpha_background_2_shape_0Constraint,alpha_background_2_shape_1Constraint,alpha_background_2_shape_2Constraint,alpha_signal_norm_uncertainty_0Constraint,alpha_signal_norm_uncertainty_1Constraint,alpha_signal_norm_uncertainty_2Constraint,alpha_signal_shape_0Constraint,alpha_signal_shape_1Constraint,alpha_signal_shape_2Constraint)
globalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
likelihoodTerms:(Region1_model)
obsAndWeight:(weightVar,obs_x_Region1)
observables:(obs_x_Region1)
observablesSet:(obs_x_Region1)
shapeList:(L_x_model0_Region1_overallSyst_x_HistSyst,L_x_model1_Region1_overallSyst_x_HistSyst,L_x_model2_Region1_overallSyst_x_HistSyst,L_x_signal_Region1_overallSyst_x_HistSyst)

generic objects
---------------
RooStats::ModelConfig::ModelConfig

Setting Parameter(s) of Interest as: SignalStrength 


-------------------
Starting to process Region2 channel with 1 observables
lumi str = [1,0,10]
lumi Error str = nominalLumi[1,0,1.2],0.02
Gaussian::alpha_background_0_shape_0Constraint(alpha_background_0_shape_0,nom_alpha_background_0_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_1Constraint(alpha_background_0_shape_1,nom_alpha_background_0_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_0_shape_2Constraint(alpha_background_0_shape_2,nom_alpha_background_0_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_0Constraint(alpha_background_1_shape_0,nom_alpha_background_1_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_1Constraint(alpha_background_1_shape_1,nom_alpha_background_1_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_1_shape_2Constraint(alpha_background_1_shape_2,nom_alpha_background_1_shape_2[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_0Constraint(alpha_background_2_shape_0,nom_alpha_background_2_shape_0[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_1Constraint(alpha_background_2_shape_1,nom_alpha_background_2_shape_1[0.,-10,10],1.)
Gaussian::alpha_background_2_shape_2Constraint(alpha_background_2_shape_2,nom_alpha_background_2_shape_2[0.,-10,10],1.)
making normFactor: SignalStrength
Gaussian::alpha_signal_shape_0Constraint(alpha_signal_shape_0,nom_alpha_signal_shape_0[0.,-10,10],1.)
Gaussian::alpha_signal_shape_1Constraint(alpha_signal_shape_1,nom_alpha_signal_shape_1[0.,-10,10],1.)
Gaussian::alpha_signal_shape_2Constraint(alpha_signal_shape_2,nom_alpha_signal_shape_2[0.,-10,10],1.)
-----------------------------------------
import model into workspace
RooDataSet::AsimovData[obs_x_Region2,weight:binWeightAsimov] = 1 entries (1089 weighted)

RooWorkspace(Region2) Region2 workspace contents

variables
---------
(Lumi,SignalStrength,alpha_background_0_shape_0,alpha_background_0_shape_1,alpha_background_0_shape_2,alpha_background_1_shape_0,alpha_background_1_shape_1,alpha_background_1_shape_2,alpha_background_2_shape_0,alpha_background_2_shape_1,alpha_background_2_shape_2,alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2,alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2,alpha_signal_shape_0,alpha_signal_shape_1,alpha_signal_shape_2,binWidth_obs_x_Region2_0,binWidth_obs_x_Region2_1,binWidth_obs_x_Region2_2,binWidth_obs_x_Region2_3,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2,nominalLumi,obs_x_Region2,weightVar)

p.d.f.s
-------
RooRealSumPdf::Region2_model[ binWidth_obs_x_Region2_0 * L_x_model0_Region2_overallSyst_x_HistSyst + binWidth_obs_x_Region2_1 * L_x_model1_Region2_overallSyst_x_HistSyst + binWidth_obs_x_Region2_2 * L_x_model2_Region2_overallSyst_x_HistSyst + binWidth_obs_x_Region2_3 * L_x_signal_Region2_overallSyst_x_HistSyst ] = 1089/1089
RooGaussian::alpha_background_0_shape_0Constraint[ x=alpha_background_0_shape_0 mean=nom_alpha_background_0_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_1Constraint[ x=alpha_background_0_shape_1 mean=nom_alpha_background_0_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_0_shape_2Constraint[ x=alpha_background_0_shape_2 mean=nom_alpha_background_0_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_0Constraint[ x=alpha_background_1_shape_0 mean=nom_alpha_background_1_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_1Constraint[ x=alpha_background_1_shape_1 mean=nom_alpha_background_1_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_1_shape_2Constraint[ x=alpha_background_1_shape_2 mean=nom_alpha_background_1_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_0Constraint[ x=alpha_background_2_shape_0 mean=nom_alpha_background_2_shape_0 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_1Constraint[ x=alpha_background_2_shape_1 mean=nom_alpha_background_2_shape_1 sigma=1 ] = 1
RooGaussian::alpha_background_2_shape_2Constraint[ x=alpha_background_2_shape_2 mean=nom_alpha_background_2_shape_2 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_0Constraint[ x=alpha_background_norm_uncertainty_0 mean=nom_alpha_background_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_1Constraint[ x=alpha_background_norm_uncertainty_1 mean=nom_alpha_background_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_background_norm_uncertainty_2Constraint[ x=alpha_background_norm_uncertainty_2 mean=nom_alpha_background_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_0Constraint[ x=alpha_signal_norm_uncertainty_0 mean=nom_alpha_signal_norm_uncertainty_0 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_1Constraint[ x=alpha_signal_norm_uncertainty_1 mean=nom_alpha_signal_norm_uncertainty_1 sigma=1 ] = 1
RooGaussian::alpha_signal_norm_uncertainty_2Constraint[ x=alpha_signal_norm_uncertainty_2 mean=nom_alpha_signal_norm_uncertainty_2 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_0Constraint[ x=alpha_signal_shape_0 mean=nom_alpha_signal_shape_0 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_1Constraint[ x=alpha_signal_shape_1 mean=nom_alpha_signal_shape_1 sigma=1 ] = 1
RooGaussian::alpha_signal_shape_2Constraint[ x=alpha_signal_shape_2 mean=nom_alpha_signal_shape_2 sigma=1 ] = 1
RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.02 ] = 1
RooProdPdf::model_Region2[ lumiConstraint * alpha_background_norm_uncertainty_0Constraint * alpha_background_norm_uncertainty_1Constraint * alpha_background_norm_uncertainty_2Constraint * alpha_background_0_shape_0Constraint * alpha_background_0_shape_1Constraint * alpha_background_0_shape_2Constraint * alpha_background_1_shape_0Constraint * alpha_background_1_shape_1Constraint * alpha_background_1_shape_2Constraint * alpha_background_2_shape_0Constraint * alpha_background_2_shape_1Constraint * alpha_background_2_shape_2Constraint * alpha_signal_norm_uncertainty_0Constraint * alpha_signal_norm_uncertainty_1Constraint * alpha_signal_norm_uncertainty_2Constraint * alpha_signal_shape_0Constraint * alpha_signal_shape_1Constraint * alpha_signal_shape_2Constraint * Region2_model(obs_x_Region2) ] = 1089

functions
--------
RooProduct::L_x_model0_Region2_overallSyst_x_HistSyst[ 1 * model0_Region2_overallSyst_x_HistSyst ] = 90
RooProduct::L_x_model1_Region2_overallSyst_x_HistSyst[ 1 * model1_Region2_overallSyst_x_HistSyst ] = 90
RooProduct::L_x_model2_Region2_overallSyst_x_HistSyst[ 1 * model2_Region2_overallSyst_x_HistSyst ] = 900
RooProduct::L_x_signal_Region2_overallSyst_x_HistSyst[ 1 * signal_Region2_overallSyst_x_HistSyst ] = 9
PiecewiseInterpolation::model0_Region2_Hist_alpha[ ] = 90
RooHistFunc::model0_Region2_Hist_alpha_0high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model0_Region2_Hist_alpha_0low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model0_Region2_Hist_alpha_1high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model0_Region2_Hist_alpha_1low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model0_Region2_Hist_alpha_2high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model0_Region2_Hist_alpha_2low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model0_Region2_Hist_alphanominal[ depList=(obs_x_Region2) ] = 90
RooStats::HistFactory::FlexibleInterpVar::model0_Region2_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model0_Region2_overallSyst_x_HistSyst[ model0_Region2_Hist_alpha * model0_Region2_epsilon ] = 90
PiecewiseInterpolation::model1_Region2_Hist_alpha[ ] = 90
RooHistFunc::model1_Region2_Hist_alpha_0high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model1_Region2_Hist_alpha_0low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model1_Region2_Hist_alpha_1high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model1_Region2_Hist_alpha_1low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model1_Region2_Hist_alpha_2high[ depList=(obs_x_Region2) ] = 99
RooHistFunc::model1_Region2_Hist_alpha_2low[ depList=(obs_x_Region2) ] = 81
RooHistFunc::model1_Region2_Hist_alphanominal[ depList=(obs_x_Region2) ] = 90
RooStats::HistFactory::FlexibleInterpVar::model1_Region2_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model1_Region2_overallSyst_x_HistSyst[ model1_Region2_Hist_alpha * model1_Region2_epsilon ] = 90
PiecewiseInterpolation::model2_Region2_Hist_alpha[ ] = 900
RooHistFunc::model2_Region2_Hist_alpha_0high[ depList=(obs_x_Region2) ] = 990
RooHistFunc::model2_Region2_Hist_alpha_0low[ depList=(obs_x_Region2) ] = 810
RooHistFunc::model2_Region2_Hist_alpha_1high[ depList=(obs_x_Region2) ] = 990
RooHistFunc::model2_Region2_Hist_alpha_1low[ depList=(obs_x_Region2) ] = 810
RooHistFunc::model2_Region2_Hist_alpha_2high[ depList=(obs_x_Region2) ] = 990
RooHistFunc::model2_Region2_Hist_alpha_2low[ depList=(obs_x_Region2) ] = 810
RooHistFunc::model2_Region2_Hist_alphanominal[ depList=(obs_x_Region2) ] = 900
RooStats::HistFactory::FlexibleInterpVar::model2_Region2_epsilon[ paramList=(alpha_background_norm_uncertainty_0,alpha_background_norm_uncertainty_1,alpha_background_norm_uncertainty_2) ] = 1
RooProduct::model2_Region2_overallSyst_x_HistSyst[ model2_Region2_Hist_alpha * model2_Region2_epsilon ] = 900
PiecewiseInterpolation::signal_Region2_Hist_alpha[ ] = 9
RooHistFunc::signal_Region2_Hist_alpha_0high[ depList=(obs_x_Region2) ] = 9.9
RooHistFunc::signal_Region2_Hist_alpha_0low[ depList=(obs_x_Region2) ] = 8.1
RooHistFunc::signal_Region2_Hist_alpha_1high[ depList=(obs_x_Region2) ] = 9.9
RooHistFunc::signal_Region2_Hist_alpha_1low[ depList=(obs_x_Region2) ] = 8.1
RooHistFunc::signal_Region2_Hist_alpha_2high[ depList=(obs_x_Region2) ] = 9.9
RooHistFunc::signal_Region2_Hist_alpha_2low[ depList=(obs_x_Region2) ] = 8.1
RooHistFunc::signal_Region2_Hist_alphanominal[ depList=(obs_x_Region2) ] = 9
RooStats::HistFactory::FlexibleInterpVar::signal_Region2_epsilon[ paramList=(alpha_signal_norm_uncertainty_0,alpha_signal_norm_uncertainty_1,alpha_signal_norm_uncertainty_2) ] = 1
RooProduct::signal_Region2_overallNorm_x_sigma_epsilon[ SignalStrength * signal_Region2_epsilon ] = 1
RooProduct::signal_Region2_overallSyst_x_HistSyst[ signal_Region2_Hist_alpha * signal_Region2_overallNorm_x_sigma_epsilon ] = 9

datasets
--------
RooDataSet::asimovData(obs_x_Region2)
RooDataSet::obsData(obs_x_Region2)

embedded datasets (in pdfs and functions)
-----------------------------------------
RooDataHist::model0_Region2_Hist_alphanominalDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_0lowDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_1lowDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_2lowDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_0highDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_1highDHist(obs_x_Region2)
RooDataHist::model0_Region2_Hist_alpha_2highDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alphanominalDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_0lowDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_1lowDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_2lowDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_0highDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_1highDHist(obs_x_Region2)
RooDataHist::model1_Region2_Hist_alpha_2highDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alphanominalDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_0lowDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_1lowDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_2lowDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_0highDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_1highDHist(obs_x_Region2)
RooDataHist::model2_Region2_Hist_alpha_2highDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alphanominalDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_0lowDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_1lowDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_2lowDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_0highDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_1highDHist(obs_x_Region2)
RooDataHist::signal_Region2_Hist_alpha_2highDHist(obs_x_Region2)

named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
ModelConfig_Observables:(obs_x_Region2)
coefList:(binWidth_obs_x_Region2_0,binWidth_obs_x_Region2_1,binWidth_obs_x_Region2_2,binWidth_obs_x_Region2_3)
constraintTerms:(lumiConstraint,alpha_background_norm_uncertainty_0Constraint,alpha_background_norm_uncertainty_1Constraint,alpha_background_norm_uncertainty_2Constraint,alpha_background_0_shape_0Constraint,alpha_background_0_shape_1Constraint,alpha_background_0_shape_2Constraint,alpha_background_1_shape_0Constraint,alpha_background_1_shape_1Constraint,alpha_background_1_shape_2Constraint,alpha_background_2_shape_0Constraint,alpha_background_2_shape_1Constraint,alpha_background_2_shape_2Constraint,alpha_signal_norm_uncertainty_0Constraint,alpha_signal_norm_uncertainty_1Constraint,alpha_signal_norm_uncertainty_2Constraint,alpha_signal_shape_0Constraint,alpha_signal_shape_1Constraint,alpha_signal_shape_2Constraint)
globalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
likelihoodTerms:(Region2_model)
obsAndWeight:(weightVar,obs_x_Region2)
observables:(obs_x_Region2)
observablesSet:(obs_x_Region2)
shapeList:(L_x_model0_Region2_overallSyst_x_HistSyst,L_x_model1_Region2_overallSyst_x_HistSyst,L_x_model2_Region2_overallSyst_x_HistSyst,L_x_signal_Region2_overallSyst_x_HistSyst)

generic objects
---------------
RooStats::ModelConfig::ModelConfig

Setting Parameter(s) of Interest as: SignalStrength 
full list of observables:
RooArgList:: = (obs_x_Region0,obs_x_Region1,obs_x_Region2)


------------------
 Entering combination
-----------------------------------------
create toy data for Region0,Region1,Region2
RooDataSet::AsimovData0[obs_x_Region0,channelCat,weight:binWeightAsimov] = 1 entries (990 weighted)
RooDataSet::AsimovData1[obs_x_Region1,channelCat,weight:binWeightAsimov] = 1 entries (1089 weighted)
RooDataSet::AsimovData2[obs_x_Region2,channelCat,weight:binWeightAsimov] = 1 entries (1089 weighted)
Merging data for channel Region0
Merging data for channel Region1
Merging data for channel Region2

RooWorkspace(combined) combined contents

variables
---------
(channelCat,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2,obs_x_Region0,obs_x_Region1,obs_x_Region2,weightVar)

datasets
--------
RooDataSet::asimovData(obs_x_Region0,obs_x_Region1,obs_x_Region2,weightVar,channelCat)
RooDataSet::obsData(channelCat,obs_x_Region0,obs_x_Region1,obs_x_Region2)

named sets
----------
ModelConfig_GlobalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
ModelConfig_Observables:(obs_x_Region0,obs_x_Region1,obs_x_Region2,weightVar,channelCat)
globalObservables:(nom_alpha_background_norm_uncertainty_0,nom_alpha_background_norm_uncertainty_1,nom_alpha_background_norm_uncertainty_2,nom_alpha_background_0_shape_0,nom_alpha_background_0_shape_1,nom_alpha_background_0_shape_2,nom_alpha_background_1_shape_0,nom_alpha_background_1_shape_1,nom_alpha_background_1_shape_2,nom_alpha_background_2_shape_0,nom_alpha_background_2_shape_1,nom_alpha_background_2_shape_2,nom_alpha_signal_norm_uncertainty_0,nom_alpha_signal_norm_uncertainty_1,nom_alpha_signal_norm_uncertainty_2,nom_alpha_signal_shape_0,nom_alpha_signal_shape_1,nom_alpha_signal_shape_2)
observables:(obs_x_Region0,obs_x_Region1,obs_x_Region2,weightVar,channelCat)



----------------
 Importing combined model
setting Lumi constant
Setting Parameter(s) of Interest as: SignalStrength

ok that seemed to work





In [10]:



    


events = 1000
chans = 1
nps = 1
for bins in [1,10,20,30,40,50,60,70,80,90,100]:
    workspace = get_workspace(nchannels = chans, events = events, nbins = bins, nnps = nps)
    workspace.SetName('BinnedWorkspace')
    workspace.writeToFile("output/workspace{}channels{}events{}bins{}nps.root".format(chans, events, bins, nps))



    


















    






this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0










    






Warning in <TROOT::Append>: Replacing existing TH1: observed0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: above_expected0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: above_expected0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: above_expected0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: expected0sample0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: expected0sample0 (Potential memory leak).
Warning in <TROOT::Append>: Replacing existing TH1: expected0sample0 (Potential memory leak).

















In [11]:



    


events = 1000
chans = 1
nps = 1
bins = 10
for events in [1000,100000,10000000]:
    workspace = get_workspace(nchannels = chans, events = events, nbins = bins, nnps = nps)
    workspace.SetName('BinnedWorkspace')
    workspace.writeToFile("output/workspace{}channels{}events{}bins{}nps.root".format(chans, events, bins, nps))



    


















    






this uncertainty is 0
this uncertainty is 0
this uncertainty is 0

















In [16]:



    


events = 1000
chans = 1
nps = 1
bins = 10
for chans in range(1,10):
    workspace = get_workspace(nchannels = chans, events = events, nbins = bins, nnps = nps)
    workspace.SetName('BinnedWorkspace')
    workspace.writeToFile("output/workspace{}channels{}events{}bins{}nps.root".format(chans, events, bins, nps))



    


















    






this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0
this uncertainty is 0

















In [15]:



    


events = 1000
chans = 1
nps = 1
bins = 10
for nps in range(10):
    workspace = get_workspace(nchannels = chans, events = events, nbins = bins, nnps = nps)
    workspace.SetName('BinnedWorkspace')
    workspace.writeToFile("output/workspace{}channels{}events{}bins{}nps.root".format(chans, events, bins, nps))



    


















    






this uncertainty is 0
this uncertainty is 0
this uncertainty is 1
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 4
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 4
this uncertainty is 5
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 4
this uncertainty is 5
this uncertainty is 6
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 4
this uncertainty is 5
this uncertainty is 6
this uncertainty is 7
this uncertainty is 0
this uncertainty is 1
this uncertainty is 2
this uncertainty is 3
this uncertainty is 4
this uncertainty is 5
this uncertainty is 6
this uncertainty is 7
this uncertainty is 8

















In [ ]:

Content source: roofit-dev/parallel-roofit-scripts

Similar notebooks:

notebook.community | gallery | about