PseudoJet FlavourKTPlugin::combine(const PseudoJet& jet1, const PseudoJet& jet2)
{
  PseudoJet newjet = PseudoJet(0., 0., 0., 0.);
  if (jet1.E() != 0. and jet2.E() != 0.) {
    newjet = jet1 + jet2;
  }
  addFlavour(newjet, combineFlavour(getFlavour(jet1), getFlavour(jet2)));
  return newjet;
}
 void init(const PseudoJet& jet) {
   if (jet.E() == 0.) { // beam representative
     pt2 = 0.;
     rap = 0.;
     phi = 0.;
   } else {
     pt2 = jet.pt2();
     rap = jet.rapidity();
     phi = jet.phi();
   }
   lhid = FlavourKTPlugin::getFlavour(jet);
 }