Process.h

Go to the documentation of this file.

#ifndef CombineTools_Process_h
#define CombineTools_Process_h
#include <memory>
#include <string>
#include "TH1.h"
#include "RooAbsPdf.h"
#include "RooAbsReal.h"
#include "RooAbsData.h"
#include "RooRealVar.h"
#include "CombineHarvester/CombineTools/interface/MakeUnique.h"
#include "CombineHarvester/CombineTools/interface/Object.h"
 
namespace ch {
 
class Process : public Object {
 public:
  Process();
  ~Process();
  Process(Process const& other);
  Process(Process&& other);
  Process& operator=(Process other);
 
 
  void set_rate(double const& rate) { rate_ = rate; }
  double rate() const {
    double base = 1.;
    if (pdf_ && !dynamic_cast<RooAbsPdf*>(pdf_) && cached_obs_) {
      if (!cached_int_) {
        cached_int_ = pdf_->createIntegral(RooArgSet(*cached_obs_));
      }
      base = cached_int_->getVal();
    }
    return norm_ ? base * norm_->getVal() * rate_ : base * rate_;
  }
 
  double no_norm_rate() const { return rate_; }
 
  void set_shape(std::unique_ptr<TH1> shape, bool set_rate);
 
  void set_shape(TH1 const& shape, bool set_rate);
 
  TH1 const* shape() const { return shape_.get(); }
 
  std::unique_ptr<TH1> ClonedShape() const;
  std::unique_ptr<TH1> ClonedScaledShape() const;
 
  TH1F ShapeAsTH1F() const;
 
  void set_pdf(RooAbsReal* pdf) { pdf_ = pdf; }
  RooAbsReal const* pdf() const { return pdf_; }
 
  void set_data(RooAbsData* data) { data_ = data; }
  RooAbsData const* data() const { return data_; }
 
  void set_norm(RooAbsReal* norm) { norm_ = norm; }
  RooAbsReal const* norm() const { return norm_; }
 
  void set_observable(RooRealVar* obs) { cached_obs_ = obs; }
  RooRealVar * observable() const { return cached_obs_; }
 
  std::string to_string() const;
  friend std::ostream& operator<< (std::ostream &out, Process const& val);
  static std::ostream& PrintHeader(std::ostream &out);
 
 private:
  double rate_;
  std::unique_ptr<TH1> shape_;
  RooAbsReal* pdf_;
  RooAbsData* data_;
  RooAbsReal* norm_;
  RooRealVar* cached_obs_;
  mutable RooAbsReal* cached_int_;
 
  friend void swap(Process& first, Process& second);
};
}
 
#endif