NumuCCIncPionTemplateFitter.h

Go to the documentation of this file.

 
#pragma once

#include "CAFAna/Analysis/Plots.h"

//ROOT Includes
#include "Utilities/rootlogon.C"
#include "TH1.h"
#include "TH2.h"
#include "THStack.h"
#include "TCanvas.h"
#include "TLegend.h"
#include "TStyle.h"
#include "TFile.h"
#include <vector> 
#include <math.h>
#include "TLine.h"
#include "TFile.h"
#include "TTree.h"
#include "TCanvas.h"
#include "TF1.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TF2.h"
#include "TH2.h"
#include "TCutG.h"
#include "TMath.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TRandom.h"
#include "TGraph.h"
#include "TColor.h"
#include "TF3.h"
#include "TH3.h"
#include "TH3F.h"
#include "TMinuit.h"
#include "TLegend.h"
#include "THStack.h"
#include "TMatrix.h"
#include "TGraphErrors.h"
#include "TGraph.h"
#include "TDecompLU.h"
#include "TDecompSVD.h"
#include "TRandom.h"

//C++
#include <iostream>
#include <fstream>
#include <iomanip> 
#include <vector>
#include <cmath>
#include <math.h> 
#include <string>
#include <stdint.h>
#include <sstream>
#include <functional>

namespace ana
{

  class NumuCCIncPionTemplateFitter{
  public:
    NumuCCIncPionTemplateFitter(TH2F* da,
                                std::vector<TH2F*> templates,
                                std::vector<TH2F*> analysis,
                                TH2F* cv,
                                std::vector<TH2F*> systs_up,
                                std::vector<TH2F*> systs_down,
                                std::vector<TH2F*> systs_oneside,
                                std::vector<TH2F*> multiverse,
                                std::vector<TH2F*> ppfx,
                                int numbins);
    virtual ~NumuCCIncPionTemplateFitter() {};

    std::pair<double,double> signal_wei;
    std::pair<double,double> other_wei;
    std::pair<double,double> cc0pi_wei;
    double fit_covariance_matrix[3][3];
    double chisquared;

    void PrintValueHolders(int binx);
    void FCNMinuit();
    std::vector<TH2F*> doFullFit();
    //Function That Can Be used after the full fit is performed
    std::vector<TH1F*> getPIDFitTemplates(int binx);
    std::vector<std::pair<TH2F*,TH2F*>> getFitNormalizationAndErrors();
    //Functions That Can Be used at any point
    std::vector<TH1F*> getPIDTemplates(int binx);
    TH2F* getCovarianceMatrixTemplateBins(int binx);
    TH2F* getCorrelationMatrixTemplateBins(int binx);
    TH2F* getCovarianceBetweenSysts(int binx);
    TH2F* getCorrelationBetweenSysts(int binx);
    void FillCovarianceMatrixExtra(int binx);


  private:

    void FillValueHolders(int binx);
    bool CheckSystsSize();
    std::vector<TH2F*> GetAnalysisTemplates();
    std::string GetAnalysisBinCaption(int binx);
    bool CheckIfBinShouldBeUsed(int binx);
    std::vector<TH1F*> CalculateOneSigmaShift(int binx);
    void FillCovarianceMatrix(int binx);
    void FillStatisticCovarianceMatrix(int binx);
    double GetCovarianceMatrixValue(int binx, int biny);
    double GetStatisticCovarianceMatrixValue(int binx, int biny);
    static void fcn3Var(Int_t &npar, Double_t *gin, 
                    Double_t &f, Double_t *par, Int_t iflag);
    double myFunction3Vars(double numu_par,double nue_par,double nc_par);
    static void fcn2Var(Int_t &npar, Double_t *gin, 
                    Double_t &f, Double_t *par, Int_t iflag);
    double myFunction2Vars(double bkgd_par,double nue_par);
    void SetBinWeightsAndErrors2D(int binx);
    void PropagateFitUncertainty3D(int binx);
    std::vector<std::pair<double,double>> GetTemplateWeightsAndErrors(int binx);
    std::vector<TH1F*> GetTemplates(int binx,std::string name);
    TH2F* FillCovarianceBetweenSysts(int binx);


    TH2F* hData2D;
    std::vector<TH2F*> hTemplates2D;
    std::vector<TH2F*> hAnalysis2D;
    TH2F* hCV2D;
    std::vector<TH2F*> hSystsUp2D;
    std::vector<TH2F*> hSystsDown2D;
    std::vector<TH2F*> hSystsOneSided;
    std::vector<TH2F*> hSystsMultiverse;
    std::vector<TH2F*> hSystsPPFX;
    const int NumberTemplateBins;
    TH1F* hSignalWeights1D;
    TH1F* hNumuCCWeights1D;
    TH1F* hNCWeights1D;
    TH1F* hSignalError1D;
    TH1F* hNumuCCError1D;
    TH1F* hNCError1D;
    TH1F* hChiSquaredResults1D;

    
    double DataValues[25];
    double SignalLike_Values[25];
    double CC0Pi_Values[25];
    double Other_Values[25];
    double CovarianceMatrix[25][25];
    double StatCovarianceMatrix[25][25];
    //*/


  };//end class def
}