Commit e7c85977 authored by Enrico Bothmann's avatar Enrico Bothmann

Add missing OrderQCD/EW implementations

- OpenLoops_Born did not overwrite the implementation of its (abstract)
  superclass Tree_ME2_Base
- This led to a wrong renormalisation scale reweighting of (at least)
  loop-induced processes using OpenLoops
- To prevent this in the future, OrderQCD/EW are now declared as pure virtual
  functions in Tree_ME2_Base
parent 0bba9098
......@@ -35,6 +35,16 @@ double OpenLoops_Born::Calc(const Vec4D_Vector& momenta)
return m_symfac*result;
}
int OpenLoops_Born::OrderQCD(const int &id)
{
return OpenLoops_Interface::GetIntParameter("coupling_qcd_0");
}
int OpenLoops_Born::OrderEW(const int &id)
{
return OpenLoops_Interface::GetIntParameter("coupling_ew_0");
}
}
using namespace OpenLoops;
......@@ -64,5 +74,3 @@ operator()(const Process_Info &pi) const
return NULL;
}
......@@ -24,6 +24,9 @@ namespace OpenLoops {
double Calc(const ATOOLS::Vec4D_Vector& momenta);
int OrderQCD(const int &id=-1);
int OrderEW(const int &id=-1);
};
}
......
......@@ -46,16 +46,6 @@ void Tree_ME2_Base::FillCombinations
{
}
int Tree_ME2_Base::OrderQCD(const int &id)
{
return -1;
}
int Tree_ME2_Base::OrderEW(const int &id)
{
return -1;
}
double Tree_ME2_Base::TR() const
{
return 0.5;
......
......@@ -51,8 +51,8 @@ namespace PHASIC {
(std::set<std::pair<size_t,size_t> > &combs,
std::map<size_t,ATOOLS::Flavour_Vector> &fls);
virtual int OrderQCD(const int &id=-1);
virtual int OrderEW(const int &id=-1);
virtual int OrderQCD(const int &id=-1) = 0;
virtual int OrderEW(const int &id=-1) = 0;
virtual double TR() const;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment