// $Id$
//
//    File: DEventProcessor_pppmp0.cc
// Created: Tue Aug  6 10:51:29 EDT 2013
// Creator: yqiang (on Linux halldw1.jlab.org 2.6.32-358.11.1.el6.x86_64 x86_64)
//

#include "DEventProcessor_pppmp0.h"

// Routine used to create our DEventProcessor

extern "C"
{
	void InitPlugin(JApplication *locApplication)
	{
		InitJANAPlugin(locApplication);
		locApplication->AddProcessor(new DEventProcessor_pppmp0()); //register this plugin
		locApplication->AddFactoryGenerator(new DFactoryGenerator_pppmp0()); //register the factory generator
	}
} // "C"

//------------------
// init
//------------------
jerror_t DEventProcessor_pppmp0::init(void)
{
	// This is called once at program startup. If you are creating
	// and filling historgrams in this plugin, you should lock the
	// ROOT mutex like this:
	//
	// japp->RootWriteLock();
	//  ... create historgrams or trees ...
	// japp->RootUnLock();
	//

	return NOERROR;
}

//------------------
// brun
//------------------
jerror_t DEventProcessor_pppmp0::brun(jana::JEventLoop* locEventLoop, int locRunNumber)
{
	// This is called whenever the run number changes

	//Recommeded: Initialize reaction-independent analysis actions (nothing should be done if already initialized)
		//Warning: If running JANA with multiple plugins, these actions should only be executed in ONE of the plugins.
			//Else you will double/triple/etc.-fill your histograms!!
//	dHistogramAction_TrackMultiplicity.Initialize(locEventLoop);
//	dHistogramAction_ThrownParticleKinematics.Initialize(locEventLoop);
//	dHistogramAction_DetectedParticleKinematics.Initialize(locEventLoop);
//	dHistogramAction_GenReconTrackComparison.Initialize(locEventLoop);

	//Recommeded: Create output ROOT TTrees (nothing is done if already created)
	const DEventWriterROOT* locEventWriterROOT = NULL;
	locEventLoop->GetSingle(locEventWriterROOT);
	locEventWriterROOT->Create_DataTrees(locEventLoop);

	return NOERROR;
}

//------------------
// evnt
//------------------
jerror_t DEventProcessor_pppmp0::evnt(jana::JEventLoop* locEventLoop, int locEventNumber)
{
	// This is called for every event. Use of common resources like writing
	// to a file or filling a histogram should be mutex protected. Using
	// locEventLoop->Get(...) to get reconstructed objects (and thereby activating the
	// reconstruction algorithm) should be done outside of any mutex lock
	// since multiple threads may call this method at the same time.
	//
	// Here's an example:
	//
	// vector<const MyDataClass*> mydataclasses;
	// locEventLoop->Get(mydataclasses);
	//
	// japp->RootWriteLock();
	//  ... fill historgrams or trees ...
	// japp->RootUnLock();

	// DOCUMENTATION:
	// ANALYSIS library: https://halldweb1.jlab.org/wiki/index.php/GlueX_Analysis_Software

	//Recommended: Execute reaction-independent actions (fill histograms)
		//Warning: If running JANA with multiple plugins, these actions should only be executed in ONE of the plugins.
			//Else you will double/triple/etc.-fill your histograms!!
//	dHistogramAction_TrackMultiplicity(locEventLoop);
//	dHistogramAction_ThrownParticleKinematics(locEventLoop);
//	dHistogramAction_DetectedParticleKinematics(locEventLoop);
//	dHistogramAction_GenReconTrackComparison(locEventLoop);

	//Recommended: Write surviving particle combinations (if any) to output ROOT TTree
		//If no cuts are performed by the analysis actions added to a DReaction, then this saves all of its particle combinations. 
		//The event writer gets the DAnalysisResults objects from JANA, performing the analysis. 
	const DEventWriterROOT* locEventWriterROOT = NULL;
	locEventLoop->GetSingle(locEventWriterROOT);
	// string is DReaction factory tag: will fill trees for all DReactions that are defined in the specified factory
	locEventWriterROOT->Fill_DataTrees(locEventLoop, "pppmp0");

	//Optional: Get all particle combinations for all DReactions. 
		//If kinematic fits were requested, these contain both the measured and kinematic-fit track parameters
		//No cuts from DAnalysisActions are placed on these combos
//	vector<const DParticleCombo*> locParticleCombos;
//	locEventLoop->Get(locParticleCombos);
//	for(size_t loc_i = 0; loc_i < locParticleCombos.size(); ++loc_i)
//	{
//		const DParticleCombo* locParticleCombo = locParticleCombos[loc_i];
//		if(locParticleCombo->Get_Reaction()->Get_ReactionName() != "pppmp0")
//			continue; // particle combination was for a different reaction
//		//perform further analysis steps here...
//	}

	//Optional: Get the analysis results for all DReactions. 
		//Getting these objects triggers the analysis, if it wasn't performed already. 
		//These objects contain the DParticleCombo objects that survived the DAnalysisAction cuts that were added to the DReactions
//	vector<const DAnalysisResults*> locAnalysisResultsVector;
//	locEventLoop->Get(locAnalysisResultsVector);
//	for(size_t loc_i = 0; loc_i < locAnalysisResultsVector.size(); ++loc_i)
//	{
//		const DAnalysisResults* locAnalysisResults = locAnalysisResultsVector[loc_i];
//		if(locAnalysisResults->Get_Reaction()->Get_ReactionName() != "pppmp0")
//			continue; // analysis results were for a different reaction
//		//get the DParticleCombo objects for this DReaction that survived all of the DAnalysisAction cuts
//		deque<const DParticleCombo*> locPassedParticleCombos;
//		locAnalysisResults->Get_PassedParticleCombos(locPassedParticleCombos);
//		for(size_t loc_j = 0; loc_j < locPassedParticleCombos.size(); ++loc_j)
//		{
//			const DParticleCombo* locPassedParticleCombo = locPassedParticleCombos[loc_j];
//			//perform further analysis steps here...
//		}
//	}

	/*
	//Optional: Save event to output REST file. Use this to create skims.
	const DEventWriterREST* locEventWriterREST = NULL;
	locEventLoop->GetSingle(locEventWriterREST);
	locEventWriterREST->Write_RESTEvent(locEventLoop, "pppmp0"); //string is part of output file name
	*/

	return NOERROR;
}

//------------------
// erun
//------------------
jerror_t DEventProcessor_pppmp0::erun(void)
{
	// This is called whenever the run number changes, before it is
	// changed to give you a chance to clean up before processing
	// events from the next run number.
	return NOERROR;
}

//------------------
// fini
//------------------
jerror_t DEventProcessor_pppmp0::fini(void)
{
	// Called before program exit after event processing is finished.
	return NOERROR;
}