// $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/DVertex.h>
#include <PID/DChargedTrack.h>
#include <BCAL/DBCALShower.h>
#include <BCAL/DBCALHit.h>
#include <FCAL/DFCALShower.h>
#include <DVector3.h>
#include <DMatrixDSym.h>
#include <TMath.h>

class DNeutralParticleHypothesis_factory:public jana::JFactory<DNeutralParticleHypothesis>{
public:
  DNeutralParticleHypothesis_factory();
  ~DNeutralParticleHypothesis_factory(){};
  
  bool BCALShowerLacksEntryHit(vector<const DBCALShower*> &shower);
  bool BCALShowerLikelyEM(const DBCALShower* shower);
  void Build_ErrorMatrix(const DVector3 &locPathVector, double locEnergy, 
			 const DMatrixDSym& locVariances, DMatrixDSym& locErrorMatrix);
  void Calc_Variances(const DNeutralShower *locNeutralShower, 
		      double locParticleEnergyUncertainty, DMatrixDSym &locVariances);
  
private:
  DParticleID *dPIDAlgorithm;
  double dTargetLength;
  double dTargetRadius;
  double dTargetCenter_Z;

  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.
  
  bool CENTER_VERTEX, REQ_BCAL_ENTRY_HIT, TEST_EM_SHOWER;
  float EM_BACK_E_THRES, EM_CENT_THRES;
};

#endif // _DNeutralParticleHypothesis_factory_