SRParticleTruth.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \file    SRParticleTruth.h
////////////////////////////////////////////////////////////////////////
#ifndef SRPARTICLETRUTH_H
#define SRPARTICLETRUTH_H

#include "StandardRecord/SRLorentzVector.h"
#include "StandardRecord/SRVector3D.h"
#include "StandardRecord/SREnums.h"

namespace caf
{
  /// \brief The truth information of reco objects within a slice
  class SRParticleTruth 
    {
    public:
      SRParticleTruth();
      ~SRParticleTruth();

      int   pdg;         ///< PDG Code of the best matched truth particle
      int   motherpdg;   ///< PDG Code of the mother of the best matched truth particle
      float eff;         ///< True deposited energy collection efficiency for the best matched particle relative to the slice
      float pur;         ///< True deposited energy purity for the best matched particle
      
      SRLorentzVector p; ///< True energy 4-vector of the best matched particle
      SRLorentzVector motherp;///< True energy 4-vector of the mother particle
      int  trkID;        ///< GEANT trackId for particle
      SRVector3D start;  ///< Start point of true particle in detector coordinates (cm)
      /// vector containing pdg of mother particles.
      /// Index 0 is the immedieate mother, the last particle in the index
      /// will be a primary to the neutrino interaction.
      /// If this particle is a primary the vector will be empty
      std::vector<int > motherlist; 
      /// vector containing pdg of the immediate daughter
      /// particles. If there are no daughters the vector will be empty.
      std::vector<int > daughterlist; 
      float           visE;           ///< Visible Energy in detector, all summed FLSHits that made CellHits [GeV]
      float           visEinslc;      ///< Visible Energy in detector, slice summed FLSHits that made CellHits [GeV]
      float           daughterVisE;           ///< Visible Energy in detector for all daughters of this particle, all summed FLSHits that made CellHits [GeV]
      float           daughterVisEinslc;      ///< Visible Energy in detector for all daughters of this particle, slice summed FLSHits that made CellHits [GeV]
      float           visEBirks;           ///< Visible Energy in detector, all summed FLSHits that made CellHits [GeV] with birks suppression
      float           visEinslcBirks;      ///< Visible Energy in detector, slice summed FLSHits that made CellHits [GeV] with birks suppression
      float           daughterVisEBirks;           ///< Visible Energy in detector for all daughters of this particle, all summed FLSHits that made CellHits [GeV] with birks suppression
      float           daughterVisEinslcBirks;      ///< Visible Energy in detector for all daughters of this particle, slice summed FLSHits that made CellHits [GeV] with birks suppression
      
      
      std::vector<float> daughterVisEnergies;    ///< Energy of each particle contributing to the prong
      gen_process_t      processMax;   ///< The process conributing the most the prong
      float              processParticleE; ///< Energy of the particle causing the process that contirbuted the most
      std::vector<float> primNeutronE;    ///< Energy of the primary neutron that is linked to the prong, if one exists
      std::vector<float> primNeutronProcessE; ///<Fls energy contributing to the prong from neutronInelastic scatters with a proton/photon 
                                              ///<in the final state and elastic with proton in final state. Order (inelastic proton, inelastic photon, elastic proton)
      
    };
  
} // end namespace

#endif // SREPARTICLETRUTH_H
//////////////////////////////////////////////////////////////////////////////