SingleTOFAnalysis.h

Go to the documentation of this file.

#ifndef _SINGLE_TOF_ANALYSIS_H_
#define _SINGLE_TOF_ANALYSIS_H_
 
 
#include "TTree.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TGraph.h"
#include "TMultiGraph.h"
 
#include <string>
#include <vector>
class SingleTOFAnalysis
{
public:
  SingleTOFAnalysis();
 
 
  SingleTOFAnalysis(double a_flightPath,
                    double a_timeCalib = 0
                    );
 
  int loadExpDataTrees(std::string a_fileName,int a_ch);
 
 
  int addResponseTree(std::string a_filename);
 
  TH1* generatePSDPlot();
  
 
  std::vector<TH1*> generateCalibrationAndPSDCutPlots(std::string a_psdBoundFileName);
 
  int loadPSDBounds(std::string a_psdBoundFileName);
 
  TH1* correctGammaPeak(TH1* a_gammaTimeHistory,
                        double a_startTime,
                        double a_stopTime
                        );
  TH3* createResponseMatrix(TH1* a_gammaTimeHistory
                           );
  TMultiGraph* getPSDBoundPlot();
 
 
  TH2F* producePH_V_TOF(TH1* a_spectrum,
                      TH3* a_responseMatrix
                      );
 
  //getters and setters
  int setTimeBinWidth(double a_timeBinWidth);
  int setTofMin(double a_tofMin);
  int setTofMax(double a_tofMax);
  int setPhRes(double a_phRes);
  int setPhMin(double a_phMin);
  int setPhMax(double a_phMax);
  //getters
  double getTimeBinWidth();
  double getTofMin();
  double getTofMax();
  double getPhRes();
  double getPhMin();
  double getPhMax();
 
 
private:
 
  //required for calculations
  double m_flightPath;
  double m_gammaFlightTime;
  double m_timeCalib;
  double m_tdcConversion;
  //histogram deffinitions
  double m_timeBinWidth;
  double m_tofMin;
  double m_tofMax;
  double m_phRes;
  double m_phMin;
  double m_phMax;
 
  /*
  tmpTreePtr->Branch("amp",&curEn);
  tmpTreePtr->Branch("shape",&curShape);
  tmpTreePtr->Branch("pileUpCliping",&piledOrClipped);
  tmpTreePtr->Branch("hitTime",&deltaTs);
  */
  TTree* m_dataTree;
  //variable for tree association
  int               m_amp;
  float             m_shape;
  bool              m_pileUpCliping;
  std::vector<int>* m_hitTime;
  //is the data tree present
  bool   m_hasExpData;
  bool   dataLoaded();
  TTree* m_respTree;
  //variables for tree association
  double m_ph;
  double m_en;
  double m_tof;
  //is the data tree present
  bool   m_hasSimData;
  bool   simLoaded();
  TH1* m_neutronFluxHistory;
 
  std::vector<TGraph*> m_psdBoundGraphs;
  int getPSDMetric(double a_amp,double a_shape);
 
 
  TH1* m_psdPlot;
 
 
};
 
#endif