Difference between revisions of "Mattione GlueX Analysis Factories"

From GlueXWiki
Jump to: navigation, search
(DTrackCombo)
(DTrackCombo)
Line 41: Line 41:
 
  const DReaction* <span style="color:#0000FF">dReaction</span>;
 
  const DReaction* <span style="color:#0000FF">dReaction</span>;
 
  deque<const DTrackComboStep*> <span style="color:#0000FF">dTrackComboSteps</span>;
 
  deque<const DTrackComboStep*> <span style="color:#0000FF">dTrackComboSteps</span>;
const DKinFitResults* <span style="color:#0000FF">dKinFitResults</span>;
+
const DKinFitResults* <span style="color:#0000FF">dKinFitResults</span>;
 
+
 
  deque<const DChargedTrack*> <span style="color:#0000FF">dUnusedPositiveTracks</span>;
 
  deque<const DChargedTrack*> <span style="color:#0000FF">dUnusedPositiveTracks</span>;
 
  deque<const DChargedTrack*> <span style="color:#0000FF">dUnusedNegativeTracks</span>;
 
  deque<const DChargedTrack*> <span style="color:#0000FF">dUnusedNegativeTracks</span>;

Revision as of 15:41, 12 August 2012

Summary

  • In his/her plugin, a user specifies the reactions he/she wants to study (DReaction objects), then just asks JANA for the results (DTrackCombo objects).
    • Creates all possible track combinations for the desired reaction; user then cuts on the likelihood that those combinations represent the hypothesized reaction.

DReaction

  • Will be located in src/libraries/ANALYSIS
  • User creates one of these in their plugin for each reaction they want to analyze (example shown later below) (can analyze more than one at once).
class DReaction : public JObject
{
	private:
		//REACTION AND ANALYSIS
		string dReactionName; //must be unique: e.g. "Y2175"
		deque<const DReactionStep*> dReactionSteps;
		deque<const DAnalysisAction*> dAnalysisActions;

		//CONTROL VARIABLES
		unsigned short int dFitTypeFlag; //0 is none, 1 is p4 only, 2 is p4 & v, 3 is p4, v, & t
		deque<size_t> dDecayingParticlesExcludedFromKinfit; //to exclude decaying particles from the kinematic fit (resonances are automatically excluded)
};
class DReactionStep
{
	private:
		Particle_t dInitialParticleID; //e.g. lambda, phi, gamma
		Particle_t dTargetParticleID; //Unknown for no target
		deque<Particle_t> dFinalParticleIDs;

		int dMissingParticleIndex; //-1 for no missing particles, else final state particle at this index is missing (0 -> x)
};

DTrackCombo

  • Will be located in src/libraries/ANALYSIS
  • DTrackCombo_factory_PreKinFit: Generates a DTrackCombo object for each possible combination of detected tracks and DBeamPhotons (configurable time cut used to eliminate way-out-of-time photons).
    • PID is NOT assumed, all possibilities are tested: For each q = +/-/0 particle in a given DReaction, all detected q = +/-/0 tracks are tested as to whether or not they are that particle.
  • DKinFitResults_factory: Performs a single kinematic fit testing whether each DTrackCombo matches the DReaction: includes all mass, vertex, timing constraints if requested.
  • DTrackCombo_factory: Create new DTrackCombo objects, using the kinematic fit track information and results.
class DTrackCombo : public JObject
{
	private:
		const DReaction* dReaction;
		deque<const DTrackComboStep*> dTrackComboSteps;
		const DKinFitResults* dKinFitResults;

		deque<const DChargedTrack*> dUnusedPositiveTracks;
		deque<const DChargedTrack*> dUnusedNegativeTracks;
		deque<const DNeutralShower*> dUnusedNeutralShowers;
};
class DTrackComboStep
{
	private:
		//PIDS:
		Particle_t dInitialParticleID; //e.g. lambda, phi, gamma
		Particle_t dTargetParticleID; //Unknown for no target
		deque<Particle_t> dFinalParticleIDs;

		//INITIAL PARTICLES:
		DKinematicData* dInitialParticle; //kinfit result, else measured
		DKinematicData* dInitialParticle_Measured; //e.g. DBeamPhoton
		DKinematicData* dTargetParticle; //NULL for no target

		//FINAL PARTICLES:
		deque<const DKinematicData*> dFinalParticles; //kinfit result, else measured //e.g. DChargedTrackHypothesis, DNeutralParticleHypothesis
		deque<const DKinematicData*> dFinalParticles_Measured; //e.g. DChargedTrackHypothesis, DNeutralParticleHypothesis
		deque<const JObject*> dFinalParticleSourceObjects; //original DChargedTrack or DNeutralShower objects

		//CONTROL VARIABLES:
		int dInitialParticleDecayFromStepIndex; //points to which DTrackComboStep represents the production of the initial particle (e.g. Lambda)
		deque<int> dDecayStepIndices; //let's you know whether a final state particle is decaying, missing, or detected.  if decaying, points to which DTrackComboStep represents its decay
};

Cutting & Histogramming: DAnalysisAction, DAnalysisResults

* DAnalysisAction: virtual base class that encapsulates the setup (e.g. histogram creation) and execution (cut, hist, etc.) of the contained action for each DTrackCombo of the given DReaction that has passed all previous cuts.
class DAnalysisAction
{
	public:
		//returns true/false if passes/fails cut/hist/etc.
		virtual bool operator()(JEventLoop* locEventLoop, map<const DTrackCombo*, bool>& locTrackCombos) const = 0;
	protected:
		DReaction* dReaction; //tells it to ignore DTrackCombo objects from other reactions
};

User 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 DAnalysisResults*> locAnalysisResultsVector;
	locEventLoop->Get(locAnalysisResultsVector);

	//SAVE OUTPUT HERE

	return NOERROR;
}

User Plugin DReaction Factory

  • DReaction.h located in DANA libraries, DReaction_factory.h located in user plugin (along with a factory generator for it).
  • Below shows how to set up DReaction for a Y(2175) analysis.
//------------------
// init
//------------------
jerror_t DReaction_factory::init(void)
{
	// Setting the PERSISTANT prevents JANA from deleting
	// the objects every event so we only create them once.
	SetFactoryFlag(PERSISTANT);

	DReaction* locReaction;
	DReactionStep* locReactionStep;

	// Make as many DReaction objects as desired
	locReaction = new DReaction();
	locReaction->Set_ReactionName("Y2175");
	locReaction->Set_FitTypeFlag(3); //kinfit p4, v, & t simultaneously (0 is no kinfit)
 
	//g, p -> phi(2175), (p)
	locReactionStep = new DReactionStep();
	locReactionStep->Set_InitialParticleID(Gamma);
	locReactionStep->Set_TargetParticleID(Proton);
	locReactionStep->Add_FinalParticleID(Unknown); //phi(2175)
	locReactionStep->Add_FinalParticleID(Proton);
	locReactionStep->Set_MissingParticleIndex(1); //proton missing
	locReaction->Add_ReactionStep(locReactionStep);
	dReactionStepPool.push_back(locReactionStep); //prevent memory leak

	//phi(2175) -> pi+, pi-, phi
	locReactionStep = new DReactionStep();
	locReactionStep->Set_InitialParticleID(Unknown); //phi(2175)
	locReactionStep->Set_TargetParticleID(Unknown); //no target for this step
	locReactionStep->Add_FinalParticleID(PiPlus);
	locReactionStep->Add_FinalParticleID(PiMinus);
	locReactionStep->Add_FinalParticleID(phiMeson);
	locReactionStep->Set_MissingParticleIndex(-1); //none missing
	locReaction->Add_ReactionStep(locReactionStep);
	dReactionStepPool.push_back(locReactionStep); //prevent memory leak

	//phi -> K+, K-
	locReactionStep = new DReactionStep();
	locReactionStep->Set_InitialParticleID(Unknown); //phi(2175)
	locReactionStep->Set_TargetParticleID(Unknown); //no target for this step
	locReactionStep->Add_FinalParticleID(KPlus);
	locReactionStep->Add_FinalParticleID(KMinus);
 	locReactionStep->Set_MissingParticleIndex(-1); //none missing
 	locReaction->Add_ReactionStep(locReactionStep);
	dReactionStepPool.push_back(locReactionStep); //prevent memory leak

	//Track reconstruction quality/multiplicity
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TrackMultiplicity());
	dReaction->Add_AnalysisAction(new DHistogramFunctor_GenReconTrackComparison());

	//Vertex reconstruction and cuts of tracks with bogus vertex-z
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TrackComboKinematics(1)); //1: initial data
	dReaction->Add_AnalysisAction(new DCutFunctor_AllVertexZ(45.0, 85.0));
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TrackVertexComparison());

	//PID
	dReaction->Add_AnalysisAction(new DHistogramFunctor_PID());
	dReaction->Add_AnalysisAction(new DCutFunctor_AllPID(0.01, 1)); //1%, all charged tracks
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TruePID());

	//Phi Invariant Mass Cut & PID Check
	dReaction->Add_AnalysisAction(new DHistogramFunctor_PhiMass(1)); //1: initial data
	dReaction->Add_AnalysisAction(new DCutFunctor_InvariantMass(phiMeson, 1, 1.004, 1.04)); //~3sigma
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TruePID());

	//Initial Proton and Y2175 Mass Distributions
	dReaction->Add_AnalysisAction(new DHistogramFunctor_ProtonMass(1)); //1: initial data
	dReaction->Add_AnalysisAction(new DHistogramFunctor_Phi2175Mass(1)); //1: initial data

	//Missing Proton Kinematic Fit
	dReaction->Add_AnalysisAction(new DHistogramFunctor_KinFitResults(0.1));
	dReaction->Add_AnalysisAction(new DCutFunctor_KinFitFOM(0.01); //1%

	//Final Proton and Y2175 Mass Distributions
	dReaction->Add_AnalysisAction(new DHistogramFunctor_ProtonMass(1)); //1: initial data
	dReaction->Add_AnalysisAction(new DHistogramFunctor_Phi2175Mass(1)); //1: initial data
	dReaction->Add_AnalysisAction(new DHistogramFunctor_Phi2175Mass(2)); //2: final data

	//Final track kinematics & PID Check
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TrackComboKinematics(2)); //2: final data
	dReaction->Add_AnalysisAction(new DHistogramFunctor_TruePID());

	_data.push_back(locReaction);

	return NOERROR;
}

Tupling

  • Haven't really gotten this far yet...