float MelaReweighterWW::weightStoH(){
//setupDaughters(_isVBF, id_l1, id_l2, id_n1, id_n2, l1, l2, n1, n2, associated, idsAssociated);
float meSpow;
_mela->setProcess(TVar::HSMHiggs, TVar::MCFM, _isVBF ? TVar::JJVBF_S : TVar::ZZGG);
// Added here -- U. Sarica
_mela->setMelaHiggsMassWidth(_mpole, _width, 0);
//
if (!_isVBF)
_mela->computeP(meSpow, false);
else
_mela->computeProdDecP(meSpow, false);
float meH;
_mela->setProcess(TVar::HSMHiggs, TVar::MCFM, _isVBF ? TVar::JJVBF : TVar::ZZGG);
//Added here -- U. Sarica
_mela->setMelaHiggsMassWidth(125., 4.07e-3, 0);
//
if (!_isVBF)
_mela->computeP(meH, false);
else
_mela->computeProdDecP(meH, false);
//_mela->resetInputEvent();
//std::cout << "### meS = " << meS << " meSpow = " << meSpow << " meB = " << meB << " meSBI = " << meSBI << std::endl;
//std::cout << "### weight = " << (meSBI-meS-meB)/meSpow << std::endl;
return (meH)/meSpow;
}
float MelaReweighterWW::weightStoSBI(){
//setupDaughters(isVBF, id_l1, id_l2, id_n1, id_n2, l1, l2, n1, n2, associated, idsAssociated);
float meS;
_mela->setProcess(TVar::HSMHiggs, TVar::MCFM, _isVBF ? TVar::JJVBF_S : TVar::ZZGG);
// Added here -- U. Sarica
_mela->setMelaHiggsMassWidth(_mpole, _width, 0);
_mela->setMelaHiggsMassWidth(125., 4.07e-3, 1);
//
if (!_isVBF)
_mela->computeP(meS, false);
else
_mela->computeProdDecP(meS, false);
float meSBI;
_mela->setProcess(TVar::bkgWW_SMHiggs, TVar::MCFM, _isVBF ? TVar::JJEW : TVar::ZZGG);
// Added here -- U. Sarica
_mela->setMelaHiggsMassWidth(_mpole, _width, 0);
_mela->setMelaHiggsMassWidth(125., 4.07e-3, 1);
//
if (!_isVBF)
_mela->computeP(meSBI, false);
else
_mela->computeProdDecP(meSBI, false);
_mela->resetInputEvent();
return meSBI/meS;
}