本文整理汇总了C++中atools::Flavour::HadMass方法的典型用法代码示例。如果您正苦于以下问题:C++ Flavour::HadMass方法的具体用法?C++ Flavour::HadMass怎么用?C++ Flavour::HadMass使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类atools::Flavour
的用法示例。
在下文中一共展示了Flavour::HadMass方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: GetXPDF
double Hadron_Remnant::GetXPDF(ATOOLS::Flavour flavour,double scale)
{
double cut, x;
cut=2.0*(flavour.HadMass()+m_hardpt.PPerp()/
sqr(m_companions.size()))/p_beam->OutMomentum()[0];
// assume heavy flavours have been pair-produced
// => scale should be approximately (2m)^2
scale=Max(scale,4.0*sqr(flavour.Mass()));
if (scale<p_pdfbase->Q2Min()) {
msg_Error()<<"Hadron_Remnant::GetXPDF("<<flavour<<","<<scale<<"): "
<<"Scale under-runs minimum given by PDF: "
<<scale<<" < "<<p_pdfbase->Q2Min()<<std::endl;
scale=1.001*p_pdfbase->Q2Min();
}
unsigned int xtrials, pdftrials=0;
while (true) {
++pdftrials;
xtrials=0;
do {
++xtrials;
x=m_xrem*ran->Get();
if (xtrials>=m_maxtrials) {
x=Min(cut,0.999999*p_pdfbase->RescaleFactor());
break;
}
if (x>p_pdfbase->RescaleFactor()) continue;
} while (x<cut);
if (x>p_pdfbase->XMin() && x<p_pdfbase->XMax()) {
p_pdfbase->Calculate(x,scale);
}
else {
m_xscheme=0; return 0.01;
}
if (pdftrials>=m_maxtrials) { m_xscheme=0; return 0.01; }
if (p_pdfbase->GetXPDF(flavour)/x>ran->Get()) return x;
}
return 0.0;
}