Difference between revisions of "Mattione GlueX Analysis Factories"

Revision as of 00:31, 16 July 2012

Plugin DReaction Factory

• DReaction.h located in DANA libraries, DReaction.h_factory located in user plugin (along with a factory generator for it).
• Below shows how to set up DReaction for a Y(2175) analysis.

jerror_t DReaction_factory::init(void)
{
	// Make as many DReaction objects as desired
	DReaction* locReaction = new DReaction();
	{
		//g, p -> Y(2175), (p)
		DReactionStep locReactionStep;
		locReactionStep.dInitialParticleID = Gamma;
		locReactionStep.dTargetParticleID = Proton;
		locReactionStep.dFinalParticleIDs.push_back(Unknown); //Y(2175)
		locReactionStep.dFinalParticleIDs.push_back(Proton);
		locReactionStep.dMissingParticleIndex = 1; //proton missing
		locReaction->dReactionSteps.push_back(locReactionStep);
	}
	{
		//Y(2175) -> pi+, pi-, phi
		DReactionStep locReactionStep;
		locReactionStep.dInitialParticleID = Unknown; //Y(2175)
		locReactionStep.dTargetParticleID = Unknown;
		locReactionStep.dFinalParticleIDs.push_back(PiPlus);
		locReactionStep.dFinalParticleIDs.push_back(PiMinus);
		locReactionStep.dFinalParticleIDs.push_back(phiMeson);
		locReactionStep.dMissingParticleIndex = -1; //none missing
		locReaction->dReactionSteps.push_back(locReactionStep);
	}
	{
		//phi -> K+, K-
		DReactionStep locReactionStep;
		locReactionStep.dInitialParticleID = phiMeson;
		locReactionStep.dTargetParticleID = Unknown;
		locReactionStep.dFinalParticleIDs.push_back(KPlus);
		locReactionStep.dFinalParticleIDs.push_back(KMinus);
		locReactionStep.dMissingParticleIndex = -1; //none missing
		locReaction->dReactionSteps.push_back(locReactionStep);
	}
	_data.push_back(locReaction);
}

Plugin DEventProcessor

extern "C" {
	void InitPlugin(JApplication *app){
		InitJANAPlugin(app);
		app->AddProcessor(new DEventProcessor_Y2175());
		app->AddFactoryGenerator(new DFactoryGenerator_DReaction());
	}
} // "C"

jerror_t DEventProcessor_Y2175::evnt(JEventLoop *locEventLoop, int eventnumber)
{
	vector<const DKinFitTrackCombo*> locKinFitTrackCombo;
	locEventLoop->Get(locKinFitTrackCombo);

	//DO HISTS/CUTS/TUPLES HERE

	return NOERROR;
}