// $Id$
//
//    File: JEventProcessor_CascadeStudies.h
// Created: Thu Jul 11 10:46:01 EDT 2019
// Creator: staylor (on Linux ifarm1401.jlab.org 3.10.0-327.el7.x86_64 x86_64)
//

#ifndef _JEventProcessor_CascadeStudies_
#define _JEventProcessor_CascadeStudies_

#include <JANA/JEventProcessor.h>
#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_CascadeStudies:public jana::JEventProcessor{
 public:
  JEventProcessor_CascadeStudies();
  ~JEventProcessor_CascadeStudies();
  const char* className(void){return "JEventProcessor_CascadeStudies";}
  
 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 CascadeAnalysis(DKinFitter *dKinFitter,
		       const DChargedTrack *proton_hyp,
		       const DChargedTrack *kplus_hyp1, 
		       const DChargedTrack *kplus_hyp2,       
		       const DChargedTrack *piminus_hyp,
		       const DChargedTrack *piminus_hyp2,
		       double weight);

  
  void FillPIDHistos(double t0_rf,
		     const DChargedTrackHypothesis *proton_hyp,
		     const DTrackTimeBased *proton_track,
		     const DChargedTrackHypothesis *kplus_hyp1,
		     const DTrackTimeBased *kplus_track1,  
		     const DChargedTrackHypothesis *kplus_hyp2,
		     const DTrackTimeBased *kplus_track2, 
		     const DChargedTrackHypothesis *piminus_hyp1,
		     const DTrackTimeBased *piminus_track1, 
		     const DChargedTrackHypothesis *piminus_hyp2, 
		     const DTrackTimeBased *piminus_track2,
		     const DParticleID *pid_algorithm);
  
  TH1F *CascadeConfidenceLevel,*CascadeChiSq;
  TH1F *P2PimMass,*PPimMass,*LambdaPimMass;
  TH2F *Xi_KpTOFdt,*Xi_KpFCALdt,*Xi_KpBCALdt,*Xi_KpdEdxCDC,*Xi_KpdEdxFDC;
  TH2F *Xi_KpEOverPFCAL,*Xi_KpEOverPBCAL; 
  TH2F *Xi_KpDenom,*Xi_KpNumer_Kp,*Xi_KpNumer_Pip,*Xi_KpNumer_P,*Xi_KpNumer_Ep;

  double CL_CUT;
  
    
};

#endif // _JEventProcessor_CascadeStudies_