[Rivet] Rivet

[Andrii Verbytskyi]

Dear Andy,
please have a look at
Rivet/src/Projections/Thrust.cc


Maybe I'm mistaken, but n is constant there.


  // Do the general case thrust calculation
  void _calcT(const vector<Vector3>& momenta, double& t, Vector3& taxis)
{
    // This function implements the iterative algorithm as described in
the
    // Pythia manual. We take eight (four) different starting vectors
    // constructed from the four (three) leading particles to make sure
that
    // we don't find a local maximum.
    vector<Vector3> p = momenta;
    assert(p.size() >= 3);                                    <---------
    unsigned int n = 3;                                       <---------
    if (p.size() == 3) n = 3;                                <---------
    vector<Vector3> tvec;
    vector<double> tval;
    std::sort(p.begin(), p.end(), mod2Cmp);
    for (int i = 0 ; i < intpow(2, n-1); ++i) {


best regards,
Andrii