TP0DReconModule.hxx

Go to the documentation of this file.

///
/// For questions or suggestions about this module please contact the
/// current responsible and CC in the eventAnalysis package manager.
///
/// 14-Feb-2012: Current responsible for this module is,
/// Ian Taylor (itaylor [*a*t*] nngroup.physics.sunysb.edu)
///
#ifndef TP0DReconModule_hxx_seen
#define TP0DReconModule_hxx_seen

#include <string>

#include <TClonesArray.h>
#include <TLorentzVector.h>
#include <TNamed.h>
#include <TPRegexp.h>
#include <TTree.h>
#include <TVector3.h>

#include <TND280Event.hxx>
#include <TReconCluster.hxx>
#include <TReconPID.hxx>
#include <TReconShower.hxx>
#include <TReconTrack.hxx>
#include <TReconVertex.hxx>

#include <TG4VHit.hxx>
#include <TMCDigit.hxx>
#include <TMCHit.hxx>

#include "TAnalysisReconModuleBase.hxx"

namespace ND {
class TP0DReconModule;
};

/// Used to fill a TTree of summary information for the P0D
/// reconstruction algorithms.
class ND::TP0DReconModule : public TAnalysisReconModuleBase {
 public:
  /// An internal Id and object type for each Reconstruction object.
  struct ObjectID {
    /// An internal Id for each Reconstruction object.
    /// These Ids are sequential from 0 for each object type within each event
    /// and correspond to their position within the respective TClonesArrays
    /// e.g.
    /// ND::TP0DReconModule::fNVertices.
    /// They are set within the various Fill<type>() methods e.g.
    /// ND::TP0DReconModule::FillVertexObject
    ///\dontinclude TP0DReconModule.cxx
    ///\skip ObjectID vertexID;
    ///\until vertexID.id
    /// and used by the, per event, reconstruction objects to keep track of
    /// which
    /// other objects relate to this one e.g.
    /// ND::TP0DReconModule::TP0DVertex::Particles.
    short id;
    enum OType_t {
      kBlank = 0,
      kVertex,
      kParticle,
      kTrack,
      kShower,
      kCluster,
      kHit
    };
    /// The reconstruction object type that this internal Id refers to.
    /// Not explicitly stored in output TTree. Used in
    /// ND::TP0DReconModule::FillReconObject to determine which Fill<type>()
    /// Method to call.
    ///\dontinclude TP0DReconModule.cxx
    ///\skip td::string reco_class_name
    ///\until }
    OType_t type;
    ObjectID() {
      id = -1;
      type = kBlank;
    }
  };

  //###################################################
  ///\short Used to summarise the various reconstruction algorithm results that
  /// p0dRecon employs.
  class TP0DAlgoRes : public TObject {
   public:
    virtual ~TP0DAlgoRes();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    /// The name of the p0dRecon algorithm that created this object.
    std::string AlgorithmName;

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< algorithm result.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this algorithm result.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< algorithm result.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< algorithm result.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this algorithm result.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DNode objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fNodes which hold Nodes related to this
    ///< algorithm result.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this algorithm result.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< algorithm result.
    short NHits;  ///< A count of the Hits associated with this object.

    //###################################################

   
    UInt_t UniqueID;///< UNUSED: A unique ID used for reconstruction objects.
    ///<\details This is not used for this class. Filled by
    ///< ND::TP0DReconModule::FillBaseObject which is not called by
    ///< ND::TP0DReconModule::FillAlgorithmResult.

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    std::string FullName;///< The name of the p0dRecon algorithm that corresponds to this
    ///< algorithm result summary.
    std::vector<short> AlgoResults;///< Holds internal Ids of daughter algorithm results to this one.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DAlgoRes objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fAlgoResults which hold algorithm
    ///< results related to this algorithm result.
    short Parent;///< The internal Id of the parent algorithm result of this one.
    ///<\details An ND::TP0DReconModule::ObjectID::id of an
    ///< ND::TP0DReconModule::TP0DAlgoRes which is equivalent to
    ///< the index of ND::TP0DReconModule::fAlgoResults which is the parent of
    ///< this one. If this is the root result defaults to '-1'.
    short UsedHitCluster;///< The internal Id of the cluster of hits which were used in a
    ///< reconstruction object by this algorithm.
    ///<\details An ND::TP0DReconModule::ObjectID::id of an
    ///< ND::TP0DReconModule::TP0DCluster which is equivalent to
    ///< the index of ND::TP0DReconModule::fClusters which is the used hit
    ///< 'selection'.
    short UnusedHitCluster;///< The internal Id of the cluster of hits which were not used in a
    ///< reconstruction object by this algorithm.
    ///<\details An ND::TP0DReconModule::ObjectID::id of an
    ///< ND::TP0DReconModule::TP0DCluster which is equivalent to
    ///< the index of ND::TP0DReconModule::fClusters which is the unused hit
    ///< 'selection'.

    ClassDef(TP0DReconModule::TP0DAlgoRes, 5);
  };

  //###################################################
  /// Contains a summary of the reconstruction information in a TReconVertex.
  class TP0DVertex : public TObject {
   public:
    virtual ~TP0DVertex();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    /// The name of the p0dRecon algorithm that created this object.
    std::string AlgorithmName;

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< reconstructed vertex.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this reconstructed vertex.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< reconstructed vertex.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< reconstructed vertex.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this reconstructed vertex.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DNode objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fNodes which hold Nodes related to this
    ///< reconstructed vertex.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this reconstructed vertex.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstructed vertex.
    short NHits;  ///< A count of the Hits associated with this object.

    //###################################################

    UInt_t UniqueID;///< The Unique ID of the corresponding ND::TReconVertex, used for  global-subdetector matching.
    ///<\details In ND::TP0DReconModule::FillBaseObject
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline basePtr->UniqueID

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    int Status;///< The reported Status of the corresponding ND::TReconVertex.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->Status
    ///< Use ND::TReconBase::StateBits to check the which status flags are set.
    float Quality;///< The reported reconstruction 'quality' of the corresponding  ND::TReconVertex.
    ///<\details Uses ND::TReconBase::GetQuality
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->Quality
    int NDOF;///< The Number of Degrees of Freedom in the reconstruction of the
    ///< corresponding ND::TReconVertex.
    ///<\details Retreived from ND::TReconBase::GetNDOF
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->NDOF

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DVertex::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the ND::THits which are constituents of this vertex.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DVertex::Truth_HitCount and TP0DVertex::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DVertex::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth  trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DVertex::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    TLorentzVector Position;///< Position of the TVertexState of the corresponding ND::TReconVertex.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TVertexState> vertexState
    ///<\skipline p0dVertex->Position
    TLorentzVector PosVariance;///< Variance on the position of the TVertexState of the corresponding  ND::TReconVertex.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TVertexState> vertexState
    ///<\skipline p0dVertex->PosVariance

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this vertex's position.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    float Fiducial;///< The distance to the outside of the fiducial volume
    ///<(negative if the vertex is outside).
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dVertex->Fiducial =
    ///<\until ;

    ClassDef(TP0DReconModule::TP0DVertex, 4);
  };

  /// A summary of the reconstruction information in a TReconPID.  Do NOT
  /// access nodes for a TP0DParticle.  Instead access the TP0DTrack or
  /// TP0DShower associated with the TP0DParticle and access the TP0DNodes
  /// associated with those objects. Currently, TP0DParticle's Nodes are not
  /// filled correctly.
  class TP0DParticle : public TObject {
   public:
    virtual ~TP0DParticle();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    std::string AlgorithmName;///< The name of the p0dRecon algorithm that created this object.

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this reconstructed particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< reconstructed particle.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this reconstructed  particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this reconstructed particle.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this reconstructed particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< reconstructed particle.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this reconstructed particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< reconstructed particle.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this reconstructed particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this reconstructed particle.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this reconstructed particle  (probably not what you want for a TP0DParticle).  This vector has the  index in the ND::TP0DReconModule::fNodes vector of the nodes  associated with this object. The nodes for a TP0DParticle don't hold  the expected information for most P0D analyses, so you probably should  be accessing the nodes of the constituent track.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this  reconstructed particle.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstructed particle.
        short NHits;  ///< A count of the Hits associated with this object.
    //###################################################

    UInt_t UniqueID;///< The Unique ID of the corresponding ND::TReconPID, used for  global-subdetector matching.
    ///<\details In ND::TP0DReconModule::FillBaseObject
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline basePtr->UniqueID

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    int Status;///< The reported Status of the corresponding ND::TReconPID.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->Status
    ///< Use ND::TReconBase::StateBits to check the which status flags are set.
    float Quality;///< The reported reconstruction 'quality' of the corresponding  ND::TReconPID.
    ///<\details Uses ND::TReconBase::GetQuality
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->Quality
    int NDOF;///< The Number of Degrees of Freedom in the reconstruction of the
    ///< corresponding ND::TReconPID.
    ///<\details Retreived from ND::TReconBase::GetNDOF
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dVertex->NDOF

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the     ND::THits which are constituents of this vertex.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DParticle::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DParticle::Truth_HitCount and TP0DParticle::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DParticle::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth  trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DParticle::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    float SideDeposit;///< The deposited charge in the side bars.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dParticle->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skip int layer  = ND::TGeomInfo::P0D().GetBar(geomid).GetLayer();
    ///<\until int bar
    ///<\skip if(bar < 4 || (layer == 0 && bar > 121)
    ///<\until }
    float EndDeposit;///< The deposited charge in the end p0dules.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dParticle->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skipline int p0dule = ND::TGeomInfo::P0D().GetBar(geomid).GetP0Dule();
    ///<\skip if(p0dule > 38){
    ///<\until }

    TLorentzVector Position;///< Position of the TPIDState of the corresponding ND::TReconPID.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TPIDState> particleState
    ///<\skipline p0dParticle->Position
    TLorentzVector PosVariance;///< Variance on the position of the TPIDState of the corresponding  ND::TReconPID.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TPIDState> particleState
    ///<\skipline p0dParticle->PosVariance

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this particle's position.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    TVector3 Direction;///< Direction of the TPIDState of the corresponding ND::TReconPID.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TPIDState> particleState
    ///<\skipline p0dParticle->Direction
    TVector3 DirVariance;///< Variance on the TP0DParticle::Direction of the TPIDState of the  corresponding ND::TReconPID.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TPIDState> particleState
    ///<\skipline p0dParticle->DirVariance

    float Momentum;///< The reconstructed momentum of the corresponding TReconPID.

    float Charge;  ///< The reconstructed charge of the corresponding TReconPID.

    std::vector<std::string> realPIDNames;///< Holds the variable names of ND::TRealDatum stored in the TReconPID.
    ///< The index of a given variable's name corresponds to the index of it's
    ///< value held in TP0DParticle::realPIDValues.
    std::vector<std::vector<float> > realPIDValues;///< Holds the PID variable values from ND::TRealDatum stored in the  TReconPID.
    ///<\details The index of a given variable's value vector corresponds to
    ///< the index of it's name in TP0DParticle::realPIDNames.

    std::vector<std::string> integerPIDNames;///< Holds the variable names of ND::TIntegerDatum stored in the TReconPID.
    ///< The index of a given variable's name corresponds to the index of it's
    ///< value held in TP0DParticle::integerPIDValues.
    std::vector<std::vector<short> > integerPIDValues;///< Holds the PID variable values from ND::TIntegerDatum stored in the TReconPID.
    ///<\details The index of a given variable's value vector corresponds to
    ///< the index of it's name in TP0DParticle::integerPIDNames.

    std::vector<short> PID;///< Potential PIDs (ND::TReconPID::ParticleId) matching with
    ///< TP0DParticle::PID_Weight values.

    std::vector<float> PID_weight;///< PID weights for each PID in TP0DParticle::PID.

    ClassDef(TP0DReconModule::TP0DParticle, 5);
  };

  /// Contains a summary of the reconstruction information in a TReconShower.
  class TP0DShower : public TObject {
   public:
    virtual ~TP0DShower();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    std::string AlgorithmName;///< The name of the p0dRecon algorithm that created this object.

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< reconstructed shower.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this reconstructed shower.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< reconstructed shower.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< reconstructed shower.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this reconstructed shower.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DNode objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fNodes which hold Nodes related to this
    ///< reconstructed shower.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this reconstructed shower.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstructed shower.
    short NHits;  ///< A count of the Hits associated with this object.

    //###################################################

    UInt_t UniqueID;///< The Unique ID of the corresponding ND::TReconShower, used for  global-subdetector matching.
    ///<\details In ND::TP0DReconModule::FillBaseObject
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline basePtr->UniqueID

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    int Status;///< The reported Status of the corresponding ND::TReconShower.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dShower->Status
    ///< Use ND::TReconBase::StateBits to check the which status flags are set.
    float Quality;///< The reported reconstruction 'quality' of the corresponding  ND::TReconShower.
    ///<\details Uses ND::TReconBase::GetQuality
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dShower->Quality
    int NDOF;///< The Number of Degrees of Freedom in the reconstruction of the corresponding ND::TReconShower.
    ///<\details Retreived from ND::TReconBase::GetNDOF
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dShower->NDOF

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the  ND::THits which are constituents of this shower.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DShower::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the ND::THits which are constituents of this shower.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DShower::Truth_HitCount and TP0DShower::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DShower::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DShower::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    float EDeposit;///< The total reconstructed energy deposit of the corresponding ND::TReconShower.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dShower->EDeposit

    float SideDeposit;///< The deposited charge in the side bars.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dShower->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skip int layer  = ND::TGeomInfo::P0D().GetBar(geomid).GetLayer();
    ///<\until int bar
    ///<\skip if(bar < 4 || (layer == 0 && bar > 121)
    ///<\until }
    float EndDeposit;///< The deposited charge in the end p0dules.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dShower->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skipline int p0dule = ND::TGeomInfo::P0D().GetBar(geomid).GetP0Dule();
    ///<\skip if(p0dule > 38){
    ///<\until }

    TLorentzVector Position;///< Position of the TShowerState of the corresponding ND::TReconShower.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TShowerState> showerState
    ///<\skipline p0dShower->Position
    TLorentzVector PosVariance;///< Variance on the position of the TShowerState of the corresponding  ND::TReconShower.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TShowerState> showerState
    ///<\skipline p0dShower->PosVariance

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this shower's position.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    TVector3 Direction;///< Direction of the TShowerState of the corresponding ND::TReconShower.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TShowerState> showerState
    ///<\skipline p0dPShower->Direction
    TVector3 DirVariance;///< Variance on the TP0DShower::Direction of the TShowerState of the corresponding ND::TReconShower.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TShowerState> showerState
    ///<\skipline p0dShower->DirVariance

    float Cone;///< The opening angle of the TShowerState cone of the corresponding  TReconShower.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TShowerState> showerState
    ///<\skipline p0dShower->Cone

    float Width;///< The 'width' of the shower, the extent perpendicular to the direction.
    ///< This is calculated here. Each of the THits, within the TReconShower,
    ///< contributes to the width based on it's, charged weighted,
    ///< perpendicular distance from the TReconShower's position. This is done
    ///< in 2D.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip for (hit = showerHits->begin();
    ///<\until {
    ///<\skipline - p0dShower->Position.Vect();
    ///<\skipline TVector3 para
    ///<\skipline TVector3 perp = diff - para;
    ///<\skip if ((*hit)->IsXHit()) {
    ///<\until p0dShower->Width
    ///<\skip p0dShower->SideDeposit
    ///<\skip }
    ///<\skip }
    ///<\skipline }
    ///<\skipline p0dShower->Width /= std::max(totalCharge,1.0);

    float Length;///< The RMS length of the shower, the extent parallel to the direction.
    ///< This is calculated here. Each of the THits, within the TReconShower,
    ///< contributes to the width based on it's, charged weighted,
    ///< perpendicular distance from the TReconShower's position. This is done
    ///< in 2D.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip for (hit = showerHits->begin();
    ///<\until {
    ///<\skipline - p0dShower->Position.Vect();
    ///<\skipline TVector3 para
    ///<\skip if ((*hit)->IsXHit()) {
    ///<\until totalCharge += q;
    ///<\skipline avgLen
    ///<\skipline avgLen2
    ///<\skip p0dShower->SideDeposit
    ///<\skip }
    ///<\skip }
    ///<\skipline }
    ///<\skip avgLen
    ///<\until p0dShower->Length

    ClassDef(TP0DReconModule::TP0DShower, 5);
  };

  /// Contains a summary of the reconstruction information in a TReconCluster.
  class TP0DCluster : public TObject {
   public:
    virtual ~TP0DCluster();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    std::string AlgorithmName;///< The name of the p0dRecon algorithm that created this object.

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< reconstructed cluster.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this reconstructed cluster.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< reconstructed cluster.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< reconstructed cluster.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this reconstructed cluster.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DNode objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fNodes which hold Nodes related to this
    ///< reconstructed cluster.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this reconstructed cluster.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstructed cluster.
    short NHits;  ///< A count of the Hits associated with this object.

    //###################################################

    UInt_t UniqueID;///< The Unique ID of the corresponding ND::TReconCluster, used for  global-subdetector matching.
    ///<\details In ND::TP0DReconModule::FillBaseObject
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline basePtr->UniqueID

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DCluster::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DCluster::Truth_HitCount and TP0DCluster::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DCluster::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth  trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DCluster::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    short NFiducialHits;///< The number of THits in the cluster which are within the fiducial volume.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip float fiducial
    ///<\until ;
    ///<\skipline p0dCluster->NFiducialHits

    float EDeposit;///< The total reconstructed energy deposit of the corresponding  ND::TReconCluster.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dCluster->EDeposit

    TLorentzVector Position;///< Position of the corresponding ND::TReconCluster.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dCluster->Position
    TLorentzVector PosVariance;///< Variance on the position of the corresponding  ND::TReconCluster.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dCluster->PosVariance

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this cluster's positon.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    static const int arraySize = 9;///< Size of the TP0DCluster::Moments array.
    float Moments[arraySize];///< Moments of the cluster stored as a flat 3x3 matrix.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline TMatrixF moments
    ///<\skip std::copy
    ///<\until ;

    ClassDef(TP0DReconModule::TP0DCluster, 7);
  };

  ///< Contains a summary of the reconstruction information in a TReconTrack.
  class TP0DTrack : public TObject {
   public:
    virtual ~TP0DTrack();

    //###################################################
    // These items are used by templated code, and must be the same
    // for all of: P0DAlgoRes, P0DVertex, P0DParticle, P0DTrack,
    // P0DShower and P0DCluster.
    //###################################################

    std::string AlgorithmName;///< The name of the p0dRecon algorithm that created this object.

    //###################################################

    short Cycle;  ///< The cycle number, based on the first hit.

    //###################################################

    std::vector<short> Vertices;///< Holds internal Ids of Vertices pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DVertex objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fVertices which hold Vertices related to this
    ///< reconstructed track.
    std::vector<short> Particles;///< Holds internal Ids of Particles pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DParticle objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fParticles which hold Particles related
    ///< to this reconstructed track.
    std::vector<short> Tracks;///< Holds internal Ids of Tracks pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DTrack objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fTracks which hold Tracks related to this
    ///< reconstructed track.
    std::vector<short> Showers;///< Holds internal Ids of Showers pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DShower objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fShowers which hold Showers related to this
    ///< reconstructed track.
    std::vector<short> Clusters;///< Holds internal Ids of Clusters pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DCluster objects which are equivalent to
    ///< indicies of ND::TP0DReconModule::fClusters which hold Clusters related
    ///< to this reconstructed track.
    std::vector<short> Nodes;///< Holds internal Ids of Nodes pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DNode objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fNodes which hold Nodes related to this
    ///< reconstructed track.
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this reconstructed track.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstructed track.
    short NHits;  ///< A count of the Hits associated with this object.

    //###################################################

    UInt_t UniqueID;///< The Unique ID of the corresponding ND::TReconTrack, used for  global-subdetector matching.
    ///<\details In ND::TP0DReconModule::FillBaseObject
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline basePtr->UniqueID

    //###################################################
    // Objects below here are unique to their individual classes.
    //###################################################

    int Status;///< The reported Status of the corresponding ND::TReconTrack.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dTrack->Status
    ///< Use ND::TReconBase::StateBits to check the which status flags are set.
    float Quality;///< The reported reconstruction 'quality' of the corresponding  ND::TReconTrack.
    ///<\details Uses ND::TReconBase::GetQuality
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dTrack->Quality
    int NDOF;///< The Number of Degrees of Freedom in the reconstruction of the  corresponding ND::TReconTrack.
    ///<\details Retreived from ND::TReconBase::GetNDOF
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dTrack->NDOF

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the ND::THits which are constituents of this track.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DTrack::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the  ND::THits which are constituents of this track.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DTrack::Truth_HitCount and TP0DTrack::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DTrack::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth  trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DTrack::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    float EDeposit;///< The total reconstructed energy deposit of the corresponding  ND::TReconTrack.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline p0dTrack->EDeposit

    float SideDeposit;///< The deposited charge in the side bars.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dTrack->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skip int layer  = ND::TGeomInfo::P0D().GetBar(geomid).GetLayer();
    ///<\until int bar
    ///<\skip if(bar < 4 || (layer == 0 && bar > 121)
    ///<\until }
    float EndDeposit;///< The deposited charge in the end p0dules.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip p0dTrack->SideDeposit
    ///<\skipline ND::TGeometryId geomid = (*hit)->GetGeomId();
    ///<\skipline int p0dule = ND::TGeomInfo::P0D().GetBar(geomid).GetP0Dule();
    ///<\skip if(p0dule > 38){
    ///<\until }

    TLorentzVector Position;///< Position of the TTrackState of the corresponding ND::TReconTrack.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TTrackState> trackState
    ///<\skipline p0dTrack->Position
    TLorentzVector PosVariance;///< Variance on the position of the TTrackState of the corresponding  ND::TReconTrack.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TTrackState> trackState
    ///<\skipline p0dTrack->PosVariance

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this track's position.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    TVector3 Direction;///< Direction of the TTrackState of the corresponding ND::TReconTrack.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TTrackState> trackState
    ///<\skipline p0dPTrack->Direction
    TVector3 DirVariance;///< Variance on the TP0DTrack::Direction of the TTrackState of the  corresponding ND::TReconTrack.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TTrackState> trackState
    ///<\skipline p0dTrack->DirVariance

    float Length;///< The length of the track calculated by summing over the distances
    ///< between the constituent TReconNodes.

    ClassDef(TP0DReconModule::TP0DTrack, 6);
  };

  ///< Contains a summary of the reconstruction information in a TReconNode.
  class TP0DNode : public TObject {
   public:
    virtual ~TP0DNode();
    std::vector<short> Hits;///< Holds internal Ids of Hits pertinent to this reconstruction node.
    ///<\details Holds ND::TP0DReconModule::ObjectID::id of
    ///< ND::TP0DReconModule::TP0DHit objects which are equivalent to indicies
    ///< of ND::TP0DReconModule::fHits which hold Hits related to this
    ///< reconstruction node.

    // Truth level hit information.  For the MC hits in an object we
    // store: the true trajectory IDs, the number of hits contributed
    // by each and the charge share.

    std::vector<int> Truth_PrimaryTrajIDs;///< The vector of primary true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details These are
    ///< 'pseudo-primaries', which may or may not be primaries for a
    ///< vertex.  It is the 'interesting' primary that is defined by
    ///< nd280Geant4Sim, and can include particles such as photons from pi0
    ///< decay and michel electrons from muon decay.  If you want to be
    ///< sure of finding a vertex, use the TP0DNode::Truth_TrajIDs list and
    ///< get primary ID from the trajectories.
    ///< See ND::TP0DReconModule::HitTruthPrimaryInfo for how they are
    ///< calculated.

    std::vector<int> Truth_TrajIDs;///< The vector of true trajectory IDs which contribute to the  ND::THits which are constituents of this vertex.
    ///<\details The index of a given trajectory id within this vector is used
    ///< as the index for the corresponding information held in
    ///< TP0DNode::Truth_HitCount and TP0DNode::Truth_ChargeShare.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<short> Truth_HitCount;///< The number of THits that each truth trajectory contributed to.
    ///<\details This doesn't have to equal the
    ///< number of hits, hits can be shared and some hits may not have
    ///< associated trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DNode::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    std::vector<float> Truth_ChargeShare;///< The weighted reconstructed charge deposited by each truth  trajectory.
    ///<\details For each THit, a trajectory will
    ///< receive a proportional share of the charge, based on the total
    ///< true charge deposited by all contributing trajectories.
    ///< Indexed by the position of the relevant trajectory ID in
    ///< TP0DNode::Truth_TrajIDs.
    ///< See ND::TP0DReconModule::HitTruthInfo for how they are
    ///< calculated.

    float EDeposit;///< The total reconstructed energy deposit of the  ND::TMEDepositState of the relevant ND::TReconNode.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skip ND::THandle<ND::TReconBase> nodeObject = node->GetObject();
    ///<\until p0dNode->EDeposit
    ///<\skipline }
    ///<\skipline }
    ///<\skip p0dNode->Truth_PrimaryTrajIDs
    ///<\skip }
    ///<\skip }
    ///<\skipline }

    TLorentzVector Position;///< Position of the TMPositionDirectionState of the corresponding  ND::TReconNode.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TReconState> s1 = node->GetState();
    ///<\skip posDirState =
    ///<\until ND::GetPointer(s1));
    ///<\skipline p0dNode->Position
    TLorentzVector PosVariance;///< Variance on the position of the TMPositionDirectionState of the  corresponding ND::TReconPID.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TReconState> s1 = node->GetState();
    ///<\skip posDirState =
    ///<\until ND::GetPointer(s1));
    ///<\skipline p0dNode->PosVar

    short ValidDimensions;///< Bit flag demarkating valid dimensions for this node's position.
    ///< values(valid dimensions) : 1(x),2(y),4(z),3(xy),5(xz),6(yz),7(xyz)

    TVector3 Direction;///< Direction of the TMPositionDirectionState of the corresponding  ND::TReconNode.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TReconState> s1 = node->GetState();
    ///<\skip posDirState =
    ///<\until );
    ///<\skipline p0dNode->Direction
    TVector3 DirVariance;///< Variance on the TP0DNode::Direction of the  TMPositionDirectionState of the corresponding ND::TReconNode.
    ///<\details
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TReconState> s1 = node->GetState();
    ///<\skip posDirState =
    ///<\until );
    ///<\skipline p0dNode->DirVar

    ClassDef(TP0DReconModule::TP0DNode, 5);
  };

  ///< Contains a summary of the reconstruction information in a TSingleHit.
  class TP0DHit : public TObject {
   public:
    virtual ~TP0DHit() {}

    UInt_t GeomID;  ///< Geometry Id of the TSingleHit.

    UInt_t ChanID;  ///< Channel Id of the TSingleHit.

    float Charge;///< Reconstructed hit charge without attenuation correction from the  corresponding TSingleHit.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TSingleHit> sh
    ///<\skipline p0dHit->Charge

    float Time;///< Reconstructed hit time from the corresponding TSingleHit.
    ///<\dontinclude TP0DReconModule.cxx
    ///<\skipline ND::THandle<ND::TSingleHit> sh
    ///<\skipline p0dHit->Time

    ClassDef(TP0DReconModule::TP0DHit, 1);
  };

 public:
  TP0DReconModule(const char *name = "P0D",
                  const char *title = "P0D Recon Module");
  virtual ~TP0DReconModule();

  virtual EType GetTreeType() const { return kRecon; }

  virtual Bool_t ProcessFirstEvent(ND::TND280Event &);

 protected:
  virtual void InitializeModule();

  virtual void InitializeBranches();
  virtual bool FillTree(ND::TND280Event &);

  virtual short FillAlgorithmResult(const ND::THandle<ND::TAlgorithmResult>,
                                    short);

  template <class T>
  void FillBaseObject(T basePtr, ND::THandle<ND::TReconBase> baseObject,
                      bool saveHits);

  virtual ObjectID FillReconObject(const ND::THandle<ND::TReconBase>,
                                   bool saveHits);

  virtual ObjectID FillVertexObject(const ND::THandle<ND::TReconVertex>,
                                    bool saveHits);

  virtual ObjectID FillParticleObject(const ND::THandle<ND::TReconPID>,
                                      bool saveHits);

  virtual ObjectID FillTrackObject(const ND::THandle<ND::TReconTrack>,
                                   bool saveHits);

  virtual ObjectID FillShowerObject(const ND::THandle<ND::TReconShower>,
                                    bool saveHits);

  virtual ObjectID FillClusterObject(const ND::THandle<ND::TReconCluster>,
                                     bool saveHits);

  virtual std::map<int, std::pair<int, float> > HitTruthInfo(
      const ND::THandle<ND::THitSelection>);
  virtual std::vector<int> HitTruthPrimaryInfo(
      const ND::THandle<ND::THitSelection>);

  virtual short FillNode(const ND::THandle<ND::TReconNode>, bool saveHits);

  short FillHit(const ND::THandle<ND::THit>);

  short GetCycle(const ND::THandle<ND::THitSelection>);

  short CountValidDimensions(TVector3 posVar);

 public:
  // Tree Entries
  Short_t fNAlgoResults;  ///<  The number of Algorithm Results
  Short_t fNVertices;     ///<  The number of added vertices
  Short_t fNParticles;    ///<  The number of particles
  Short_t fNShowers;      ///<  The number of showers
  Short_t fNTracks;       ///<  The number of tracks
  Short_t fNNodes;        ///<  The number of track nodes
  Short_t fNClusters;     ///<  The number of clusters
  Short_t fNHits;         ///<  The number of hits that are saved

  TClonesArray
      *fAlgoResults;         ///<  The TP0DAlgoRes vector of Algorithm Results.
  TClonesArray *fVertices;   ///<  The TP0DVertex vector of vertices.
  TClonesArray *fParticles;  ///<  The TP0DParticle vector of particles.
  TClonesArray *fShowers;    ///<  The TP0DShower vector of showers.
  TClonesArray *fTracks;     ///<  The TP0DTrack vector of tracks.
  TClonesArray *fNodes;      ///<  The TP0DNode vector of nodes.
  TClonesArray *fHits;       ///<  The TP0DHit vector of hits.
  TClonesArray *fClusters;   ///<  The TP0DCluster vector of clusters

 private:
  ///\short If an incoming TAlgorithmResult's name matches a pattern within this
  /// vector it will not be summarised and stored.
  std::vector<TPRegexp> fRejectAlgoResultList;

  /// Scratch space to keep track of already summarised THits.
  std::map<UInt_t, short> fTempHitMap;

};

#endif