// $Id$
//
//    File: DFCALDigiHit_factory_DAVECUT.cc
// Created: Mon Nov 24 22:13:17 EST 2014
// Creator: davidl (on Linux gluon47.jlab.org 2.6.32-358.23.2.el6.x86_64 x86_64)
//


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

#include <stdint.h>

#include <DAQ/Df250PulsePedestal.h>

#include "DFCALDigiHit_factory_DAVECUT.h"
#include "JFactoryGenerator_DFCALDigiHit_DAVECUT.h"
using namespace jana;

// Routine used to create our JEventProcessor
#include <JANA/JApplication.h>
//extern "C"{
//void InitPlugin(JApplication *app){
//	InitJANAPlugin(app);
//	app->AddFactoryGenerator(new JFactoryGenerator_DFCALDigiHit_DAVECUT());
//}
//} // "C"

//------------------
// init
//------------------
jerror_t DFCALDigiHit_factory_DAVECUT::init(void)
{
	flags = NOT_OBJECT_OWNER;

	return NOERROR;
}

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

//------------------
// evnt
//------------------
jerror_t DFCALDigiHit_factory_DAVECUT::evnt(JEventLoop *loop, int eventnumber)
{
	vector<const DFCALDigiHit*> fcaldigihits;
	loop->Get(fcaldigihits, "", false);

	for(uint32_t i=0; i<fcaldigihits.size(); i++){
		const DFCALDigiHit *fcaldigihit = fcaldigihits[i];
		
		double pulse_integral = (double)fcaldigihit->pulse_integral;
		if( pulse_integral < 5000.0 ) continue;
		
		const Df250PulsePedestal *pp = NULL;
		fcaldigihit->GetSingle(pp);
		if(!pp) continue;
		double pulse_peak = (double)pp->pulse_peak;
		
		double cut = 5000.0 + 40000.0/3500.0*pulse_peak;
		if( pulse_integral > cut ) continue;
		
		_data.push_back((DFCALDigiHit*)fcaldigihit);
	}
//_DBG_<<"fcaldigihits.size()="<<fcaldigihits.size()<<"  _data.size()="<<_data.size()<<endl;
	return NOERROR;
}

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

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