GenesisAnalysis.h

Go to the documentation of this file.

#include "MesyDAQPostProcessing.h"
#include "GENESISConfig.h"
//project includes
#include "MDPP16DataStructures.h"
#include "GenesisEventStructures.h"
#include "MesyDAQEvent.h"
#include "ConfigClasses.h"
#include "ConfigBuilder.h"
#include "PostProcBase.h"
#include "NSDPhysicsCalcs.h"
//root includes
#include "TChain.h"
#include "TFile.h"
#include "TH1.h"
//c++ includes
#include <vector>
#include <string>
#include <deque>
#include <utility>
#include <algorithm>
using std::vector;
using std::string;
using std::deque;
using std::sort;
using std::pair;
 
class GenesisAnalysis : public PostProcBase<GenesisEvent> 
{
public:
  GenesisAnalysis(string a_configFileName);
  void readConfigFile(string a_configFileName);
  void accumulateBasicQuantities(string a_outFileName 
                                      = "genEventBasicQuantities.root");
 
  //this function builds a tree of information for the scintillator detectors
  //that applies the timing calibrations before producing time differences
  //with the rf. It should be used in conjunction with an expectation of the 
  //arrival times of the beam produced gammas to develop a RF time offset.
  //this requires that your config file is complete with timing and gain 
  //calibrations for the scintillators.
  //the resultant tree contains the following variables
  // int detID
  // uint32_t wallTime
  // int isNeutron
  // double amp
  // double shape
  // std::vector<double> rfDiffs
  void buildRFCalibrationTree(string a_outFileName 
                                      = "rfCalibrationResults.root",
                              double a_dsFraction = 0.1
                              );
  void buildGammaSpectra(string a_outFileName = "spec1.root", 
                          bool a_requireScint=false,
                          bool a_integralCal = false,
                          int a_maxE = 8192
                          );
  void buildTOFHists(string a_outFileName = "tofHists.root");
  void buildPSDFHists(string a_outFileName = "psdHists.root");
  void findGammaGammaTimeOffset(
      string a_outFileNameRoot = "gammaTimeOffset.root",
      string a_outFileNameTxt = "gammaTimeOffset.txt",
      int a_doFit=1,
      std::array<double,3> a_offset = {0,0,0}
      );
  void buildOutgoingTOFHist(string a_outFileName = "outgoingTOF.root");
  void reduceToCalibratedTrees(string a_outFileName = "calibrateEvents.root",
                               int a_detOption = 0);
  void reduceToCloverTree(string a_outFileName = "calibratedGamma.root");
  void buildEfficiencyAnaTree(bool a_buildForNeutrons,
                              string a_outFileName="",
                              double a_enThreshold=.2);
  void buildEfficiencyAnaTreeScintScint(string a_outFileName,
                              double a_enThreshold=.2,
                              bool a_gammaGamma=false,
                              std::array<double,3> a_offset = {0,0,0});
  void makeRandomTOFHists(string a_outFileName="",
                              double a_enThreshold=.2);
  void  buildGammaGammaCoins(string a_fileName="gammaGammaCoinHists.root");
  void buildNeutronSinglesTree(string a_fileName="neutronSinglesTree.root");
  virtual std::ostream& printEvent( std::ostream& a_stream = std::cout,
                                      int a_eventToPrint = -1
                                     );
private:
  GenesisAnalysis();
  GenesisConfig m_conf;
  //overload of pure virtual functions to enable inheritance functionality
  void setTreeName();
  void setBranchName();
  vector<MDPPScintEvent> m_scintEvents;
  map<int, MDPPCloverEvent> m_cloverEvents;
  vector<TH1*> m_coinTreeTestHists;
  map<int, map<int ,vector<TH1*>>> m_ch_hist;
  map<DetType, map<int, vector<TH1*>>> m_basicHists;
 
  int timeDiff(string a_outFile);
  int m_lastPrintedEvent;
  double m_TDCRes;
  bool sortTrig(vector<double>& v1, vector<double>& v2);
 
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
};