void QGraph::keyPressEvent ( QKeyEvent * e ){
switch(e->key()){
case Qt::Key_Right:
viewport.setLeft(viewport.left() +1); viewport.setRight(viewport.right() +1);
break;
case Qt::Key_Left:
viewport.setLeft(viewport.left() -1); viewport.setRight(viewport.right() -1);
break;
case Qt::Key_Down:
viewport.setTop(viewport.top() -1); viewport.setBottom(viewport.bottom() -1);
break;
case Qt::Key_Up:
viewport.setTop(viewport.top() +1); viewport.setBottom(viewport.bottom() +1);
break;
case Qt::Key_Minus:
resolucio=(resolucio*viewport.width())/(viewport.width()+2);
viewport.setCoords(viewport.left() - 1, viewport.top() +1, viewport.right() + 1, viewport.bottom() -1);
update_scale();
break;
case Qt::Key_Plus:
if(viewport.height() < -3 && viewport.width() > 3){
resolucio=(resolucio*viewport.width())/(viewport.width()-2);
viewport.setCoords(viewport.left() + 1, viewport.top() -1, viewport.right() -1, viewport.bottom() +1);
update_scale();
} else return;
break;
default:
return;
}
valid=false;
update_points();
this->repaint(false);
}
void QGraph::mouseReleaseEvent(QMouseEvent *e){
pushed &= ~e->button();
this->setCursor(QCursor(Qt::CrossCursor));
if(e->button() == Qt::LeftButton){
if((toViewport(press) - toViewport(e->pos())).isNull())
return;
QPoint p=toViewport(e->pos())+viewport.topLeft();
QPoint p1=toViewport(press)+viewport.topLeft();
viewport = QRect(p1, p);
if(viewport.bottom() > viewport.top()){
int a=viewport.bottom();
viewport.setBottom(viewport.top());
viewport.setTop(a);
}
if(viewport.left() > viewport.right()){
int a=viewport.left();
viewport.setLeft(viewport.right());
viewport.setRight(a);
}
update_scale();
update_points();
sendStatus(QString("(%1, %2)-(%3, %4)").arg(viewport.left()).arg(viewport.top()).arg(viewport.right()).arg(viewport.bottom()));
}
this->repaint(false);
}
void SNmssm_susy_scale_constraint::apply()
{
assert(snmssm && "Error: pointer to SNmssm<Two_scale> cannot be zero");
snmssm->calcDrBarPars();
update_scale();
const double mtrun = snmssm->displayDrBarPars().mt;
snmssm->rewsb(pp.signMu, mtrun);
}
void QGraph::wheelEvent(QWheelEvent *e){
int d = e->delta()>0 ? -1 : 1;
if(viewport.left()-d < 1 && viewport.top()+d > 1 && viewport.right()+d > 1 && viewport.bottom()-d < 1) {
viewport.setLeft(viewport.left() - d);
viewport.setTop(viewport.top() + d);
viewport.setRight(viewport.right() + d);
viewport.setBottom(viewport.bottom() - d);
update_scale();
update_points();
}
sendStatus(QString("(%1, %2)-(%3, %4)").arg(viewport.left()).arg(viewport.top()).arg(viewport.right()).arg(viewport.bottom()));
}
void CMSSM_susy_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: CMSSM_susy_scale_constraint::apply():"
" model pointer must not be zero");
model->calculate_DRbar_masses();
update_scale();
// apply user-defined susy scale constraints
// the parameters, which are fixed by the EWSB eqs., will now be
// defined at this scale (at the EWSB loop level defined in the
// model)
model->solve_ewsb();
}
static void process_tuple(Tuple* t) {
// Store incoming information
static char conditions_buffer[16];
switch(t->key) {
case KEY_SCALE:
if (strcmp(t->value->cstring, "0") == 0)
scale = 0;
else if (strcmp(t->value->cstring, "1") == 0)
scale = 1;
else if (strcmp(t->value->cstring, "2") == 0)
scale = 2;
else if (strcmp(t->value->cstring, "3") == 0)
scale = 3;
write_settings_to_memory();
handle_weather();
break;
case KEY_CLOCK:
if (strcmp(t->value->cstring, "0") == 0)
clock_format = 0;
else if (strcmp(t->value->cstring, "1") == 0)
clock_format = 1;
write_settings_to_memory();
update_time();
break;
case KEY_DATE:
if (strcmp(t->value->cstring, "0") == 0)
date_format = 0;
else if (strcmp(t->value->cstring, "1") == 0)
date_format = 1;
else if (strcmp(t->value->cstring, "2") == 0)
date_format = 2;
write_settings_to_memory();
handle_date();
break;
case KEY_TEMPERATURE:
update_scale((int)t->value->int32);
break;
case KEY_CONDITIONS:
snprintf(conditions_buffer, sizeof(conditions_buffer), "%s", t->value->cstring);
text_layer_set_text(condition_layer, conditions_buffer);
break;
default:
APP_LOG(APP_LOG_LEVEL_ERROR, "Key %d not recognized!", (int)t->key);
break;
}
}
void genericE6SSM_susy_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: genericE6SSM_susy_scale_constraint:"
" model pointer must not be zero");
model->calculate_DRbar_parameters();
update_scale();
// apply user-defined susy scale constraints
const auto BmuPrimeInput = INPUTPARAMETER(BmuPrimeInput);
const auto muPrimeInput = INPUTPARAMETER(muPrimeInput);
const auto vSInput = INPUTPARAMETER(vSInput);
MODEL->set_BMuPr(BmuPrimeInput);
MODEL->set_MuPr(muPrimeInput);
MODEL->set_vs(vSInput);
// the parameters, which are fixed by the EWSB eqs., will now be
// defined at this scale (at the EWSB loop level defined in the
// model)
model->solve_ewsb();
}
void QGraph::resizeEvent(QResizeEvent *){
buffer.resize(this->size());
update_scale();
pintafunc(NULL);
}
void lowNE6SSM_susy_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: lowNE6SSM_susy_scale_constraint:"
" model pointer must not be zero");
model->calculate_DRbar_parameters();
update_scale();
// apply user-defined susy scale constraints
const auto BmuPrimeInput = INPUTPARAMETER(BmuPrimeInput);
const auto muPrimeInput = INPUTPARAMETER(muPrimeInput);
const auto ssumInput = INPUTPARAMETER(ssumInput);
const auto TanTheta = INPUTPARAMETER(TanTheta);
const auto hEInput = INPUTPARAMETER(hEInput);
const auto SigmaLInput = INPUTPARAMETER(SigmaLInput);
const auto KappaPrInput = INPUTPARAMETER(KappaPrInput);
const auto SigmaxInput = INPUTPARAMETER(SigmaxInput);
const auto gDInput = INPUTPARAMETER(gDInput);
const auto KappaInput = INPUTPARAMETER(KappaInput);
const auto Lambda12Input = INPUTPARAMETER(Lambda12Input);
const auto LambdaxInput = INPUTPARAMETER(LambdaxInput);
const auto fuInput = INPUTPARAMETER(fuInput);
const auto fdInput = INPUTPARAMETER(fdInput);
const auto MuPhiInput = INPUTPARAMETER(MuPhiInput);
const auto XiFInput = INPUTPARAMETER(XiFInput);
const auto g1pInput = INPUTPARAMETER(g1pInput);
const auto vphiInput = INPUTPARAMETER(vphiInput);
const auto TYdInput = INPUTPARAMETER(TYdInput);
const auto ThEInput = INPUTPARAMETER(ThEInput);
const auto TYeInput = INPUTPARAMETER(TYeInput);
const auto TSigmaLInput = INPUTPARAMETER(TSigmaLInput);
const auto TKappaPrInput = INPUTPARAMETER(TKappaPrInput);
const auto TSigmaxInput = INPUTPARAMETER(TSigmaxInput);
const auto TgDInput = INPUTPARAMETER(TgDInput);
const auto TKappaInput = INPUTPARAMETER(TKappaInput);
const auto TLambda12Input = INPUTPARAMETER(TLambda12Input);
const auto TLambdaxInput = INPUTPARAMETER(TLambdaxInput);
const auto TfuInput = INPUTPARAMETER(TfuInput);
const auto TfdInput = INPUTPARAMETER(TfdInput);
const auto TYuInput = INPUTPARAMETER(TYuInput);
const auto BMuPhiInput = INPUTPARAMETER(BMuPhiInput);
const auto LXiFInput = INPUTPARAMETER(LXiFInput);
const auto mq2Input = INPUTPARAMETER(mq2Input);
const auto ml2Input = INPUTPARAMETER(ml2Input);
const auto md2Input = INPUTPARAMETER(md2Input);
const auto mu2Input = INPUTPARAMETER(mu2Input);
const auto me2Input = INPUTPARAMETER(me2Input);
const auto mH1I2Input = INPUTPARAMETER(mH1I2Input);
const auto mH2I2Input = INPUTPARAMETER(mH2I2Input);
const auto mSI2Input = INPUTPARAMETER(mSI2Input);
const auto mDx2Input = INPUTPARAMETER(mDx2Input);
const auto mDxbar2Input = INPUTPARAMETER(mDxbar2Input);
const auto mHp2Input = INPUTPARAMETER(mHp2Input);
const auto mHpbar2Input = INPUTPARAMETER(mHpbar2Input);
const auto MassBInput = INPUTPARAMETER(MassBInput);
const auto MassWBInput = INPUTPARAMETER(MassWBInput);
const auto MassGInput = INPUTPARAMETER(MassGInput);
const auto MassBpInput = INPUTPARAMETER(MassBpInput);
MODEL->set_BMuPr(BmuPrimeInput);
MODEL->set_MuPr(muPrimeInput);
MODEL->set_vs(ssumInput/Sqrt(1 + Sqr(TanTheta)));
MODEL->set_vsb((ssumInput*TanTheta)/Sqrt(1 + Sqr(TanTheta)));
MODEL->set_hE(hEInput);
MODEL->set_SigmaL(SigmaLInput);
MODEL->set_KappaPr(KappaPrInput);
MODEL->set_Sigmax(SigmaxInput);
MODEL->set_gD(gDInput);
MODEL->set_Kappa(KappaInput);
MODEL->set_Lambda12(Lambda12Input);
MODEL->set_Lambdax(LambdaxInput);
MODEL->set_fu(fuInput);
MODEL->set_fd(fdInput);
MODEL->set_MuPhi(MuPhiInput);
MODEL->set_XiF(XiFInput);
MODEL->set_g1p(g1pInput);
MODEL->set_vphi(vphiInput);
MODEL->set_TYd(TYdInput);
MODEL->set_ThE(ThEInput);
MODEL->set_TYe(TYeInput);
MODEL->set_TSigmaL(TSigmaLInput);
MODEL->set_TKappaPr(TKappaPrInput);
MODEL->set_TSigmax(TSigmaxInput);
MODEL->set_TgD(TgDInput);
MODEL->set_TKappa(TKappaInput);
MODEL->set_TLambda12(TLambda12Input);
MODEL->set_TLambdax(TLambdaxInput);
MODEL->set_Tfu(TfuInput);
MODEL->set_Tfd(TfdInput);
MODEL->set_TYu(TYuInput);
MODEL->set_BMuPhi(BMuPhiInput);
MODEL->set_LXiF(LXiFInput);
MODEL->set_mq2(mq2Input);
MODEL->set_ml2(ml2Input);
MODEL->set_md2(md2Input);
MODEL->set_mu2(mu2Input);
MODEL->set_me2(me2Input);
MODEL->set_mH1I2(mH1I2Input);
MODEL->set_mH2I2(mH2I2Input);
MODEL->set_mSI2(mSI2Input);
MODEL->set_mDx2(mDx2Input);
MODEL->set_mDxbar2(mDxbar2Input);
MODEL->set_mHp2(mHp2Input);
MODEL->set_mHpbar2(mHpbar2Input);
MODEL->set_MassB(MassBInput);
MODEL->set_MassWB(MassWBInput);
MODEL->set_MassG(MassGInput);
MODEL->set_MassBp(MassBpInput);
// the parameters, which are fixed by the EWSB eqs., will now be
// defined at this scale (at the EWSB loop level defined in the
// model)
model->solve_ewsb();
}
void CNE6SSMSusy_high_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: CNE6SSMSusy_high_scale_constraint::apply():"
" model pointer must not be zero");
if (std::fabs(model->get_g1()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSMSusy_high_scale_constraint: Non-perturbative gauge "
"coupling g1 = " << model->get_g1());
#endif
model->set_g1(3.54491);
}
if (std::fabs(model->get_g2()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSMSusy_high_scale_constraint: Non-perturbative gauge "
"coupling g2 = " << model->get_g2());
#endif
model->set_g2(3.54491);
}
if (std::fabs(model->get_g3()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSMSusy_high_scale_constraint: Non-perturbative gauge "
"coupling g3 = " << model->get_g3());
#endif
model->set_g3(3.54491);
}
update_scale();
const auto MuPrInput = INPUTPARAMETER(MuPrInput);
const auto MuPhiInput = INPUTPARAMETER(MuPhiInput);
const auto SigmaLInput = INPUTPARAMETER(SigmaLInput);
const auto KappaPrInput = INPUTPARAMETER(KappaPrInput);
const auto SigmaxInput = INPUTPARAMETER(SigmaxInput);
const auto gDInput = INPUTPARAMETER(gDInput);
const auto hEInput = INPUTPARAMETER(hEInput);
const auto KappaInput = INPUTPARAMETER(KappaInput);
const auto Lambda12Input = INPUTPARAMETER(Lambda12Input);
const auto LambdaxInput = INPUTPARAMETER(LambdaxInput);
const auto fuInput = INPUTPARAMETER(fuInput);
const auto fdInput = INPUTPARAMETER(fdInput);
const auto g1 = MODELPARAMETER(g1);
MODEL->set_g1p(g1);
MODEL->set_MuPr(MuPrInput);
MODEL->set_MuPhi(MuPhiInput);
MODEL->set_SigmaL(SigmaLInput);
MODEL->set_KappaPr(KappaPrInput);
MODEL->set_Sigmax(SigmaxInput);
MODEL->set_gD(gDInput);
MODEL->set_hE(hEInput);
MODEL->set_Kappa(KappaInput);
MODEL->set_Lambda12(Lambda12Input);
MODEL->set_Lambdax(LambdaxInput);
MODEL->set_fu(fuInput);
MODEL->set_fd(fdInput);
{
const auto g1 = MODELPARAMETER(g1);
const auto g2 = MODELPARAMETER(g2);
const auto g3 = MODELPARAMETER(g3);
const auto g1p = MODELPARAMETER(g1p);
const auto Yd = MODELPARAMETER(Yd);
const auto hE = MODELPARAMETER(hE);
const auto Ye = MODELPARAMETER(Ye);
const auto SigmaL = MODELPARAMETER(SigmaL);
const auto KappaPr = MODELPARAMETER(KappaPr);
const auto Sigmax = MODELPARAMETER(Sigmax);
const auto gD = MODELPARAMETER(gD);
const auto Kappa = MODELPARAMETER(Kappa);
const auto Lambda12 = MODELPARAMETER(Lambda12);
const auto Lambdax = MODELPARAMETER(Lambdax);
const auto fu = MODELPARAMETER(fu);
const auto fd = MODELPARAMETER(fd);
const auto Yu = MODELPARAMETER(Yu);
if (MaxAbsValue(g1) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g1", MaxAbsValue(g1), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g1");
if (MaxAbsValue(g2) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g2", MaxAbsValue(g2), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g2");
if (MaxAbsValue(g3) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g3", MaxAbsValue(g3), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g3");
if (MaxAbsValue(g1p) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g1p", MaxAbsValue(g1p), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g1p");
if (MaxAbsValue(Yd) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yd", MaxAbsValue(Yd), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yd");
if (MaxAbsValue(hE) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("hE", MaxAbsValue(hE), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("hE");
if (MaxAbsValue(Ye) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Ye", MaxAbsValue(Ye), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Ye");
if (MaxAbsValue(SigmaL) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("SigmaL", MaxAbsValue(SigmaL), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("SigmaL");
if (MaxAbsValue(KappaPr) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("KappaPr", MaxAbsValue(KappaPr), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("KappaPr");
if (MaxAbsValue(Sigmax) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Sigmax", MaxAbsValue(Sigmax), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Sigmax");
if (MaxAbsValue(gD) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("gD", MaxAbsValue(gD), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("gD");
if (MaxAbsValue(Kappa) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Kappa", MaxAbsValue(Kappa), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Kappa");
if (MaxAbsValue(Lambda12) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Lambda12", MaxAbsValue(Lambda12), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Lambda12");
if (MaxAbsValue(Lambdax) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Lambdax", MaxAbsValue(Lambdax), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Lambdax");
if (MaxAbsValue(fu) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("fu", MaxAbsValue(fu), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("fu");
if (MaxAbsValue(fd) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("fd", MaxAbsValue(fd), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("fd");
if (MaxAbsValue(Yu) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yu", MaxAbsValue(Yu), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yu");
}
}
void CNE6SSM_high_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: CNE6SSM_high_scale_constraint::apply():"
" model pointer must not be zero");
if (std::fabs(model->get_g1()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSM_high_scale_constraint: Non-perturbative gauge "
"coupling g1 = " << model->get_g1());
#endif
model->set_g1(3.54491);
}
if (std::fabs(model->get_g2()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSM_high_scale_constraint: Non-perturbative gauge "
"coupling g2 = " << model->get_g2());
#endif
model->set_g2(3.54491);
}
if (std::fabs(model->get_g3()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CNE6SSM_high_scale_constraint: Non-perturbative gauge "
"coupling g3 = " << model->get_g3());
#endif
model->set_g3(3.54491);
}
update_scale();
const auto MuPrInput = INPUTPARAMETER(MuPrInput);
const auto BMuPrInput = INPUTPARAMETER(BMuPrInput);
const auto MuPhiInput = INPUTPARAMETER(MuPhiInput);
const auto BMuPhiInput = INPUTPARAMETER(BMuPhiInput);
const auto SigmaLInput = INPUTPARAMETER(SigmaLInput);
const auto KappaPrInput = INPUTPARAMETER(KappaPrInput);
const auto SigmaxInput = INPUTPARAMETER(SigmaxInput);
const auto gDInput = INPUTPARAMETER(gDInput);
const auto hEInput = INPUTPARAMETER(hEInput);
const auto KappaInput = INPUTPARAMETER(KappaInput);
const auto Lambda12Input = INPUTPARAMETER(Lambda12Input);
const auto fuInput = INPUTPARAMETER(fuInput);
const auto fdInput = INPUTPARAMETER(fdInput);
const auto Azero = INPUTPARAMETER(Azero);
const auto m0 = INPUTPARAMETER(m0);
const auto m12 = INPUTPARAMETER(m12);
const auto g1 = MODELPARAMETER(g1);
const auto Ye = MODELPARAMETER(Ye);
const auto Yd = MODELPARAMETER(Yd);
const auto Yu = MODELPARAMETER(Yu);
const auto Lambdax = MODELPARAMETER(Lambdax);
MODEL->set_g1p(Re(g1));
MODEL->set_MuPr(Re(MuPrInput));
MODEL->set_BMuPr(Re(BMuPrInput));
MODEL->set_MuPhi(Re(MuPhiInput));
MODEL->set_BMuPhi(Re(BMuPhiInput));
MODEL->set_SigmaL(Re(SigmaLInput));
MODEL->set_KappaPr(Re(KappaPrInput));
MODEL->set_Sigmax(Re(SigmaxInput));
MODEL->set_gD((gDInput).real());
MODEL->set_hE((hEInput).real());
MODEL->set_Kappa((KappaInput).real());
MODEL->set_Lambda12((Lambda12Input).real());
MODEL->set_fu((fuInput).real());
MODEL->set_fd((fdInput).real());
MODEL->set_TYe((Azero*Ye).real());
MODEL->set_TYd((Azero*Yd).real());
MODEL->set_TYu((Azero*Yu).real());
MODEL->set_TKappaPr(Re(Azero*KappaPrInput));
MODEL->set_TSigmax(Re(Azero*SigmaxInput));
MODEL->set_ThE((Azero*hEInput).real());
MODEL->set_TSigmaL(Re(Azero*SigmaLInput));
MODEL->set_TgD((Azero*gDInput).real());
MODEL->set_Tfu((Azero*fuInput).real());
MODEL->set_Tfd((Azero*fdInput).real());
MODEL->set_TKappa((Azero*KappaInput).real());
MODEL->set_TLambda12((Azero*Lambda12Input).real());
MODEL->set_TLambdax(Re(Azero*Lambdax));
MODEL->set_mHd2(Re(Sqr(m0)));
MODEL->set_mHu2(Re(Sqr(m0)));
MODEL->set_ms2(Re(Sqr(m0)));
MODEL->set_msbar2(Re(Sqr(m0)));
MODEL->set_mphi2(Re(Sqr(m0)));
MODEL->set_mHp2(Re(Sqr(m0)));
MODEL->set_mHpbar2(Re(Sqr(m0)));
MODEL->set_mH1I2((Sqr(m0)*UNITMATRIX(2)).real());
MODEL->set_mH2I2((Sqr(m0)*UNITMATRIX(2)).real());
MODEL->set_mSI2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_mq2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_ml2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_md2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_mu2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_me2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_mDx2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_mDxbar2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_MassB(Re(m12));
MODEL->set_MassWB(Re(m12));
MODEL->set_MassG(Re(m12));
MODEL->set_MassBp(Re(m12));
{
const auto g1 = MODELPARAMETER(g1);
const auto g2 = MODELPARAMETER(g2);
const auto g3 = MODELPARAMETER(g3);
const auto g1p = MODELPARAMETER(g1p);
const auto Yd = MODELPARAMETER(Yd);
const auto hE = MODELPARAMETER(hE);
const auto Ye = MODELPARAMETER(Ye);
const auto SigmaL = MODELPARAMETER(SigmaL);
const auto KappaPr = MODELPARAMETER(KappaPr);
const auto Sigmax = MODELPARAMETER(Sigmax);
const auto gD = MODELPARAMETER(gD);
const auto Kappa = MODELPARAMETER(Kappa);
const auto Lambda12 = MODELPARAMETER(Lambda12);
const auto Lambdax = MODELPARAMETER(Lambdax);
const auto fu = MODELPARAMETER(fu);
const auto fd = MODELPARAMETER(fd);
const auto Yu = MODELPARAMETER(Yu);
if (MaxAbsValue(g1) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g1", MaxAbsValue(g1), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g1");
if (MaxAbsValue(g2) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g2", MaxAbsValue(g2), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g2");
if (MaxAbsValue(g3) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g3", MaxAbsValue(g3), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g3");
if (MaxAbsValue(g1p) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g1p", MaxAbsValue(g1p), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g1p");
if (MaxAbsValue(Yd) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yd", MaxAbsValue(Yd), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yd");
if (MaxAbsValue(hE) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("hE", MaxAbsValue(hE), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("hE");
if (MaxAbsValue(Ye) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Ye", MaxAbsValue(Ye), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Ye");
if (MaxAbsValue(SigmaL) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("SigmaL", MaxAbsValue(SigmaL), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("SigmaL");
if (MaxAbsValue(KappaPr) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("KappaPr", MaxAbsValue(KappaPr), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("KappaPr");
if (MaxAbsValue(Sigmax) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Sigmax", MaxAbsValue(Sigmax), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Sigmax");
if (MaxAbsValue(gD) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("gD", MaxAbsValue(gD), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("gD");
if (MaxAbsValue(Kappa) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Kappa", MaxAbsValue(Kappa), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Kappa");
if (MaxAbsValue(Lambda12) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Lambda12", MaxAbsValue(Lambda12), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Lambda12");
if (MaxAbsValue(Lambdax) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Lambdax", MaxAbsValue(Lambdax), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Lambdax");
if (MaxAbsValue(fu) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("fu", MaxAbsValue(fu), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("fu");
if (MaxAbsValue(fd) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("fd", MaxAbsValue(fd), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("fd");
if (MaxAbsValue(Yu) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yu", MaxAbsValue(Yu), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yu");
}
}
void CMSSM_high_scale_constraint<Two_scale>::apply()
{
assert(model && "Error: CMSSM_high_scale_constraint::apply():"
" model pointer must not be zero");
if (std::fabs(model->get_g1()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CMSSM_high_scale_constraint: Non-perturbative gauge "
"coupling g1 = " << model->get_g1());
#endif
model->set_g1(3.54491);
}
if (std::fabs(model->get_g2()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CMSSM_high_scale_constraint: Non-perturbative gauge "
"coupling g2 = " << model->get_g2());
#endif
model->set_g2(3.54491);
}
if (std::fabs(model->get_g3()) > 3.54491) {
#ifdef ENABLE_VERBOSE
ERROR("CMSSM_high_scale_constraint: Non-perturbative gauge "
"coupling g3 = " << model->get_g3());
#endif
model->set_g3(3.54491);
}
update_scale();
const auto Azero = INPUTPARAMETER(Azero);
const auto m0 = INPUTPARAMETER(m0);
const auto m12 = INPUTPARAMETER(m12);
const auto Ye = MODELPARAMETER(Ye);
const auto Yd = MODELPARAMETER(Yd);
const auto Yu = MODELPARAMETER(Yu);
MODEL->set_TYe((Azero*Ye).real());
MODEL->set_TYd((Azero*Yd).real());
MODEL->set_TYu((Azero*Yu).real());
MODEL->set_mHd2(Re(Sqr(m0)));
MODEL->set_mHu2(Re(Sqr(m0)));
MODEL->set_mq2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_ml2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_md2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_mu2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_me2((Sqr(m0)*UNITMATRIX(3)).real());
MODEL->set_MassB(Re(m12));
MODEL->set_MassWB(Re(m12));
MODEL->set_MassG(Re(m12));
{
const auto g1 = MODELPARAMETER(g1);
const auto g2 = MODELPARAMETER(g2);
const auto g3 = MODELPARAMETER(g3);
const auto Yd = MODELPARAMETER(Yd);
const auto Ye = MODELPARAMETER(Ye);
const auto Yu = MODELPARAMETER(Yu);
if (MaxAbsValue(g1) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g1", MaxAbsValue(g1), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g1");
if (MaxAbsValue(g2) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g2", MaxAbsValue(g2), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g2");
if (MaxAbsValue(g3) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("g3", MaxAbsValue(g3), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("g3");
if (MaxAbsValue(Yd) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yd", MaxAbsValue(Yd), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yd");
if (MaxAbsValue(Ye) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Ye", MaxAbsValue(Ye), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Ye");
if (MaxAbsValue(Yu) > 3.5449077018110318)
model->get_problems().flag_non_perturbative_parameter("Yu", MaxAbsValue(Yu), model->get_scale(), 3.5449077018110318);
else
model->get_problems().unflag_non_perturbative_parameter("Yu");
}
}
