// $Id$
//
//    File: JEventProcessor_EpEmStudies.h
// Created: Fri Jun 28 10:04:30 EDT 2019
// Creator: staylor (on Linux ifarm1401.jlab.org 3.10.0-327.el7.x86_64 x86_64)
//

#ifndef _JEventProcessor_EpEmStudies_
#define _JEventProcessor_EpEmStudies_

#include <JANA/JEventProcessor.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <KINFITTER/DKinFitter.h>
#include <ANALYSIS/DKinFitUtils_GlueX.h>
#include <PID/DChargedTrack.h>
#include <PID/DParticleID.h>
#include <ANALYSIS/DAnalysisUtilities.h>

class JEventProcessor_EpEmStudies:public jana::JEventProcessor{
 public:
  JEventProcessor_EpEmStudies();
  ~JEventProcessor_EpEmStudies();
  const char* className(void){return "JEventProcessor_EpEmStudies";}
  
 private:
  jerror_t init(void);						///< Called once at program start.
  jerror_t brun(jana::JEventLoop *eventLoop, int32_t runnumber);	///< Called everytime a new run number is detected.
  jerror_t evnt(jana::JEventLoop *eventLoop, uint64_t eventnumber);	///< Called every event.
  jerror_t erun(void);						///< Called everytime run number changes, provided brun has been called.
  jerror_t fini(void);						///< Called after last event of last event source has been processed.

  bool EpEmGammaAnalysis(DKinFitter *dKinFitter,
			 const DChargedTrack *proton_hyp,
			 const DChargedTrack *positron_hyp,
			 const DChargedTrack *electron_hyp,double weight); 
  bool EpEmAnalysis(DKinFitter *dKinFitter,
		    const DChargedTrack *proton_hyp,
		    const DChargedTrack *positron_hyp,
		    const DChargedTrack *electron_hyp,double weight);
  void FillPIDHistos(double t0_rf,
		     const DChargedTrackHypothesis *proton_hyp,
		     const DTrackTimeBased *proton_track,
		     const DChargedTrackHypothesis *positron_hyp,
		     const DTrackTimeBased *positron_track,
		     const DChargedTrackHypothesis *electron_hyp,
		     const DTrackTimeBased *electron_track,
		     const DParticleID *pid_algorithm);
    
  TH1F *EpEmGammaConfidenceLevel,*EpEmGammaChiSq,*EpEmGammaMass,*EpEmMass; 
  TH1F *EpEmConfidenceLevel,*EpEmChiSq,*EpEmMassNoGamma;
  
  
  TH2F *ElectronTOFdt,*ElectronEOverPFCAL,*ElectronFCALdt,*ElectronBCALdt;
  TH2F *ElectronEOverPBCAL,*ElectrondEdxCDC,*ElectrondEdxFDC;
  TH2F *ElectronFCALdtPull,*ElectronTOFdtPull,*ElectronBCALdtPull;
  TH1F *ElectronTOFdtCL,*ElectronFCALdtCL,*ElectronBCALdtCL;
  TH1F *ElectrondEdxCDCCL,*ElectrondEdxFDCCL;
  TH2F *ElectrondEdxCDCPull,*ElectrondEdxFDCPull;
  TH2F *ElectronEOverPPullBCAL,*ElectronEOverPPullFCAL;
  TH1F *ElectronEOverPBCALCL,*ElectronEOverPFCALCL;

  TH2F *EmDenom,*EmNumer_Km,*EmNumer_Pim,*EmNumer_aP,*EmNumer_Em;
  TH2F *EpDenom,*EpNumer_Kp,*EpNumer_Pip,*EpNumer_P,*EpNumer_Ep;

  double CL_CUT;
};

#endif // _JEventProcessor_EpEmStudies_