Here, we set/reset the Momentum Error
Resolution could be run-dependent!
Bool_t simuflag;
simuflag=gHades->getCurrentEvent()->getCategory(catGeantKine);
if (!pParams->getCut("MomentumError_SEC_0", mom_err_sec_0))
{
if(simuflag) mom_err_sec_0=0.02*1.26; else mom_err_sec_0=0.06;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_0 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_0" << std::endl;
}
if (!pParams->getCut("MomentumError_SEC_1", mom_err_sec_1))
{
if(simuflag) mom_err_sec_1=0.02*1.26; else mom_err_sec_1=0.046;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_1 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_1" << std::endl;
}
if (!pParams->getCut("MomentumError_SEC_2", mom_err_sec_2))
{
if(simuflag) mom_err_sec_2=0.025*1.26; else mom_err_sec_2=0.17;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_2 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_2" << std::endl;
}
if (!pParams->getCut("MomentumError_SEC_3", mom_err_sec_3))
{
if(simuflag) mom_err_sec_3=0.02*1.26; else mom_err_sec_3=0.044;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_3 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_3" << std::endl;
}
if (!pParams->getCut("MomentumError_SEC_4", mom_err_sec_4))
{
if(simuflag) mom_err_sec_4=0.02*1.26; else mom_err_sec_4=0.04;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_4 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_4" << std::endl;
}
if (!pParams->getCut("MomentumError_SEC_5", mom_err_sec_5))
{
if(simuflag) mom_err_sec_5=0.025*1.26; else mom_err_sec_5=0.071;
std::cout << "HHypPPPipPimPi0Kine::init: MomentumError_SEC_5 not found" << std::endl;
std::cout << "using hardcoded MomentumError_SEC_5" << std::endl;
}
std::cout << "MomentumError_SEC_0 is: " << mom_err_sec_0 << std::endl;
std::cout << "MomentumError_SEC_1 is: " << mom_err_sec_1 << std::endl;
std::cout << "MomentumError_SEC_2 is: " << mom_err_sec_2 << std::endl;
std::cout << "MomentumError_SEC_3 is: " << mom_err_sec_3 << std::endl;
std::cout << "MomentumError_SEC_4 is: " << mom_err_sec_4 << std::endl;
std::cout << "MomentumError_SEC_5 is: " << mom_err_sec_5 << std::endl;