Changes in external/fastjet/ClusterSequence.cc [1d208a2:10e33bc] in git

File:

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

external/fastjet/ClusterSequence.cc

@@ -1 +1 @@
 //FJSTARTHEADER
-// $Id: ClusterSequence.cc 4154 2016-07-20 16:20:48Z soyez $
+// $Id: ClusterSequence.cc 3809 2015-02-20 13:05:13Z soyez $
 //
 // Copyright (c) 2005-2014, Matteo Cacciari, Gavin P. Salam and Gregory Soyez
@@ -152 +152 @@
   // set the pointer in the wrapper to this object to NULL to say that
   // we're going out of scope
-  if (_structure_shared_ptr){
-    ClusterSequenceStructure* csi = dynamic_cast<ClusterSequenceStructure*>(_structure_shared_ptr.get()); 
+  if (_structure_shared_ptr()){
+    ClusterSequenceStructure* csi = dynamic_cast<ClusterSequenceStructure*>(_structure_shared_ptr()); 
     // normally the csi is purely internal so it really should not be
     // NULL i.e assert should be OK
@@ -724 +724 @@
   }
   
-  //bool code_should_never_reach_here = false;
-  //assert(code_should_never_reach_here);
-
-  assert(0 && "Code should never reach here");
-
+  bool code_should_never_reach_here = false;
+  assert(code_should_never_reach_here); 
   return N2MHTLazy9;
 
@@ -778 +775 @@
 
 
-ClusterSequence & ClusterSequence::operator=(const ClusterSequence & cs) {
-  // self assignment is trivial
-  if (&cs != this) {
-    _deletes_self_when_unused = false;
-    transfer_from_sequence(cs);
-  }
-  return *this;
-}
-
 //----------------------------------------------------------------------
 // transfer the sequence contained in other_seq into our own;
@@ -827 +815 @@
 
   // clean up existing structure
-  if (_structure_shared_ptr) {
+  if (_structure_shared_ptr()) {
     // If there are jets associated with an old version of the CS and
     // a new one, keeping track of when to delete the CS becomes more
@@ -835 +823 @@
     // anything that is currently associated with the cluster sequence
     // should be told that its cluster sequence no longer exists
-    ClusterSequenceStructure* csi = dynamic_cast<ClusterSequenceStructure*>(_structure_shared_ptr.get()); 
+    ClusterSequenceStructure* csi = dynamic_cast<ClusterSequenceStructure*>(_structure_shared_ptr()); 
     assert(csi != NULL);
     csi->set_associated_cs(NULL);
@@ -1467 +1455 @@
 // initialise the history in a standard way
 void ClusterSequence::_add_step_to_history (
-               //NO_LONGER_USED: const int step_number,
-               const int parent1, 
+               const int step_number, const int parent1, 
                const int parent2, const int jetp_index,
                const double dij) {
@@ -1482 +1469 @@
 
   int local_step = _history.size()-1;
-  //#ifndef __NO_ASSERTS__
-  //assert(local_step == step_number);
-  //#endif
+  assert(local_step == step_number);
 
   // sanity check: make sure the particles have not already been recombined
@@ -1684 +1669 @@
   int hist_j = _jets[jet_j].cluster_hist_index();
 
-  _add_step_to_history(min(hist_i, hist_j), max(hist_i,hist_j),
+  _add_step_to_history(newstep_k, min(hist_i, hist_j), max(hist_i,hist_j),
                        newjet_k, dij);
-
-  //  _add_step_to_history(newstep_k, min(hist_i, hist_j), max(hist_i,hist_j),
-  //                   newjet_k, dij);
-
 
 }
@@ -1699 +1680 @@
 void ClusterSequence::_do_iB_recombination_step(
                                   const int jet_i, const double diB) {
+  // get history index
+  int newstep_k = _history.size();
+
   // recombine the jet with the beam
-  _add_step_to_history(_jets[jet_i].cluster_hist_index(),BeamJet,
+  _add_step_to_history(newstep_k,_jets[jet_i].cluster_hist_index(),BeamJet,
                        Invalid, diB);
-
-  // // get history index
-  // int newstep_k = _history.size();
-  // 
-  // _add_step_to_history(newstep_k,_jets[jet_i].cluster_hist_index(),BeamJet,
-  //                   Invalid, diB);
 
 }