MadChi2Calc.h

Go to the documentation of this file.

#ifndef madchi2calc_h
#define madchi2calc_h
#include <iostream>
#include "TH1.h"
#include "TH2.h"
#include "TMatrix.h"

class MadChi2Calc {

 public:
  
  MadChi2Calc(int ct=0);
  ~MadChi2Calc();
  
  void SetCalcType(int ct) {calcType = ct;}
  int GetCalcType() {return calcType;}

  double GetChi2(TH1F *,TH1F *,TH1F *,double);
  double GetChi2(TH2F *,TH2F *,TH2F *,double);
  double GetChi2(int,double*,double*,double*,double,double Min = 0,
                 double Max = 20);
  
  void SetFluxMap(TH2F *hist = 0);
  void SetSysErrors(double,double,double*,double);

  void SetNormalisationError(double err) {NormErr = err;}
  double GetNormalisationError() {return NormErr;}

  void SetBinToBinError(double err) {BinErr = err;}
  double GetBinToBinError() {return BinErr;}

  void SetBinHeightError(double *err) {HeightErr = err;}
  double *GetBinHeightError() {return HeightErr;}
  
  void SetBkdXSecError(double err) {BkdXSecErr = err;}
  double GetBkdXSecError() {return BkdXSecErr;}

 protected:

  double SumChi2(int,double*,double*,double*,double,double,double);
  double BasicChi2(double sig,double mc,double bkd=0);
  double PoisDistChi2(double sig,double mc,double bkd=0);
  double WithSysErrChi2(int,double*,double*,double*,double,double,double);
  double WithSysErrLoopsChi2(int,double *,double *,double *,double);

  int calcType; //determines which method to use to calculate chi2
                // 1 = BasicChi2; 2 = WithSysErrChi2; 3 = WithSysErrLoopsChi2;
                // default = PoisDistChi2

  TH2F *FluxMapHist; //No flux map used by default

  //bin to bin errors:
  double BinErr; //2% by default

  //Energy dependent bin height errors:
  double *HeightErr; //2% by default - assumes binning is same as histo binning

  //Overall normalisation error:
  double NormErr; //2% by default
  
  //Uncertainty in cross-section of Background:
  double BkdXSecErr; //5% by default

  // Combination of NormErr and BkdXSecErr allow uncertainties 
  // in signal and bkd cross-sections to be accounted for
  // Note!! Not all errors are included in every Chi2Calc method

};
#endif // #ifdef madchi2calc_h

---