Uploaded image for project: 'ROOT'
  1. ROOT
  2. ROOT-5311

HistFactory output is too verbose and cannot be controlled

    XMLWordPrintable

    Details

    • Type: Improvement
    • Status: Closed (View Workflow)
    • Priority: Medium
    • Resolution: Fixed
    • Affects Version/s: 5.34/00
    • Fix Version/s: 6.22/00
    • Component/s: RooFit/RooStats
    • Labels:
      None
    • Environment:

      N/A

      Description

      As the subject states, I think the default output of HistFactory's HistoToWorkspaceFactoryFast (etc..) is rather verbose and prints out loads of junk I don't need to see. I would like an option to turn this off. Yes, some of this output is fine when debugging but I don't want to see it fill up my terminal or log files when I don't need it. Can a verbosity flag be implemented? Even better would be to send all output through ROOT's logging system at the appropriate level (most of what it prints now would be better off at the DEBUG level) instead of printing on stdout directly. Also the poor formatting of the output could be improved. Looks rather hastily slapped together at the moment.

      Sample output from a fake simple analysis setup below. Actual analysis configs print out MUCH more than this (too long to paste here).

      -------------------
      Starting to process VBF channel with 1 observables
      lumi str = [1,0,10]
      lumi Error str = nominalLumi[1,0,2],0.1
      making normFactor: x
      making normFactor: y
      making normFactor: z
      Gaussian::alpha_xConstraint(alpha_x,nom_alpha_x[0.,-10,10],1.)
      Gaussian::alpha_yConstraint(alpha_y,nom_alpha_y[0.,-10,10],1.)
      Gaussian::alpha_zConstraint(alpha_z,nom_alpha_z[0.,-10,10],1.)
      -----------------------------------------
      import model into workspace
      WARNING::Detected a bin with zero expected events- skip it
      RooDataSet::AsimovData[obs_x_VBF,weight:binWeightAsimov] = 9 entries (10000 weighted)

      RooWorkspace(VBF) VBF workspace contents

      variables
      ---------
      (Lumi,alpha_x,alpha_y,alpha_z,binWidth_obs_x_VBF_0,nom_alpha_x,nom_alpha_y,nom_alpha_z,nominalLumi,obs_x_VBF,weightVar,x,y,z)

      p.d.f.s
      -------
      RooRealSumPdf::VBF_model[ binWidth_obs_x_VBF_0 * L_x_QCD_VBF_overallSyst_x_HistSyst ] = 1/10000
      RooGaussian::alpha_xConstraint[ x=alpha_x mean=nom_alpha_x sigma=1 ] = 1
      RooGaussian::alpha_yConstraint[ x=alpha_y mean=nom_alpha_y sigma=1 ] = 1
      RooGaussian::alpha_zConstraint[ x=alpha_z mean=nom_alpha_z sigma=1 ] = 1
      RooGaussian::lumiConstraint[ x=Lumi mean=nominalLumi sigma=0.1 ] = 1
      RooProdPdf::model_VBF[ lumiConstraint * alpha_xConstraint * alpha_yConstraint * alpha_zConstraint * VBF_model(obs_x_VBF) ] = 1

      functions
      --------
      RooProduct::L_x_QCD_VBF_overallSyst_x_HistSyst[ Lumi * QCD_VBF_overallSyst_x_HistSyst ] = 1
      PiecewiseInterpolation::QCD_VBF_Hist_alpha[ ] = 1
      RooHistFunc::QCD_VBF_Hist_alpha_0high[ depList=(obs_x_VBF) ] = 0
      RooHistFunc::QCD_VBF_Hist_alpha_0low[ depList=(obs_x_VBF) ] = 1
      RooHistFunc::QCD_VBF_Hist_alpha_1high[ depList=(obs_x_VBF) ] = 0
      RooHistFunc::QCD_VBF_Hist_alpha_1low[ depList=(obs_x_VBF) ] = 0
      RooHistFunc::QCD_VBF_Hist_alpha_2high[ depList=(obs_x_VBF) ] = 0
      RooHistFunc::QCD_VBF_Hist_alpha_2low[ depList=(obs_x_VBF) ] = 0
      RooHistFunc::QCD_VBF_Hist_alphanominal[ depList=(obs_x_VBF) ] = 1
      RooProduct::QCD_VBF_overallNorm_x_sigma_epsilon[ x * y * z * QCD_VBF_epsilon ] = 1
      RooProduct::QCD_VBF_overallSyst_x_HistSyst[ QCD_VBF_Hist_alpha * QCD_VBF_overallNorm_x_sigma_epsilon ] = 1

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

      named sets
      ----------
      ModelConfig_GlobalObservables:(nominalLumi,nom_alpha_x,nom_alpha_y,nom_alpha_z)
      ModelConfig_Observables:(obs_x_VBF)
      coefList:(binWidth_obs_x_VBF_0)
      constraintTerms:(lumiConstraint,alpha_xConstraint,alpha_yConstraint,alpha_zConstraint)
      globalObservables:(nominalLumi,nom_alpha_x,nom_alpha_y,nom_alpha_z)
      likelihoodTerms:(VBF_model)
      obsAndWeight:(weightVar,obs_x_VBF)
      observables:(obs_x_VBF)
      observablesSet:(obs_x_VBF)
      shapeList:(L_x_QCD_VBF_overallSyst_x_HistSyst)

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

      Warining: No Parametetrs of interest are set
      Setting Parameter(s) of Interest as:
      full list of observables:
      RooArgList:: = (obs_x_VBF)

      ------------------
      Entering combination
      -----------------------------------------
      create toy data for VBF
      WARNING::Detected a bin with zero expected events- skip it
      RooDataSet::AsimovData0[obs_x_VBF,channelCat,weight:binWeightAsimov] = 9 entries (10000 weighted)
      Merging data for channel VBF

      RooWorkspace(combined) combined contents

      variables
      ---------
      (channelCat,nom_alpha_x,nom_alpha_y,nom_alpha_z,nominalLumi,obs_x_VBF,weightVar)

      datasets
      --------
      RooDataSet::asimovData(obs_x_VBF,weightVar,channelCat)
      RooDataSet::obsData(channelCat,obs_x_VBF)

      named sets
      ----------
      ModelConfig_GlobalObservables:(nominalLumi,nom_alpha_x,nom_alpha_y,nom_alpha_z)
      ModelConfig_Observables:(obs_x_VBF,weightVar,channelCat)
      globalObservables:(nominalLumi,nom_alpha_x,nom_alpha_y,nom_alpha_z)
      observables:(obs_x_VBF,weightVar,channelCat)

      ----------------
      Importing combined model
      Warining: No Parametetrs of interest are set
      Setting Parameter(s) of Interest as:

        Attachments

          Activity

            People

            • Assignee:
              shageboe Stephan Hageboeck
              Reporter:
              end Noel Dawe (Inactive)
            • Votes:
              0 Vote for this issue
              Watchers:
              1 Start watching this issue

              Dates

              • Created:
                Updated:
                Resolved:
                PlannedEnd:
                PlannedStart: