No OneTemporary
Actions

Size

1 KB

Subscribers

None

View Options

	diff --git a/src/ScaleFunction.cc b/src/ScaleFunction.cc
	index c802bac..dfe80bc 100644
	--- a/src/ScaleFunction.cc
	+++ b/src/ScaleFunction.cc
	@@ -1,37 +1,37 @@
	#include "RHEJ/ScaleFunction.hh"
	#include <numeric>

	namespace RHEJ{
	namespace{
	- double ht(Event const & ev){
	+ double ht(ClusteredEvent const & ev){
	double result = 0.;
	- for(auto const & parton: ev.outgoing) result += parton.p.perp();
	+ for(auto const & parton: ev.outgoing()) result += parton.p.perp();
	return result;
	}
	}

	EventParameters Ht::operator()(ClusteredEvent const & ev) const{
	const double mu = ht(ev);
	return {mu, mu, ev.central().weight};
	}

	EventParameters MaxJetPperp::operator()(ClusteredEvent const & ev) const{
	const double mu = sorted_by_pt(ev.jets()).front().pt();
	return {mu, mu, ev.central().weight};
	}

	EventParameters JetInvariantMass::operator()(ClusteredEvent const & ev) const{
	const double mu = std::accumulate(
	begin(ev.jets()), end(ev.jets()), fastjet::PseudoJet{}
	).m();
	return {mu, mu, ev.central().weight};
	}

	EventParameters Product::operator()(ClusteredEvent const &ev) const{
	EventParameters s = (*base_scale_)(ev);
	s.muf *= factor_;
	s.mur *= factor_;
	return s;
	}

	}