double HiggsKvKfgen::computeGammaTotalRatio() const
{
double gtt_SM = trueSM.computeGammaHgg_tt();
double gbb_SM = trueSM.computeGammaHgg_bb();
double gtb_SM = trueSM.computeGammaHgg_tb();
double Ggg_over_GggSM = (computeKt() * computeKt() * gtt_SM
+ computeKb() * computeKb() * gbb_SM
+ computeKt() * computeKb() * gtb_SM) / (gtt_SM + gbb_SM + gtb_SM);
return ((Ggg_over_GggSM * trueSM.computeBrHtogg()
+ computeKW() * computeKW() * trueSM.computeBrHtoWW()
+ computeKZ() * computeKZ() * trueSM.computeBrHtoZZ()
+ computeKZga() * computeKZga() * trueSM.computeBrHtoZga()
+ computeKgaga() * computeKgaga() * trueSM.computeBrHtogaga()
+ computeKmu() * computeKmu() * trueSM.computeBrHtomumu()
+ computeKtau() * computeKtau() * trueSM.computeBrHtotautau()
+ computeKc() * computeKc() * trueSM.computeBrHtocc()
+ computeKb() * computeKb() * trueSM.computeBrHtobb())
/ (trueSM.computeBrHtogg()
+ trueSM.computeBrHtoWW()
+ trueSM.computeBrHtoZZ()
+ trueSM.computeBrHtoZga()
+ trueSM.computeBrHtogaga()
+ trueSM.computeBrHtomumu()
+ trueSM.computeBrHtotautau()
+ trueSM.computeBrHtocc()
+ trueSM.computeBrHtobb() - BrHinv));
}
double HiggsKvKfgen::BrHggRatio() const
{
double gtt_SM = trueSM.computeGammaHgg_tt();
double gbb_SM = trueSM.computeGammaHgg_bb();
double gtb_SM = trueSM.computeGammaHgg_tb();
double Ggg_over_GggSM = (computeKt() * computeKt() * gtt_SM
+ computeKb() * computeKb() * gbb_SM
+ computeKt() * computeKb() * gtb_SM) / (gtt_SM + gbb_SM + gtb_SM);
return (Ggg_over_GggSM / computeGammaTotalRatio());
}
double HiggsKvKfgen::computeKg(const double sqrt_s) const
{
if (sqrt_s == 1.96)
return computeKt();
double sigmatt_SM = trueSM.computeSigmaggH_tt(sqrt_s);
double sigmabb_SM = trueSM.computeSigmaggH_bb(sqrt_s);
double sigmatb_SM = trueSM.computeSigmaggH_tb(sqrt_s);
return (sqrt((computeKt() * computeKt() * sigmatt_SM
+ computeKb() * computeKb() * sigmabb_SM
+ computeKt() * computeKb() * sigmatb_SM)
/ (sigmatt_SM + sigmabb_SM + sigmatb_SM)));
}
double HiggsKvKf::computeGammaTotalRatio() const
{
return ((computeKg() * computeKg() * trueSM.computeBrHtogg()
+ computeKW() * computeKW() * trueSM.computeBrHtoWW()
+ computeKZ() * computeKZ() * trueSM.computeBrHtoZZ()
+ computeKZga() * computeKZga() * trueSM.computeBrHtoZga()
+ computeKgaga() * computeKgaga() * trueSM.computeBrHtogaga()
+ computeKmu() * computeKmu() * trueSM.computeBrHtomumu()
+ computeKtau() * computeKtau() * trueSM.computeBrHtotautau()
+ computeKc() * computeKc() * trueSM.computeBrHtocc()
+ computeKb() * computeKb() * trueSM.computeBrHtobb())
/ (trueSM.computeBrHtogg()
+ trueSM.computeBrHtoWW()
+ trueSM.computeBrHtoZZ()
+ trueSM.computeBrHtoZga()
+ trueSM.computeBrHtogaga()
+ trueSM.computeBrHtomumu()
+ trueSM.computeBrHtotautau()
+ trueSM.computeBrHtocc()
+ trueSM.computeBrHtobb() - BrHinv));
}
double HiggsKvKfgen::BrHbbRatio() const
{
return (computeKb() * computeKb() / computeGammaTotalRatio());
}