// $Id$
//
//    File: JEventProcessor_pip_pim.h
// Created: Sat Apr 14 13:18:21 EDT 2012
// Creator: davidl (on Darwin genmacbook.local 11.3.0 i386)
//

#ifndef _JEventProcessor_pip_pim_
#define _JEventProcessor_pip_pim_

#include <pthread.h>

#include <TH2.h>

#include <JANA/JEventProcessor.h>
#include <DLorentzVector.h>
#include <TLorentzVector.h>
#include <TTree.h>

class JEventProcessor_pip_pim:public jana::JEventProcessor{
	public:
		JEventProcessor_pip_pim();
		~JEventProcessor_pip_pim();
		const char* className(void){return "JEventProcessor_pip_pim";}

	private:
		jerror_t init(void);						///< Called once at program start.
		jerror_t brun(jana::JEventLoop *eventLoop, int runnumber);	///< Called everytime a new run number is detected.
		jerror_t evnt(jana::JEventLoop *eventLoop, int 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.

		TH1D *mass_pip_pim;
		TH1D *mass_pip_pim_gen;
	
		TLorentzVector *v4_pip, *v4_pim;
		TLorentzVector *v4_pip_gen, *v4_pim_gen;
		TTree *pip_tree, *pim_tree;
		TTree *pip_tree_gen, *pim_tree_gen;
	
		TH1D *prob_pip, *prob_pim;
		TH1D *prob_track_pip, *prob_track_pim;
		
		pthread_rwlock_t root_lock;
	
		void FillInvariantMassHisto(TH1D *h, vector<DLorentzVector> &pip, vector<DLorentzVector> &pim);
};

#endif // _JEventProcessor_pip_pim_