mono::Track3D Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-26/Monopole/Track3D.h"

Public Member Functions
	Track3D ()
	Track3D (Track2D const &xz_track, Track2D const &yz_track)
TVector3	start () const
TVector3	velocity () const
double	dxdz () const
double	dydz () const
double	dzdt (geo::View_t view) const
double	r2_xt () const
double	r2_yt () const
double	r2_zt (geo::View_t view) const
double	beta () const
double	max_time_gap_xz () const
double	max_time_gap_yz () const
rb::Cluster	cluster () const

Private Member Functions
double	get_covariance (double const &a, double const &mean_a, double const &b, double const &mean_b) const
double	get_variance (double const &a, double const &mean_a) const
bool	time_fit (Track2D const &track, double &dvdt, double &dzdt, double &r2_vt, double &r2_zt) const
double	time_gap (art::PtrVector< rb::CellHit > const &hits) const

Private Attributes
double	default_value_
double	x0_
double	y0_
double	z0_
double	dxdz_
double	dydz_
double	dxdt_
double	dydt_
double	dzdt_
double	dzdt_x_track_
double	dzdt_y_track_
double	r2_xt_
double	r2_yt_
double	r2_zt_x_track_
double	r2_zt_y_track_
double	max_time_gap_xz_
double	max_time_gap_yz_
double	beta_
Track2D	x_track_
Track2D	y_track_
rb::Cluster	cluster_

Detailed Description

Definition at line 21 of file Track3D.h.

Constructor & Destructor Documentation

mono::Track3D::Track3D

(

)

Definition at line 18 of file Track3D.cxx.

{
}

mono::Track3D::Track3D	(	Track2D const &	xz_track,
		Track2D const &	yz_track
	)

Definition at line 24 of file Track3D.cxx.

References abs(), rb::Cluster::Add(), rb::Cluster::AllCells(), beta_, mono::Track2D::cluster, cluster_, dxdt_, dxdz_, dydt_, dydz_, dzdt_, dzdt_x_track_, dzdt_y_track_, mono::Track2D::intercept, max_time_gap_xz_, max_time_gap_yz_, util::pythag(), r2_xt_, r2_yt_, r2_zt_x_track_, r2_zt_y_track_, mono::Track2D::slope, mono::Constants::SPEED_OF_LIGHT, time_fit(), time_gap(), velocity(), x0_, rb::Cluster::XCells(), y0_, rb::Cluster::YCells(), and z0_.

                                                                   :
  default_value_  (-9999),
  dxdt_           (default_value_),
  dydt_           (default_value_),
  dzdt_           (default_value_),
  dzdt_x_track_   (default_value_),
  dzdt_y_track_   (default_value_),
  r2_xt_          (default_value_),
  r2_yt_          (default_value_),
  r2_zt_x_track_  (default_value_),
  r2_zt_y_track_  (default_value_),
  max_time_gap_xz_(default_value_),
  max_time_gap_yz_(default_value_),
  beta_           (default_value_),
  x_track_        (x_track),
  y_track_        (y_track)
{
  x0_ = x_track.intercept;
  y0_ = y_track.intercept;
  z0_ = 0.0;

  dxdz_ = x_track.slope;
  dydz_ = y_track.slope;

  bool xt_fit_valid =
    time_fit(x_track, dxdt_, dzdt_x_track_, r2_xt_, r2_zt_x_track_);
  bool yt_fit_valid =
    time_fit(y_track, dydt_, dzdt_y_track_, r2_yt_, r2_zt_y_track_);

  if (xt_fit_valid && yt_fit_valid)
  {
    dzdt_ = dzdt_x_track_;
    double velocity = util::pythag(dxdt_, dydt_, dzdt_);
    beta_ = velocity / Constants::SPEED_OF_LIGHT;
  }


  /*
    Hit Matching

    Note that this can be done here or at the beginning of the constructor,
    but it does not affect the x_track and y_track themselves, which are used
    for the fit.
  */
  for (auto const& x_hit : x_track.cluster.AllCells())
  {
    bool has_match = false;
    for (auto const& y_hit : y_track.cluster.AllCells())
      if (abs(int(x_hit->Plane()) - int(y_hit->Plane())) <= 1)
        has_match = true;
    
    if (has_match)
      cluster_.Add(x_hit);
  }

  for (auto const& y_hit : y_track.cluster.AllCells())
  {
    bool has_match = false;
    for (auto const& x_hit : x_track.cluster.AllCells())
      if (abs(int(x_hit->Plane()) - int(y_hit->Plane())) <= 1)
        has_match = true;

    if (has_match)
      cluster_.Add(y_hit);
  }

  max_time_gap_xz_ = time_gap(cluster_.XCells());
  max_time_gap_yz_ = time_gap(cluster_.YCells());
}

Member Function Documentation

double mono::Track3D::beta

(

)

const

Definition at line 178 of file Track3D.cxx.

References beta_.

Referenced by mono::MonopoleTrack::filter().

{
  return beta_;
}

rb::Cluster mono::Track3D::cluster

(

)

const

Definition at line 185 of file Track3D.cxx.

References cluster_, and time_fit().

Referenced by mono::MonopoleTrack::filter(), mono::Monopole::n_mc_hits(), and mono::operator<<().

{
  return cluster_;
}

double mono::Track3D::dxdz

(

)

const

Definition at line 110 of file Track3D.cxx.

References dxdz_.

Referenced by mono::operator<<().

{
  return dxdz_;
}

double mono::Track3D::dydz

(

)

const

Definition at line 117 of file Track3D.cxx.

References dydz_.

Referenced by mono::operator<<().

{
  return dydz_;
}

double mono::Track3D::dzdt

(

geo::View_t

view

)

const

Definition at line 124 of file Track3D.cxx.

References default_value_, dzdt_x_track_, dzdt_y_track_, kX, and kY.

{
  if (view == geo::View_t::kX)
    return dzdt_x_track_;

  if (view == geo::View_t::kY)
    return dzdt_y_track_;

  return default_value_;
}

double mono::Track3D::get_covariance ( double const &  a, double const &  mean_a, double const &  b, double const &  mean_b  ) const

private

Definition at line 299 of file Track3D.cxx.

Referenced by get_variance(), and time_fit().

{
  return (a - mean_a) * (b - mean_b);
}

double mono::Track3D::get_variance ( double const &  a, double const &  mean_a  ) const

private

Definition at line 291 of file Track3D.cxx.

References get_covariance().

Referenced by time_fit(), and time_gap().

{
  return (a - mean_a) * (a - mean_a);
}

double mono::Track3D::max_time_gap_xz

(

)

const

Definition at line 164 of file Track3D.cxx.

References max_time_gap_xz_.

{
  return max_time_gap_xz_;
}

double mono::Track3D::max_time_gap_yz

(

)

const

Definition at line 171 of file Track3D.cxx.

References max_time_gap_yz_.

{
  return max_time_gap_yz_;
}

double mono::Track3D::r2_xt

(

)

const

Definition at line 137 of file Track3D.cxx.

References r2_xt_.

{
  return r2_xt_;
}

double mono::Track3D::r2_yt

(

)

const

Definition at line 144 of file Track3D.cxx.

References r2_yt_.

{
  return r2_yt_;
}

double mono::Track3D::r2_zt

(

geo::View_t

view

)

const

Definition at line 151 of file Track3D.cxx.

References default_value_, kX, kY, r2_zt_x_track_, and r2_zt_y_track_.

{
  if (view == geo::View_t::kX)
    return r2_zt_x_track_;

  if (view == geo::View_t::kY)
    return r2_zt_y_track_;

  return default_value_;
}

TVector3 mono::Track3D::start

(

)

const

Definition at line 96 of file Track3D.cxx.

References x0_, y0_, and z0_.

{
  return TVector3(x0_, y0_, z0_);
}

bool mono::Track3D::time_fit ( Track2D const &  track, double &  dvdt, double &  dzdt, double &  r2_vt, double &  r2_zt  ) const

private

Definition at line 193 of file Track3D.cxx.

References rb::Cluster::AllCells(), geo::GeometryBase::CellInfo(), geo::GeometryBase::CellTpos(), mono::Track2D::cluster, default_value_, get_covariance(), get_variance(), extractScale::mean, rb::Cluster::NCell(), cet::pow(), confusionMatrixTree::t, registry_explorer::v, variance(), and test::z.

Referenced by cluster(), and Track3D().

{
  dvdt = default_value_;
  dzdt = default_value_;
  double N = track.cluster.NCell();
  // we need at least two hits to get a non-zero variance
  if (N <= 2)
    return false;
  
  art::ServiceHandle<geo::Geometry> geometry;

  std::map<std::string, double> mean = { {"v", 0.0}, {"z", 0.0}, {"t", 0.0} };
  for (auto const& hit : track.cluster.AllCells())
  {
    double xyz[3];
    geometry->CellInfo(hit->Plane(), hit->Cell(), NULL, xyz);
    TVector3 hit_vector(xyz);
    
    mean.at("v") += geometry->CellTpos(hit->Plane(), hit->Cell()) / N;
    mean.at("z") += hit_vector.z() / N;
    mean.at("t") += hit->TNS() / N;
  }

  std::map<std::string, double> variance =
    { {"t", 0.0}, {"v", 0.0}, {"z", 0.0}, {"vt", 0.0}, {"zt", 0.0} };
  for (auto const& hit : track.cluster.AllCells())
  {
    double xyz[3];
    geometry->CellInfo(hit->Plane(), hit->Cell(), NULL, xyz);
    TVector3 hit_vector(xyz);
    double v = geometry->CellTpos(hit->Plane(), hit->Cell());
    double z = hit_vector.z();
    double t = hit->TNS();

    variance.at("t")  += get_variance(t, mean.at("t")) / (N - 1);
    variance.at("v")  += get_variance(v, mean.at("v")) / (N - 1);
    variance.at("z")  += get_variance(z, mean.at("z")) / (N - 1);
    variance.at("vt") +=
      get_covariance(v, mean.at("v"), t, mean.at("t")) / (N - 1);
    variance.at("zt") +=
      get_covariance(z, mean.at("z"), t, mean.at("t")) / (N - 1);
  }
  
  if (variance.at("t") == 0  || variance.at("v") == 0 || variance.at("z") == 0)
    return false;
  // throw art::Exception(art::errors::InvalidNumber)
  //   << "Monopole/Track3D.cxx\n"
  //   << "Variance is zero.  Aborting.\n" << std::endl;


  //
  // Set all of the results:
  //
  dvdt = variance.at("vt") / variance.at("t");
  dzdt = variance.at("zt") / variance.at("t");

  // The intercept is currently not used.
  // double v0 = mean.at("v") - dvdt * mean.at("t");
  // double z0 = mean.at("z") - dzdt * mean.at("t");

  // Calculate the squared correlation coefficienct (r^2):
  r2_vt =
    std::pow(variance.at("vt"), 2) / (variance.at("v") * variance.at("t"));
  r2_zt =
    std::pow(variance.at("zt"), 2) / (variance.at("z") * variance.at("t"));

  return true;
}

double mono::Track3D::time_gap ( art::PtrVector< rb::CellHit > const &  hits ) const

private

Definition at line 266 of file Track3D.cxx.

References default_value_, get_variance(), std::max(), fillBadChanDBTables::result, and plot_timing_data::times.

Referenced by Track3D().

{
  double result = default_value_;

  std::set<double> times;
  for (auto const& hit : hits)
    times.insert(hit->TNS());

  std::adjacent_find
  (times.begin(), times.end(), [&](double object, double next_object)
  {
    double gap = next_object - object;
    result = std::max(gap, result);

    // This return is meaningless, but a boolean return is required by
    // adjacent find.  We choose false, so the loop goes until the end.
    return false;
  });

  return result;
}

TVector3 mono::Track3D::velocity

(

)

const

Definition at line 103 of file Track3D.cxx.

References dxdt_, dydt_, and dzdt_.

Referenced by Track3D().

{
  return TVector3(dxdt_, dydt_, dzdt_);
}

Member Data Documentation

double mono::Track3D::beta_

private

Definition at line 68 of file Track3D.h.

Referenced by beta(), and Track3D().

rb::Cluster mono::Track3D::cluster_

private

Definition at line 71 of file Track3D.h.

Referenced by cluster(), and Track3D().

double mono::Track3D::default_value_

private

Definition at line 59 of file Track3D.h.

Referenced by dzdt(), r2_zt(), time_fit(), and time_gap().

double mono::Track3D::dxdt_

private

Definition at line 63 of file Track3D.h.

Referenced by Track3D(), and velocity().

double mono::Track3D::dxdz_

private

Definition at line 62 of file Track3D.h.

Referenced by dxdz(), and Track3D().

double mono::Track3D::dydt_

private

Definition at line 63 of file Track3D.h.

Referenced by Track3D(), and velocity().

double mono::Track3D::dydz_

private

Definition at line 62 of file Track3D.h.

Referenced by dydz(), and Track3D().

double mono::Track3D::dzdt_

private

Definition at line 63 of file Track3D.h.

Referenced by Track3D(), and velocity().

double mono::Track3D::dzdt_x_track_

private

Definition at line 64 of file Track3D.h.

Referenced by dzdt(), and Track3D().

double mono::Track3D::dzdt_y_track_

private

Definition at line 64 of file Track3D.h.

Referenced by dzdt(), and Track3D().

double mono::Track3D::max_time_gap_xz_

private

Definition at line 67 of file Track3D.h.

Referenced by max_time_gap_xz(), and Track3D().

double mono::Track3D::max_time_gap_yz_

private

Definition at line 67 of file Track3D.h.

Referenced by max_time_gap_yz(), and Track3D().

double mono::Track3D::r2_xt_

private

Definition at line 65 of file Track3D.h.

Referenced by r2_xt(), and Track3D().

double mono::Track3D::r2_yt_

private

Definition at line 65 of file Track3D.h.

Referenced by r2_yt(), and Track3D().

double mono::Track3D::r2_zt_x_track_

private

Definition at line 66 of file Track3D.h.

Referenced by r2_zt(), and Track3D().

double mono::Track3D::r2_zt_y_track_

private

Definition at line 66 of file Track3D.h.

Referenced by r2_zt(), and Track3D().

double mono::Track3D::x0_

private

Definition at line 61 of file Track3D.h.

Referenced by start(), and Track3D().

Track2D mono::Track3D::x_track_

private

Definition at line 70 of file Track3D.h.

double mono::Track3D::y0_

private

Definition at line 61 of file Track3D.h.

Referenced by start(), and Track3D().

Track2D mono::Track3D::y_track_

private

Definition at line 70 of file Track3D.h.

double mono::Track3D::z0_

private

Definition at line 61 of file Track3D.h.

Referenced by start(), and Track3D().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-26/Monopole/Track3D.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-26/Monopole/Track3D.cxx