Commit ada1386b authored by Frank Siegert's avatar Frank Siegert
Browse files

Limit MENLOPS BBar/B weight to 10.0 by default to avoid potential spikes in unweighted events.

parent 3f67af38
......@@ -34,6 +34,7 @@ Perturbative_Interface::Perturbative_Interface
m_cmode=ToType<int>(rpa->gen.Variable("METS_CLUSTER_MODE"));
m_bbarmode=read.GetValue<int>("METS_BBAR_MODE",1);
m_globalkfac=read.GetValue<double>("GLOBAL_KFAC",0.);
m_maxkfac=read.GetValue<double>("MENLOPS_MAX_KFAC",10.0);
}
Perturbative_Interface::Perturbative_Interface
......@@ -211,7 +212,7 @@ bool Perturbative_Interface::LocalKFactor(ATOOLS::Cluster_Amplitude* ampl)
MCatNLO_Process* mcnloproc=dynamic_cast<MCatNLO_Process*>(procs[i]);
if (mcnloproc) {
double K(mcnloproc->LocalKFactor(*ampl));
if (K==0.0) continue;
if (K==0.0 || K>m_maxkfac) continue;
m_weight*=K;
return true;
}
......
......@@ -31,7 +31,7 @@ namespace SHERPA {
ATOOLS::Cluster_Amplitude *p_ampl;
double m_weight, m_globalkfac;
double m_weight, m_globalkfac, m_maxkfac;
int m_cmode, m_bbarmode;
bool SetColours(ATOOLS::Cluster_Amplitude *ampl,ATOOLS::Blob *const bl);
......
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