HistFactory output is too verbose and cannot be controlled

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: