17 TVector3 end1 = track1.
Start() + track1.
Dir();
18 TVector3 end2 = track2.
Start() + track2.
Dir();
25 TVector3 F1, TVector3 F2)
28 double sigma1, sigma2;
42 sigma1 = -(alpha * anyV) / (alpha *
delta2_);
43 sigma2 = (beta * anyV) / (beta *
delta1_);
This is the version 2, works for multi-track event.
void DoCalculation(TVector3 I1, TVector3 I2, TVector3 F1, TVector3 F2)
fvar< T > fabs(const fvar< T > &x)
A rb::Prong with full reconstructed trajectory.
def ratio(spec1, spec2, nbins, pot, binrange=(0, 1))
virtual TVector3 Start() const
DOCAInfo(rb::Track track1, rb::Track track2)
virtual TVector3 Dir() const
Unit vector describing prong direction.
TVector3 VertexReco() const