diff --git a/src/lowe.cpp b/src/lowe.cpp
index 52521d86a..4d44da7ae 100644
--- a/src/lowe.cpp
+++ b/src/lowe.cpp
@@ -65,9 +65,9 @@ const std::array<std::string, NUMBER_OF_LOW_ENERGY_INPUT_PARAMETERS> QedQcd_inpu
    "GFermi",
    "MZ_pole", "MW_pole",
    "Mv1_pole", "Mv2_pole", "Mv3_pole",
-   "MElectron_pole", "MMuon_pole", "MTau_pole",
-   "MU_2GeV", "MS_2GeV", "MT_pole",
-   "MD_2GeV", "mc_mc", "mb_mb",
+   "Me_pole", "Mm_pole", "Mtau_pole",
+   "mu_2GeV", "ms_2GeV", "Mt_pole",
+   "md_2GeV", "mc_mc", "mb_mb",
    "CKM_theta_12", "CKM_theta_13", "CKM_theta_23", "CKM_delta",
    "PMNS_theta_12", "PMNS_theta_13", "PMNS_theta_23", "PMNS_delta", "PMNS_alpha_1", "PMNS_alpha_2"
 };
@@ -129,18 +129,18 @@ QedQcd::QedQcd()
   a(1) = flexiblesusy::Electroweak_constants::alpha3;
   input(alpha_em_MSbar_at_MZ) = flexiblesusy::Electroweak_constants::aem;
   input(alpha_s_MSbar_at_MZ) = flexiblesusy::Electroweak_constants::alpha3;
-  input(MT_pole) = flexiblesusy::Electroweak_constants::PMTOP;
+  input(Mt_pole) = flexiblesusy::Electroweak_constants::PMTOP;
   input(mb_mb) = flexiblesusy::Electroweak_constants::MBOTTOM;
-  input(MTau_pole) = flexiblesusy::Electroweak_constants::MTAU;
-  input(MMuon_pole) = flexiblesusy::Electroweak_constants::MMUON;
-  input(MElectron_pole) = flexiblesusy::Electroweak_constants::MELECTRON;
+  input(Mtau_pole) = flexiblesusy::Electroweak_constants::MTAU;
+  input(Mm_pole) = flexiblesusy::Electroweak_constants::MMUON;
+  input(Me_pole) = flexiblesusy::Electroweak_constants::MELECTRON;
   input(MW_pole) = flexiblesusy::Electroweak_constants::MW;
   input(MZ_pole) = flexiblesusy::Electroweak_constants::MZ;
   input(GFermi) = flexiblesusy::Electroweak_constants::gfermi;
   input(mc_mc) = flexiblesusy::Electroweak_constants::MCHARM;
-  input(MU_2GeV) = flexiblesusy::Electroweak_constants::MUP;
-  input(MD_2GeV) = flexiblesusy::Electroweak_constants::MDOWN;
-  input(MS_2GeV) = flexiblesusy::Electroweak_constants::MSTRANGE;
+  input(mu_2GeV) = flexiblesusy::Electroweak_constants::MUP;
+  input(md_2GeV) = flexiblesusy::Electroweak_constants::MDOWN;
+  input(ms_2GeV) = flexiblesusy::Electroweak_constants::MSTRANGE;
   set_scale(flexiblesusy::Electroweak_constants::MZ);
   set_loops(3);
   set_thresholds(1);
diff --git a/src/lowe.h b/src/lowe.h
index da0f18754..02f67469e 100644
--- a/src/lowe.h
+++ b/src/lowe.h
@@ -49,9 +49,9 @@ enum QedQcd_input_parmeters : int {
    GFermi,
    MZ_pole, MW_pole,
    Mv1_pole, Mv2_pole, Mv3_pole,
-   MElectron_pole, MMuon_pole, MTau_pole,
-   MU_2GeV, MS_2GeV, MT_pole,
-   MD_2GeV, mc_mc, mb_mb,
+   Me_pole, Mm_pole, Mtau_pole,
+   mu_2GeV, ms_2GeV, Mt_pole,
+   md_2GeV, mc_mc, mb_mb,
    CKM_theta_12, CKM_theta_13, CKM_theta_23, CKM_delta,
    PMNS_theta_12, PMNS_theta_13, PMNS_theta_23, PMNS_delta, PMNS_alpha_1, PMNS_alpha_2,
    NUMBER_OF_LOW_ENERGY_INPUT_PARAMETERS
@@ -99,16 +99,16 @@ class QedQcd: public flexiblesusy::Beta_function
   virtual void set(const Eigen::ArrayXd&) override;
   virtual Eigen::ArrayXd beta() const override;
 
-  void setPoleMt(double mt) { input(MT_pole) = mt; }; ///< set pole top mass
+  void setPoleMt(double mt) { input(Mt_pole) = mt; }; ///< set pole top mass
   void setPoleMb(double mb) { mbPole = mb; }; ///< set pole bottom mass
-  void setPoleMtau(double mtau) { input(MTau_pole) = mtau; }; ///< set pole tau mass
-  void setPoleMmuon(double m) { input(MMuon_pole) = m; } ///< set pole muon mass
-  void setPoleMel(double m) { input(MElectron_pole) = m; } ///< set pole electron mass
+  void setPoleMtau(double mtau) { input(Mtau_pole) = mtau; }; ///< set pole tau mass
+  void setPoleMmuon(double m) { input(Mm_pole) = m; } ///< set pole muon mass
+  void setPoleMel(double m) { input(Me_pole) = m; } ///< set pole electron mass
   void setMbMb(double mb)   { input(mb_mb) = mb;   }; ///< set mb(mb)
   void setMcMc(double mc)   { input(mc_mc) = mc;   }  ///< set mc(mc)
-  void setMu2GeV(double mu) { input(MU_2GeV) = mu; } ///< set mu(2 GeV)
-  void setMd2GeV(double md) { input(MD_2GeV) = md; } ///< set md(2 GeV)
-  void setMs2GeV(double ms) { input(MS_2GeV) = ms; } ///< set ms(2 GeV)
+  void setMu2GeV(double mu) { input(mu_2GeV) = mu; } ///< set mu(2 GeV)
+  void setMd2GeV(double md) { input(md_2GeV) = md; } ///< set md(2 GeV)
+  void setMs2GeV(double ms) { input(ms_2GeV) = ms; } ///< set ms(2 GeV)
   void setPoleMW(double mw) { input(MW_pole) = mw; } ///< set W boson pole mass
   void setPoleMZ(double mz) { input(MZ_pole) = mz; } ///< set Z boson pole mass
   /// sets a running quark mass
@@ -133,13 +133,13 @@ class QedQcd: public flexiblesusy::Beta_function
   /// Displays input parameters
   Input_t displayInput() const { return input; }
   /// Display pole top mass
-  double displayPoleMt() const { return input(MT_pole); };
+  double displayPoleMt() const { return input(Mt_pole); };
   /// Display pole tau mass
-  double displayPoleMtau() const { return input(MTau_pole); };
+  double displayPoleMtau() const { return input(Mtau_pole); };
   /// Display pole muon mass
-  double displayPoleMmuon() const { return input(MMuon_pole); };
+  double displayPoleMmuon() const { return input(Mm_pole); };
   /// Display pole electron mass
-  double displayPoleMel() const { return input(MElectron_pole); };
+  double displayPoleMel() const { return input(Me_pole); };
   /// Returns bottom "pole" mass
   double displayPoleMb() const { return mbPole; };
   /// Returns W boson pole mass
@@ -171,11 +171,11 @@ class QedQcd: public flexiblesusy::Beta_function
   /// Returns mc(mc) MSbar
   double displayMcMc() const { return input(mc_mc); }
   /// Returns mu(2 GeV)
-  double displayMu2GeV() const { return input(MU_2GeV); }
+  double displayMu2GeV() const { return input(mu_2GeV); }
   /// Returns md(2 GeV)
-  double displayMd2GeV() const { return input(MD_2GeV); }
+  double displayMd2GeV() const { return input(md_2GeV); }
   /// Returns ms(2 GeV)
-  double displayMs2GeV() const { return input(MS_2GeV); }
+  double displayMs2GeV() const { return input(ms_2GeV); }
   /// returns CKM parameters
   flexiblesusy::CKM_parameters displayCKM() const;
   /// Returns real CKM matrix