/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus
*
* MAUS is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* MAUS is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with MAUS. If not, see .
*
*/
/** @class SciFiBasePRTrack Base class for scifi pattern recogntion tracks. */
#ifndef _SRC_COMMON_CPP_DATASTRUCTURE_SCIFIBASEPRTRACK_HH_
#define _SRC_COMMON_CPP_DATASTRUCTURE_SCIFIBASEPRTRACK_HH_
// C++ headers
#include
// ROOT headers
#include "TRef.h"
#include "TRefArray.h"
#include "TObject.h"
// MAUS headers
#include "src/common_cpp/DataStructure/SciFiSpacePoint.hh"
#include "src/common_cpp/Utils/VersionNumber.hh"
namespace MAUS {
typedef std::vector DoubleArray;
class SciFiBasePRTrack : public TObject {
public:
/** Default constructor */
SciFiBasePRTrack();
/** Construct some variables */
explicit SciFiBasePRTrack(DoubleArray cov);
/** Construct some variables */
SciFiBasePRTrack(DoubleArray cov, SciFiSpacePointPArray spoints);
/** Copy constructor */
SciFiBasePRTrack(const SciFiBasePRTrack& track);
/** Default destructor */
virtual ~SciFiBasePRTrack();
/** Assignment operator */
SciFiBasePRTrack& operator=(const SciFiBasePRTrack& track);
/** Get the covariance matrix elements of the least square fit parameters */
std::vector get_covariance() const { return _covariance; }
/** Set the covariance matrix elements of the least square fit parameters */
void set_covariance(const std::vector& covariance) { _covariance = covariance; }
/** Get the vector holding pointers to the spacepoints used by the track */
TRefArray* get_spacepoints() const { return _spoints; }
/** Set the vector holding pointers to the spacepoints
* used by the track from a TRefArray */
void set_spacepoints(TRefArray* spoints);
/** Return a vector holding pointers to the spacepoints
* (requires extracting them from the TRefArray holder
*/
SciFiSpacePointPArray get_spacepoints_pointers() const;
/** Set the vector holding pointers to the spacepoints
* used by the track from a vector of spacepoint pointers */
void set_spacepoints_pointers(const SciFiSpacePointPArray &spoints);
/** Set the reconstructed reference frame position */
void set_reference_position(ThreeVector pos) { _position = pos; }
/** Get the reconstructed reference frame position */
ThreeVector get_reference_position() const { return _position; }
/** Set the reconstructed reference frame momentum */
void set_reference_momentum(ThreeVector mom) { _momentum = mom; }
/** Get the reconstructed reference frame momentum */
ThreeVector get_reference_momentum() const { return _momentum; }
/** Get the global chi-squared value */
double get_chi_squared() const { return _chi_sq; }
/** Set the global chi-squared value */
void set_chi_squared(double chi_sq) { _chi_sq = chi_sq; }
/** Get the Number of Degrees of Freedom */
int get_ndf() const { return _ndf; }
/** Set the Number of Degrees of Freedom */
void set_ndf(int ndf) { _ndf = ndf; }
private:
TRefArray* _spoints;
/**
* Vector of covariances from the least squared fits.
* In the helical track case these are the covariances of:
* alpha-alpha, alpha-beta, alpha-gamma, beta-alpha, beta-beta, beta-gamma,
* gamma-alpha, gamma-beta, gamma-gamma, c_sz-c_sz, c_sz-m_sz, m_sz-c_sz, m_sz-m_sz.
* In the straight track case these are the covariances of:
* c_x-c_x, c_x-m_x, m_x-c_x, m_x-m_x, c_y-c_y, c_y-m_y, m_y-c_y, m_y-m_y.
*/
DoubleArray _covariance; // The covariance matrix derived from the fit
ThreeVector _position; // Reconstructed reference plane position
ThreeVector _momentum; // Reconstructed reference plane momentum
double _chi_sq;
int _ndf;
MAUS_VERSIONED_CLASS_DEF(SciFiBasePRTrack)
};
typedef std::vector SciFiBasePRTrackPArray;
typedef std::vector SciFiBasePRTrackArray;
} // ~namespace MAUS
#endif