TTOFAna.h

Go to the documentation of this file.

#ifndef _TTOF_ANA_
#define _TTOF_ANA_
//project 
#include "PostProcBase.h"
#include "TTOFEvent.h"
#include "TTOFConfig.h"
#include "CutManager.h"
#include "CutBuilder.h"
#include <Eigen/Dense>
//root
#include "TH1.h"
#include "TH2.h"
//c++
#include <string>
#include <map>
#include <random>
class TTOFAna : public PostProcBase<TTOFEvent>
{
public:
  TTOFAna();
  TTOFAna(std::string a_configFile);
  void buildThetaEstimationPlots(std::string a_outFile = "detThetatHist.root",
                                 double a_startLight = 1.4,
                                 double a_stopLight = 14,
                                 double a_step = 0.8
                                 );
  void testTransform();
  void addGraphicalPIDCuts(string a_cutFile);
  vector<TH2F*> buildRotatedPIDPlot(int a_intP1,
                             int a_intP2,
                             int a_intP3,
                             int a_intP4,
                             double a_lightThreshold=2.0,
                             double a_maxLight=10.0,
                             double a_theta=-18.3);
  double calcFOM(TH2F* a_histX1, TH2F* a_histX2);
  void plotFOMVsEner(TH2F* a_hist);
  void getPIDSlice(double a_eLow, double a_eHigh);
  double fitPIDSlice(double a_ener,
                     double a_mean1=-.038,
                     double a_mean2=-.03,
                     double a_amp1=100, 
                     double a_amp2=100,
                     double a_sig1=.001,
                     double a_sig2=.001);
  void makeFOMVsLight();
 
  void procFileToCLYCTOTWFTree(std::string a_outFileName);
 
  std::vector<TH2F*> buildResponseFunction(double a_deltaT=4.);
 
  void testUnfoldingWithReg(TH2F* a_response,
                            TH1F* a_lightYield,
                            TH1F* a_guessFlux
                            );
 
  TH2F* m_pidVsLight;
  TH1D* m_currLightSlice;
  std::map<double, vector<double>> m_fomFits;
  std::map<double, Eigen::MatrixXf> m_fomFitCovs;
  std::map<string, bool> m_cutState;
private:
  TTOFConfig m_config;
  CutManager m_cutManager;
  bool passesCuts(TTOFEvent* a_event);
  //these are used to functionalize the templated base class
  void setTreeName();
  void setBranchName();
  //used to store histograms from various routines
  std::map<std::string , TH1*> m_hists;
  //used to control simulation smearing and histogram generation
  //set up random generator
  std::default_random_engine generator;
};
 
 
#endif