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"); } }