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

#ifndef _JEventProcessor_LambdaStudies_
#define _JEventProcessor_LambdaStudies_

#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_LambdaStudies:public jana::JEventProcessor{
 public:
  JEventProcessor_LambdaStudies();
  ~JEventProcessor_LambdaStudies();
  const char* className(void){return "JEventProcessor_LambdaStudies";}
  
 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 LambdaAnalysis(DKinFitter *dKinFitter,
		      const DChargedTrack *proton_hyp,
		      const DChargedTrack *kplus_hyp,
		      const DChargedTrack *piminus_hyp,double weight);
  void FillPIDHistos(double t0_rf,
		     const DChargedTrackHypothesis *proton_hyp,
		     const DTrackTimeBased *proton_track,
		     const DChargedTrackHypothesis *kplus_hyp,
		     const DTrackTimeBased *kplus_track,
		     const DChargedTrackHypothesis *piminus_hyp,
		     const DTrackTimeBased *piminus_track,
		     const DParticleID *pid_algorithm);
  
  TH1F *LambdaConfidenceLevel,*LambdaChiSq,*PimPMass;
  TH2F *KpDenom_Lambda,*KpNumer_Kp_Lambda,*KpNumer_P_Lambda,*KpNumer_Ep_Lambda,*KpNumer_Pip_Lambda; 
  TH2F *PimDenom_Lambda,*PimNumer_Km_Lambda,*PimNumer_aP_Lambda,*PimNumer_Em_Lambda,*PimNumer_Pim_Lambda;
  TH2F *PDenom_Lambda,*PNumer_Pip_Lambda,*PNumer_Kp_Lambda,*PNumer_P_Lambda,*PNumer_Ep_Lambda; 
  TH2F *KpdEdxCDCPull_Lambda,*KpdEdxFDCPull_Lambda;
  TH1F *KpdEdxCDCCL_Lambda,*KpdEdxFDCCL_Lambda; 
  TH2F *PTOFdt_Lambda,*PBCALdt_Lambda,*PFCALdt_Lambda,*PdEdxCDC_Lambda,*PdEdxFDC_Lambda;  
  TH2F *KpTOFdt_Lambda,*KpBCALdt_Lambda,*KpFCALdt_Lambda,*KpdEdxCDC_Lambda,*KpdEdxFDC_Lambda;
  TH2F *KpBCALdtPull,*KpTOFdtPull,*KpFCALdtPull;
  TH1F *KpBCALdtCL,*KpFCALdtCL,*KpTOFdtCL;

  double CL_CUT;
};

#endif // _JEventProcessor_LambdaStudies_