void gaugegravityBcs1( MssmSoftsusy & m, const DoubleVector & inputParameters ) { double M_moduli_local, M_gauge_local, M_mess_local ; M_moduli_local = inputParameters.display(1) ; double m1, m2, m3 ; m.setGaugeCoupling( 1, global_g1 ); m.setGaugeCoupling( 2, global_g2 ); m.setGaugeCoupling( 3, global_g3 ); m1 = l1 * M_moduli_local ; m2 = l2 * M_moduli_local ; m3 = l3 * M_moduli_local ; m.setGauginoMass( 1, m1 ) ; m.setGauginoMass( 2, m2 ) ; m.setGauginoMass( 3, m3 ) ; double mqlsq , mllsq, mursq, mdrsq, mersq ; double mhusq , mhdsq; double M_moduli_sqr_local; M_moduli_sqr_local = M_moduli_local*M_moduli_local ; mqlsq = ( 1 - nQ ) * M_moduli_sqr_local ; mllsq = ( 1 - nL ) * M_moduli_sqr_local ; mursq = ( 1 - nU ) * M_moduli_sqr_local ; mdrsq = ( 1 - nD ) * M_moduli_sqr_local ; mersq = ( 1 - nE ) * M_moduli_sqr_local ; mhusq = ( 1 - nHu) * M_moduli_sqr_local ; mhdsq = ( 1 - nHd) * M_moduli_sqr_local ; DoubleMatrix id(3, 3); id(1, 1) = 1.0; id(2, 2) = 1.0; id(3, 3) = 1.0; m.setSoftMassMatrix(mQl, mqlsq * id); m.setSoftMassMatrix(mUr, mursq * id); m.setSoftMassMatrix(mDr, mdrsq * id); m.setSoftMassMatrix(mLl, mllsq * id); m.setSoftMassMatrix(mEr, mersq * id); m.setMh1Squared(mhdsq); m.setMh2Squared(mhusq); double A_HuQU , A_HdQD, A_HdLE ; A_HuQU = ( 3 - nHu - nQ - nU ) * M_moduli_local ; A_HdQD = ( 3 - nHd - nQ - nD ) * M_moduli_local ; A_HdLE = ( 3 - nHd - nL - nE ) * M_moduli_local ; m.setTrilinearElement(UA, 1, 1, m.displayYukawaElement(YU, 1, 1) * A_HuQU); m.setTrilinearElement(UA, 2, 2, m.displayYukawaElement(YU, 2, 2) * A_HuQU); m.setTrilinearElement(UA, 3, 3, m.displayYukawaElement(YU, 3, 3) * A_HuQU); m.setTrilinearElement(DA, 1, 1, m.displayYukawaElement(YD, 1, 1) * A_HdQD); m.setTrilinearElement(DA, 2, 2, m.displayYukawaElement(YD, 2, 2) * A_HdQD); m.setTrilinearElement(DA, 3, 3, m.displayYukawaElement(YD, 3, 3) * A_HdQD); m.setTrilinearElement(EA, 1, 1, m.displayYukawaElement(YE, 1, 1) * A_HdLE); m.setTrilinearElement(EA, 2, 2, m.displayYukawaElement(YE, 2, 2) * A_HdLE); m.setTrilinearElement(EA, 3, 3, m.displayYukawaElement(YE, 3, 3) * A_HdLE); }
void gaugegravityBcs2( MssmSoftsusy & m, const DoubleVector & inputParameters ) { double alpha1, alpha2 , alpha3 ; alpha1 = sqr(m.displayGaugeCoupling(1)) / ( 4.0 * PI ) ; alpha2 = sqr(m.displayGaugeCoupling(2)) / ( 4.0 * PI ) ; alpha3 = sqr(m.displayGaugeCoupling(3)) / ( 4.0 * PI ) ; double M_gauge_local ; M_gauge_local = inputParameters.display(2); m.setGaugeCoupling( 1, global_g1 ); m.setGaugeCoupling( 2, global_g2 ); m.setGaugeCoupling( 3, global_g3 ); double m1 , m2 , m3 ; m1 = inter_gaugino1 - N * alpha1 / (4.0*PI) * M_gauge_local ; m2 = inter_gaugino2 - N * alpha2 / (4.0*PI) * M_gauge_local ; m3 = inter_gaugino3 - N * alpha3 / (4.0*PI) * M_gauge_local ; m.setGauginoMass(1, m1) ; m.setGauginoMass(2, m2) ; m.setGauginoMass(3, m3) ; double mqlsq , mllsq, mursq, mdrsq, mersq ; double mhusq , mhdsq; double M_gauge_sqr_local; M_gauge_sqr_local = sqr( M_gauge_local ) ; mursq = 2.0/sqr( 4.0 * PI )*N * M_gauge_sqr_local * ( 4.0/3.0 * sqr(alpha3) + 0.6 * 4.0 / 9.0 * sqr(alpha1) ) ; mdrsq = 2.0/sqr( 4.0 * PI )*N * M_gauge_sqr_local * ( 4.0/3.0 * sqr(alpha3) + 0.6 * 1.0 / 9.0 * sqr(alpha1) ) ; mersq = 2.0/sqr( 4.0 * PI )*N * M_gauge_sqr_local * ( 0.6 * sqr(alpha1) ) ; mqlsq = 2.0/sqr( 4.0 * PI )*N * M_gauge_sqr_local * ( 4.0/3.0 * sqr(alpha3) + 0.75 * sqr(alpha2) + 0.6 /36.0 * sqr(alpha1)) ; mllsq = 2.0/sqr( 4.0 * PI )*N * M_gauge_sqr_local * ( 0.75 * sqr(alpha2) + 0.6*0.25 * sqr(alpha1) ) ; mhusq = mllsq; mhdsq = mllsq; DoubleMatrix id(3, 3); id(1, 1) = 1.0; id(2, 2) = 1.0; id(3, 3) = 1.0; m.setSoftMassMatrix(mQl, inter_massmQl + mqlsq * id); m.setSoftMassMatrix(mUr, inter_massmUr + mursq * id); m.setSoftMassMatrix(mDr, inter_massmDr + mdrsq * id); m.setSoftMassMatrix(mLl, inter_massmLl + mllsq * id); m.setSoftMassMatrix(mEr, inter_massmEr + mersq * id); m.setMh2Squared(inter_massmHu+ mhusq); m.setMh1Squared(inter_massmHd+ mhdsq); m.setTrilinearElement(UA, 1, 1, m.displayYukawaElement(YU, 1, 1) * inter_A_HuQU(1,1)); m.setTrilinearElement(UA, 2, 2, m.displayYukawaElement(YU, 2, 2) * inter_A_HuQU(2,2)); m.setTrilinearElement(UA, 3, 3, m.displayYukawaElement(YU, 3, 3) * inter_A_HuQU(3,3)); m.setTrilinearElement(DA, 1, 1, m.displayYukawaElement(YD, 1, 1) * inter_A_HdQD(1,1)); m.setTrilinearElement(DA, 2, 2, m.displayYukawaElement(YD, 2, 2) * inter_A_HdQD(2,2)); m.setTrilinearElement(DA, 3, 3, m.displayYukawaElement(YD, 3, 3) * inter_A_HdQD(3,3)); m.setTrilinearElement(EA, 1, 1, m.displayYukawaElement(YE, 1, 1) * inter_A_HdLE(1,1)); m.setTrilinearElement(EA, 2, 2, m.displayYukawaElement(YE, 2, 2) * inter_A_HdLE(2,2)); m.setTrilinearElement(EA, 3, 3, m.displayYukawaElement(YE, 3, 3) * inter_A_HdLE(3,3)); // cout << "In gaugegravityBcs2" << endl; // cout << "gaugino(1) = " << inter_gaugino1 << " " << m1 << endl // << "gaugino(2) = " << inter_gaugino2 << " " << m2 << endl // << "gaugino(3) = " << inter_gaugino3 << " " << m3 << endl; // cout << "inter_massmQl + mqlsq * id (3,3) = " << inter_massmQl(3,3) + mqlsq // << endl; }
void printLineOut(double m0, double m12, double a0, double tanb, MssmSoftsusy twoLoop, MssmSoftsusy twoLoopAs, MssmSoftsusy & twoLoopMt, MssmSoftsusy & twoLoopMb, MssmSoftsusy & threeLoop, double omega2, double omegaAs, double omegaMt, double omegaMb, double omega3, double msqAv2, double msqAvAs, double msqAvMt, double msqAvMb, double msqAv3, double cs, double csAs, double csMt, double csMb, double cs3, double dA, double dAAs, double dAMt, double dAMb, double dA3, double dY, double dYAs, double dYMt, double dYMb, double dY3, double dYp, double dYpAs, double dYpMt, double dYpMb, double dYp3) { /// Any problems and get gnuplot to ignore the point unless it's only musq /// wrong sign... if (twoLoop.displayProblem().test() || threeLoop.displayProblem().test()) if (!testMuSqWrongSignOnlyProblem(twoLoop.displayProblem()) || !testMuSqWrongSignOnlyProblem(threeLoop.displayProblem())) cout << "# "; cout << m0 ///< 1 << " " << m12 ///< 2 << " " << a0 ///< 3 << " " << tanb ///< 4 << " " << twoLoop.displayPhys().mh0(1) ///< 5 << " " << twoLoopAs.displayPhys().mh0(1) ///< 6 << " " << twoLoopMt.displayPhys().mh0(1) ///< 7 << " " << twoLoopMb.displayPhys().mh0(1) ///< 8 << " " << threeLoop.displayPhys().mh0(1) ///< 9 << " " << twoLoop.displayPhys().mA0(1) ///< 10 << " " << twoLoopAs.displayPhys().mA0(1) ///< 11 << " " << twoLoopMt.displayPhys().mA0(1) ///< 12 << " " << twoLoopMb.displayPhys().mA0(1) ///< 13 << " " << threeLoop.displayPhys().mA0(1) ///< 14 << " " << twoLoop.displayPhys().mh0(2) ///< 15 << " " << twoLoopAs.displayPhys().mh0(2) ///< 16 << " " << twoLoopMt.displayPhys().mh0(2) ///< 17 << " " << twoLoopMb.displayPhys().mh0(2) ///< 18 << " " << threeLoop.displayPhys().mh0(2) ///< 19 << " " << twoLoop.displayPhys().mHpm ///< 20 << " " << twoLoopAs.displayPhys().mHpm ///< 21 << " " << twoLoopMt.displayPhys().mHpm ///< 22 << " " << twoLoopMb.displayPhys().mHpm ///< 23 << " " << threeLoop.displayPhys().mHpm ///< 24 << " " << twoLoop.displayPhys().mGluino ///< 25 << " " << twoLoopAs.displayPhys().mGluino ///< 26 << " " << twoLoopMt.displayPhys().mGluino ///< 27 << " " << twoLoopMb.displayPhys().mGluino ///< 28 << " " << threeLoop.displayPhys().mGluino ///< 29 << " " << msqAv2 ///< 30 << " " << msqAvAs ///< 31 << " " << msqAvMt ///< 32 << " " << msqAvMb ///< 33 << " " << msqAv3 ///< 34 << " " << twoLoop.displayPhys().me(1, 1) ///< 35 << " " << twoLoopAs.displayPhys().me(1, 1) ///< 36 << " " << twoLoopMt.displayPhys().me(1, 1) ///< 37 << " " << twoLoopMb.displayPhys().me(1, 1) ///< 38 << " " << threeLoop.displayPhys().me(1, 1) ///< 39 << " " << twoLoop.displayPhys().me(2, 1) ///< 40 << " " << twoLoopAs.displayPhys().me(2, 1) ///< 41 << " " << twoLoopMt.displayPhys().me(2, 1) ///< 42 << " " << twoLoopMb.displayPhys().me(2, 1) ///< 43 << " " << threeLoop.displayPhys().me(2, 1) ///< 44 << " " << fabs(twoLoop.displayPhys().mneut(1)) ///< 45 << " " << fabs(twoLoopAs.displayPhys().mneut(1)) ///< 46 << " " << fabs(twoLoopMt.displayPhys().mneut(1)) ///< 47 << " " << fabs(twoLoopMb.displayPhys().mneut(1)) ///< 48 << " " << fabs(threeLoop.displayPhys().mneut(1)) ///< 49 << " " << fabs(twoLoop.displayPhys().mneut(2)) ///< 50 << " " << fabs(twoLoopAs.displayPhys().mneut(2)) ///< 51 << " " << fabs(twoLoopMt.displayPhys().mneut(2)) ///< 52 << " " << fabs(twoLoopMb.displayPhys().mneut(2)) ///< 53 << " " << fabs(threeLoop.displayPhys().mneut(2)) ///< 54 << " " << fabs(twoLoop.displayPhys().mneut(3)) ///< 55 << " " << fabs(twoLoopAs.displayPhys().mneut(3)) ///< 56 << " " << fabs(twoLoopMt.displayPhys().mneut(3)) ///< 57 << " " << fabs(twoLoopMb.displayPhys().mneut(3)) ///< 58 << " " << fabs(threeLoop.displayPhys().mneut(3)) ///< 59 << " " << fabs(twoLoop.displayPhys().mneut(4)) ///< 60 << " " << fabs(twoLoopAs.displayPhys().mneut(4)) ///< 61 << " " << fabs(twoLoopMt.displayPhys().mneut(4)) ///< 62 << " " << fabs(twoLoopMb.displayPhys().mneut(4)) ///< 63 << " " << fabs(threeLoop.displayPhys().mneut(4)) ///< 64 << " " << twoLoop.displayPhys().mu(1, 3) ///< 65 << " " << twoLoopAs.displayPhys().mu(1, 3) ///< 66 << " " << twoLoopMt.displayPhys().mu(1, 3) ///< 67 << " " << twoLoopMb.displayPhys().mu(1, 3) ///< 68 << " " << threeLoop.displayPhys().mu(1, 3) ///< 69 << " " << twoLoop.displayPhys().mu(2, 3) ///< 70 << " " << twoLoopAs.displayPhys().mu(2, 3) ///< 71 << " " << twoLoopMt.displayPhys().mu(2, 3) ///< 72 << " " << twoLoopMb.displayPhys().mu(2, 3) ///< 73 << " " << threeLoop.displayPhys().mu(2, 3) ///< 74 << " " << twoLoop.displayPhys().md(1, 3) ///< 75 << " " << twoLoopAs.displayPhys().md(1, 3) ///< 76 << " " << twoLoopMt.displayPhys().md(1, 3) ///< 77 << " " << twoLoopMb.displayPhys().md(1, 3) ///< 78 << " " << threeLoop.displayPhys().md(1, 3) ///< 79 << " " << twoLoop.displayPhys().md(2, 3) ///< 80 << " " << twoLoopAs.displayPhys().md(2, 3) ///< 81 << " " << twoLoopMt.displayPhys().md(2, 3) ///< 82 << " " << twoLoopMb.displayPhys().md(2, 3) ///< 83 << " " << threeLoop.displayPhys().md(2, 3) ///< 84 << " " << twoLoop.displayPhys().me(1, 3) ///< 85 << " " << twoLoopAs.displayPhys().me(1, 3) ///< 86 << " " << twoLoopMt.displayPhys().me(1, 3) ///< 87 << " " << twoLoopMb.displayPhys().me(1, 3) ///< 88 << " " << threeLoop.displayPhys().me(1, 3) ///< 89 << " " << twoLoop.displayPhys().me(2, 3) ///< 90 << " " << twoLoopAs.displayPhys().me(2, 3) ///< 91 << " " << twoLoopMt.displayPhys().me(2, 3) ///< 92 << " " << twoLoopMb.displayPhys().me(2, 3) ///< 93 << " " << threeLoop.displayPhys().me(2, 3) ///< 94 << " " << twoLoop.displayYukawaElement(YU, 3, 3) ///< 95 << " " << twoLoopAs.displayYukawaElement(YU, 3, 3) ///< 96 << " " << twoLoopMt.displayYukawaElement(YU, 3, 3) ///< 97 << " " << twoLoopMb.displayYukawaElement(YU, 3, 3) ///< 98 << " " << threeLoop.displayYukawaElement(YU, 3, 3) ///< 99 << " " << twoLoop.displayYukawaElement(YD, 3, 3) ///< 100 << " " << twoLoopAs.displayYukawaElement(YD, 3, 3) ///< 101 << " " << twoLoopMt.displayYukawaElement(YD, 3, 3) ///< 102 << " " << twoLoopMb.displayYukawaElement(YD, 3, 3) ///< 103 << " " << threeLoop.displayYukawaElement(YD, 3, 3) ///< 104 << " " << twoLoop.displayYukawaElement(YE, 3, 3) ///< 105 << " " << twoLoopAs.displayYukawaElement(YE, 3, 3) ///< 106 << " " << twoLoopMt.displayYukawaElement(YE, 3, 3) ///< 107 << " " << twoLoopMb.displayYukawaElement(YE, 3, 3) ///< 108 << " " << threeLoop.displayYukawaElement(YE, 3, 3); ///< 109 int facMusq3= 1., facMusq2 = 1., facMusq2As = 1., facMusq2Mt = 1., facMusq2Mb = 1.; if (twoLoop.displayProblem().muSqWrongSign) facMusq2 = -1. ; if (twoLoopAs.displayProblem().muSqWrongSign) facMusq2As = -1. ; if (twoLoopMt.displayProblem().muSqWrongSign) facMusq2Mt = -1. ; if (twoLoopMb.displayProblem().muSqWrongSign) facMusq2Mb = -1. ; if (threeLoop.displayProblem().muSqWrongSign) facMusq3 = -1. ; cout << " " << twoLoop.displaySusyMu() * facMusq2 ///< 110 << " " << twoLoopAs.displaySusyMu() * facMusq2As ///< 111 << " " << twoLoopMt.displaySusyMu() * facMusq2Mt ///< 112 << " " << twoLoopMb.displaySusyMu() * facMusq2Mb ///< 113 << " " << threeLoop.displaySusyMu() * facMusq3 ///< 114 << " " << fabs(twoLoop.displayPhys().mch(1)) ///< 115 << " " << fabs(twoLoopAs.displayPhys().mch(1)) ///< 116 << " " << fabs(twoLoopMt.displayPhys().mch(1)) ///< 117 << " " << fabs(twoLoopMb.displayPhys().mch(1)) ///< 118 << " " << fabs(threeLoop.displayPhys().mch(1)) ///< 119 << " " << fabs(twoLoop.displayPhys().mch(2)) ///< 120 << " " << fabs(twoLoopAs.displayPhys().mch(2)) ///< 121 << " " << fabs(twoLoopMt.displayPhys().mch(2)) ///< 122 << " " << fabs(twoLoopMb.displayPhys().mch(2)) ///< 123 << " " << fabs(threeLoop.displayPhys().mch(2)) ///< 124 << " " << omega2 ///< 125 << " " << omegaAs ///< 126 << " " << omegaMt ///< 127 << " " << omegaMb ///< 128 << " " << omega3 ///< 129 << " " << twoLoop.displayGaugeCoupling(3) ///< 130 << " " << twoLoopAs.displayGaugeCoupling(3) ///< 131 << " " << twoLoopMt.displayGaugeCoupling(3) ///< 132 << " " << twoLoopMb.displayGaugeCoupling(3) ///< 133 << " " << threeLoop.displayGaugeCoupling(3) ///< 134 << " " << cs ///< 135 << " " << csAs ///< 136 << " " << csMt ///< 137 << " " << csMb ///< 138 << " " << cs3 ///< 139 << " " << QEWSB ///< 140 << " " << dA ///< 141 << " " << dAAs ///< 142 << " " << dAMt ///< 143 << " " << dAMb ///< 144 << " " << dA3 ///< 145 << " " << dY ///< 146 << " " << dYAs ///< 147 << " " << dYMt ///< 148 << " " << dYMb ///< 149 << " " << dY3 ///< 150 << " " << dYp ///< 151 << " " << dYpAs ///< 152 << " " << dYpMt ///< 153 << " " << dYpMb ///< 154 << " " << dYp3; ///< 155 if (twoLoop.displayProblem().test()) cout << "# 2-loop problem: " << twoLoop.displayProblem(); if (threeLoop.displayProblem().test()) cout << "# 3-loop problem " << threeLoop.displayProblem(); cout << endl; }