Changeset 10e33bc in git for external/fastjet/ClusterSequence.cc

Timestamp:

Jun 26, 2015, 9:03:11 PM (9 years ago)

Author:

Pavel Demin <pavel.demin@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

1bc87da

Parents:

8497ac6

Message:

update FastJet library to 3.1.2

File:

• external/fastjet/ClusterSequence.cc (modified) (3 diffs)

external/fastjet/ClusterSequence.cc

@@ -1 +1 @@
 //FJSTARTHEADER
-// $Id: ClusterSequence.cc 3685 2014-09-11 20:15:00Z salam $
+// $Id: ClusterSequence.cc 3809 2015-02-20 13:05:13Z soyez $
 //
 // Copyright (c) 2005-2014, Matteo Cacciari, Gavin P. Salam and Gregory Soyez
@@ -372 +372 @@
     _plugin_activated = false;
 
+  } else if (_strategy == N2MHTLazy9AntiKtSeparateGhosts) {
 #ifndef __FJCORE__
-  } else if (_strategy == N2MHTLazy9AntiKtSeparateGhosts) {
     // attempt to use an external tiling routine -- it manipulates
     // the CS history via the plugin mechanism
@@ -1471 +1471 @@
   assert(local_step == step_number);
 
+  // sanity check: make sure the particles have not already been recombined
+  //
+  // Note that good practice would make this an assert (since this is
+  // a serious internal issue). However, we decided to throw an
+  // InternalError so that the end user can decide to catch it and
+  // retry the clustering with a different strategy.
+
   assert(parent1 >= 0);
+  if (_history[parent1].child != Invalid){
+    throw InternalError("trying to recomine an object that has previsously been recombined");
+  }
   _history[parent1].child = local_step;
-  if (parent2 >= 0) {_history[parent2].child = local_step;}
+  if (parent2 >= 0) {
+    if (_history[parent2].child != Invalid){
+      throw InternalError("trying to recomine an object that has previsously been recombined");
+    }
+    _history[parent2].child = local_step;
+  }
 
   // get cross-referencing right from PseudoJets