#include "Calorimetry.h"

Public Member Functions
	Calorimetry ()

	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)

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

Public Attributes
float	fKineticEnergy
	determined kinetic energy More...

std::vector< float >	fdEdx
	dE/dx, should be same size as fResidualRange More...

std::vector< float >	fdQdx
	dQ/dx More...

std::vector< float >	fResidualRange
	range from end of track More...

std::vector< float >	fDeadWireResR
	dead wire residual range, collection plane More...

float	fRange
	total range of track More...

std::vector< float >	fTrkPitch
	track pitch on collection plane More...

std::vector< Point_t >	fXYZ
	coordinates of space points; for a discussion on the object type for coordinates see recob::tracking::Coord_t. More...

std::vector< size_t >	fTpIndices
	indices of original trajectory points on track More...

Private Attributes
geo::PlaneID	fPlaneID

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

Detailed Description

Definition at line 22 of file Calorimetry.h.

Constructor & Destructor Documentation

anab::Calorimetry::Calorimetry ( )

Definition at line 17 of file Calorimetry.cxx.

References fDeadWireResR, fdEdx, fdQdx, fResidualRange, fTpIndices, fTrkPitch, and fXYZ.

                            : fKineticEnergy(0.), fRange(0.)
   {
     fdEdx.clear();
     fdQdx.clear();
     fResidualRange.clear();
     fDeadWireResR.clear();
     fTrkPitch.clear();
     fXYZ.clear();
     fTpIndices.clear();
   }

anab::Calorimetry::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
	)

Definition at line 29 of file Calorimetry.cxx.

References fDeadWireResR, fdEdx, fdQdx, fKineticEnergy, fPlaneID, fRange, fResidualRange, fTrkPitch, fXYZ, KineticEnergy(), and Range().

   {
 
     fKineticEnergy = KineticEnergy;
     fRange = Range;
     for (size_t i = 0; i != dQdx.size(); ++i) {
       fTrkPitch.push_back(TrkPitch);
       fXYZ.push_back({-999., -999., -999.});
     }
     if (dEdx.size() != resRange.size())
       throw cet::exception("anab::Calorimetry") << "dE/dx and residual range vectors "
                                                 << "have different sizes, this is a problem.\n";
     fdEdx.resize(dEdx.size());
     fdQdx.resize(dQdx.size());
     fResidualRange.resize(resRange.size());
     for (size_t i = 0; i < dEdx.size(); ++i) {
       fdEdx[i] = dEdx[i];
       fdQdx[i] = dQdx[i];
       fResidualRange[i] = resRange[i];
     }
 
     fDeadWireResR.resize(deadwire.size());
     for (size_t i = 0; i < deadwire.size(); ++i) {
       fDeadWireResR[i] = deadwire[i];
     }
 
     fPlaneID = planeID;
   }

anab::Calorimetry::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
	)

Definition at line 66 of file Calorimetry.cxx.

References fDeadWireResR, fdEdx, fdQdx, fKineticEnergy, fPlaneID, fRange, fResidualRange, fTrkPitch, fXYZ, KineticEnergy(), and Range().

   {
 
     fPlaneID = planeID;
     fKineticEnergy = KineticEnergy;
     fRange = Range;
     fTrkPitch = TrkPitch;
     if (dEdx.size() != resRange.size())
       throw cet::exception("anab::Calorimetry") << "dE/dx and residual range vectors "
                                                 << "have different sizes, this is a problem.\n";
     for (size_t i = 0; i != dQdx.size(); ++i) {
       fXYZ.push_back({-999., -999., -999.});
     }
     fdEdx.resize(dEdx.size());
     fdQdx.resize(dQdx.size());
     fResidualRange.resize(resRange.size());
     for (size_t i = 0; i < dEdx.size(); ++i) {
       fdEdx[i] = dEdx[i];
       fdQdx[i] = dQdx[i];
       fResidualRange[i] = resRange[i];
     }
 
     fDeadWireResR.resize(deadwire.size());
     for (size_t i = 0; i < deadwire.size(); ++i) {
       fDeadWireResR[i] = deadwire[i];
     }
   }

anab::Calorimetry::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
	)

Definition at line 102 of file Calorimetry.cxx.

     : Calorimetry(KineticEnergy,
                   dEdx,
                   dQdx,
                   resRange,
                   deadwire,
                   Range,
                   TrkPitch,
                   XYZ,
                   std::vector<size_t>(),
                   planeID)
   {}

anab::Calorimetry::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
	)

Definition at line 123 of file Calorimetry.cxx.

References dEdx(), dQdx(), fDeadWireResR, fdEdx, fdQdx, fKineticEnergy, fPlaneID, fRange, fResidualRange, fTpIndices, fTrkPitch, fXYZ, KineticEnergy(), Range(), TpIndices(), and XYZ().

   {
 
     if (dEdx.size() != resRange.size() || dEdx.size() != dQdx.size() ||
         dEdx.size() != TrkPitch.size() || dEdx.size() != XYZ.size() ||
         (TpIndices.size() > 0 && dEdx.size() != TpIndices.size()))
       throw cet::exception("anab::Calorimetry") << "Input vectors "
                                                 << "have different sizes, this is a problem.\n";
     fPlaneID = planeID;
     fKineticEnergy = KineticEnergy;
     fRange = Range;
     fTrkPitch = TrkPitch;
     fdEdx = dEdx;
     fdQdx = dQdx;
     fResidualRange = resRange;
     fXYZ = XYZ;
     fTpIndices = TpIndices;
     fDeadWireResR = deadwire;
   }

Member Function Documentation

const std::vector< float > & anab::Calorimetry::DeadWireResRC ( ) const

inline

Definition at line 109 of file Calorimetry.h.

References fDeadWireResR.

Referenced by pid::Chi2PIDAlg::DoParticleID().

 {
   return fDeadWireResR;
 }

const std::vector< float > & anab::Calorimetry::dEdx ( ) const

inline

Definition at line 97 of file Calorimetry.h.

References fdEdx.

Referenced by Calorimetry(), pid::Chi2PIDAlg::DoParticleID(), and pid::PIDAAlg::RunPIDAAlg().

 {
   return fdEdx;
 }

const std::vector< float > & anab::Calorimetry::dQdx ( ) const

inline

Definition at line 101 of file Calorimetry.h.

References fdQdx.

Referenced by Calorimetry().

 {
   return fdQdx;
 }

const float & anab::Calorimetry::KineticEnergy ( ) const

inline

Definition at line 113 of file Calorimetry.h.

References fKineticEnergy.

Referenced by Calorimetry(), and pid::PIDAAlg::FillPIDAProperties().

 {
   return fKineticEnergy;
 }

const geo::PlaneID & anab::Calorimetry::PlaneID ( ) const

inline

Definition at line 140 of file Calorimetry.h.

References fPlaneID.

Referenced by pid::Chi2PIDAlg::DoParticleID(), and pid::PIDAAlg::FillPIDAProperties().

 {
   return fPlaneID;
 }

const float & anab::Calorimetry::Range ( ) const

inline

Definition at line 117 of file Calorimetry.h.

References fRange.

Referenced by Calorimetry(), and pid::PIDAAlg::FillPIDAProperties().

 {
   return fRange;
 }

const std::vector< float > & anab::Calorimetry::ResidualRange ( ) const

inline

Definition at line 105 of file Calorimetry.h.

References fResidualRange.

Referenced by pid::Chi2PIDAlg::DoParticleID(), and pid::PIDAAlg::RunPIDAAlg().

 {
   return fResidualRange;
 }

const std::vector< size_t > & anab::Calorimetry::TpIndices ( ) const

inline

Definition at line 136 of file Calorimetry.h.

References fTpIndices.

Referenced by Calorimetry().

 {
   return fTpIndices;
 }

float anab::Calorimetry::TrkPitchC ( ) const

inline

Definition at line 125 of file Calorimetry.h.

References fTrkPitch.

 {
   if (fTrkPitch.size())
     return fTrkPitch[0];
   else
     return 0;
 }

const std::vector< float > & anab::Calorimetry::TrkPitchVec ( ) const

inline

Definition at line 121 of file Calorimetry.h.

References fTrkPitch.

 {
   return fTrkPitch;
 }

const std::vector< anab::Point_t > & anab::Calorimetry::XYZ ( ) const

inline

Definition at line 132 of file Calorimetry.h.

References fXYZ.

Referenced by Calorimetry().

 {
   return fXYZ;
 }

Friends And Related Function Documentation

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

friend

Definition at line 155 of file Calorimetry.cxx.

   {
     o << "Kinetic Energy: " << a.fKineticEnergy << "\n Range: " << a.fRange << std::endl;
 
     for (size_t n = 0; n < a.fdEdx.size(); ++n)
       o << "dE/dx=" << a.fdEdx[n] << " Residual range=" << a.fResidualRange[n]
         << " dQ/dx=" << a.fdQdx[n] << " (x,y,z)=(" << a.fXYZ[n].X() << "," << a.fXYZ[n].Y() << ","
         << a.fXYZ[n].Z() << ")"
         << " pitch=" << a.fTrkPitch[n] << " planeID=(" << a.fPlaneID.Cryostat << ","
         << a.fPlaneID.TPC << "," << a.fPlaneID.Plane << ")" << std::endl;
 
     return o;
   }