33 : _minBinToPrint(-1), _maxBinToPrint(-1),
34 _label_print_width(20), _bin_print_width(12) {
37 std::cout <<
"Error: The supplied ModelConfig is nullptr " << std::endl;
44 std::cout <<
"Error: The pdf found in the ModelConfig: " <<
mc->GetName()
45 <<
" is nullptr" << std::endl;
55 std::cout <<
"Error: Observable set in the ModelConfig: " <<
mc->GetName()
56 <<
" is nullptr" << std::endl;
61 <<
" found in ModelConfig: " <<
mc->GetName()
62 <<
" has no entries." << std::endl;
79 _minBinToPrint(-1), _maxBinToPrint(-1),
80 _label_print_width(20), _bin_print_width(12) {
83 auto file = std::make_unique<TFile>(FileName.c_str());
85 std::cout <<
"Error: Failed to open file: " << FileName << std::endl;
92 std::cout <<
"Error: Failed to get workspace: " <<
WorkspaceName
93 <<
" from file: " << FileName << std::endl;
102 <<
" in file: " << FileName << std::endl;
109 std::cout <<
"Error: The pdf found in the ModelConfig: " <<
ModelConfigName
110 <<
" is nullptr" << std::endl;
120 std::cout <<
"Error: Observable set in the ModelConfig: " <<
ModelConfigName
121 <<
" is nullptr" << std::endl;
127 <<
" in file: " << FileName
128 <<
" has no entries." << std::endl;
143 _minBinToPrint(-1), _maxBinToPrint(-1),
144 _label_print_width(20), _bin_print_width(12) {
148 std::cout <<
"Error: The supplied pdf is nullptr" << std::endl;
158 std::cout <<
"Error: Supplied Observable set is nullptr" << std::endl;
161 if( observables->
empty() ) {
162 std::cout <<
"Error: Observable list: " << observables->
GetName()
163 <<
" has no entries." << std::endl;
190 std::cout << std::endl;
198 std::map< std::string, RooAbsPdf* >::iterator
itr;
202 std::cout <<
"Warning: Could not find channel: " << channel
208 if( pdf ==
nullptr ) {
209 std::cout <<
"Warning: Pdf associated with channel: " << channel
210 <<
" is nullptr" << std::endl;
222 std::cout << std::endl << channel <<
":" << std::endl;
264 std::string
tmp_name = channel +
"_pretty_tmp";
349 std::cout << std::endl;
359 std::cout << std::endl;
362 std::cout << std::setw(30) <<
"Parameter";
363 std::cout << std::setw(15) <<
"Value"
364 << std::setw(15) <<
"Error Low"
365 << std::setw(15) <<
"Error High"
372 std::cout << std::setw(30) << param->GetName();
373 std::cout << std::setw(15) << param->getVal();
374 if( !param->isConstant() ) {
375 std::cout << std::setw(15) << param->getErrorLo() << std::setw(15) << param->getErrorHi();
377 std::cout<< std::endl;
379 std::cout << std::endl;
392 std::cout << std::endl;
395 std::cout << std::setw(30) <<
"Parameter";
396 std::cout << std::setw(15) <<
"Value"
397 << std::setw(15) <<
"Error Low"
398 << std::setw(15) <<
"Error High"
405 std::cout << std::setw(30) << param->GetName();
406 std::cout << std::setw(15) << param->getVal();
407 if( !param->isConstant() ) {
408 std::cout << std::setw(15) << param->getErrorLo() << std::setw(15) << param->getErrorHi();
410 std::cout<< std::endl;
412 std::cout << std::endl;
417 const std::string&
sample,
427 std::cout << std::endl;
430 std::cout << std::setw(30) <<
"Parameter";
431 std::cout << std::setw(15) <<
"Value"
432 << std::setw(15) <<
"Error Low"
433 << std::setw(15) <<
"Error High"
440 std::cout << std::setw(30) << param->GetName();
441 std::cout << std::setw(15) << param->getVal();
442 if( !param->isConstant() ) {
443 std::cout << std::setw(15) << param->getErrorLo() << std::setw(15) << param->getErrorHi();
445 std::cout<< std::endl;
447 std::cout << std::endl;
478 std::map< std::string, std::map< std::string, RooAbsReal*> >::iterator
channel_itr;
481 std::cout <<
"Error: Channel: " << channel <<
" not found in Navigation" << std::endl;
493 std::map< std::string, std::map< std::string, RooAbsReal*> >::iterator
channel_itr;
496 std::cout <<
"Error: Channel: " << channel <<
" not found in Navigation" << std::endl;
501 std::map< std::string, RooAbsReal*>::iterator
sample_itr;
504 std::cout <<
"Error: Sample: " <<
sample <<
" not found in Navigation" << std::endl;
516 std::map< std::string, RooArgSet*>::iterator
channel_itr;
519 std::cout <<
"Error: Channel: " << channel <<
" not found in Navigation" << std::endl;
585 else total_hist->SetName( (channel +
"_hist").c_str() );
608 const std::string&
name) {
615 for(
unsigned int i=0; i < samples.size(); ++i) {
629 const std::string&
name) {
647 return channel == item->GetName();
655 int dim = vars.
size();
660 std::cout <<
"Error: To Create Histogram from RooDataSet" << std::endl;
662 }
else if( dim==1 ) {
681 std::cout <<
"Error: To Create Histogram from RooDataSet, Dimension must be 1, 2, or 3" << std::endl;
682 std::cout <<
"Observables: " << std::endl;
699 if(
data!=
nullptr ) {
719 std::string ClassName = arg->
ClassName();
720 if( ClassName ==
"RooProduct" ) {
790 for (
auto *arg : *
comps) {
791 std::string ClassName = arg->
ClassName();
792 if( ClassName ==
"RooRealSumPdf" ) {
822 if(
SampleName.find(
"L_x_") != std::string::npos ) {
852 for (
auto arg : *
comps) {
853 std::string
ArgName = arg->GetName();
866 std::unique_ptr<RooArgSet> parameters{parent->
getParameters(&args)};
867 for (
auto *param : *parameters) {
868 std::string
ParamName = param->GetName();
899 if( param==
nullptr ) {
900 std::cout <<
"Error: Couldn't Find parameter: " <<
parameter <<
" in model."
906 if(
parameter.find(
"gamma_stat_") != std::string::npos ) {
913 if(
term==
nullptr ) {
914 std::cout <<
"Error: Couldn't Find constraint term for parameter: " <<
parameter
927 if( constraintTerm==
nullptr ) {
928 std::cout <<
"Error: Cannot get uncertainty because parameter: " <<
parameter
929 <<
" has no constraint term" << std::endl;
940 std::cout <<
"Error: Constraint type is an empty string."
941 <<
" This simply should not be." << std::endl;
954 if(
parameter.find(
"alpha_")!=std::string::npos ) {
957 else if(
parameter.find(
"gamma_stat_")!=std::string::npos ) {
964 std::cout <<
"Error: Failed to find the 'sigma' node: " <<
sigmaName
965 <<
" in the RooGaussian: " << constraintTerm->
GetName() << std::endl;
976 std::cout <<
"Error: Failed to find the nominal 'tau' node: " <<
tauName
977 <<
" for the RooPoisson: " << constraintTerm->
GetName() << std::endl;
984 std::cout <<
"Error: Encountered unknown constraint type for Stat Uncertainties: "
998 std::cout <<
"Error: Cannot replace node: " <<
ToReplace
999 <<
" because this node wasn't found in: " <<
fModel->
GetName()
1024 std::cout <<
"Replaced: " <<
ToReplace <<
" with: " << ReplaceWith->
GetName()
1025 <<
" in node: " << client->GetName() << std::endl;
1038 const std::string&
sample) {
1073 std::cout << std::endl;
1074 std::cout <<
"Channel: " << channel <<
" Sample: " <<
sample << std::endl;
1077 for(
unsigned int i=0; i <
num_bins; ++i) {
1080 std::stringstream
sstr;
1084 std::cout << std::endl;
1094 std::unique_ptr<TH1> hist;
1121 std::string
name ) {
1136 int dim = vars.
size();
1159 std::cout <<
"Error: To Create Histogram from RooAbsReal function, Dimension must be 1, 2, or 3" << std::endl;
1169 const RooArgSet* observables =
mc->GetObservables();
1188 std::cout << std::endl;
1191 std::cout << std::setw(30) <<
"Parameter";
1192 std::cout << std::setw(15) <<
"Value"
1193 << std::setw(15) <<
"Error Low"
1194 << std::setw(15) <<
"Error High"
1200 std::string
ParamName = param->GetName();
1209 std::cout <<
"Setting param: " <<
ParamName <<
" constant"
1210 <<
" (matches regex: " <<
regExpr <<
")" << std::endl;
1217 if( !arg )
return nullptr;
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void data
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t index
const_iterator begin() const
const_iterator end() const
Common abstract base class for objects that represent a value and a "shape" in RooFit.
RooFit::OwningPtr< RooArgSet > getParameters(const RooAbsData *data, bool stripDisconnected=true) const
Create a list of leaf nodes in the arg tree starting with ourself as top node that don't match any of...
bool redirectServers(const RooAbsCollection &newServerList, bool mustReplaceAll=false, bool nameChange=false, bool isRecursionStep=false)
Replace all direct servers of this object with the new servers in newServerList.
RooFit::OwningPtr< RooArgSet > getComponents() const
Create a RooArgSet with all components (branch nodes) of the expression tree headed by this object.
bool getAttribute(const Text_t *name) const
Check if a named attribute is set. By default, all attributes are unset.
void setAttribute(const Text_t *name, bool value=true)
Set (default) or clear a named boolean attribute of this object.
RooAbsArg * findServer(const char *name) const
Return server of this with name name. Returns nullptr if not found.
const char * GetName() const override
Returns name of object.
virtual bool add(const RooAbsArg &var, bool silent=false)
Add the specified argument to list.
Storage_t::size_type size() const
void Print(Option_t *options=nullptr) const override
This method must be overridden when a class wants to print itself.
Abstract interface for all probability density functions.
Abstract base class for objects that represent a real value and implements functionality common to al...
TH1 * createHistogram(RooStringView varNameList, Int_t xbins=0, Int_t ybins=0, Int_t zbins=0) const
Create and fill a ROOT histogram TH1, TH2 or TH3 with the values of this function for the variables w...
RooArgList is a container object that can hold multiple RooAbsArg objects.
RooAbsArg * at(Int_t idx) const
Return object at given index, or nullptr if index is out of range.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Object to represent discrete states.
Container class to hold unbinned data.
static RooMsgService & instance()
Return reference to singleton instance.
Represents the product of a given set of RooAbsReal objects.
Implements a PDF constructed from a sum of functions:
const RooArgList & funcList() const
Variable that can be changed from the outside.
Facilitates simultaneous fitting of multiple PDFs to subsets of a given dataset.
void PrintParameters(bool IncludeConstantParams=false)
Print the current values and errors of pdf parameters.
std::vector< std::string > fChannelNameVec
The list of channels.
RooArgSet * GetObservableSet(const std::string &channel)
Get the set of observables for a given channel.
TH1 * GetDataHist(RooDataSet *data, const std::string &channel, const std::string &name="")
Get a histogram from the dataset for this channel.
RooArgSet * fObservables
The observables.
void _GetNodes(ModelConfig *mc)
Fetch the node information for the pdf in question, and save it in the various collections in this cl...
std::map< std::string, RooAbsPdf * > fChannelSumNodeMap
Map of channel names to pdf without constraint.
void ReplaceNode(const std::string &ToReplace, RooAbsArg *ReplaceWith)
Find a node in the pdf and replace it with a new node These nodes can be functions,...
void PrintState()
Should pretty print all channels and the current values

void PrintSampleComponents(const std::string &channel, const std::string &sample)
Print the different components that make up a sample (NormFactors, Statistical Uncertainties,...
void SetPrintWidths(const std::string &channel)
Set the title and bin widths.
RooAbsPdf * fModel
The HistFactory Pdf Pointer.
TH1 * GetChannelHist(const std::string &channel, const std::string &name="")
Get the total channel histogram for this channel.
TH1 * GetSampleHist(const std::string &channel, const std::string &sample, const std::string &name="")
The (current) histogram for that sample This includes all parameters and interpolation.
TH1 * MakeHistFromRooFunction(RooAbsReal *func, RooArgList vars, std::string name="Hist")
Make a histogram from a function Edit so it can take a RooArgSet of parameters.
void PrintMultiDimHist(TH1 *hist, int bin_print_width)
Print a histogram's contents to the screen void PrettyPrintHistogram(TH1* hist);.
void PrintSampleParameters(const std::string &channel, const std::string &sample, bool IncludeConstantParams=false)
Print parameters that effect a particular sample.
std::map< std::string, RooAbsReal * > GetSampleFunctionMap(const std::string &channel)
Get a map of sample names to their functions for a particular channel.
RooAbsReal * SampleFunction(const std::string &channel, const std::string &sample)
Get the RooAbsReal function for a given sample in a given channel.
HistFactoryNavigation(ModelConfig *mc)
Initialize based on an already-created HistFactory Model.
double GetBinValue(int bin, const std::string &channel)
The value of the ith bin for the total in that channel.
RooAbsArg * findChild(const std::string &name, RooAbsReal *parent) const
Internal method implementation of finding a daughter node from a parent node (looping over all genera...
std::map< std::string, std::map< std::string, RooAbsReal * > > fChannelSampleFunctionMap
Map of Map of Channel, Sample names to Function Nodes Used by doing: fChannelSampleFunctionMap["MyCha...
std::map< std::string, RooAbsPdf * > fChannelPdfMap
Map of channel names to their full pdf's.
std::vector< std::string > GetChannelSampleList(const std::string &channel)
void PrintModelAndData(RooDataSet *data)
Print the model and the data, comparing channel by channel.
std::map< std::string, RooArgSet * > fChannelObservMap
Map of channel names to their set of ovservables.
RooAbsPdf * GetChannelPdf(const std::string &channel)
void PrintDataSet(RooDataSet *data, const std::string &channel="")
Print a "HistFactory style" RooDataSet in a readable way.
void PrintChannelParameters(const std::string &channel, bool IncludeConstantParams=false)
Print parameters that effect a particular channel.
void SetConstant(const std::string ®Expr=".*", bool constant=true)
RooArgSet _GetAllProducts(RooProduct *node)
Recursively get all products of products.
THStack * GetChannelStack(const std::string &channel, const std::string &name="")
Get a stack of all samples in a channel.
RooAbsReal * GetConstraintTerm(const std::string ¶meter)
Get the constraint term for a given systematic (alpha or gamma)
void DrawChannel(const std::string &channel, RooDataSet *data=nullptr)
Draw a stack of the channel, and include data if the pointer is supplied.
RooRealVar * var(const std::string &varName) const
double GetConstraintUncertainty(const std::string ¶meter)
Get the uncertainty based on the constraint term for a given systematic.
< A class that holds configuration information for a model using a workspace as a store
Persistable container for RooFit projects.
TObject * obj(RooStringView name) const
Return any type of object (RooAbsArg, RooAbsData or generic object) with given name)
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
TH1 is the base class of all histogram classes in ROOT.
virtual Int_t GetNbinsY() const
virtual Int_t GetNbinsZ() const
virtual Int_t GetNbinsX() const
Bool_t IsBinUnderflow(Int_t bin, Int_t axis=0) const
Return true if the bin is underflow.
Bool_t IsBinOverflow(Int_t bin, Int_t axis=0) const
Return true if the bin is overflow.
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
The Histogram stack class.
virtual void Add(TH1 *h, Option_t *option="")
Add a new histogram to the list.
void Draw(Option_t *chopt="") override
Draw this stack with its current attributes.
const char * GetName() const override
Returns name of object.
virtual const char * ClassName() const
Returns name of class to which the object belongs.
Regular expression class.
RooCmdArg YVar(const RooAbsRealLValue &var, const RooCmdArg &arg={})
RooCmdArg Scaling(bool flag)
RooCmdArg Binning(const RooAbsBinning &binning)
MsgLevel
Verbosity level for RooMsgService::StreamConfig in RooMsgService.
Namespace for the RooStats classes.