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"

namespace anab {

  class Calorimetry{
  public:
    
    Calorimetry();

    double              fKineticEnergy;   
    std::vector<double> fdEdx;            
    std::vector<double> fdQdx;            
    std::vector<double> fResidualRange;   
    std::vector<double> fDeadWireResR;    
    double              fRange;           
    std::vector<double> fTrkPitch;        
    std::vector<TVector3> fXYZ;           

  private:
    geo::PlaneID        fPlaneID;

  public:

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

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

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

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

    const std::vector<double>& dEdx()          const; 
    const std::vector<double>& dQdx()          const; 
    const std::vector<double>& ResidualRange() const; 
    const std::vector<double>& DeadWireResRC() const; 
    const double&              KineticEnergy() const; 
    const double&              Range()         const; 
    double                     TrkPitchC()     const; 
    const std::vector<double>& TrkPitchVec()   const;
    const std::vector<TVector3>& XYZ()         const;
    const geo::PlaneID&        PlaneID()       const;

    
  };

}


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

#endif //ANAB_CALORIMETRY_H