MultiHistoChiSquare.h

Go to the documentation of this file.

#ifndef PhysicsTools_Utilities_MultiHistoChiSquare_h
#define PhysicsTools_Utilities_MultiHistoChiSquare_h
#include "PhysicsTools/Utilities/interface/RootMinuitResultPrinter.h"
#include "TMath.h"
#include "TH1.h"

namespace fit {
  namespace helper {
    struct MultiHistoChiSquareNoArg {};
  }  // namespace helper

  template <typename T1,
            typename T2 = helper::MultiHistoChiSquareNoArg,
            typename T3 = helper::MultiHistoChiSquareNoArg,
            typename T4 = helper::MultiHistoChiSquareNoArg,
            typename T5 = helper::MultiHistoChiSquareNoArg,
            typename T6 = helper::MultiHistoChiSquareNoArg>
  class MultiHistoChiSquare {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1, typename TT2, typename TT3, typename TT4, typename TT5, typename TT6>
    MultiHistoChiSquare(TT1 &t1,
                        TH1 *histo1,
                        TT2 &t2,
                        TH1 *histo2,
                        TT3 &t3,
                        TH1 *histo3,
                        TT4 &t4,
                        TH1 *histo4,
                        TT5 &t5,
                        TH1 *histo5,
                        TT6 &t6,
                        TH1 *histo6,
                        double rangeMin,
                        double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax),
          chi2_(t2, histo2, rangeMin, rangeMax),
          chi3_(t3, histo3, rangeMin, rangeMax),
          chi4_(t4, histo4, rangeMin, rangeMax),
          chi5_(t5, histo5, rangeMin, rangeMax),
          chi6_(t6, histo6, rangeMin, rangeMax) {}
    double operator()() const {
      double chi2 = chi1_() + chi2_() + chi3_() + chi4_() + chi5_() + chi6_();
      static size_t count = 0;
      ++count;
      if (count % 10 == 0)
        return chi2;
    }
    void setHistos(TH1 *histo1, TH1 *histo2, TH1 *histo3, TH1 *histo4, TH1 *histo5, TH1 *histo6) {
      chi1_.setHistos(histo1);
      chi2_.setHistos(histo2);
      chi3_.setHistos(histo3);
      chi4_.setHistos(histo4);
      chi5_.setHistos(histo5);
      chi6_.setHistos(histo6);
    }
    size_t numberOfBins() const {
      return chi1_.numberOfBins() + chi2_.numberOfBins() + chi3_.numberOfBins() + chi4_.numberOfBins() +
             chi5_.numberOfBins() + chi6_.numberOfBins();
    }
    T1 &function1() { return chi1_.function(); }
    const T1 &function1() const { return chi1_.function(); }
    T2 &function2() { return chi2_.function(); }
    const T2 &function2() const { return chi2_.function(); }
    T3 &function3() { return chi3_.function(); }
    const T3 &function3() const { return chi3_.function(); }
    T4 &function4() { return chi4_.function(); }
    const T4 &function4() const { return chi4_.function(); }
    T5 &function5() { return chi5_.function(); }
    const T5 &function5() const { return chi5_.function(); }
    T6 &function6() { return chi6_.function(); }
    const T6 &function6() const { return chi6_.function(); }

  private:
    T1 chi1_;
    T2 chi2_;
    T3 chi3_;
    T4 chi4_;
    T5 chi5_;
    T6 chi6_;
  };

  template <typename T1, typename T2, typename T3, typename T4, typename T5>
  class MultiHistoChiSquare<T1, T2, T3, T4, T5, helper::MultiHistoChiSquareNoArg> {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1, typename TT2, typename TT3, typename TT4, typename TT5>
    MultiHistoChiSquare(TT1 &t1,
                        TH1 *histo1,
                        TT2 &t2,
                        TH1 *histo2,
                        TT3 &t3,
                        TH1 *histo3,
                        TT4 &t4,
                        TH1 *histo4,
                        TT5 &t5,
                        TH1 *histo5,
                        double rangeMin,
                        double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax),
          chi2_(t2, histo2, rangeMin, rangeMax),
          chi3_(t3, histo3, rangeMin, rangeMax),
          chi4_(t4, histo4, rangeMin, rangeMax),
          chi5_(t5, histo5, rangeMin, rangeMax) {}
    double operator()() const {
      double chi2 = chi1_() + chi2_() + chi3_() + chi4_() + chi5_();
      static size_t count = 0;
      ++count;
      return chi2;
    }
    void setHistos(TH1 *histo1, TH1 *histo2, TH1 *histo3, TH1 *histo4, TH1 *histo5) {
      chi1_.setHistos(histo1);
      chi2_.setHistos(histo2);
      chi3_.setHistos(histo3);
      chi4_.setHistos(histo4);
      chi5_.setHistos(histo5);
    }
    size_t numberOfBins() const {
      return chi1_.numberOfBins() + chi2_.numberOfBins() + chi3_.numberOfBins() + chi4_.numberOfBins() +
             chi5_.numberOfBins();
    }
    T1 &function1() { return chi1_.function(); }
    const T1 &function1() const { return chi1_.function(); }
    T2 &function2() { return chi2_.function(); }
    const T2 &function2() const { return chi2_.function(); }
    T3 &function3() { return chi3_.function(); }
    const T3 &function3() const { return chi3_.function(); }
    T4 &function4() { return chi4_.function(); }
    const T4 &function4() const { return chi4_.function(); }
    T5 &function5() { return chi5_.function(); }
    const T5 &function5() const { return chi5_.function(); }

  private:
    T1 chi1_;
    T2 chi2_;
    T3 chi3_;
    T4 chi4_;
    T5 chi5_;
  };

  template <typename T1, typename T2, typename T3, typename T4>
  class MultiHistoChiSquare<T1, T2, T3, T4, helper::MultiHistoChiSquareNoArg, helper::MultiHistoChiSquareNoArg> {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1, typename TT2, typename TT3, typename TT4>
    MultiHistoChiSquare(TT1 &t1,
                        TH1 *histo1,
                        TT2 &t2,
                        TH1 *histo2,
                        TT3 &t3,
                        TH1 *histo3,
                        TT4 &t4,
                        TH1 *histo4,
                        double rangeMin,
                        double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax),
          chi2_(t2, histo2, rangeMin, rangeMax),
          chi3_(t3, histo3, rangeMin, rangeMax),
          chi4_(t4, histo4, rangeMin, rangeMax) {}
    double operator()() const {
      double chi2 = chi1_() + chi2_() + chi3_() + chi4_();
      static size_t count = 0;
      ++count;
      return chi2;
    }
    void setHistos(TH1 *histo1, TH1 *histo2, TH1 *histo3, TH1 *histo4) {
      chi1_.setHistos(histo1);
      chi2_.setHistos(histo2);
      chi3_.setHistos(histo3);
      chi4_.setHistos(histo4);
    }
    size_t numberOfBins() const {
      return chi1_.numberOfBins() + chi2_.numberOfBins() + chi3_.numberOfBins() + chi4_.numberOfBins();
    }
    T1 &function1() { return chi1_.function(); }
    const T1 &function1() const { return chi1_.function(); }
    T2 &function2() { return chi2_.function(); }
    const T2 &function2() const { return chi2_.function(); }
    T3 &function3() { return chi3_.function(); }
    const T3 &function3() const { return chi3_.function(); }
    T4 &function4() { return chi4_.function(); }
    const T4 &function4() const { return chi4_.function(); }

  private:
    T1 chi1_;
    T2 chi2_;
    T3 chi3_;
    T4 chi4_;
  };

  template <typename T1, typename T2, typename T3>
  class MultiHistoChiSquare<T1,
                            T2,
                            T3,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg> {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1, typename TT2, typename TT3>
    MultiHistoChiSquare(
        TT1 &t1, TH1 *histo1, TT2 &t2, TH1 *histo2, TT3 &t3, TH1 *histo3, double rangeMin, double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax),
          chi2_(t2, histo2, rangeMin, rangeMax),
          chi3_(t3, histo3, rangeMin, rangeMax) {}
    double operator()() const { return chi1_() + chi2_() + chi3_(); }
    void setHistos(TH1 *histo1, TH1 *histo2, TH1 *histo3) {
      chi1_.setHistos(histo1);
      chi2_.setHistos(histo2);
      chi3_.setHistos(histo3);
    }
    size_t numberOfBins() const { return chi1_.numberOfBins() + chi2_.numberOfBins() + chi3_.numberOfBins(); }
    T1 &function1() { return chi1_.function(); }
    const T1 &function1() const { return chi1_.function(); }
    T2 &function2() { return chi2_.function(); }
    const T2 &function2() const { return chi2_.function(); }
    T3 &function3() { return chi3_.function(); }
    const T3 &function3() const { return chi3_.function(); }

  private:
    T1 chi1_;
    T2 chi2_;
    T3 chi3_;
  };

  template <typename T1, typename T2>
  class MultiHistoChiSquare<T1,
                            T2,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg> {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1, typename TT2>
    MultiHistoChiSquare(TT1 &t1, TH1 *histo1, TT2 &t2, TH1 *histo2, double rangeMin, double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax), chi2_(t2, histo2, rangeMin, rangeMax) {}
    double operator()() const { return chi1_() + chi2_(); }
    void setHistos(TH1 *histo1, TH1 *histo2) {
      chi1_.setHistos(histo1);
      chi2_.setHistos(histo2);
    }
    size_t numberOfBins() const { return chi1_.numberOfBins() + chi2_.numberOfBins(); }

  private:
    T1 chi1_;
    T2 chi2_;
  };

  template <typename T1>
  class MultiHistoChiSquare<T1,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg,
                            helper::MultiHistoChiSquareNoArg> {
  public:
    MultiHistoChiSquare() {}
    template <typename TT1>
    MultiHistoChiSquare(TT1 &t1, TH1 *histo1, double rangeMin, double rangeMax)
        : chi1_(t1, histo1, rangeMin, rangeMax) {}
    double operator()() const { return chi1_(); }
    void setHistos(TH1 *histo1) { chi1_.setHistos(histo1); }
    size_t numberOfBins() const { return chi1_.numberOfBins(); }

  private:
    T1 chi1_;
  };

  template <typename T1, typename T2, typename T3, typename T4, typename T5, typename T6>
  struct RootMinuitResultPrinter<MultiHistoChiSquare<T1, T2, T3, T4, T5, T6> > {
    static void print(double amin,
                      unsigned int numberOfFreeParameters,
                      const MultiHistoChiSquare<T1, T2, T3, T4, T5, T6> &f) {
      unsigned int ndof = f.numberOfBins() - numberOfFreeParameters;
      std::cout << "chi-squared/n.d.o.f. = " << amin << "/" << ndof << " = " << amin / ndof
                << "; prob: " << TMath::Prob(amin, ndof) << std::endl;
    }
  };
}  // namespace fit

#endif