Changes in external/fastjet/contribs/Nsubjettiness/Njettiness.cc [1d208a2:973b92a] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
external/fastjet/contribs/Nsubjettiness/Njettiness.cc
r1d208a2 r973b92a 5 5 // Jesse Thaler, Ken Van Tilburg, Christopher K. Vermilion, and TJ Wilkason 6 6 // 7 // $Id: Njettiness.cc 933 2016-04-04 22:23:32Z jthaler $7 // $Id: Njettiness.cc 821 2015-06-15 18:50:53Z jthaler $ 8 8 //---------------------------------------------------------------------- 9 9 // This file is part of FastJet contrib. … … 46 46 // setAxes for Manual mode 47 47 void Njettiness::setAxes(const std::vector<fastjet::PseudoJet> & myAxes) { 48 if (_axes_def ->needsManualAxes()) {48 if (_axes_def()->needsManualAxes()) { 49 49 _currentAxes = myAxes; 50 50 } else { … … 59 59 _currentAxes = inputJets; 60 60 _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(); 71 62 _seedAxes = _currentAxes; 72 63 _currentPartition = TauPartition(n_jets); // empty partition 73 64 } else { 74 assert(_axes_def ); // this should never fail.65 assert(_axes_def()); // this should never fail. 75 66 76 if (_axes_def ->needsManualAxes()) { // if manual mode67 if (_axes_def()->needsManualAxes()) { // if manual mode 77 68 // take current axes as seeds 78 69 _seedAxes = _currentAxes; 79 70 80 71 // 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()); 82 73 } else { // non-manual axes 83 74 84 75 //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()); 86 77 87 78 // 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()); 89 80 90 81 // NOTE: The above two function calls are combined in "AxesDefinition::get_axes"
Note:
See TracChangeset
for help on using the changeset viewer.