// $Id$
//
//    File: DNeutralParticleHypothesis_factory.h
// Created: Tue Aug  9 14:29:24 EST 2011
// Creator: pmatt (on Linux ifarml6 2.6.18-128.el5 x86_64)
//

#ifndef _DNeutralParticleHypothesis_factory_
#define _DNeutralParticleHypothesis_factory_

#include <JANA/JFactory.h>
#include <PID/DNeutralParticleHypothesis.h>
#include <PID/DParticleID.h>
#include <PID/DNeutralShower.h>
#include <PID/DChargedTrackHypothesis.h>
#include <PID/DChargedTrack.h>
#include <BCAL/DBCALShower.h>
#include <FCAL/DFCALShower.h>
#include <DVector3.h>
#include <DMatrixDSym.h>
#include <DMatrix.h>
#include <TMath.h>

class DNeutralParticleHypothesis_factory:public jana::JFactory<DNeutralParticleHypothesis>{
	public:
		DNeutralParticleHypothesis_factory(){};
		~DNeutralParticleHypothesis_factory(){};

		void Build_ErrorMatrix(const DVector3 &locPathVector, double locEnergy, const DMatrixDSym& locVariances, DMatrixDSym& locErrorMatrix);
		void Calc_Variances(const DNeutralShower *locNeutralShower, double locTimeDifferenceVariance, DMatrixDSym &locVariances);

	private:
		DParticleID *dPIDAlgorithm;
		double dTargetLength;
		double dTargetRadius;
		DVector3 dTargetCenter;

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

#endif // _DNeutralParticleHypothesis_factory_