Merge pull request #229 from JeffersonLab/flatten_working

Added particle vertex info to FlattenForFSRoot

FlattenForFSRoot/Documentation/GlueXFSRootFormat.tex

@@ -70,27 +70,36 @@ \section{Event Information}
   NumBeam:          number of beam hypotheses from tagger
   NumCombos:        number of combos
 
-  ProdVx:       production vertex parameters
-  ProdVy:         "
+  ProdVx:       production vertex parameters (from X4_Production)
+  ProdVy:         "                       (NOT particle vertices)
   ProdVz:         "
   ProdVt:         "
   PxPB:         kinematically fit beam parameters
-  PyPB:           "
-  PzPB:           "
+  PyPB:           "             (four-momentum -- Px,Py,Pz,En
+  PzPB:           "                 and vertex -- Vx,Vy,Vz)
   EnPB:           "
+  VxPB:           "
+  VyPB:           "
+  VzPB:           "
   RPxPB:        measured beam parameters
-  RPyPB:          "
-  RPzPB:          "
+  RPyPB:          "             (four-momentum -- Px,Py,Pz,En
+  RPzPB:          "                 and vertex -- Vx,Vy,Vz)
   REnPB:          "
+  RVxPB:          "
+  RVyPB:          "
+  RVzPB:          "
   MCPxPB:       thrown beam parameters (MC only)
-  MCPyPB:         "
-  MCPzPB:         "
+  MCPyPB:         "             (four-momentum -- Px,Py,Pz,En
+  MCPzPB:         "                 and vertex -- Vx,Vy,Vz)
   MCEnPB:         "
+  MCVxPB:         "
+  MCVyPB:         "
+  MCVzPB:         "
 \end{verbatim}
 
 
 
-\section{Particle Four-Momenta}
+\section{Particle Four-Momenta and Vertices}
 \label{sec:nt:momentum}
 
 The ``Particle Four-Momenta'' variables contain the four-momentum for each particle in the final state:
@@ -102,9 +111,17 @@ \section{Particle Four-Momenta}
   (prefix)EnP(n):  energy of particle (n)
 \end{verbatim}
 
-Different types of four-momenta are distinguished using prefixes.  Raw four-momenta have a prefix~{\tt R}; the final four-momenta (the fully-constrained four-momenta resulting from the kinematic fit) have no prefix; MC truth information have a prefix~{\tt MC}.
+The vertices for all non-decaying final state particles (including $\pi^\pm$, $\gamma$, etc., but not including decaying particles like $K_S$ or $\pi^0$) are also kept:
 
-Different particles are differentiated using the postfix {\tt P(n)}, where {\tt (n)} is the number of the particle in the ordered list.  Four-momenta for secondaries originating from particle~{\tt (n)}, such as the two $\gamma$'s from a $\pi^0$, are recorded using {\tt P(n)a} and {\tt P(n)b}, where the ordering follows the same conventions as above, or, in the case of identical daughter particles, no ordering is assumed.  As two examples: in the process $\gamma p \to\pi^+\pi^-J/\psi p; J/\psi\to\mu^+\mu^-$, the raw energy of the $\pi^+$ is given by {\tt REnP4}; and in the process $\gamma p \to \pi^+\pi^-J/\psi p; J/\psi\to\pi^+\pi^-\pi^0$, the y-momentum of a photon from the $\pi^0$ decay, after the kinematic fit, is given by {\tt PyP6b}.
+\begin{verbatim}
+  (prefix)VxP(n):  x vertex of particle (n)
+  (prefix)VyP(n):  y vertex of particle (n)
+  (prefix)VzP(n):  z vertex of particle (n)
+\end{verbatim}
+
+Different types of four-momenta and vertices are distinguished using prefixes.  Raw (i.e. measured) values have a prefix~{\tt R}; final (i.e. kinematically fit) values have no prefix; MC truth information have a prefix~{\tt MC}.
+
+Different particles are differentiated using the postfix {\tt P(n)}, where {\tt (n)} is the number of the particle in the ordered list.  Four-momenta and vertices for secondaries originating from particle~{\tt (n)}, such as the two $\gamma$'s from a $\pi^0$, are recorded using {\tt P(n)a} and {\tt P(n)b}, where the ordering follows the same conventions as above, or, in the case of identical daughter particles, no ordering is assumed.  As two examples: in the process $\gamma p \to\pi^+\pi^-J/\psi p; J/\psi\to\mu^+\mu^-$, the raw energy of the $\pi^+$ is given by {\tt REnP4}; and in the process $\gamma p \to \pi^+\pi^-J/\psi p; J/\psi\to\pi^+\pi^-\pi^0$, the y-momentum of a photon from the $\pi^0$ decay, after the kinematic fit, is given by {\tt PyP6b}.