// $Id$
//
//    File: DVertex_factory_THROWN.cc
// Created: Tue Apr  6 17:01:54 EDT 2010
// Creator: pmatt (on Darwin Amelia.local 9.8.0 i386)
//


#include <iostream>
#include <iomanip>
using namespace std;

#include <TROOT.h>
#include <TMath.h>
#include "DVertex_factory_THROWN.h"
using namespace jana;

//------------------
// init
//------------------
jerror_t DVertex_factory_THROWN::init(void)
{
	dTargetCenter = 65.0;
	return NOERROR;
}

//------------------
// brun
//------------------
jerror_t DVertex_factory_THROWN::brun(jana::JEventLoop *loop, int runnumber)
{
	return NOERROR;
}

//------------------
// evnt
//------------------
jerror_t DVertex_factory_THROWN::evnt(JEventLoop *loop, int eventnumber)
{
	vector<const DChargedTrack*> locChargedTracks;
	loop->Get(locChargedTracks);
	vector<const DMCThrown*> locThrownTracks;
	loop->Get(locThrownTracks);

	if(locThrownTracks.size() == 0)
		return RESOURCE_UNAVAILABLE;

	DVertex* locVertex = new DVertex;
	locVertex->locCovarianceMatrix.ResizeTo(3,3);
	locVertex->dSpacetimeVertex.SetVect(locThrownTracks[0]->position());
	locVertex->dSpacetimeVertex.SetT(0.0 + (locThrownTracks[0]->position().Z() - dTargetCenter)/SPEED_OF_LIGHT);
	locVertex->dTimeUncertainty = 0.;

	// Add list of tracks used to create this vertex
	for(unsigned int loc_j = 0; loc_j < locChargedTracks.size(); loc_j++)
		locVertex->dChargedTracks.push_back(locChargedTracks[loc_j]);
	_data.push_back(locVertex);	

	return NOERROR;
}

//------------------
// erun
//------------------
jerror_t DVertex_factory_THROWN::erun(void)
{
	return NOERROR;
}

//------------------
// fini
//------------------
jerror_t DVertex_factory_THROWN::fini(void)
{
	return NOERROR;
}