GlueX Analysis Software
From GlueXWiki
Summary
- In his/her plugin, a user specifies the reaction(s) he/she wants to study (DReaction), along with the histograms and cuts to perform (if any) (DAnalysisAction), then just asks JANA for the results (DAnalysisResults).
- Through this process, the program creates all possible track combinations for the desired reaction (DParticleCombo, tag "PreKinFit"), kinematic fits the event to the reaction (DKinFitResults), and updates the tracks with the new kinfit information (DParticleCombo, no tag).
- The DAnalysisAction objects encapsulate the setup and execution of an action (e.g. cut, histogram) into a single object.
- Many common actions are predefined in DANA (see libraries/ANALYSIS/DHistogramActions.h and libraries/ANALYSIS/DCutActions.h), but the user can write custom ones for their analysis in their plugin.
- The DAnalysisResults objects indicate which DParticleCombo objects have passed/failed the DAnalysisAction cuts for each DReaction.
Quick Start
- NOTE: For an example of the below steps, see one of the example plugins below (e.g. b1pi).
- Create a new plugin containing a DEventProcessor, DReaction_factory, and DFactoryGenerator (for the DReaction_factory).
- It is recommended that you create these objects with unique names (and a unique tag for the DReaction_factory) so that this plugin can be executed with other analysis plugins simultaneously.
- Create desired custom DAnalysisAction objects for your analysis (if any).
- See existing analysis actions in sim-recon/src/libraries/ANALYSIS/DHistogramActions.h and sim-recon/src/libraries/ANALYSIS/DCutActions.h for examples.
- In DReaction_factory::init(void), create a DReaction object (with a unique name!) for each analysis you want to perform, specify the kinematic fit type for it, specify the DAnalysisActions for each, and add them to _data.
- The DAnalysisAction objects are executed sequentially for each possible particle combination (DParticleCombo) matching the given DReaction, until it fails a cut.
- If more than one DAnalysisAction object of a given type is used in a DReaction, all of their constructors must be called with a unique identifier string (including the default "").
- More details on analysis actions in sim-recon/src/libraries/ANALYSIS/: DAnalysisAction.*, DHistogramActions.*, DCutActions.*
- Setup the DEventProcessor.
- In DEventProcessor, make sure to grab the DAnalysisResults objects from the JEventLoop in the evnt() method, and to add the DFactoryGenerator to the application in the InitPlugin() method.
- Compile and run hd_root with your plugin (preferably on REST data).
Example Plugins
- b1pi: sim-recon/src/programs/Analysis/plugins/b1pi_hists/
- monitoring_hists: sim-recon/src/programs/Analysis/plugins/monitoring_hists/
Known Issues
- Kinematic fitting with vertex constraints that contain only one detected particle do not converge very often.
- e.g. The Λ production vertex in γ, p → K+, Λ; Λ → p, π-
- Normally can include the beam photon in the vertex fit, but not if the uncertainties are all zero (MC).
- Technically the kinematic fitter works, but the analysis framework currently does not give an accurate enough initial guess for the vertex position. You can fit this case manually if your input guess is accurate enough.
- e.g. The Λ production vertex in γ, p → K+, Λ; Λ → p, π-
- Spacetime kinematic fits don't work yet.
