WCTrack.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////////
/// \file    WCTrack.h
/// \brief   Encapsulation of reconstructed Wire Chamber track.
///          Part of beamline reconstruction for NOvA test beam.
/// \author  Mike Wallbank (University of Cincinnati) <wallbank@fnal.gov>
/// \date    November 2018
////////////////////////////////////////////////////////////////////////////

#ifndef BRBWCTRACK_H
#define BRBWCTRACK_H

// framework

// stl
#include <iostream>

// root
#include "TVector2.h"
#include "TVector3.h"

// -----------------------------------------------------------------------
namespace brb {

  // -----------------------------------------------------------------------
  class WCTrack {

  public:

    WCTrack() {}
    WCTrack(float momentum, float y_kink, TVector3 delta_dist,
            TVector3 mag_entry, float dist_to_mag_axis,
            TVector2 xy_face, std::vector<TVector3> wcHits,
            TVector3 dir, float theta, float phi, float residual);

    float Momentum() const;
    float YKink() const;
    TVector3 DeltaDist() const;
    TVector3 MagnetEntryPoint() const;
    float TransDistToMagAxis() const;
    TVector2 XYFace() const;
    TVector3 Dir() const;
    float Theta() const;
    float Phi() const;
    float Residual() const;
    TVector3 WCHit(unsigned int wc) const;

    int WC(size_t hit) const;
    float HitWire(size_t hit) const;
    float HitPosition(int wc, int wire) const;
    size_t NHits() const;

  private:

    float fMomentum;               ///< Reconstructed momentum in the XZ plane (coord system origin is at secondary target)
    float fYKink;                  ///< Angle difference between upstream and downstream tracks
    TVector3 fDeltaDist;           ///< Distance between upstream and downstream track ends
    TVector3 fMagnetEntryPoint;    ///< X,Y,Z intersect of upstream WC track with front face of magnet
    float fTransDistToMagAxis;     ///< Transverse distance from WC track intercept to central axis of magnet
    TVector2 fXYFace;              ///< X and Y position of the track on the upstream face of the detector
    std::vector<TVector3> fWCHit;  ///< Hits on each chamber used to make the track
    TVector3 fDir;                 ///< Unit vector describing direction
    float fTheta;                  ///< Theta defined from the Z axis
    float fPhi;                    ///< Phi defined counterclockwise from the X axis
    float fResidual;               ///< Returns the goodness of fit to a linear regression for points used in track.

    // These are indexed by hit: each hit is represented by the same index in all three.
    std::vector<int> fWC;        ///< Defined from 1 to 4, like the Wire Chambers
    std::vector<float> fHitWire;
    float fHitPosition[4][3];    ///< WC is first index, dimension (x,y,z) as the second index. A [4][3] object.  

  };

} // namespace

#endif