Calorimetry.h

Go to the documentation of this file.

// \version 
//
// \brief Definition of data product to hold Calorimetry information
//
// \author brebel@fnal.gov, tjyang@fnal.gov
//
#ifndef ANAB_CALORIMETRY_H
#define ANAB_CALORIMETRY_H

#include <vector>
#include <iosfwd>
#include <iostream>
#include <iomanip>
#include <TVector3.h>

#include "larcoreobj/SimpleTypesAndConstants/geo_types.h"
#include "lardataobj/RecoBase/TrackingTypes.h"

namespace anab {

  using Point_t = recob::tracking::Point_t;

  class Calorimetry{
  public:
    
    Calorimetry();

    float                fKineticEnergy;   
    std::vector<float>   fdEdx;            
    std::vector<float>   fdQdx;            
    std::vector<float>   fResidualRange;   
    std::vector<float>   fDeadWireResR;    
    float                fRange;           
    std::vector<float>   fTrkPitch;        
    std::vector<Point_t> fXYZ;             
    std::vector<size_t>  fTpIndices;       

  private:
    geo::PlaneID        fPlaneID;

  public:

    Calorimetry(float KinematicEnergy,
                std::vector<float> const& dEdx,
                std::vector<float> const& dQdx,
                std::vector<float> const& resRange,
                std::vector<float> const& deadwire,
                float Range,
                float TrkPitch,
                geo::PlaneID planeID);

    Calorimetry(float KineticEnergy,
                std::vector<float> const& dEdx,
                std::vector<float> const& dQdx,
                std::vector<float> const& resRange,
                std::vector<float> const& deadwire,
                float Range,
                std::vector<float> const& TrkPitch,
                geo::PlaneID planeID);

    Calorimetry(float KineticEnergy,
                std::vector<float> const& dEdx,
                std::vector<float> const& dQdx,
                std::vector<float> const& resRange,
                std::vector<float> const& deadwire,
                float Range,
                std::vector<float> const& TrkPitch,
                std::vector<Point_t> const& XYZ,
                geo::PlaneID planeID);

    Calorimetry(float KineticEnergy,
                std::vector<float> const& dEdx,
                std::vector<float> const& dQdx,
                std::vector<float> const& resRange,
                std::vector<float> const& deadwire,
                float Range,
                std::vector<float> const& TrkPitch,
                std::vector<Point_t> const& XYZ,
                std::vector<size_t> const& TpIndices,
                geo::PlaneID planeID);

    friend std::ostream& operator << (std::ostream &o, Calorimetry const& a);

    const std::vector<float>& dEdx()          const;
    const std::vector<float>& dQdx()          const;
    const std::vector<float>& ResidualRange() const;
    const std::vector<float>& DeadWireResRC() const;
    const float&              KineticEnergy() const;
    const float&              Range()         const;
    float                     TrkPitchC()     const;
    const std::vector<float>& TrkPitchVec()   const;
    const std::vector<Point_t>& XYZ()         const;
    const std::vector<size_t>& TpIndices()    const;
    const geo::PlaneID&        PlaneID()       const;

    
  };

}


inline const std::vector<float>& anab::Calorimetry::dEdx()          const { return fdEdx;          }
inline const std::vector<float>& anab::Calorimetry::dQdx()          const { return fdQdx;          }
inline const std::vector<float>& anab::Calorimetry::ResidualRange() const { return fResidualRange; }
inline const std::vector<float>& anab::Calorimetry::DeadWireResRC() const { return fDeadWireResR;  }
inline const float&              anab::Calorimetry::KineticEnergy() const { return fKineticEnergy; }
inline const float&              anab::Calorimetry::Range()         const { return fRange;         }
inline const std::vector<float>& anab::Calorimetry::TrkPitchVec()   const { return fTrkPitch;      }
inline float                     anab::Calorimetry::TrkPitchC()     const
{
  if (fTrkPitch.size())
    return fTrkPitch[0];
  else return 0;
}
inline const std::vector<anab::Point_t>& anab::Calorimetry::XYZ()         const { return fXYZ;      }
inline const std::vector<size_t>&        anab::Calorimetry::TpIndices()   const { return fTpIndices;}
inline const geo::PlaneID& anab::Calorimetry::PlaneID()              const { return fPlaneID; }

#endif //ANAB_CALORIMETRY_H