// $Id$
//
//    File: DPiMinus_factory.cc
// Created: Sat Apr 14 12:13:05 EDT 2012
// Creator: davidl (on Darwin genmacbook.local 11.3.0 i386)
//


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

#include <PID/DChargedTrack.h>

#include "DPiMinus_factory.h"
using namespace jana;

//------------------
// init
//------------------
jerror_t DPiMinus_factory::init(void)
{
	return NOERROR;
}

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

//------------------
// evnt
//------------------
jerror_t DPiMinus_factory::evnt(JEventLoop *loop, int eventnumber)
{
	// Get list of all charged tracks
	vector<const DChargedTrack*> chargedtracks;
	loop->Get(chargedtracks);
	
	// Loop over charged tracks looking for ones where the most
	// probable hypothesis is a pi+
	for(unsigned int i=0; i<chargedtracks.size(); i++){
		const DChargedTrack *chargedtrack = chargedtracks[i];
		if(chargedtrack->dChargedTrackHypotheses.size()<1)continue;
		
		const DChargedTrackHypothesis* hypothesis = chargedtrack->dChargedTrackHypotheses[0];
		if(hypothesis->PID() == PiMinus){
			// Most probable hypothesis is a pi+. Make copy and store it in factory
			DPiMinus *pim = new DPiMinus(hypothesis);
			_data.push_back(pim);
		}
	}
	
	return NOERROR;
}

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

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