STOFPostProcessing.h

Go to the documentation of this file.

#ifndef _STOF_POST_PROCESSING_H_
#define _STOF_POST_PROCESSING_H_
//project
#include "STOFConfig.h"
#include "STOFEvent.h"
#include "COMPASSFileManagement.h"
//root
#include "TGraph.h"
class STOFPostProcessing : public CompassPostProcBase
{
public:
  STOFPostProcessing();
 
  STOFPostProcessing(std::string a_writeMode);
 
  void loadConfig(std::string a_jsonConfigName);
 
  void genTimeDiffMultiPMTTarg(std::string a_outFile);
 
 
  void genAmpCalibHists(std::string a_outFile);
 
  void genPSDCalibHists(std::string a_outFile);
 
  void genCoinTimingCalibrationHists(std::string a_outFile,
                                     bool        a_targDTCut = false
                                    );
 
 
  void genIncTimingCalibrationHists(std::string a_outFile);
  void reduceToSTOFEvents(std::string a_outfile = "stofEvents.root");
 
  int buildCoincidences(double a_coinWindow,
                        int    a_minMult = 2
                        );
 
 
  int addCoinsFromFile(std::string a_fileName, 
                       double      a_coinWindow,
                       int a_minMult = 2
                       );
  void checkPSDCuts(int a_ch,
                    string a_cutFile
                    );
 
  TGraph* investigateCoincidenceWindow(int a_startWindow,
                                       int a_stopWindow,
                                       int a_step
                                       );
  void setRunStart(double a_runStart);
private:
  STOFConfig m_config;
 
  bool characterizeCoin(compassCoin&           a_coinEv,
                        std::map<int,DetType>& a_targMap,
                        std::map<int,DetType>& a_scatMap,
                        std::vector<int>&      a_tarChannels,
                        std::vector<int>&      a_scatChannels,
                        bool&                  a_hasRF,
                        double&                a_incTime,
                        double&                a_coinTime,
                        double&                a_evWallTime
                        );
  bool reduceCoinToSTOFEvent(compassCoin&  a_coinEv,
                             STOFEvent* a_stofEv
                             );
 
  double getMPMTDT(compassCoin& a_ev,
                     int          a_detID
                     );
 
  int m_maxRawPH;
  double m_runStartTime;
 
 
};
 
#endif