// $Id$
//
//    File: JEventProcessor_PipPimPi0Studies.h
// Created: Fri Jan 17 16:34:34 EST 2020
// Creator: staylor (on Linux ifarm1802.jlab.org 3.10.0-1062.4.1.el7.x86_64 x86_64)
//

#ifndef _JEventProcessor_PipPimPi0Studies_
#define _JEventProcessor_PipPimPi0Studies_

#include <JANA/JEventProcessor.h>

#include <PID/DBeamPhoton.h>
#include <TTree.h>
#include <TH1F.h>
#include <TH2F.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>
#include <PID/DNeutralShower.h>
#include <BCAL/DBCALShower.h>
#include <TAGGER/DTAGMGeometry.h>
#include <TAGGER/DTAGHGeometry.h>

class JEventProcessor_PipPimPi0Studies:public jana::JEventProcessor{
 public:
  JEventProcessor_PipPimPi0Studies();
  ~JEventProcessor_PipPimPi0Studies();
  const char* className(void){return "JEventProcessor_PipPimPi0Studies";}
  
 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.

  void FillParticleVectors(vector<const DChargedTrack *>&tracks,
			   map<Particle_t, vector<const DTrackTimeBased *> >&particles
			   ) const;

  static bool SortParticleProbability(const pair<double,Particle_t>a,
				     const pair<double,Particle_t>b){
   return a.first>b.first;
  }
  TH2F *AcceptedP,*ThrownP;

  TTree *threePiQTree,*threePiQMCTree;
  float treeval[100],mctreeval[100];

  double FCAL_THRESHOLD,BCAL_THRESHOLD,SPLIT_CUT,FCAL_RADIAL_CUT,BCAL_Z_CUT;
  double TIME_CUT;
};

#endif // _JEventProcessor_PipPimPi0Studies_