diff --git a/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.cxx b/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.cxx
index 55b6e6db29748b5f391eb905f166cdd04f4ff656..660e3a504a05dc4a60e53c41ec44986dab6989b9 100644
--- a/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.cxx
+++ b/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.cxx
@@ -192,7 +192,7 @@ std::vector<double> iFatras::MultipleScatteringSamplerGeneralMixture::getSemigau
   return scattering_params;
 }
 
-double iFatras::MultipleScatteringSamplerGeneralMixture::simGaussmix(std::vector<double> scattering_params) const{
+double iFatras::MultipleScatteringSamplerGeneralMixture::simGaussmix(const std::vector<double>& scattering_params) const{
   double sigma_tot = scattering_params[0];
   double var1 = scattering_params[1];
   double var2 = scattering_params[2];
@@ -205,7 +205,7 @@ double iFatras::MultipleScatteringSamplerGeneralMixture::simGaussmix(std::vector
   return sqrt(var2)*sqrt(-2*log(u))*sigma_tot;
 }
 
-double iFatras::MultipleScatteringSamplerGeneralMixture::simSemigauss(std::vector<double> scattering_params) const{
+double iFatras::MultipleScatteringSamplerGeneralMixture::simSemigauss(const std::vector<double>& scattering_params) const{
   double a = scattering_params[0];
   double b = scattering_params[1];
   double var1 = scattering_params[2];
diff --git a/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.h b/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.h
index 4e307daa199f9c853af178f58756e65979ef9d4d..9e77fc601d2253a0b16a12561a9eef0b88062dc7 100644
--- a/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.h
+++ b/Simulation/ISF/ISF_Fatras/ISF_FatrasTools/src/MultipleScatteringSamplerGeneralMixture.h
@@ -75,9 +75,9 @@ namespace iFatras {
       //!< General mixture model: get parameters for semi-gaussian mixture
       static std::vector<double> 	getSemigauss(double beta,double p,double dOverX0,double Z, double scale) ;
       //!< General mixture model: simulate semi-gaussian mixture
-      double 	simGaussmix(std::vector<double> scattering_params) const;
+      double 	simGaussmix(const std::vector<double>& scattering_params) const;
       //!< General mixture model: simulate gaussian mixture
-      double 	simSemigauss(std::vector<double> scattering_params) const;
+      double 	simSemigauss(const std::vector<double>& scattering_params) const;
     };