// $Id$
//
//    File: JEventProcessor_AntiProtonStudies.h
// Created: Mon Jul  8 13:08:29 EDT 2019
// Creator: staylor (on Linux ifarm1801 3.10.0-327.el7.x86_64 x86_64)
//

#ifndef _JEventProcessor_AntiProtonStudies_
#define _JEventProcessor_AntiProtonStudies_

#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 <TAGGER/DTAGMGeometry.h>
#include <TAGGER/DTAGHGeometry.h>
#include <PID/DMCReaction.h>
#include <TRACKING/DMCThrown.h>

class JEventProcessor_AntiProtonStudies:public jana::JEventProcessor{
 public:
  JEventProcessor_AntiProtonStudies();
  ~JEventProcessor_AntiProtonStudies();
  const char* className(void){return "JEventProcessor_AntiProtonStudies";}
  
 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 LambdaLambdaBarAnalysis(DKinFitter *dKinFitter,
			       const DChargedTrack *piplus_hyp,
			       const DChargedTrack *proton_hyp1, 
			       const DChargedTrack *proton_hyp2,	       
			       const DChargedTrack *piminus_hyp,
			       const DChargedTrack *antiproton_hyp,
			       double weight);
  void FillPIDHistos(double t0_rf,const DChargedTrackHypothesis *proton_hyp1,
		     const DTrackTimeBased *proton_track1,
		     const DChargedTrackHypothesis *proton_hyp2,
		     const DTrackTimeBased *proton_track2,  
		     const DChargedTrackHypothesis *piplus_hyp,
		     const DTrackTimeBased *piplus_track, 
		     const DChargedTrackHypothesis *piminus_hyp,
		     const DTrackTimeBased *piminus_track, 
		     const DChargedTrackHypothesis *antiproton_hyp, 
		     const DTrackTimeBased *antiproton_track,
		     const DParticleID *pid_algorithm);

  TH1F *LambdaLambdaBarConfidenceLevel,*LambdaLambdaBarChiSq;
  TH2F *AntiPPipVsPPim;
  TH2F *AntiPTOFdt,*AntiPFCALdt,*AntiPBCALdt,*AntiPdEdxCDC,*AntiPdEdxFDC;
  TH2F *AntiPEOverPFCAL,*AntiPEOverPBCAL; 
  TH2F *PTOFdt,*PFCALdt,*PBCALdt,*PdEdxCDC,*PdEdxFDC;
  TH2F *PEOverPFCAL,*PEOverPBCAL;

  TH2F *PDenom,*PNumer_Kp,*PNumer_Pip,*PNumer_P,*PNumer_Ep;
  TH2F *AntiPDenom,*AntiPNumer_Km,*AntiPNumer_Pim,*AntiPNumer_aP,*AntiPNumer_Em;

  double CL_CUT;
};

#endif // _JEventProcessor_AntiProtonStudies_