Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • external/fastjet/contribs/Nsubjettiness/Njettiness.cc

    r1d208a2 r973b92a  
    55//  Jesse Thaler, Ken Van Tilburg, Christopher K. Vermilion, and TJ Wilkason
    66//
    7 //  $Id: Njettiness.cc 933 2016-04-04 22:23:32Z jthaler $
     7//  $Id: Njettiness.cc 821 2015-06-15 18:50:53Z jthaler $
    88//----------------------------------------------------------------------
    99// This file is part of FastJet contrib.
     
    4646// setAxes for Manual mode
    4747void Njettiness::setAxes(const std::vector<fastjet::PseudoJet> & myAxes) {
    48    if (_axes_def->needsManualAxes()) {
     48   if (_axes_def()->needsManualAxes()) {
    4949      _currentAxes = myAxes;
    5050   } else {
     
    5959      _currentAxes = inputJets;
    6060      _currentAxes.resize(n_jets,fastjet::PseudoJet(0.0,0.0,0.0,0.0));
    61      
    62      // Put in empty tau components
    63      std::vector<double> dummy_jet_pieces;
    64      _current_tau_components = TauComponents(UNDEFINED_SHAPE,
    65                                              dummy_jet_pieces,
    66                                              0.0,
    67                                              1.0,
    68                                              _currentAxes,
    69                                              _currentAxes
    70                                              );
     61      _current_tau_components = TauComponents();
    7162      _seedAxes = _currentAxes;
    7263      _currentPartition = TauPartition(n_jets); // empty partition
    7364   } else {
    74       assert(_axes_def); // this should never fail.
     65      assert(_axes_def()); // this should never fail.
    7566     
    76       if (_axes_def->needsManualAxes()) { // if manual mode
     67      if (_axes_def()->needsManualAxes()) { // if manual mode
    7768         // take current axes as seeds
    7869         _seedAxes = _currentAxes;
    7970         
    8071         // refine axes if requested
    81          _currentAxes = _axes_def->get_refined_axes(n_jets,inputJets,_seedAxes, _measure_def.get());
     72         _currentAxes = _axes_def->get_refined_axes(n_jets,inputJets,_seedAxes, _measure_def());
    8273      } else { // non-manual axes
    8374         
    8475          //set starting point for minimization
    85          _seedAxes = _axes_def->get_starting_axes(n_jets,inputJets,_measure_def.get());
     76         _seedAxes = _axes_def->get_starting_axes(n_jets,inputJets,_measure_def());
    8677         
    8778         // refine axes as needed
    88          _currentAxes = _axes_def->get_refined_axes(n_jets,inputJets,_seedAxes, _measure_def.get());
     79         _currentAxes = _axes_def->get_refined_axes(n_jets,inputJets,_seedAxes, _measure_def());
    8980         
    9081         // NOTE:  The above two function calls are combined in "AxesDefinition::get_axes"
Note: See TracChangeset for help on using the changeset viewer.